Quantcast
Channel: Intel® oneAPI Math Kernel Library & Intel® Math Kernel Library
Viewing all articles
Browse latest Browse all 2652

How to avoid redundant factorizations in Pardiso for a transient analysis

$
0
0

Hello,

For a transient problem, [A]{x}t={b}t ,the RHS {b}is time-varying and depends on the previous {b}t-dt, but [A] 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/...

Thx!


Viewing all articles
Browse latest Browse all 2652

Latest Images

Trending Articles



Latest Images

<script src="https://jsc.adskeeper.com/r/s/rssing.com.1596347.js" async> </script>