Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
7243 Discussions

Pardiso format Vector Manipulation Library in MKL

Ngu_Soon_Hui
Beginner
714 Views

I just wonder whether there is any vector/matrix manipulation library in MKL if my data is stored in pardiso format? The library shall include, but not limited to factorization, linear system solver, matrix vector product.

This is because for band matrix, lapack has readily available matrix/vector manipulation library (see for example, sgbtrf and sgbmv), which comes in handy when one needs to do ARPACK eigen solvers ( see an example here). However, a lot of times we don't use band matrix storage format and instead, we opt for pardiso format ( is it CSR or CSC or a variant of any of the two? I'm not very sure).

A similar library for Pardiso would be great aid for everyone.

0 Kudos
1 Solution
ArturGuzik
Valued Contributor I
714 Views

I'm not entirely sure what you have in mind, however, you may take a look at the Sparse BLAS level 2 and 3 routines (in MKL), which provide (most probably) what you want, including all sort of multiplications, conversions etc.

A.

View solution in original post

0 Kudos
2 Replies
ArturGuzik
Valued Contributor I
715 Views

I'm not entirely sure what you have in mind, however, you may take a look at the Sparse BLAS level 2 and 3 routines (in MKL), which provide (most probably) what you want, including all sort of multiplications, conversions etc.

A.

0 Kudos
Gennady_F_Intel
Moderator
714 Views
Quoting Ngu Soon Hui

I just wonder whether there is any vector/matrix manipulation library in MKL if my data is stored in pardiso format? The library shall include, but not limited to factorization, linear system solver, matrix vector product.

This is because for band matrix, lapack has readily available matrix/vector manipulation library (see for example, sgbtrf and sgbmv), which comes in handy when one needs to do ARPACK eigen solvers ( see an example here). However, a lot of times we don't use band matrix storage format and instead, we opt for pardiso format ( is it CSR or CSC or a variant of any of the two? I'm not very sure).

A similar library for Pardiso would be great aid for everyone.

one more comment on this - Pardiso supports only one CSR format ( 1-based indexes).

--Gennady

0 Kudos
Reply