1 /* { dg-do compile } */
2 /* { dg-options "-Ofast -floop-nest-optimize" } */
4 typedef struct _IO_FILE
FILE;
5 extern struct _IO_FILE
*stderr
;
9 void ps_box (int, int);
10 void plot_phi(char *fn
,rvec box
,int natoms
,rvec x
[],real phi
[])
12 real phi_max
,rr
,gg
,bb
,fac
,dx
,x0
,y0
;
14 for(i
=0; (i
<natoms
); i
++)
15 phi_max
=((phi_max
> __builtin_fabs(phi
[i
]))
16 ? phi_max
: __builtin_fabs(phi
[i
]));
17 if (__builtin_fabs(phi_max
)<1.2e-38)
18 __builtin_fprintf(stderr
, "X");
19 ps_box((real
)(fac
*box
[0]-1),(real
)(fac
*box
[1]-1));
20 for(i
=0; (i
<natoms
); i
++)
24 gg
=bb
=(1.0+(phi
[i
]/phi_max
));
26 rr
=gg
=(1.0-(phi
[i
]/phi_max
));