1  /* { dg-do compile } */
       2  /* { dg-additional-options "-ffast-math" } */
       3  
       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  }