adding check to prevent negative numbers

PMEAL · Jan 26, 2024 · 51c3715 · 51c3715
1 parent 93b5bb5
commit 51c3715
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/openpnm/models/geometry/pore_surface_area/_funcs.py b/openpnm/models/geometry/pore_surface_area/_funcs.py
@@ -42,6 +42,7 @@ def sphere(
     Tca = network[throat_cross_sectional_area]
     _np.subtract.at(value, network.conns[:, 0], Tca)
     _np.subtract.at(value, network.conns[:, 1], Tca)
+    value = _np.clip(value, 1e-12, _np.inf)
     return value
 
 
@@ -72,6 +73,7 @@ def circle(
     Tca = network[throat_cross_sectional_area]
     _np.subtract.at(value, network.conns[:, 0], Tca)
     _np.subtract.at(value, network.conns[:, 1], Tca)
+    value = _np.clip(value, 1e-12, _np.inf)
     return value
 
 
@@ -99,6 +101,7 @@ def cube(
     Tca = network[throat_cross_sectional_area]
     _np.subtract.at(value, network.conns[:, 0], Tca)
     _np.subtract.at(value, network.conns[:, 1], Tca)
+    value = _np.clip(value, 1e-12, _np.inf)
     return value
 
 
@@ -126,4 +129,5 @@ def square(
     Tca = network[throat_cross_sectional_area]
     _np.subtract.at(value, network.conns[:, 0], Tca)
     _np.subtract.at(value, network.conns[:, 1], Tca)
+    value = _np.clip(value, 1e-12, _np.inf)
     return value