rastr 0.7.0__tar.gz → 0.7.1__tar.gz

{rastr-0.7.0 → rastr-0.7.1}/PKG-INFO

@@ -1,11 +1,11 @@
 Metadata-Version: 2.4
 Name: rastr
-Version: 0.7.0
+Version: 0.7.1
 Summary: Geospatial Raster datatype library for Python.
 Project-URL: Source Code, https://github.com/tonkintaylor/rastr
 Project-URL: Bug Tracker, https://github.com/tonkintaylor/rastr/issues
 Project-URL: Releases, https://github.com/tonkintaylor/rastr/releases
-Project-URL: Source Archive, https://github.com/tonkintaylor/rastr/archive/2b485cc676121c82f468dca7733e444c3033abbe.zip
+Project-URL: Source Archive, https://github.com/tonkintaylor/rastr/archive/f84ab299005d87e37339d53759e005e85c9e5ad2.zip
 Author-email: Tonkin & Taylor Limited <Sub-DisciplineData+AnalyticsStaff@tonkintaylor.co.nz>, Nathan McDougall <nmcdougall@tonkintaylor.co.nz>, Ben Karl <bkarl@tonkintaylor.co.nz>
 License-Expression: MIT
 License-File: LICENSE

{rastr-0.7.0 → rastr-0.7.1}/src/rastr/_version.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '0.7.0'
-__version_tuple__ = version_tuple = (0, 7, 0)
+__version__ = version = '0.7.1'
+__version_tuple__ = version_tuple = (0, 7, 1)
 
 __commit_id__ = commit_id = None

{rastr-0.7.0 → rastr-0.7.1}/src/rastr/gis/smooth.py

@@ -132,10 +132,14 @@ def _recursive_eval(slice4: NDArray, tangents: NDArray, ts: NDArray) -> NDArray:
     for r in range(1, 4):
         idx = max(r - 2, 0)
         denom = tangents[r - idx : 4 - idx] - tangents[idx : 4 - r + idx]
-        denom = np.where(denom == 0, np.finfo(float).eps, denom)  # avoid div 0
 
         # Compute weights for all parameter values at once
-        left_w = (tangents[r - idx : 4 - idx][None, :] - ts[:, None]) / denom
+        with np.errstate(divide="ignore", invalid="ignore"):
+            left_w = (tangents[r - idx : 4 - idx][None, :] - ts[:, None]) / denom
+
+        left_w[:, denom == 0] = (
+            0.5  # Use 0.5 (midpoint) when points are identical (zero denominator)
+        )
         right_w = 1 - left_w
 
         # Weighted sums between consecutive points

{rastr-0.7.0 → rastr-0.7.1}/tests/rastr/gis/test_smooth.py

@@ -174,3 +174,18 @@ class TestRecursiveEval:
         # Assert
         assert result.shape == (1, 2)
         assert np.allclose(result, np.array([[0.0, 0.0]]))
+
+    def test_denom_zero_non_edge_point(self):
+        # To address regression in https://github.com/tonkintaylor/rastr/issues/276
+
+        # Arrange
+        coords = np.array([[0.0, 0.0], [1.0, 1.0], [1.0, 1.0], [2.0, 0.0]])
+        # This will force denom == 0 in _recursive_eval for the middle segment
+
+        # Act
+        result = _recursive_eval(
+            coords, np.asarray([0.0, 1.0, 1.0, 2.0]), np.array([0.5])
+        )
+        # Assert
+        assert result.shape == (1, 2)
+        assert np.allclose(result, np.array([[1.0, 0.75]]))