diff options
Diffstat (limited to 'ML_CCZ4/src/ML_CCZ4_convertFromADMBase.cc')
-rw-r--r-- | ML_CCZ4/src/ML_CCZ4_convertFromADMBase.cc | 20 |
1 files changed, 10 insertions, 10 deletions
diff --git a/ML_CCZ4/src/ML_CCZ4_convertFromADMBase.cc b/ML_CCZ4/src/ML_CCZ4_convertFromADMBase.cc index 8995328..1e0e4c3 100644 --- a/ML_CCZ4/src/ML_CCZ4_convertFromADMBase.cc +++ b/ML_CCZ4/src/ML_CCZ4_convertFromADMBase.cc @@ -117,12 +117,12 @@ static void ML_CCZ4_convertFromADMBase_Body(cGH const * restrict const cctkGH, i 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); @@ -257,7 +257,7 @@ static void ML_CCZ4_convertFromADMBase_Body(cGH const * restrict const cctkGH, i CCTK_REAL_VEC g33 = gzzL; CCTK_REAL_VEC detg = - knmsub(g22,kmul(g13,g13),knmsub(g11,kmul(g23,g23),kmadd(g33,kmsub(g11,g22,kmul(g12,g12)),kmul(g12,kmul(g13,kmul(g23,ToReal(2.))))))); + knmsub(g22,kmul(g13,g13),knmsub(g11,kmul(g23,g23),kmadd(g33,kmsub(g11,g22,kmul(g12,g12)),kmul(g12,kmul(g13,kmul(g23,ToReal(2))))))); CCTK_REAL_VEC gu11 = kdiv(kmsub(g22,g33,kmul(g23,g23)),detg); @@ -283,7 +283,7 @@ static void ML_CCZ4_convertFromADMBase_Body(cGH const * restrict const cctkGH, i { phiL = kmul(klog(detg),ToReal(0.0833333333333333333333333333333)); - em4phi = kexp(kmul(phiL,ToReal(-4.))); + em4phi = kexp(kmul(phiL,ToReal(-4))); } CCTK_REAL_VEC gt11L = kmul(em4phi,g11); @@ -299,7 +299,7 @@ static void ML_CCZ4_convertFromADMBase_Body(cGH const * restrict const cctkGH, i CCTK_REAL_VEC gt33L = kmul(em4phi,g33); trKL = - kmadd(kxxL,gu11,kmadd(kyyL,gu22,kmadd(kzzL,gu33,kmul(kmadd(kxyL,gu12,kmadd(kxzL,gu13,kmul(kyzL,gu23))),ToReal(2.))))); + kmadd(kxxL,gu11,kmadd(kyyL,gu22,kmadd(kzzL,gu33,kmul(kmadd(kxyL,gu12,kmadd(kxzL,gu13,kmul(kyzL,gu23))),ToReal(2))))); CCTK_REAL_VEC At11L = kmul(em4phi,kmadd(trKL,kmul(g11,ToReal(-0.333333333333333333333333333333)),kxxL)); @@ -327,7 +327,7 @@ static void ML_CCZ4_convertFromADMBase_Body(cGH const * restrict const cctkGH, i CCTK_REAL_VEC beta3L = betazL; - CCTK_REAL_VEC ThetaL = ToReal(0.); + CCTK_REAL_VEC ThetaL = ToReal(0); /* Copy local copies back to grid functions */ vec_store_partial_prepare(i,lc_imin,lc_imax); |