Using iterative-direct solver in pardiso

Intel® oneAPI Math Kernel Library

Ask questions and share information with other developers who use Intel® Math Kernel Library.

Using iterative-direct solver in pardiso

447 Views

I'm trying to use the iterative-direct solver in Pardiso by changing the default value of iparm[3]. However that didn't work.

I have a nonsymmetric matrix, so I set iparm[3] = 61. However, when I print out iparm[19], which contains the CG/CGS diagnostics information, it always says 0. That implies CG wasn't run at all. In my code, when I first solve the linear system, I set phase to be 13. After that, I set the phase to be 23 (so that it doesn't have to factorize every time). The input for my matrix type mtype is 11.

I then tried a symmetric positive definite matrix and set iparm[3] = 62. I set mtype to be 2. That didn't work either: iparm[19] returns 0.

Did I miss anything?

Link Copied

0 Replies

Community support is provided Monday to Friday. Other contact methods are available here.

Intel does not verify all solutions, including but not limited to any file transfers that may appear in this community. Accordingly, Intel disclaims all express and implied warranties, including without limitation, the implied warranties of merchantability, fitness for a particular purpose, and non-infringement, as well as any warranty arising from course of performance, course of dealing, or usage in trade.

For more complete information about compiler optimizations, see our Optimization Notice.