Hi!
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!
/Anna