Address review: plain-int count, loosen cross-algorithm tolerance (#3597)

Author: brendancol

xrspatial/geotiff/_symbology.py

@@ -268,7 +268,9 @@ def update(self, buf):
             vals = buf.ravel() if bool(mask.all()) else buf[mask]
         else:
             vals = buf.ravel()
-        n_b = vals.size
+        # Plain ``int`` so the running count (and the moment arithmetic
+        # built on it) stays in unbounded Python ints, not numpy scalars.
+        n_b = int(vals.size)
         if n_b == 0:
             return
         # Reductions with float64 accumulators on the original buffer:

xrspatial/geotiff/tests/write/test_color_ramp_single_pass_3597.py

@@ -203,14 +203,14 @@ def test_streaming_stats_matches_finite_stats_uneven_slabs():
     for r0, r1 in [(0, 5), (5, 6), (6, 40), (40, 64)]:
         acc.update(_BASE[r0:r1])
     da = xr.DataArray(_BASE, dims=("y", "x"))
-    assert acc.result() == pytest.approx(_finite_stats(da, None), rel=1e-12)
+    assert acc.result() == pytest.approx(_finite_stats(da, None), rel=1e-9)
 
 
 def test_streaming_stats_nan_nodata_treated_as_unset():
     acc = StreamingStats(nodata=float("nan"))
     acc.update(_BASE)
     da = xr.DataArray(_BASE, dims=("y", "x"))
-    assert acc.result() == pytest.approx(_finite_stats(da, None), rel=1e-12)
+    assert acc.result() == pytest.approx(_finite_stats(da, None), rel=1e-9)
 
 
 def test_streaming_stats_empty_returns_none():