Merge branch 'master' into resize_rot90

Author: vtavana

CHANGELOG.md

@@ -4,6 +4,15 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [0.17.0] - MM/DD/2025
+
+### Added
+
+### Change
+
+### Fixed
+
+
 ## [0.16.0] - 09/DD/2024
 
 This release reaches an important milestone by making offloading fully asynchronous. Calls to `dpnp` submit tasks for execution to DPC++ runtime and return without waiting for execution of these tasks to finish. The sequential semantics a user comes to expect from execution of Python script is preserved though.

CMakeLists.txt

@@ -1,7 +1,7 @@
 cmake_minimum_required(VERSION 3.21...3.27 FATAL_ERROR)
 
 project(dpnp
-  VERSION 0.16
+  VERSION 0.17
   LANGUAGES CXX
   DESCRIPTION "NumPy-like API accelerated by SYCL."
 )

dpnp/dpnp_iface_mathematical.py

@@ -224,8 +224,7 @@ def _gradient_num_diff_2nd_order_interior(
         # fix the shape for broadcasting
         shape = [1] * ndim
         shape[axis] = -1
-        # TODO: use shape.setter once dpctl#1699 is resolved
-        # a.shape = b.shape = c.shape = shape
+
         a = a.reshape(shape)
         b = b.reshape(shape)
         c = c.reshape(shape)

tests/third_party/cupy/core_tests/test_core.py

@@ -0,0 +1,124 @@
+import sys
+import unittest
+
+import numpy
+import pytest
+
+import dpnp as cupy
+from tests.third_party.cupy import testing
+
+
+class TestSize(unittest.TestCase):
+    # def tearDown(self):
+    #     # Free huge memory for slow test
+    #     cupy.get_default_memory_pool().free_all_blocks()
+
+    @testing.for_all_dtypes()
+    @testing.numpy_cupy_equal()
+    def test_size(self, xp, dtype):
+        a = xp.ndarray((2, 3), dtype=dtype)
+        return xp.size(a)
+
+    @testing.for_all_dtypes()
+    @testing.numpy_cupy_equal()
+    def test_size_axis(self, xp, dtype):
+        a = xp.ndarray((2, 3), dtype=dtype)
+        return xp.size(a, axis=1)
+
+    @testing.for_all_dtypes()
+    def test_size_axis_error(self, dtype):
+        for xp in (numpy, cupy):
+            a = xp.ndarray((2, 3), dtype=dtype)
+            with pytest.raises(IndexError):
+                return xp.size(a, axis=3)
+
+    @testing.numpy_cupy_equal()
+    @testing.slow
+    def test_size_huge(self, xp):
+        a = xp.ndarray(2**32, "b")  # 4 GiB
+        return xp.size(a)
+
+
+_orders = {
+    order_arg: order_expect
+    for order_expect, order_args in [
+        ("C", ["C", "c", "CONTIGUOUS", "", None]),
+        ("F", ["F", "f", "FORTRAN"]),
+    ]
+    for order_arg in order_args
+}
+
+
+@pytest.mark.skip("no cupy._core submodule")
+class TestOrder(unittest.TestCase):
+    @testing.for_orders(_orders.keys())
+    def test_ndarray(self, order):
+        order_expect = _orders[order]
+        a = core.ndarray((2, 3), order=order)
+        expect_c = order_expect == "C"
+        expect_f = order_expect == "F"
+        assert a.flags.c_contiguous == expect_c
+        assert a.flags.f_contiguous == expect_f
+
+
+@pytest.mark.skip("min_scalar_type() is not supported")
+class TestMinScalarType:
+    def test_scalar(self):
+        for v in (-129, -128, 0, 1.2, numpy.inf):
+            assert cupy.min_scalar_type(v) is numpy.min_scalar_type(v)
+
+    @testing.for_all_dtypes()
+    def test_numpy_scalar(self, dtype):
+        sc = dtype(1)
+        for v in (sc, [sc, sc]):
+            assert cupy.min_scalar_type(v) is numpy.min_scalar_type(v)
+
+    @testing.for_all_dtypes()
+    def test_cupy_scalar(self, dtype):
+        sc = cupy.array(-1).astype(dtype)
+        for v in (sc, [sc, sc]):
+            assert cupy.min_scalar_type(v) is sc.dtype
+
+    @testing.for_all_dtypes()
+    def test_numpy_ndarray(self, dtype):
+        arr = numpy.array([[-1, 1]]).astype(dtype)
+        for v in (arr, (arr, arr)):
+            assert cupy.min_scalar_type(v) is numpy.min_scalar_type(v)
+
+    @testing.for_all_dtypes()
+    def test_cupy_ndarray(self, dtype):
+        arr = cupy.array([[-1, 1]]).astype(dtype)
+        for v in (arr, (arr, arr)):
+            assert cupy.min_scalar_type(v) is arr.dtype
+
+
+@testing.parameterize(
+    *testing.product(
+        {
+            "cxx": (None, "--std=c++11"),
+        }
+    )
+)
+@pytest.mark.skip("compiling cupy headers are not supported")
+class TestCuPyHeaders(unittest.TestCase):
+    def setUp(self):
+        self.temporary_cache_dir_context = test_raw.use_temporary_cache_dir()
+        self.cache_dir = self.temporary_cache_dir_context.__enter__()
+        self.header = "\n".join(
+            ["#include <" + h + ">" for h in core._cupy_header_list]
+        )
+
+    def tearDown(self):
+        self.temporary_cache_dir_context.__exit__(*sys.exc_info())
+
+    def test_compiling_core_header(self):
+        code = r"""
+        extern "C" __global__ void _test_ker_() { }
+        """
+        code = self.header + code
+        options = () if self.cxx is None else (self.cxx,)
+        ker = cupy.RawKernel(
+            code, "_test_ker_", options=options, backend="nvrtc"
+        )
+        ker((1,), (1,), ())
+        cupy.cuda.Device().synchronize()

tests/third_party/cupy/core_tests/test_ndarray.py

@@ -154,11 +154,18 @@ def test_unsupported_type(self):
         with pytest.raises(TypeError):
             cupy.array(arr)
 
+    @pytest.mark.skip("no ndim limit")
+    @testing.with_requires("numpy>=2.0")
+    @testing.numpy_cupy_array_equal()
+    def test_upper_limit_ndim(self, xp):
+        shape = [1 for i in range(64)]
+        return xp.zeros(shape, dtype=xp.int8)
+
     @pytest.mark.skip("no ndim limit")
     def test_excessive_ndim(self):
         for xp in (numpy, cupy):
             with pytest.raises(ValueError):
-                xp.ndarray(shape=[1 for i in range(33)], dtype=xp.int8)
+                xp.ndarray(shape=[1 for i in range(65)], dtype=xp.int8)
 
 
 @testing.parameterize(
@@ -265,7 +272,8 @@ def test_shape_set(self, xp):
         return xp.array(arr.shape)
 
     @pytest.mark.skip(
-        "dpctl-1699: shape setter does not work with negative shape"
+        "dpctl-1699: shape setter does not work with negative shape "
+        "(no plan to support that)"
     )
     @testing.numpy_cupy_array_equal()
     def test_shape_set_infer(self, xp):
@@ -588,7 +596,6 @@ def test_output_type_mismatch(self):
                 wrap_take(a, i)
 
 
-@pytest.mark.skip("size() is not supported")
 class TestSize(unittest.TestCase):
     @testing.numpy_cupy_equal()
     def test_size_without_axis(self, xp):