diff options
Diffstat (limited to 'ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc')
| -rw-r--r-- | ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc | 9 |
1 files changed, 8 insertions, 1 deletions
diff --git a/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc b/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc index ba11ca7..72a6512 100644 --- a/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc +++ b/ML_BSSN/src/ML_BSSN_RHSStaticBoundary.cc @@ -53,6 +53,9 @@ extern "C" void ML_BSSN_RHSStaticBoundary_SelectBCs(CCTK_ARGUMENTS) ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_shiftrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_shiftrhs."); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_Thetarhs","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_Thetarhs."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN::ML_trace_curvrhs","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN::ML_trace_curvrhs."); @@ -292,6 +295,8 @@ static void ML_BSSN_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in CCTK_REAL_VEC Xt3rhsL = ToReal(0); + CCTK_REAL_VEC ThetarhsL = ToReal(0); + CCTK_REAL_VEC alpharhsL = ToReal(0); CCTK_REAL_VEC ArhsL = ToReal(0); @@ -331,6 +336,7 @@ static void ML_BSSN_RHSStaticBoundary_Body(cGH const * restrict const cctkGH, in 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); @@ -364,8 +370,9 @@ extern "C" void ML_BSSN_RHSStaticBoundary(CCTK_ARGUMENTS) "ML_BSSN::ML_log_confacrhs", "ML_BSSN::ML_metricrhs", "ML_BSSN::ML_shiftrhs", + "ML_BSSN::ML_Thetarhs", "ML_BSSN::ML_trace_curvrhs"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_RHSStaticBoundary", 9, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_RHSStaticBoundary", 10, groups); switch(fdOrder) { |