diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_Minkowski.cc')
-rw-r--r-- | ML_BSSN/src/ML_BSSN_Minkowski.cc | 8 |
1 files changed, 6 insertions, 2 deletions
diff --git a/ML_BSSN/src/ML_BSSN_Minkowski.cc b/ML_BSSN/src/ML_BSSN_Minkowski.cc index f1033d0..1aacb7f 100644 --- a/ML_BSSN/src/ML_BSSN_Minkowski.cc +++ b/ML_BSSN/src/ML_BSSN_Minkowski.cc @@ -222,7 +222,7 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const } /* Calculate temporaries and grid functions */ - CCTK_REAL_VEC phiL = IfThen(conformalMethod,ToReal(1),ToReal(0)); + CCTK_REAL_VEC phiL = IfThen(conformalMethod == 1,ToReal(1),ToReal(0)); CCTK_REAL_VEC gt11L = ToReal(1); @@ -256,6 +256,8 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC Xt3L = ToReal(0); + CCTK_REAL_VEC ThetaL = ToReal(0); + CCTK_REAL_VEC alphaL = ToReal(1); CCTK_REAL_VEC AL = ToReal(0); @@ -295,6 +297,7 @@ static void ML_BSSN_Minkowski_Body(cGH const * restrict const cctkGH, int const 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); vec_store_nta_partial(Xt1[index],Xt1L); vec_store_nta_partial(Xt2[index],Xt2L); @@ -328,8 +331,9 @@ extern "C" void ML_BSSN_Minkowski(CCTK_ARGUMENTS) "ML_BSSN::ML_log_confac", "ML_BSSN::ML_metric", "ML_BSSN::ML_shift", + "ML_BSSN::ML_Theta", "ML_BSSN::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Minkowski", 9, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_Minkowski", 10, groups); switch(fdOrder) { |