diff options
Diffstat (limited to 'ML_BSSN_O2/src/ML_BSSN_O2_InitRHS.cc')
-rw-r--r-- | ML_BSSN_O2/src/ML_BSSN_O2_InitRHS.cc | 6 |
1 files changed, 5 insertions, 1 deletions
diff --git a/ML_BSSN_O2/src/ML_BSSN_O2_InitRHS.cc b/ML_BSSN_O2/src/ML_BSSN_O2_InitRHS.cc index 9e86ac6..fef5c39 100644 --- a/ML_BSSN_O2/src/ML_BSSN_O2_InitRHS.cc +++ b/ML_BSSN_O2/src/ML_BSSN_O2_InitRHS.cc @@ -256,6 +256,8 @@ static void ML_BSSN_O2_InitRHS_Body(cGH const * restrict const cctkGH, int const CCTK_REAL_VEC Xt3rhsL = ToReal(0); + CCTK_REAL_VEC ThetarhsL = ToReal(0); + CCTK_REAL_VEC alpharhsL = ToReal(0); CCTK_REAL_VEC ArhsL = ToReal(0); @@ -295,6 +297,7 @@ static void ML_BSSN_O2_InitRHS_Body(cGH const * restrict const cctkGH, int const vec_store_nta_partial(gt23rhs[index],gt23rhsL); vec_store_nta_partial(gt33rhs[index],gt33rhsL); vec_store_nta_partial(phirhs[index],phirhsL); + vec_store_nta_partial(Thetarhs[index],ThetarhsL); vec_store_nta_partial(trKrhs[index],trKrhsL); vec_store_nta_partial(Xt1rhs[index],Xt1rhsL); vec_store_nta_partial(Xt2rhs[index],Xt2rhsL); @@ -328,8 +331,9 @@ extern "C" void ML_BSSN_O2_InitRHS(CCTK_ARGUMENTS) "ML_BSSN_O2::ML_log_confacrhs", "ML_BSSN_O2::ML_metricrhs", "ML_BSSN_O2::ML_shiftrhs", + "ML_BSSN_O2::ML_Thetarhs", "ML_BSSN_O2::ML_trace_curvrhs"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_O2_InitRHS", 9, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_O2_InitRHS", 10, groups); switch(fdOrder) { |