CarpetInterp: bugfix in routine which maps points onto components

darcs-hash:20050829151032-776a0-95a8c84bd10d434d1dee4c47894ab2f4276836ae.gz

1 files changed, 3 insertions, 4 deletions

diff --git a/Carpet/CarpetInterp/src/interp.cc b/Carpet/CarpetInterp/src/interp.cc
index 7b8646dee..a078f6a5b 100644
--- a/Carpet/CarpetInterp/src/interp.cc
+++ b/Carpet/CarpetInterp/src/interp.cc
@@ -758,14 +758,13 @@ namespace CarpetInterp {
       // Point could not be mapped onto any processor's domain.
       // Map the point (arbitrarily) to the first component of the
       // coarsest grid
+      rl = minrl;
+      c = 0;
+      // only warn once
       if (map_onto_processors) {
         CCTK_VWarn (CCTK_WARN_PICKY, __LINE__, __FILE__, CCTK_THORNSTRING,
                     "Interpolation point #%d at [%g,%g,%g] is not on "
                     "any grid patch", n, pos[0], pos[1], pos[2]);
-        rl = minrl;
-        c = 0;
-      } else {
-        assert (0 and "point is not mapped");
       }
     found:
       assert (rl >= minrl and rl < maxrl);