PARDISO's CGS iterative refinement

Intel® oneAPI Math Kernel Library

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

PARDISO's CGS iterative refinement

377 Views

Dear all,

I found the strange output when I use the PARDISO of Intel MKL 10.0.3.021 for Windows to solve matrix. The main parameters I set is as follows:

mtype=13
iparm(1)=1
iparm(2)=2
iparm(3)=1
iparm(4)=0
iparm(5:6)=0
iparm(8)=0
iparm(10)=8
iparm(11)=1
iparm(13)=1
iparm(18:19)=0

However, after solution, the output iparm(20)=-34, I read the manual of this function and found the parameter is related to the CGS iterative refinement, but I didn't ask for CGS iterative refinement at all!

On the other hand, if I ask for CGS iterative refinement:

phase=33
iparm(4)=61

The program will crash. Could anyone explain the output of iparm(20) and tell me how to let the PARDISO to do some CGS iterative refinements if the solution is not satisfied?

Thanks,
Zhanghong Tang

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.