Added inf check in VolumeCalculation. (#2634)

Co-authored-by: Rosie Barker <[email protected]>
Co-authored-by: Paul Romano <[email protected]>

openmc/volume.py

@@ -122,6 +122,10 @@ def __init__(self, domains, samples, lower_left=None, upper_right=None):
                 raise ValueError('Could not automatically determine bounding box '
                                  'for stochastic volume calculation.')
 
+        if np.isinf(self.lower_left).any() or np.isinf(self.upper_right).any():
+            raise ValueError('Lower-left and upper-right bounding box '
+                             'coordinates must be finite.')
+
     @property
     def ids(self):
         return self._ids

tests/unit_tests/test_volume.py

@@ -0,0 +1,15 @@
+import numpy as np
+import pytest
+
+import openmc
+
+
+def test_infinity_handling():
+    surf1 = openmc.Sphere(boundary_type="vacuum")
+    cell1 = openmc.Cell(region=-surf1)
+
+    lower_left = (-2, -np.inf, -2)
+    upper_right = (np.inf, 2, 2)
+
+    with pytest.raises(ValueError, match="must be finite"):
+        openmc.VolumeCalculation([cell1], 100, lower_left, upper_right)