#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):
Though looking at the code that is used for the MultiPatch_GetBoundarySpecification case, this seems incorrect as well. If called in LEVEL mode it will set all boundary widths to 0. And when called in SINGLEMAP mode (so that map is defined) you would end up selecting boundary conditions more than once which makes Boundary abort.
Something similar to what is done when going from GetBoundaryWidths to GetBoundaryWidth seems required, ie.:
if (CCTK_IsFunctionAliased ("MultiPatch_GetBoundarySpecification")) { int const maps = MultiPatch_GetMaps (cctkGH); for (int i = 0; i < 6; i++) nboundaryzones[i] = -1; // collect all boundary information from all maps, and get // outer boundary information only for (int m = 0; m < maps--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2288/leans-boundary-condition-setting-for-bssn