#2288: Lean's boundary condition setting for BSSN constraints only fills in a single ghostzone
Reporter:Zach Etienne
Status:new
Milestone:
Version:
Type:bug
Priority:major
Component:EinsteinToolkit thorn

Comment (by Roland Haas):

Since Baikal does not support Llama (it does not, does it?) all you should have to do is use GetBoundarySpecification instead of GetBoundarySizesAndTypes. Lean_public which does, or at least there were plans to do so, support Llama has to do more work.

Ie this should work (and does for me using repos/wvuthorns/BaikalVacuum/test/BaikalVacuum_EE_O8_sgw3d.par even when increasing resolution to dx=0.0.1 so that not all MPI ranks have a boundary face all the time) :

diff --git a/BaikalVacuum/src/BoundaryConditions.c b/BaikalVacuum/src/BoundaryConditions.c
index 7240e5b..f89a9d2 100644
--- a/BaikalVacuum/src/BoundaryConditions.c
+++ b/BaikalVacuum/src/BoundaryConditions.c
@@ -96,11 +96,11 @@ void BaikalVacuum_BoundaryConditions_aux_gfs(CCTK_ARGUMENTS) {
   CCTK_INT ierr CCTK_ATTRIBUTE_UNUSED = 0;

   CCTK_INT bndsize[6];
-  CCTK_INT is_ghostbnd[6];
-  CCTK_INT is_symbnd[6];
-  CCTK_INT is_physbnd[6];
-  GetBoundarySizesAndTypes(cctkGH, 6, bndsize, is_ghostbnd, is_symbnd, is_physbnd);
+  CCTK_INT is_internal[6];
+  CCTK_INT is_staggered[6];
+  CCTK_INT shiftout[6];

+  GetBoundarySpecification(6, bndsize, is_internal, is_staggered, shiftout);

   ierr = Boundary_SelectVarForBC(cctkGH, CCTK_ALL_FACES, bndsize[0], -1, "BaikalVacuum::HGF", "flat");
   if (ierr < 0) CCTK_ERROR("Failed to register BC for BaikalVacuum::HGF!");

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2288/leans-boundary-condition-setting-for-bssn