diff options
Diffstat (limited to 'Carpet/CarpetInterp2/src/fasterp.hh')
-rw-r--r-- | Carpet/CarpetInterp2/src/fasterp.hh | 14 |
1 files changed, 10 insertions, 4 deletions
diff --git a/Carpet/CarpetInterp2/src/fasterp.hh b/Carpet/CarpetInterp2/src/fasterp.hh index fb76c0d24..cad64a9df 100644 --- a/Carpet/CarpetInterp2/src/fasterp.hh +++ b/Carpet/CarpetInterp2/src/fasterp.hh @@ -15,6 +15,12 @@ +// Define this at all times, because otherwise out-of-bounds +// interpolations may not be detected early enough +#define CARPETINTER2_CHECK + + + namespace CarpetInterp2 { using namespace std; @@ -138,7 +144,7 @@ namespace CarpetInterp2 { struct fasterp_iloc_t { mrc_t mrc; // map, refinement level, component -#ifdef CARPET_DEBUG +#ifdef CARPETINTER2_CHECK pn_t pn; // origin of this point ivect ipos; // closest grid point (Carpet indexing) ivect ind; // closest grid point (local indexing) @@ -166,7 +172,7 @@ namespace CarpetInterp2 { CCTK_REAL coeffs[dim][max_order+1]; // interpolation coefficients bvect exact; -#ifdef CARPET_DEBUG +#ifdef CARPETINTER2_CHECK public: pn_t pn; // origin of this point mrc_t mrc; // map, refinement level, component @@ -176,14 +182,14 @@ namespace CarpetInterp2 { #endif int ind3d; // source grid point offset -#ifdef CARPET_DEBUG +#ifdef CARPETINTER2_CHECK public: ivect saved_lsh; // copy of lsh private: #endif public: - void + int calc_stencil (fasterp_iloc_t const & iloc, ivect const & lsh, int order); |