[ET Trac] [Einstein Toolkit] #181: AEILocalInterp should not off-centre the interpolation stencil by default
Einstein Toolkit
trac-noreply at einsteintoolkit.org
Mon Oct 24 08:48:22 CDT 2016
#181: AEILocalInterp should not off-centre the interpolation stencil by default
---------------------+------------------------------------------------------
Reporter: hinder | Owner:
Type: defect | Status: confirmed
Priority: major | Milestone:
Component: Cactus | Version:
Resolution: | Keywords: AEILocalInterp
---------------------+------------------------------------------------------
Comment (by knarf):
Replying to [comment:9 eschnett]:
> Let's assume you have 10 points on your grid, and your domain is [0; 1).
Your points' coordinates are thus [0.0, ..., 0.9]. If you try to
interpolate at 0.95, then this is fine, as it's inside the domain.
However, all AEILocalInterp sees is that it has grid points in the range
[0.0; 0.9], and it will thus complain.
Why should it complain? With 3 symmetry points, the actual grid should go
from -0.3 ... 1.2 (1.3?), and with an interpolation stencil of 2 (0.2) in
each direction, shouldn't 0.95 be fine? Why should AEILocalInterp care
about the size of the evolved grid, and not only on the size of the grid
it has to interpolate on?
--
Ticket URL: <https://trac.einsteintoolkit.org/ticket/181#comment:10>
Einstein Toolkit <http://einsteintoolkit.org>
The Einstein Toolkit
More information about the Trac
mailing list