fix: update test setup.

Author: paquiteau

tests/operators/test_orc.py

@@ -84,7 +84,7 @@ def test_orc_forward(
     ksp = orc_nufft.op(image_data)
     ksp_ideal = ext_mat @ image_data.ravel()
 
-    assert_correlate(ksp.squeeze(), ksp_ideal)
+    assert_correlate(ksp.squeeze(), ksp_ideal, slope_err=0.005)
 
 
 def test_orc_adjoint(orc_info, kspace_data):
@@ -93,4 +93,4 @@ def test_orc_adjoint(orc_info, kspace_data):
     img = orc_nufft.adj_op(kspace_data)
     img_ideal = ext_mat.conj().T @ kspace_data.ravel()
     img_ideal = img_ideal.reshape(orc_nufft.shape)
-    assert_correlate(img.squeeze(), img_ideal)
+    assert_correlate(img.squeeze(), img_ideal, slope_err=0.005)

tests/test_offres_exp_approx.py

@@ -26,10 +26,10 @@ def case_real2D(self, N=64, b0_range=(-300, 300)):
         b0_map, mask = make_b0map(2 * [N])
         return b0_map, None, mask
 
-    # def case_real3D(self, N=32, b0range=(-300, 300)):
-    #     """Create a real (B0 only) 3D field map."""
-    #     b0_map,  mask = make_b0map(3 * [N], b0range)
-    #     return b0_map, None, mask
+    def case_real3D(self, N=32, b0range=(-300, 300)):
+        """Create a real (B0 only) 3D field map."""
+        b0_map, mask = make_b0map(3 * [N], b0range)
+        return b0_map, None, mask
 
     def case_complex2D(self, N=64, b0range=(-300, 300), t2svalue=15.0):
         """Create a complex (R2* + 1j * B0) 2D field map."""
@@ -44,21 +44,23 @@ def case_complex2D(self, N=64, b0range=(-300, 300), t2svalue=15.0):
 
         return b0_map, r2s_map, mask
 
-    # def case_complex3D(self, N=32, b0range=(-300, 300), t2svalue=15.0):
-    #     """Create a complex (R2* + 1j * B0) 3D field map."""
-    #     # Generate real and imaginary parts
-    #     t2s_map, _ = make_t2smap(3 * [N], t2svalue)
-    #     b0_map, mask = make_b0map(3 * [N], b0range)
+    def case_complex3D(self, N=32, b0range=(-300, 300), t2svalue=15.0):
+        """Create a complex (R2* + 1j * B0) 3D field map."""
+        # Generate real and imaginary parts
+        t2s_map, _ = make_t2smap(3 * [N], t2svalue)
+        b0_map, mask = make_b0map(3 * [N], b0range)
 
-    #     # Convert T2* map to R2* map
-    #     t2s_map = t2s_map * 1e-3  # ms -> s
-    #     r2s_map = 1.0 / (t2s_map + 1e-9)  # Hz
-    #     r2s_map = mask * r2s_map
-    #     return b0_map, r2s_map, mask
+        # Convert T2* map to R2* map
+        t2s_map = t2s_map * 1e-3  # ms -> s
+        r2s_map = 1.0 / (t2s_map + 1e-9)  # Hz
+        r2s_map = mask * r2s_map
+        return b0_map, r2s_map, mask
 
 
 @_param_array_interface_np_cp
-@parametrize_with_cases("b0_map, r2s_map, mask", cases=CasesB0maps)
+@parametrize_with_cases(
+    "b0_map, r2s_map, mask", cases=[CasesB0maps.case_real2D, CasesB0maps.case_complex2D]
+)
 @parametrize("method", ["svd-full", "mti", "mfi"])
 @parametrize("L, lazy", [(40, True), (-1, True), (40, False)])
 def test_b0map_coeff(b0_map, r2s_map, mask, method, L, lazy, array_interface):