Alter elements in diagonal of D-matrix during the LDL' factorization

[Dalklint__Anna]

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

[Gennady_F_Intel]

There are a couple of C examples which you find out into mklroot/examples/solverc/source folder.

[Dalklint__Anna]

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?