unmute tests

Author: vtavana

tests/test_fft.py

@@ -13,14 +13,8 @@
     get_all_dtypes,
     get_complex_dtypes,
     get_float_dtypes,
-    is_cpu_device,
 )
 
-# aspects of default device:
-_def_device = dpctl.SyclQueue().sycl_device
-_def_dev_has_fp64 = _def_device.has_aspect_fp64
-is_gpu_with_fp64 = not is_cpu_device() and _def_dev_has_fp64
-
 
 # TODO: `assert_dtype_allclose` calls in this file have `check_only_type_kind=True`
 # since stock NumPy is currently used in public CI for code coverege which
@@ -396,6 +390,10 @@ def test_hfft_1D(self, dtype, n, norm):
     def test_hfft_1D_complex(self, dtype, n, norm):
         x = dpnp.linspace(-1, 1, 11)
         a = dpnp.sin(x) + 1j * dpnp.cos(x)
+        # input should be Hermitian
+        a[0].imag = 0
+        if n in [None, 20]:  # Nyquist mode
+            a[-1].imag = 0
         a = dpnp.asarray(a, dtype=dtype)
         a_np = dpnp.asnumpy(a)
 
@@ -446,13 +444,16 @@ def test_fft_1D(self, dtype, n, norm):
         # but dpnp return float32 if input is float32
         assert_dtype_allclose(result, expected, check_only_type_kind=True)
 
-    @pytest.mark.skipif(is_gpu_with_fp64, reason="MKLD17702")
     @pytest.mark.parametrize("dtype", get_complex_dtypes())
     @pytest.mark.parametrize("n", [None, 5, 20])
     @pytest.mark.parametrize("norm", ["forward", "backward", "ortho"])
     def test_fft_1D_complex(self, dtype, n, norm):
         x = dpnp.linspace(-1, 1, 11)
         a = dpnp.sin(x) + 1j * dpnp.cos(x)
+        # input should be Hermitian
+        a[0].imag = 0
+        if n in [None, 20]:  # Nyquist mode
+            a[-1].imag = 0
         a = dpnp.asarray(a, dtype=dtype)
         a_np = dpnp.asnumpy(a)
 
@@ -473,29 +474,39 @@ def test_fft_1D_on_2D_array(self, dtype, n, axis, norm, order):
         expected = numpy.fft.irfft(a_np, n=n, axis=axis, norm=norm)
         assert_dtype_allclose(result, expected, check_only_type_kind=True)
 
-    @pytest.mark.skipif(is_gpu_with_fp64, reason="MKLD17702")
     @pytest.mark.parametrize("dtype", get_complex_dtypes())
     @pytest.mark.parametrize("n", [None, 5, 8])
     @pytest.mark.parametrize("axis", [0, 1, 2])
     @pytest.mark.parametrize("norm", ["forward", "backward", "ortho"])
     @pytest.mark.parametrize("order", ["C", "F"])
     def test_fft_1D_on_3D_array(self, dtype, n, axis, norm, order):
-        x1 = numpy.random.uniform(-10, 10, 24)
-        x2 = numpy.random.uniform(-10, 10, 24)
+        x1 = numpy.random.uniform(-10, 10, 120)
+        x2 = numpy.random.uniform(-10, 10, 120)
         a_np = numpy.array(x1 + 1j * x2, dtype=dtype).reshape(
-            2, 3, 4, order=order
+            4, 5, 6, order=order
         )
+        # each 1-D array of input should be Hermitian
+        a_np[0].imag = 0
+        a_np[:, 0, :].imag = 0
+        a_np[..., 0].imag = 0
+        # If output size is even, this is Nyquist mode
+        a_np[-1].imag = 0
+        a_np[:, -1, :].imag = 0
+        a_np[..., -1].imag = 0
         a = dpnp.asarray(a_np)
 
         result = dpnp.fft.irfft(a, n=n, axis=axis, norm=norm)
         expected = numpy.fft.irfft(a_np, n=n, axis=axis, norm=norm)
         assert_dtype_allclose(result, expected, check_only_type_kind=True)
 
-    @pytest.mark.skipif(is_gpu_with_fp64, reason="MKLD17702")
     @pytest.mark.parametrize("n", [None, 5, 20])
     def test_fft_usm_ndarray(self, n):
         x = dpt.linspace(-1, 1, 11)
         a = dpt.sin(x) + 1j * dpt.cos(x)
+        # input should be Hermitian
+        a[0] = dpt.sin(x[0])
+        if n in [None, 20]:  # Nyquist mode
+            a[-1] = dpt.sin(x[-1])
         a_usm = dpt.asarray(a, dtype=dpt.complex64)
         a_np = dpt.asnumpy(a_usm)
         out_shape = n if n is not None else 2 * (a_usm.shape[0] - 1)
@@ -506,13 +517,16 @@ def test_fft_usm_ndarray(self, n):
         expected = numpy.fft.irfft(a_np, n=n)
         assert_dtype_allclose(result, expected, check_only_type_kind=True)
 
-    @pytest.mark.skipif(is_gpu_with_fp64, reason="MKLD17702")
     @pytest.mark.parametrize("dtype", get_complex_dtypes())
     @pytest.mark.parametrize("n", [None, 5, 20])
     @pytest.mark.parametrize("norm", ["forward", "backward", "ortho"])
     def test_fft_1D_out(self, dtype, n, norm):
         x = dpnp.linspace(-1, 1, 11)
         a = dpnp.sin(x) + 1j * dpnp.cos(x)
+        # input should be Hermitian
+        a[0].imag = 0
+        if n in [None, 20]:  # Nyquist mode
+            a[-1].imag = 0
         a = dpnp.asarray(a, dtype=dtype)
         a_np = dpnp.asnumpy(a)
 

tests/test_sycl_queue.py

@@ -1239,7 +1239,7 @@ def test_fft(func, device):
 
     expected = getattr(numpy.fft, func)(data)
     result = getattr(dpnp.fft, func)(dpnp_data)
-    assert_dtype_allclose(result, expected)
+    assert_dtype_allclose(result, expected, factor=16)
 
     expected_queue = dpnp_data.get_array().sycl_queue
     result_queue = result.get_array().sycl_queue

tests/test_umath.py

@@ -313,7 +313,7 @@ def test_large_values(self, dtype):
         assert_dtype_allclose(result, expected)
 
 
-class TestLogaddexp:
+class TestLogAddExp:
     @pytest.mark.parametrize("dtype", get_all_dtypes(no_complex=True))
     def test_logaddexp(self, dtype):
         np_array1, np_array2, expected = _get_numpy_arrays_2in_1out(