File: pr77286.c

package info (click to toggle)
gcc-arm-none-eabi 15%3A8-2019-q3-1
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 571,828 kB
  • sloc: ansic: 2,937,651; cpp: 881,644; ada: 597,189; makefile: 65,528; asm: 56,499; xml: 46,621; exp: 24,747; sh: 19,684; python: 7,256; pascal: 4,370; awk: 3,497; perl: 2,695; yacc: 316; ml: 285; f90: 234; lex: 198; objc: 194; haskell: 119
file content (43 lines) | stat: -rw-r--r-- 939 bytes parent folder | download | duplicates (6) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
 
/* { dg-do compile } */
/* { dg-additional-options "-mavx2" { target x86_64-*-* i?86-*-* } } */

typedef float real;
typedef struct
{
  int ngtc;
  real *ref_t;
  real *tau_t;
} t_grpopts;
typedef struct
{
  real T;
  real xi;
} t_grp_tcstat;
typedef struct
{
  t_grp_tcstat *tcstat;
} t_groups;
extern real *save_calloc ();
void
nosehoover_tcoupl (t_grpopts * opts, t_groups * grps, real dt, real SAfactor)
{
  static real *Qinv = ((void *) 0);
  int i;
  real reft = 0, xit, oldxi;
  if (Qinv == ((void *) 0))
    {
      (Qinv) =
	save_calloc ("Qinv", "coupling.c", 372, (opts->ngtc),
		     sizeof (*(Qinv)));
      for (i = 0; i < opts->ngtc; i++)
	if ((opts->tau_t[i] > 0))
	  Qinv[i] = 1.0 / opts->tau_t[i];
    }
  for (i = 0; (i < opts->ngtc); i++)
    {
      reft =
	(((0.0) >
	  (opts->ref_t[i] * SAfactor)) ? (0.0) : (opts->ref_t[i] * SAfactor));
      grps->tcstat[i].xi += dt * Qinv[i] * (grps->tcstat[i].T - reft);
    }
}