#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):

The function to use is actually GetBoundarySpecification not GetBoundarySizesAndTypes. Easiest is most likely to copy the code from ML_ADMConstraint’s Kranc.cc and ML_ADMConstraints_evaluate.cc files, ie:

/*********************************************************************
 * GetBoundaryWidths
 *********************************************************************/

void GetBoundaryWidths(cGH const * restrict const cctkGH, CCTK_INT nboundaryzones[6])
{
  CCTK_INT is_internal[6];
  CCTK_INT is_staggered[6];
  CCTK_INT shiftout[6];
  int ierr = -1;

  if (CCTK_IsFunctionAliased ("MultiPatch_GetBoundarySpecification")) {
    int const map = MultiPatch_GetMap (cctkGH--

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