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

## PARDISO-CSR format

Beginner
1,098 Views

Hi Everyone,

I wanted to use the PARDISO solver for my finite element Fortran code. I found the sample code (pardiso_sym_f90) from MKL library examples, and I was able to run the sample code. However, I could not find the routine to convert my matrix (unsymmetric banded matrix) to CSR format to use in PARDISO. Would you please advise me about that? Thank you.  Reza

1 Solution
Employee
1,077 Views

Hello Reza!

To answer directly your question: AFAIK, we don't have such a routine in oneMKL (but one can write some relatively simple code to do that).

I'd like to ask a question though and I think you should have a good answer before proceeding with PARDISO:
If you already have your matrix in the banded format, why don't you use a linear solver which is specialized for banded matrices?
E.g. you can find one here in LAPACK:
if you choose general + band storage matrix type.

The power of PARDISO is that it is a direct sparse linear solver which works for sparse matrices with general sparsity patterns. While it can still perform quite well on matrices with simple regular patterns (like banded), there are specialized solvers for banded matrices which should work better (because they are specialized for this kind of matrices).

Best,
Kirill

4 Replies
Employee
1,078 Views

Hello Reza!

To answer directly your question: AFAIK, we don't have such a routine in oneMKL (but one can write some relatively simple code to do that).

I'd like to ask a question though and I think you should have a good answer before proceeding with PARDISO:
If you already have your matrix in the banded format, why don't you use a linear solver which is specialized for banded matrices?
E.g. you can find one here in LAPACK:
if you choose general + band storage matrix type.

The power of PARDISO is that it is a direct sparse linear solver which works for sparse matrices with general sparsity patterns. While it can still perform quite well on matrices with simple regular patterns (like banded), there are specialized solvers for banded matrices which should work better (because they are specialized for this kind of matrices).

Best,
Kirill

Beginner
1,071 Views

Hi Kirill,

That’s right. I am using LAPACK for my code. I just wanted to try PARDISO if I can reduce the computational time on my program.

Best Regards,

Reza

Moderator
1,045 Views

Hi Reza,

If you don't have any other issues let us know if we can close the thread.

Regards,

Vidya

Moderator
1,024 Views

Hi Reza,

Thanks for the confirmation!

As this issue has been resolved, we will no longer respond to this thread.