diff --git a/.claude/sweep-test-coverage-state.csv b/.claude/sweep-test-coverage-state.csv index 041911701..f8d3d3da1 100644 --- a/.claude/sweep-test-coverage-state.csv +++ b/.claude/sweep-test-coverage-state.csv @@ -5,6 +5,7 @@ contour,2026-06-08,2704;2710;3044,MEDIUM,1;2,"Pass 2 (2026-06-08, deep-sweep tes cost_distance,2026-06-16,3367,MEDIUM,1;2,"Pass (2026-06-16 deep-sweep test-coverage, CUDA host). cost_distance is heavily tested: 1122 test lines for 1354 src lines, all 4 backends parametrized + regression tests for #1191/#880/#1252/#1262/#3340/#3341/#3343/#3344. Found one MEDIUM Cat 1+2 gap: _cost_distance_dask f_min<=0 early return (all-impassable friction, finite max_cost -> da.full NaN preserving chunks) was unreached -- numpy equiv covered by test_source_on_impassable_cell, iterative dask by test_iterative_narrow_corridor, but the bounded map_overlap wrapper shortcut was not. Filed #3367, added test_dask_all_impassable_friction_returns_nan (all-zero friction, dask+numpy chunks(3,3), max_cost=5; asserts all-NaN, dask-backed, npartitions>1). RAN + PASSED; -W error::UserWarning confirms early return taken (no iterative warning). Full file 85 passed on CUDA host. LOW (documented, not fixed): non-square cellsize numeric correctness untested (_make_meta_raster uses res=(2,3) but test_metadata_preserved checks metadata only)." dasymetric,2026-06-20,3407;3406,HIGH,2;3;4;5,"deep-sweep test-coverage on a CUDA host (CUDA available, GPU tests ran). Module is well covered (813 test loc / 834 src): 4-backend equivalence for disaggregate weighted+binary, conservation, NaN/nodata/negative-weight, limiting_variable + cupy/dask NotImplemented guards, pycnophylactic numpy+cupy+dask-raises, validate_disaggregation all backends, memory guards (#1261). Filed #3407 (test-only) for real gaps and added 4 new test classes (11 passed, 2 xfailed). Cat5 HIGH: metadata (attrs res/crs + coords) never asserted -> TestMetadataPreservation (numpy/dask). Cat3 HIGH: true 1x1 raster untested (only 1x2 strip) -> TestSinglePixel for disaggregate weighted/binary + pycnophylactic (degenerate no-shift smoothing) + dask parity. Cat2 MEDIUM: Inf weight collapses zone total to 0 (silent conservation break) -> TestInfWeight pins current behaviour. Cat4 MEDIUM: 3-class limiting_variable (multi-break + per-class caps) untested despite docstring -> TestLimitingVariableThreeClass. SOURCE BUG found (filed #3406, NOT fixed - test-only sweep): pycnophylactic raises ValueError (np.nanmax on zero-size array) when no pixel is valid for smoothing (all-NaN zones or no zone id in values); disaggregate handles same input gracefully (all-NaN). Pinned with TestPycnophylacticEmptyValid xfail(strict, raises=ValueError) -> flips red when #3406 fixed. LOW (documented, not fixed): non-square cellsize never exercised (all tests use res 0.5/0.5); disaggregate cupy/dask+cupy 1x1 + metadata not separately added (eager numpy gap was the real one, GPU dispatch already covered by TestCrossBackend)." diffusion,2026-06-20,3422,HIGH,1;2;3;4,"Pass 1 (2026-06-20, deep-sweep test-coverage, CUDA host). diffuse() dispatch table registers all 4 backends but test_diffusion.py only exercised numpy + dask+numpy. Cat 1 HIGH: cupy (_diffuse_cupy/_diffuse_step_gpu) and dask+cupy (_diffuse_dask_cupy/_diffuse_chunk_cupy) registered but never invoked -- no test ran them. Cat 4 HIGH: boundary accepts nan/nearest/reflect/wrap; only nearest+wrap tested, reflect had none. Cat 3 HIGH: 1x1 single-pixel and Nx1/1xN strip rasters never tested. Cat 2 MEDIUM: NaN tested numpy-only; Inf and all-NaN inputs untested. Filed #3422, added 14 tests (PR #3424, test-only, source untouched): cupy/dask+cupy parity vs numpy (incl. spatially-varying alpha + NaN propagation), reflect boundary across all 4 backends, 1x1 + Nx1 + 1xN (numpy + chunked dask strip), all-NaN stays NaN, Inf contamination smoke test. All 14 RAN+PASSED on a CUDA host; the 4 cupy/dask+cupy tests genuinely executed (not skipped); full file 39 passed. All paths verified correct before the tests were added -- coverage gap, not a bug. LOW (documented, not fixed): non-square cellsize (res[0]!=res[1]) never exercised -- diffuse uses res[0] as dx and assumes square cells; empty 0-row/0-col raster untested; asv benchmark absent; 'nan' boundary-mode edge=NaN behaviour not directly asserted on diffuse (covered indirectly via wrap/nearest)." +fire,2026-06-19,,HIGH,1;5,"Cat1 HIGH: dask+cupy dispatch registered but untested for rdnbr/burn_severity_class/fireline_intensity/flame_length/rate_of_spread/kbdi (only dnbr had it). Cat5 MEDIUM: only dnbr asserted attrs/coords/dims preservation. Added 6 test_numpy_equals_dask_cupy + 6 test_output_preserves_metadata; all 66 fire tests pass on a CUDA host. Cat2 LOW (not fixed): Inf inputs untested (pure per-cell math, low risk)." focal,2026-06-10,3220;3219;3225,HIGH,1;2;3;4,"Deep-sweep 2026-06-10 on CUDA host, all 4 backends executed. Filed #3220 (coverage) and added 36 tests in PR branch: Inf inputs for mean/focal_stats (HIGH Cat2 - no Inf test existed anywhere), mean NaN input (HIGH Cat2 - default excludes=[nan] semantics never asserted), 1x1 + 1xN/Nx1 strips (HIGH Cat3), empty 0-row raster numpy-only (MEDIUM Cat3), mean passes=2 == mean(mean) and excludes sentinel -9999 behavioral tests (MEDIUM Cat4), dask+cupy non-default boundary modes for mean/apply/focal_stats (MEDIUM Cat1/4). Bugs surfaced, filed separately (NOT fixed here): #3219 hotspots silently returns all zeros on Inf input (nan global std passes the std==0 guard, all 4 backends); #3225 empty raster works on numpy but crashes cupy (raw CudaAPIError) and dask (map_overlap depth ValueError). hotspots+Inf and non-numpy empty behavior left unpinned until those are fixed. Backend matrix for the 4 public funcs was already solid (all 4 backends + parity); boundary modes covered except dask+cupy. Siblings filed #3214-3217 same day (dtype/docstring/apply-default-func) - no overlap." geotiff,2026-06-12,3266,MEDIUM,1,"Pass 22 (2026-06-12, deep-sweep test-coverage): delta audit of the ~20 commits since pass 21 (06-09..06-12, mostly pack/unpack fixes + #3241 GPU streaming writer + coregister #3254/#3248). Filed #3266 (tests). Cat 1 MEDIUM: pack=True gained working gpu/dask+gpu support in #3240, but three pack features were tested numpy+dask only: float32 width preservation (#3080, test_pack_float_width_3080.py), nodata kwarg fill (#3168, test_pack_nodata_kwarg_3168.py), band-subset per-band SCALE/OFFSET rewrite (#3161, test_pack_band_subset_3161.py). Live probe on this CUDA host: all six gpu/dask+gpu legs pass today (no source bug, pure coverage gap). Added one gpu/dask-gpu parametrized round-trip test per file (6 tests, requires_gpu, RUN+passing locally) and fixed two stale docstrings claiming unpack/pack is CPU-only (wrong since #3075/#3240). Verified NOT gaps this pass: #3128 int64 sentinel tests cover eager+dask+gpu; #3241 streaming writer landed with byte-identical band-first/band-last/BytesIO/small-buffer tests; #3104 scale-zero rejection has gpu legs; #3169 revived the dead compression-corpus oracle gate. Out of scope: coregister=True lives in accessor.py (excluded module); its multi-band + polar gaps are documented as experimental caveats in docs/source/reference/geotiff.rst (#3248). LOW (carried, documented not fixed): Inf as the declared nodata sentinel never tested. || PREVIOUS: Pass 21 (2026-06-09, deep-sweep test-coverage): filed #3114 (tests) + #3112 (source bug). Cat 1 HIGH: to_geotiff(pack=True) round-trip was tested only on numpy and dask+numpy (write/test_pack_3064.py); #3075 made unpack=True work on gpu and dask+gpu reads, but no test packed a GPU-read array back. Live probe on this CUDA host: BOTH GPU legs crash today -- eager gpu raises AttributeError (cupy has no astype, the known cupy 13.6/xarray 2025.12 where/astype incompat) and dask+gpu raises TypeError (numpy fill value inside cupy.where) -- both from _pack's out.fillna(nodata) in _attrs.py; _writers/gpu.py says the pre-dispatch re-pack is supposed to make every write path work. Source bug filed as #3112; test-only PR adds test_pack_round_trip_gpu (gpu + dask-gpu params, requires_gpu, xfail(strict=True) on #3112 so the fix flips them loudly) and fixes the stale module docstring claiming GPU rejects mask_and_scale. Ran on CUDA host: 13 passed, 2 xfailed. Verified NOT gaps this pass (probed before flagging): empty/zero-band writer guard is covered (test_basic.py 2075/2095 blocks incl. gpu + dask + streaming entry points); degenerate shapes covered on all 4 read backends (read/test_degenerate_shapes.py); overview_resampling all 7 modes parametrized; missing_sources raise/warn + invalid, band_nodata first/invalid, unpack on all 4 backends, masked/parse_coordinates/lock/cache/default_name/name-deprecation all exercised; attrs contract per-backend (attrs/test_contract.py). LOW (documented, not fixed): Inf as the declared nodata sentinel is never tested (only one nan+inf data round-trip in test_edge_cases.py). || PREVIOUS: Pass 20 (2026-06-06, deep-sweep test-coverage): filed #2984 and added test_writer.py degenerate-shape GPU write coverage (Cat 1 backend + Cat 3 geometric edge). Read side already covers 1x1/1xN/Nx1 on all 4 backends (read/test_degenerate_shapes.py) and the dask streaming writer covers them (integration/test_dask_pipeline.py); the GPU write path was the gap (smallest shape in gpu/test_writer.py was 2x2). Added test_write_geotiff_gpu_degenerate_round_trip (1x1/1xN/Nx1 x none/deflate) + test_to_geotiff_dask_gpu_degenerate_round_trip (dask+cupy via gpu=True). 9 new tests RUN+passing on a CUDA host. Verified paths work first (not a source bug); transform supplied explicitly via attrs. Wider tree audit (~92k test LOC vs ~33k source): rioxarray-compat (#2961), bbox NaN/Inf/rotated, 8-backend parity matrix, codec round-trips already covered -- no other real gaps. | Pass (2026-06-05 test-coverage sweep): mature module (~31k src / ~124k test LOC, 9 test dirs). Exhaustive existing coverage -- parity/test_backend_matrix.py runs all 4 backends + VRT + HTTP + fsspec; golden_corpus full-manifest parity; read_rioxarray_compat_2961 covers masked/mask_and_scale/parse_coordinates/default_name on eager+dask. Cat1+Cat3 gap found (MEDIUM): degenerate-shape READS (1x1/1xN/Nx1) were tested only on the eager numpy reader (test_edge_cases.py) and the dask streaming WRITE path (integration/test_dask_pipeline.py); the windowed dask READ (chunks=) and GPU READ (gpu=True) on a single-pixel dimension were never exercised (smallest dask-read source in read/test_tiling is 8x8/2x32, parity fixtures 32x32/64x64). Probed: paths work today, no source bug -- pure coverage gap. Added read/test_degenerate_shapes.py (18 tests): dask read x{chunks 1,3,4} x{1x1,1xN,Nx1} + coord/transform/crs parity + GPU read + dask+gpu read. GPU cells RAN and PASSED on this CUDA host (grid-size-1 launch validated). Fixture supplies explicit attrs['transform'] (writer cannot infer pixel size from a 1-element coord axis). Branch deep-sweep-test-coverage-geotiff-degenerate-read-01. NOTE: pre-existing union-merge CRLF/duplicate-record corruption in this CSV left untouched -- appended one clean record; DictReader last-write-wins picks this one." idw,2026-06-04,2919,HIGH,1;4,"cupy/dask+cupy backends untested (Cat1 HIGH); GPU k-reject error path untested (Cat4 MED). Added 6 GPU tests, validated on CUDA host. Inf-in-points (Cat2) and attrs-preservation (Cat5) are LOW, documented not fixed." diff --git a/xrspatial/tests/test_fire.py b/xrspatial/tests/test_fire.py index 9e3f2025b..383627b10 100644 --- a/xrspatial/tests/test_fire.py +++ b/xrspatial/tests/test_fire.py @@ -209,6 +209,26 @@ def test_numpy_equals_cupy(self): np.testing.assert_allclose(cp_r.data.get(), np_r.data, rtol=1e-5, equal_nan=True) + @dask_array_available + @cuda_and_cupy_available + def test_numpy_equals_dask_cupy(self): + rng = np.random.default_rng(42) + d = rng.uniform(-1, 1, (6, 8)).astype('f4') + p = rng.uniform(0.01, 1.0, (6, 8)).astype('f4') * 1000 + np_r = rdnbr(create_test_raster(d, 'numpy'), + create_test_raster(p, 'numpy')) + dc_r = rdnbr(create_test_raster(d, 'dask+cupy'), + create_test_raster(p, 'dask+cupy')) + np.testing.assert_allclose(dc_r.data.compute().get(), np_r.data, + rtol=1e-5, equal_nan=True) + + def test_output_preserves_metadata(self): + d = np.array([[0.4, 0.3], [0.5, 0.1]], dtype=np.float32) + p = np.array([[500.0, 200.0], [800.0, 100.0]], dtype=np.float32) + d_agg = create_test_raster(d) + result = rdnbr(d_agg, create_test_raster(p)) + general_output_checks(d_agg, result, verify_dtype=False) + # --------------------------------------------------------------------------- # Burn severity class @@ -269,6 +289,21 @@ def test_numpy_equals_cupy(self): cp_r = burn_severity_class(create_test_raster(data, 'cupy')) np.testing.assert_array_equal(cp_r.data.get(), np_r.data) + @dask_array_available + @cuda_and_cupy_available + def test_numpy_equals_dask_cupy(self): + rng = np.random.default_rng(99) + data = rng.uniform(-1, 1, (8, 10)).astype('f4') + np_r = burn_severity_class(create_test_raster(data, 'numpy')) + dc_r = burn_severity_class(create_test_raster(data, 'dask+cupy')) + np.testing.assert_array_equal(dc_r.data.compute().get(), np_r.data) + + def test_output_preserves_metadata(self): + vals = np.array([[-0.5, 0.3], [0.7, 0.1]], dtype=np.float32) + agg = create_test_raster(vals) + result = burn_severity_class(agg) + general_output_checks(agg, result, verify_dtype=False) + # --------------------------------------------------------------------------- # Fireline intensity @@ -326,6 +361,26 @@ def test_numpy_equals_cupy(self): np.testing.assert_allclose(cp_r.data.get(), np_r.data, rtol=1e-5, equal_nan=True) + @dask_array_available + @cuda_and_cupy_available + def test_numpy_equals_dask_cupy(self): + rng = np.random.default_rng(7) + f = rng.uniform(0, 5, (6, 8)).astype('f4') + s = rng.uniform(0, 1, (6, 8)).astype('f4') + np_r = fireline_intensity(create_test_raster(f, 'numpy'), + create_test_raster(s, 'numpy')) + dc_r = fireline_intensity(create_test_raster(f, 'dask+cupy'), + create_test_raster(s, 'dask+cupy')) + np.testing.assert_allclose(dc_r.data.compute().get(), np_r.data, + rtol=1e-5, equal_nan=True) + + def test_output_preserves_metadata(self): + f = np.array([[2.0, 0.5], [1.0, 0.3]], dtype=np.float32) + s = np.array([[0.1, 0.2], [0.3, 0.4]], dtype=np.float32) + f_agg = create_test_raster(f) + result = fireline_intensity(f_agg, create_test_raster(s)) + general_output_checks(f_agg, result, verify_dtype=False) + # --------------------------------------------------------------------------- # Flame length @@ -382,6 +437,19 @@ def test_numpy_equals_cupy(self, intensity_data): np.testing.assert_allclose(cp_r.data.get(), np_r.data, rtol=1e-5, equal_nan=True) + @dask_array_available + @cuda_and_cupy_available + def test_numpy_equals_dask_cupy(self, intensity_data): + np_r = flame_length(create_test_raster(intensity_data, 'numpy')) + dc_r = flame_length(create_test_raster(intensity_data, 'dask+cupy')) + np.testing.assert_allclose(dc_r.data.compute().get(), np_r.data, + rtol=1e-5, equal_nan=True) + + def test_output_preserves_metadata(self, intensity_data): + agg = create_test_raster(intensity_data) + result = flame_length(agg) + general_output_checks(agg, result, verify_dtype=False) + # --------------------------------------------------------------------------- # Rate of spread (Rothermel) @@ -514,6 +582,31 @@ def test_numpy_equals_cupy(self, ros_inputs): np.testing.assert_allclose(cp_r.data.get(), np_r.data, rtol=1e-4, equal_nan=True) + @dask_array_available + @cuda_and_cupy_available + def test_numpy_equals_dask_cupy(self, ros_inputs): + slope, wind, moisture = ros_inputs + np_r = rate_of_spread( + create_test_raster(slope, 'numpy'), + create_test_raster(wind, 'numpy'), + create_test_raster(moisture, 'numpy'), + ) + dc_r = rate_of_spread( + create_test_raster(slope, 'dask+cupy'), + create_test_raster(wind, 'dask+cupy'), + create_test_raster(moisture, 'dask+cupy'), + ) + np.testing.assert_allclose(dc_r.data.compute().get(), np_r.data, + rtol=1e-4, equal_nan=True) + + def test_output_preserves_metadata(self, ros_inputs): + slope, wind, moisture = ros_inputs + slope_agg = create_test_raster(slope) + result = rate_of_spread(slope_agg, + create_test_raster(wind), + create_test_raster(moisture)) + general_output_checks(slope_agg, result, verify_dtype=False) + # --------------------------------------------------------------------------- # KBDI @@ -629,6 +722,35 @@ def test_numpy_equals_cupy(self): np.testing.assert_allclose(cp_r.data.get(), np_r.data, rtol=1e-5, equal_nan=True) + @dask_array_available + @cuda_and_cupy_available + def test_numpy_equals_dask_cupy(self): + rng = np.random.default_rng(123) + prev = rng.uniform(0, 400, (6, 8)).astype('f4') + temp = rng.uniform(15, 40, (6, 8)).astype('f4') + precip = rng.uniform(0, 20, (6, 8)).astype('f4') + np_r = kbdi(create_test_raster(prev, 'numpy'), + create_test_raster(temp, 'numpy'), + create_test_raster(precip, 'numpy'), + annual_precip=1200.0) + dc_r = kbdi(create_test_raster(prev, 'dask+cupy'), + create_test_raster(temp, 'dask+cupy'), + create_test_raster(precip, 'dask+cupy'), + annual_precip=1200.0) + np.testing.assert_allclose(dc_r.data.compute().get(), np_r.data, + rtol=1e-5, equal_nan=True) + + def test_output_preserves_metadata(self): + prev = np.array([[100.0, 200.0], [300.0, 50.0]], dtype=np.float32) + temp = np.array([[30.0, 25.0], [35.0, 20.0]], dtype=np.float32) + precip = np.array([[0.0, 5.0], [2.0, 10.0]], dtype=np.float32) + prev_agg = create_test_raster(prev) + result = kbdi(prev_agg, + create_test_raster(temp), + create_test_raster(precip), + annual_precip=1500.0) + general_output_checks(prev_agg, result, verify_dtype=False) + # --------------------------------------------------------------------------- # Accessor smoke tests