restore OutOfBounds check. assert if neighbors array not fully filled.

acts-project · Jun 14, 2023 · f596fd2 · f596fd2
1 parent e8ec0ca
commit f596fd2
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/Core/include/Acts/MagneticField/InterpolatedBFieldMap.hpp b/Core/include/Acts/MagneticField/InterpolatedBFieldMap.hpp
@@ -202,11 +202,17 @@ class InterpolatedBFieldMap : public InterpolatedMagneticField {
     std::array<Vector3, nCorners> neighbors;
     const auto& cornerIndices = m_cfg.grid.closestPointsIndices(gridPosition);
 
+    if (!isInsideLocal(gridPosition)) {
+      return MagneticFieldError::OutOfBounds;
+    }
+
     size_t i = 0;
     for (size_t index : cornerIndices) {
       neighbors.at(i++) = m_cfg.transformBField(m_cfg.grid.at(index), position);
     }
 
+    assert(i == nCorners);
+
     return FieldCell(lowerLeft, upperRight, std::move(neighbors));
   }