Adjust tolerances to accommodate results on NumPy 2 as well as 1:

  Tiny numbers of pixels differing by ~0.1% after the extraction step (depending on the NumPy 2.x version and whether my NDData #18392 patches are applied).
  Small diffs in VAR only after NDData arithmetic.

Author: jehturner

geminidr/ghost/recipes/sq/tests/test_reduce_arc.py

@@ -62,8 +62,12 @@ def test_reduce_arc(input_filename, caldict, arm, path_to_inputs, path_to_refs):
     adout = astrodata.open(os.path.join("calibrations", "processed_arc", output_filename))
     adref = astrodata.open(os.path.join(path_to_refs, output_filename))
     # Changed timestamp kw from STCKARCS -> STACKARC and don't have time to
-    # re-upload reference, so just add these to the "ignore" list
-    assert ad_compare(adref, adout, ignore_kw=['PROCARC', 'STACKARC', 'STCKARCS'])
+    # re-upload reference, so just add these to the "ignore" list.
+    # The max_miss arg accommodates just a few pixels with large (~0.1%) diffs
+    # between NumPy 1 & 2 after the extraction step, while rtol allows for
+    # ~1e-7 rounding differences in VAR after NDData arithmetic on NumPy 2.
+    assert ad_compare(adref, adout, max_miss=12, rtol=1e-6,
+                      ignore_kw=['PROCARC', 'STACKARC', 'STCKARCS'])
 
     # Need to evaluate WFIT
     arm = GhostArm(arm=adout.arm(), mode=adout.res_mode())

geminidr/ghost/recipes/sq/tests/test_reduce_flat.py

@@ -57,7 +57,7 @@ def test_reduce_flat(change_working_dir, input_filename, bias, arm,
         output_filename = p.streams['main'][0].filename
         adout = astrodata.open(os.path.join("calibrations", "processed_flat", output_filename))
         adref = astrodata.open(os.path.join(path_to_refs, output_filename))
-        assert ad_compare(adref, adout, ignore_kw=['PROCFLAT'])
+        assert ad_compare(adref, adout, rtol=1e-6, ignore_kw=['PROCFLAT'])
 
         # Comparison doesn't include "exotic" extensions
         assert hasattr(adout[0], "BLAZE")

geminidr/ghost/recipes/sq/tests/test_reduce_slit.py

@@ -88,4 +88,4 @@ def test_reduce_slit_science(input_filename, caldict, path_to_inputs,
         for output_filename in [ad.filename for ad in p.streams['main']]:
             adout = astrodata.open(os.path.join("calibrations", "processed_slit", output_filename))
             adref = astrodata.open(os.path.join(path_to_refs, output_filename))
-            assert ad_compare(adref, adout, ignore_kw=['PRSLITIM'])
+            assert ad_compare(adref, adout, rtol=1e-6, ignore_kw=['PRSLITIM'])