diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc')
-rw-r--r-- | ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc b/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc index c50a732..c4965b9 100644 --- a/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc +++ b/ML_CCZ4/src/ML_CCZ4_convertToADMBase.cc @@ -117,12 +117,12 @@ static void ML_CCZ4_convertToADMBase_Body(cGH const * restrict const cctkGH, int CCTK_REAL_VEC const p1o840dx = kdiv(ToReal(0.00119047619047619047619047619048),dx); CCTK_REAL_VEC const p1o840dy = kdiv(ToReal(0.00119047619047619047619047619048),dy); CCTK_REAL_VEC const p1o840dz = kdiv(ToReal(0.00119047619047619047619047619048),dz); - CCTK_REAL_VEC const p1odx = kdiv(ToReal(1),dx); - CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1),kmul(dx,dx)); - CCTK_REAL_VEC const p1ody = kdiv(ToReal(1),dy); - CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1),kmul(dy,dy)); - CCTK_REAL_VEC const p1odz = kdiv(ToReal(1),dz); - CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1),kmul(dz,dz)); + CCTK_REAL_VEC const p1odx = kdiv(ToReal(1.),dx); + CCTK_REAL_VEC const p1odx2 = kdiv(ToReal(1.),kmul(dx,dx)); + CCTK_REAL_VEC const p1ody = kdiv(ToReal(1.),dy); + CCTK_REAL_VEC const p1ody2 = kdiv(ToReal(1.),kmul(dy,dy)); + CCTK_REAL_VEC const p1odz = kdiv(ToReal(1.),dz); + CCTK_REAL_VEC const p1odz2 = kdiv(ToReal(1.),kmul(dz,dz)); CCTK_REAL_VEC const pm1o120dx = kdiv(ToReal(-0.00833333333333333333333333333333),dx); CCTK_REAL_VEC const pm1o120dy = kdiv(ToReal(-0.00833333333333333333333333333333),dy); CCTK_REAL_VEC const pm1o120dz = kdiv(ToReal(-0.00833333333333333333333333333333),dz); @@ -251,7 +251,7 @@ static void ML_CCZ4_convertToADMBase_Body(cGH const * restrict const cctkGH, int /* Calculate temporaries and grid functions */ CCTK_REAL_VEC e4phi = IfThen(conformalMethod == - 1,kdiv(ToReal(1),kmul(phiL,phiL)),kexp(kmul(phiL,ToReal(4)))); + 1.,kdiv(ToReal(1.),kmul(phiL,phiL)),kexp(kmul(phiL,ToReal(4.)))); gxxL = kmul(gt11L,e4phi); |