Extend getrs with trans_code argument

Author: vlad-perevezentsev

dpnp/backend/extensions/lapack/getrs.cpp

@@ -166,6 +166,7 @@ std::pair<sycl::event, sycl::event>
           const dpctl::tensor::usm_ndarray &a_array,
           const dpctl::tensor::usm_ndarray &ipiv_array,
           const dpctl::tensor::usm_ndarray &b_array,
+          const int trans_code,
           const std::vector<sycl::event> &depends)
 {
     const int a_array_nd = a_array.get_ndim();
@@ -264,11 +265,20 @@ std::pair<sycl::event, sycl::event>
     const std::int64_t lda = std::max<size_t>(1UL, n);
     const std::int64_t ldb = std::max<size_t>(1UL, n);
 
-    // Use transpose::T if the LU-factorized array is passed as C-contiguous.
-    // For F-contiguous we use transpose::N.
-    oneapi::mkl::transpose trans = is_a_array_c_contig
-                                       ? oneapi::mkl::transpose::T
-                                       : oneapi::mkl::transpose::N;
+    oneapi::mkl::transpose trans;
+    switch (trans_code) {
+    case 0:
+        trans = oneapi::mkl::transpose::N;
+        break;
+    case 1:
+        trans = oneapi::mkl::transpose::T;
+        break;
+    case 2:
+        trans = oneapi::mkl::transpose::C;
+        break;
+    default:
+        throw py::value_error("`trans_code` must be 0 (N), 1 (T), or 2 (C)");
+    }
 
     char *a_array_data = a_array.get_data();
     char *b_array_data = b_array.get_data();

dpnp/backend/extensions/lapack/getrs.hpp

@@ -37,6 +37,7 @@ extern std::pair<sycl::event, sycl::event>
           const dpctl::tensor::usm_ndarray &a_array,
           const dpctl::tensor::usm_ndarray &ipiv_array,
           const dpctl::tensor::usm_ndarray &b_array,
+          const int trans_code,
           const std::vector<sycl::event> &depends = {});
 
 extern void init_getrs_dispatch_vector(void);

dpnp/backend/extensions/lapack/lapack_py.cpp

@@ -160,7 +160,8 @@ PYBIND11_MODULE(_lapack_impl, m)
           "the solves of linear equations with an LU-factored "
           "square coefficient matrix, with multiple right-hand sides",
           py::arg("sycl_queue"), py::arg("a_array"), py::arg("ipiv_array"),
-          py::arg("b_array"), py::arg("depends") = py::list());
+          py::arg("b_array"), py::arg("trans_code"),
+          py::arg("depends") = py::list());
 
     m.def("_orgqr_batch", &lapack_ext::orgqr_batch,
           "Call `_orgqr_batch` from OneMKL LAPACK library to return "