1 /* { dg-do compile } */
2 /* { dg-additional-options "-ffast-math" } */
3
metric_carttosphere(int * cctk_lsh,double txz,double tyz,double txx,double tzz,double sint,double cosp,double cost,double tyy,double sinp,double txy,double * grp,double * grq,double * r)4 void metric_carttosphere(int *cctk_lsh, double txz, double tyz, double txx,
5 double tzz, double sint, double cosp, double cost,
6 double tyy, double sinp, double txy, double *grp,
7 double *grq, double *r)
8 {
9 int i;
10 for(i=0; i<cctk_lsh[0]*cctk_lsh[1]*cctk_lsh[2]; i++)
11 {
12 grq[i] = (cost*tyy*((sinp)*(sinp))*sint+
13 2*cosp*cost*txy*sinp*sint-
14 cost*tzz*sint+
15 ((cosp)*(cosp))*cost*txx*sint+
16 2*((cost)*(cost))*tyz*sinp-
17 tyz*sinp+
18 2*cosp*((cost)*(cost))*txz-
19 cosp*txz)*r[i];
20 grp[i] = ((-txy*((sinp)*(sinp))+
21 (cosp*tyy-cosp*txx)*sinp+
22 ((cosp)*(cosp))*txy)*sint-
23 cost*txz*sinp+cosp*cost*tyz);
24 }
25 }
26