diff options
Diffstat (limited to 'ML_BSSN_MP/src/ML_BSSN_MP_boundary.cc')
-rw-r--r-- | ML_BSSN_MP/src/ML_BSSN_MP_boundary.cc | 11 |
1 files changed, 9 insertions, 2 deletions
diff --git a/ML_BSSN_MP/src/ML_BSSN_MP_boundary.cc b/ML_BSSN_MP/src/ML_BSSN_MP_boundary.cc index f7985cc..299a0d0 100644 --- a/ML_BSSN_MP/src/ML_BSSN_MP_boundary.cc +++ b/ML_BSSN_MP/src/ML_BSSN_MP_boundary.cc @@ -53,6 +53,9 @@ extern "C" void ML_BSSN_MP_boundary_SelectBCs(CCTK_ARGUMENTS) ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_MP::ML_shift","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_MP::ML_shift."); + ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_MP::ML_Theta","flat"); + if (ierr < 0) + CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_MP::ML_Theta."); ierr = Boundary_SelectGroupForBC(cctkGH, CCTK_ALL_FACES, GenericFD_GetBoundaryWidth(cctkGH), -1 /* no table */, "ML_BSSN_MP::ML_trace_curv","flat"); if (ierr < 0) CCTK_WARN(1, "Failed to register flat BC for ML_BSSN_MP::ML_trace_curv."); @@ -258,7 +261,7 @@ static void ML_BSSN_MP_boundary_Body(cGH const * restrict const cctkGH, int cons } /* 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); @@ -292,6 +295,8 @@ static void ML_BSSN_MP_boundary_Body(cGH const * restrict const cctkGH, int cons CCTK_REAL_VEC Xt3L = ToReal(0); + CCTK_REAL_VEC ThetaL = ToReal(0); + CCTK_REAL_VEC alphaL = ToReal(1); CCTK_REAL_VEC AL = ToReal(0); @@ -331,6 +336,7 @@ static void ML_BSSN_MP_boundary_Body(cGH const * restrict const cctkGH, int cons 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); @@ -364,8 +370,9 @@ extern "C" void ML_BSSN_MP_boundary(CCTK_ARGUMENTS) "ML_BSSN_MP::ML_log_confac", "ML_BSSN_MP::ML_metric", "ML_BSSN_MP::ML_shift", + "ML_BSSN_MP::ML_Theta", "ML_BSSN_MP::ML_trace_curv"}; - GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_MP_boundary", 9, groups); + GenericFD_AssertGroupStorage(cctkGH, "ML_BSSN_MP_boundary", 10, groups); switch(fdOrder) { |