diff options
Diffstat (limited to 'ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc')
-rw-r--r-- | ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc b/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc index 5203f5b..677f542 100644 --- a/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc +++ b/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc @@ -269,7 +269,7 @@ static void ML_BSSN_O8_convertFromADMBase_Body(cGH const * restrict const cctkGH CCTK_REAL_VEC em4phi; - if (conformalMethod) + if (conformalMethod == 1) { phiL = kpow(detg,-0.166666666666666666666666666667); @@ -315,6 +315,8 @@ static void ML_BSSN_O8_convertFromADMBase_Body(cGH const * restrict const cctkGH CCTK_REAL_VEC At33L = kmul(em4phi,kmadd(trKL,kmul(g33,ToReal(-0.333333333333333333333333333333)),kzzL)); + CCTK_REAL_VEC ThetaL = ToReal(0); + CCTK_REAL_VEC alphaL = alpL; CCTK_REAL_VEC beta1L = betaxL; @@ -342,6 +344,7 @@ static void ML_BSSN_O8_convertFromADMBase_Body(cGH const * restrict const cctkGH vec_store_nta_partial(gt23[index],gt23L); vec_store_nta_partial(gt33[index],gt33L); vec_store_nta_partial(phi[index],phiL); + vec_store_nta_partial(Theta[index],ThetaL); vec_store_nta_partial(trK[index],trKL); } LC_ENDLOOP3VEC(ML_BSSN_O8_convertFromADMBase); @@ -373,8 +376,9 @@ extern "C" void ML_BSSN_O8_convertFromADMBase(CCTK_ARGUMENTS) "ML_BSSN_O8::ML_log_confac", "ML_BSSN_O8::ML_metric", "ML_BSSN_O8::ML_shift", + "ML_BSSN_O8::ML_Theta", "ML_BSSN_O8::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_O8_convertFromADMBase", 10, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_O8_convertFromADMBase", 11, groups); switch(fdOrder) { |