NewRad will extrapolate into any boundary region where at least one face is marked as a physical output boundary.

This makes it impossible to use with Cartoon2D (and Cartoon2D’s tests fail if one tries) since Cartoon2D runs have only a single interior point that could be used as a source for extrapolation.

Pull request:

https://bitbucket.org/einsteintoolkit/einsteinevolve/pull-requests/14/newrad-do-not-extrapolate-into-symmetry

changes this so that NewRad does not interpolate into any region where at least one face is a symmetry boundary.

This assumes:

• symmetry boundaries are applied afterwards (likely since NewRad is triggerd by user code only in CCTK_INITIAL and not a boundary condition)
• the symmetry boundary conditions can fill in edges and corners where physical and symmetry boundaries overlap (true for the ones in the ET)

‌

--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2449/newrad-cannot-be-used-with-cartoon2d-due