diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc | 8 |
1 files changed, 2 insertions, 6 deletions
diff --git a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc index a8e97d0..5944c01 100644 --- a/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc +++ b/ML_BSSN/src/ML_BSSN_convertToADMBaseFakeDtLapseShift.cc @@ -208,7 +208,6 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(cGH const * restrict c CCTK_REAL_VEC beta2L = vec_load(beta2[index]); CCTK_REAL_VEC beta3L = vec_load(beta3[index]); CCTK_REAL_VEC rL = vec_load(r[index]); - CCTK_REAL_VEC ThetaL = vec_load(Theta[index]); CCTK_REAL_VEC trKL = vec_load(trK[index]); CCTK_REAL_VEC Xt1L = vec_load(Xt1[index]); CCTK_REAL_VEC Xt2L = vec_load(Xt2[index]); @@ -243,9 +242,7 @@ static void ML_BSSN_convertToADMBaseFakeDtLapseShift_Body(cGH const * restrict c kfmin(ToReal(1),kexp(knmsub(rL,INV(ToReal(SpatialShiftGammaCoeffRadius)),ToReal(1)))); CCTK_REAL_VEC dtalpL = - kneg(kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kadd(trKL,kmadd(ksub(AL,trKL),ToReal(LapseACoeff),kmul(IfThen(conformalMethod - == - 2,kmul(ThetaL,ToReal(2)),ToReal(0)),kadd(ToReal(-1),ToReal(LapseACoeff)))))))); + kneg(kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kmadd(ksub(AL,trKL),ToReal(LapseACoeff),trKL)))); CCTK_REAL_VEC dtbetaxL; CCTK_REAL_VEC dtbetayL; @@ -307,9 +304,8 @@ extern "C" void ML_BSSN_convertToADMBaseFakeDtLapseShift(CCTK_ARGUMENTS) "ML_BSSN::ML_Gamma", "ML_BSSN::ML_lapse", "ML_BSSN::ML_shift", - "ML_BSSN::ML_Theta", "ML_BSSN::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_convertToADMBaseFakeDtLapseShift", 11, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_convertToADMBaseFakeDtLapseShift", 10, groups); switch(fdOrder) { |