diff options
author | Erik Schnetter <schnetter@cct.lsu.edu> | 2010-04-27 10:50:31 -0500 |
---|---|---|
committer | Barry Wardell <barry.wardell@gmail.com> | 2011-12-14 16:45:48 +0000 |
commit | 31aa0a30ac88e00db2e149c33e4e3c54237ce9d3 (patch) | |
tree | d9fd7d092c9677dd0c9afc1f455e4ea617eabfa6 /Carpet/CarpetInterp2/src | |
parent | 8352d38d2e56fb581b41e1997c536904a3594505 (diff) |
CarpetInterp2: Add debugging statements
Diffstat (limited to 'Carpet/CarpetInterp2/src')
-rw-r--r-- | Carpet/CarpetInterp2/src/fasterp.cc | 14 |
1 files changed, 14 insertions, 0 deletions
diff --git a/Carpet/CarpetInterp2/src/fasterp.cc b/Carpet/CarpetInterp2/src/fasterp.cc index 9bf990650..c59330556 100644 --- a/Carpet/CarpetInterp2/src/fasterp.cc +++ b/Carpet/CarpetInterp2/src/fasterp.cc @@ -643,6 +643,15 @@ namespace CarpetInterp2 { ostringstream msg; msg << "Interpolation point " << n << " on map " << m << " " << "at " << pos << " is outside of the grid hierarchy"; + msg << "\n" + << "rl=" << rl << " c=" << c << "\n" + << "rpos=" << rpos << "\n" + << "ipos=" << ipos << "\n" + << "lower=" << lower << "\n" + << "upper=" << upper << "\n" + << "delta=" << delta << "\n" + << "idelta=" << idelta << "\n" + << "hh=" << *hh << "\n"; CCTK_WARN (CCTK_WARN_ABORT, msg.str().c_str()); } } @@ -656,6 +665,11 @@ namespace CarpetInterp2 { assert (all (ipos % ext.stride() == ivect(0))); ipos /= ext.stride(); dpos /= rvect(ext.stride()); + if (not (all (abs(dpos) <= rvect(0.5)))) { + cout << "fasterp.cc:659\n" + << " dpos=" << dpos << "\n" + << " ext=" << ext << "\n"; + } assert (all (abs(dpos) <= rvect(0.5))); ivect const ind = ipos - ext.lower() / ext.stride(); |