I want to solve a system with mtype = -2, i.e. PAP' = LDL', using pardiso MKL.
I have understood that the function pardiso_getdiag provides the diagonal of the initial and factorized system. However, I wonder if it is possible to obtain the diagonal of the D-matrix during the factorization? More precisely I want to be able to alter the diagonal elements of D which are close to zero (replace these with ones instead).
The function mkl_pardiso_pivot is supposed to "allows you to handle diagonal elements which arise during numerical factorization that are zero or near zero". However, I have problems understanding what I am supposed to provide as input to the function. Is there any examples that can be found?
I am using the fortran version of pardiso.
Thanks in advance for the help!
Yes I have heard about that. However, I cannot seem to find these examples. I can only find a folder named "solverf", not "solverc". Is it possible for you to send them to me?