topic Any routines for sparse symmetric matrix-matrix production? in Intel® oneAPI Math Kernel Library

Any routines for sparse symmetric matrix-matrix production?

Wang__Hanrui — Sun, 07 Jul 2019 01:41:09 GMT

Hi,

I notice that MKL has special support for symmetric dense matrix production. Also MKL inspector-executor sparse BLAS routine supports sparse matrix-matrix production. My task is to conduct matrix production of two sparse symmetric matrix. I am wondering if there is any special routines supporting for this kind of special case?

Thank you very much!

Hello, we don't have (or,

Kirill_V_Intel — Tue, 09 Jul 2019 01:12:21 GMT

Hello, we don't have (or, more precisely, don't need) special routines for this case. Optimizations are done internally inside our functionality for sparse matrix-matrix product, like mkl_sparse_spmm (https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-spmm) and mkl_sparse_sp2m (https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-sp2m). Also, if you actually want to compute A * A^T or a triple product, you can use mkl_sparse_syrk(https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-syrk) or mkl_sparse_sypr (https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-sypr).

You can also find the list of all routines present in IE Sparse BLAS in the documentation.

Hope this helps!

Best,
Kirill

Quote:Voronin, Kirill (Intel)

Wang__Hanrui — Tue, 09 Jul 2019 03:26:45 GMT

Voronin, Kirill (Intel) wrote:
Hello, we don't have (or, more precisely, don't need) special routines for this case. Optimizations are done internally inside our functionality for sparse matrix-matrix product, like mkl_sparse_spmm (https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-spmm) and mkl_sparse_sp2m (https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-sp2m). Also, if you actually want to compute A * A^T or a triple product, you can use mkl_sparse_syrk(https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-syrk) or mkl_sparse_sypr (https://software.intel.com/en-us/mkl-developer-reference-c-mkl-sparse-sypr).
You can also find the list of all routines present in IE Sparse BLAS in the documentation.
Hope this helps!
Best,
Kirill

Hi Kirill,

OK, I understand. Thank you very much!

Best,

Hanrui

Re: Hello, we don't have (or,

miao0109 — Thu, 15 Jul 2021 01:57:20 GMT

hello , I learned that mkl provides the mkl_sparse_syrk interface, but it does not seem to support C:=(A^T)*B*A When the matrix is a complex . I have this requirement. When can I support it?

Re: Hello, we don't have (or,

Kirill_V_Intel — Thu, 15 Jul 2021 20:44:55 GMT

Hi!

Could you please be a bit more precise with describing your case?

First of all, the routine for triple product is mkl_sparse_sypr https://software.intel.com/content/www/us/en/develop/documentation/onemkl-developer-reference-c/top/blas-and-sparse-blas-routines/inspector-executor-sparse-blas-routines/inspector-executor-sparse-blas-execution-routines/mkl-sparse-sypr.html.

Second, it supports complex data types. However, there is a set of possible input options coming from descrB argument. E.g. the routine supports only symmetric/hermitian matrices B. Maybe, your matrix is not hermitian? Or anything else which does not meet the requirements listed in the documentation?

Thanks,
Kirill