Added new geo_toolz.

jejjohnson · Mar 20, 2024 · d139ba1 · d139ba1
1 parent cfeff88
commit d139ba1
Show file tree

Hide file tree

Showing 6 changed files with 13,161 additions and 0 deletions.
diff --git a/bayesevt/_src/geoprocessing/math.py b/bayesevt/_src/geoprocessing/math.py
@@ -0,0 +1,11 @@
+import numpy as np
+
+
+# SECOND AXIS
+def safe_reciprocal(x):
+    """Vectorized 1/x, treating x==0 manually"""
+    x = np.array(x, float)
+    near_zero = np.isclose(x, 0)
+    x[near_zero] = np.inf
+    x[~near_zero] = np.reciprocal(x[~near_zero])
+    return x
diff --git a/notebooks/evt/1.1_extremes_geotoolz.ipynb b/notebooks/evt/1.1_extremes_geotoolz.ipynb