- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

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

Link Copied

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

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

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

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?

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page