aboutsummaryrefslogtreecommitdiff
path: root/ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc
diff options
context:
space:
mode:
Diffstat (limited to 'ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc')
-rw-r--r--ML_BSSN_O8/src/ML_BSSN_O8_convertFromADMBase.cc8
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)
{