Avoid copying array when strides are default (#1031)

Co-authored-by: Alexander-Makaryev <alexander.makaryev@gmail.com>

Author: densmirn

dpnp/backend/kernels/dpnp_krnl_common.cpp

@@ -567,9 +567,9 @@ void dpnp_matmul_c(void* result_out,
     }
 
     cl::sycl::event event;
-    DPNPC_ptr_adapter<_DataType> input1_ptr(input1_in, size_m * size_k, true);
-    DPNPC_ptr_adapter<_DataType> input2_ptr(input2_in, size_k * size_n, true);
-    DPNPC_ptr_adapter<_DataType> result_ptr(result_out, size_m * size_n, true, true);
+    DPNPC_ptr_adapter<_DataType> input1_ptr(input1_in, size_m * size_k);
+    DPNPC_ptr_adapter<_DataType> input2_ptr(input2_in, size_k * size_n);
+    DPNPC_ptr_adapter<_DataType> result_ptr(result_out, size_m * size_n, false, true);
     _DataType* array_1 = input1_ptr.get_ptr();
     _DataType* array_2 = input2_ptr.get_ptr();
     _DataType* result = result_ptr.get_ptr();

dpnp/dpnp_iface.py

@@ -211,7 +211,9 @@ def get_dpnp_descriptor(ext_obj, copy_when_strides=True):
     # then this behavior can be disabled with setting "copy_when_strides"
     if copy_when_strides and getattr(ext_obj, "strides", None) is not None:
         # TODO: replace this workaround when usm_ndarray will provide such functionality
-        ext_obj = array(ext_obj)
+        shape_offsets = tuple(numpy.prod(ext_obj.shape[i+1:], dtype=numpy.int64) for i in range(ext_obj.ndim))
+        if ext_obj.strides != shape_offsets:
+            ext_obj = array(ext_obj)
 
     dpnp_desc = dpnp_descriptor(ext_obj)
     if dpnp_desc.is_valid: