Merge pull request #120 from lanl/jmm/fix-out-of-bounds

Jmm/fix out of bounds

Author: Yurlungur

spiner/piecewise_grid_1d.hpp

@@ -56,12 +56,14 @@ class PiecewiseGrid1D {
       grids_[i] = grids[i];
       pointTotals_[i] = point_tot;
       point_tot += grids[i].nPoints();
-      const Real diff = std::abs(grids_[i].min() - grids_[i - 1].max());
-      const Real avg =
-          0.5 * (std::abs(grids_[i].min()) + std::abs(grids_[i - 1].max()));
-      if ((i > 0) && ratio_(diff, avg) >= EPS_()) {
-        PORTABLE_ALWAYS_THROW_OR_ABORT(
-            "Grids must be ordered and intersect at exactly one point.");
+      if (i > 0) {
+        const Real diff = std::abs(grids_[i].min() - grids_[i - 1].max());
+        const Real avg =
+            0.5 * (std::abs(grids_[i].min()) + std::abs(grids_[i - 1].max()));
+        if (ratio_(diff, avg) >= EPS_()) {
+          PORTABLE_ALWAYS_THROW_OR_ABORT(
+              "Grids must be ordered and intersect at exactly one point.");
+        }
       }
     }
   }