is there a way to access to access multiple factorizations of sparse matrices between subsequent PARDISO calls?
More precisely, I solve the usual problem Ax=B where B is a multi-column rhs for several subsystems in different time steps of my calculation.
Since in the manual entry for the PARDISO input parameter MAXFCT it is written:
"PARDISO can process several matrices with identical matrix sparsity pattern and store the factors of these matrices at the same time. Matrices with a different sparsity structure can be kept in memory with different memory address pointers pt."
I wonder if I would be able to avoid the factorizations after the first time step.
In other words, is it possible to retrieve (from memory or whatever) the multiple factorizations without asking PARDISO to actually redo the calculation? I tried i.e. declaring the data address pointers PT as global variables, but it failed.
Thank you in advance for any advice and regards.
p.s. I use Intel Fortran 11.1.048 with MKL 10.33