topic Hi! in Intel® oneAPI Math Kernel Library

Pardiso matrix size change over the simulation

Tabrizian__Amin — Thu, 04 Jul 2019 20:52:34 GMT

Dear all,

I have just started learning PARDISO and I face a question:

How can we use PARDISO when our coefficient matrix size changes over the simulation?

I know that we can simply solve this issue by reinitializing PARDISO but maybe a special order of steps needs to be followed to avoid this inefficient using of PARDISO.

Thanks.

Hi!

Kirill_V_Intel — Tue, 09 Jul 2019 02:03:00 GMT

Hi!

First of all, please clarify what do you mean by changing the coefficient matrix size? And how do you want to re-use the original matrix on the high-level? For example do you have an extended system where the original matrix is a sub-block (maybe with slightly changed matrix values but with the same pattern)? Or do you have a completely different new larger matrix (for example when you change the discretization step when solving numerically some PDEs)?

Best,
Kirill

Hi Kirill,

Tabrizian__Amin — Tue, 09 Jul 2019 10:21:28 GMT

Hi Kirill,

Thank you for your reply. I mean how we can use PARDISO when size of our system of equations changes over the simulation? And about re-using the original matrix, our original matrix may change by adding some rows(columns) or omitting some.

Thanks,

Amin

Hi Amin,

Kirill_V_Intel — Tue, 09 Jul 2019 15:35:48 GMT

Hi Amin,

You need to formulate on the level of mathematics how you would like to re-use the matrices. PARDISO is a way to factor the matrix and then solve the system. For example you can use one of the matrices as a preconditioner. Or you can use the Schur complement technique.

Best,
Kirill