How to avoid redundant factorizations in Pardiso for a transient analysis

Xudong_L_ — Thu, 30 Nov 2017 12:15:14 GMT

Hello,

For a transient problem, {x}_t={b}_t ,the RHS {b}_tis time-varying and depends on the previous {b}_t-dt, but is constant. So during a complete analysis, we need to do factorization only once, and use this factorization to solve successive RHS.

In this situation, multi-rhs is not suitable, and I tried to call the complete Pardiso (phase=11, 22, 33 and not execute phase=-1) at the 1st time step, then only executed phase=33 at later time steps. This way did not work however.

Do you have any suggestion?

I know there is a similar question, but the answers below don't seem to work either.

https://software.intel.com/en-us/forums/intel-math-kernel-library/topic/297075#

Thx!

Hello

Ying_H_Intel — Tue, 02 Jan 2018 08:22:56 GMT

Hello

As i understand, the Pardiso phase are designed for the cases , call phase 12 one time , then call phase =33 with different RHS in loop. Do you have test case to show the problem?

Best Regards,

Ying

topic How to avoid redundant factorizations in Pardiso for a transient analysis in Intel® oneAPI Math Kernel Library

How to avoid redundant factorizations in Pardiso for a transient analysis

Hello