diff options
author | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2006-08-13 19:54:28 +0000 |
---|---|---|
committer | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2006-08-13 19:54:28 +0000 |
commit | 8bf581ab72264463f2da5136e5385011d431ffe2 (patch) | |
tree | cf375fb7bf2d842f0ca68c87fce9b22dd4d39465 /src/template.c | |
parent | 57f2e0c8839578e57770cac1beca01880e8ae607 (diff) |
fix a nasty bug where, if the interpolation coordinate was close to
a grid boundary, the actual interpolation was sometimes done as if
the interpolation coordinate was displaced by one grid spacing :( :(
also add more debugging code (introduced to help track down this bug,
but there's no harm in leaving it in for possible future use)
also add explicit checks for a couple of error cases:
* interpolation coordinate is not a finite number (eg NaN or infinity)
* grid origin/delta not a finite number (eg NaN or infinity)
* grid spacing is 0 (we need to divide by this!)
also add explicit note about the test driver in this directory
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/AEILocalInterp/trunk@35 0f49ee68-0e4f-0410-9b9c-b2c123ded7ef
Diffstat (limited to 'src/template.c')
-rw-r--r-- | src/template.c | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/src/template.c b/src/template.c index c1ad45f..09ed7b3 100644 --- a/src/template.c +++ b/src/template.c @@ -974,6 +974,8 @@ if (debug > 0) case MOLECULE_POSN_ERROR_X_GT_MAX: this_axis_status = CCTK_ERROR_INTERP_POINT_OUTSIDE; break; + case MOLECULE_POSN_ERROR_NAN: + this_axis_status = CCTK_ERROR_INTERP_COORD_NAN; default: CCTK_VWarn(BUG_MSG_SEVERITY_LEVEL, __LINE__, __FILE__, CCTK_THORNSTRING, |