diff options
Diffstat (limited to 'ML_BSSN_UPW/src/ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary.cc')
| -rw-r--r-- | ML_BSSN_UPW/src/ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary.cc | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/ML_BSSN_UPW/src/ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary.cc b/ML_BSSN_UPW/src/ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary.cc index 723f9b0..9a56020 100644 --- a/ML_BSSN_UPW/src/ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary.cc +++ b/ML_BSSN_UPW/src/ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary.cc @@ -223,6 +223,7 @@ static void ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * re 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]); @@ -257,7 +258,9 @@ static void ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary_Body(cGH const * re kfmin(ToReal(1),kexp(knmsub(rL,INV(ToReal(SpatialShiftGammaCoeffRadius)),ToReal(1)))); CCTK_REAL_VEC dtalpL = - kneg(kmul(kpow(alphaL,harmonicN),kmul(ToReal(harmonicF),kmadd(ksub(AL,trKL),ToReal(LapseACoeff),trKL)))); + 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)))))))); CCTK_REAL_VEC dtbetaxL; CCTK_REAL_VEC dtbetayL; @@ -319,8 +322,9 @@ extern "C" void ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary(CCTK_ARGUMENTS) "ML_BSSN_UPW::ML_Gamma", "ML_BSSN_UPW::ML_lapse", "ML_BSSN_UPW::ML_shift", + "ML_BSSN_UPW::ML_Theta", "ML_BSSN_UPW::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary", 10, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_UPW_convertToADMBaseDtLapseShiftBoundary", 11, groups); switch(fdOrder) { |