[petsc-users] MatMkl_CPardisoSetCntl

[Matthew Overholt]

Hi All,

 

I am using the Intel MKL CPardiso library as a PC direct solver, and I am
trying to figure out how to properly set options (the Pardiso and CPardiso
"iparm" parameter values in the Intel docs).

 

Q1: To set the Pivoting perturbation, for example, is the correct call:

ierr = MatMkl_CPardisoSetCntl( K, 10, 13 );

where the KSP matrix to be inverted is "K" and the value that I want is
"13"?  This is the "-mat_mkl_pardiso_10" option listed on the
MATSOLVERMKL_PARDISO page, and the "iparm[9]" option in the Intel docs.

 

Q2: Where do I make this call, before "MatSetUp( K )", at some point during
the creation and setup of KSP, or after calling
"PCFactorSetMatSolverPackage( pc, MATSOLVERMKL_CPARDISO )"?  I've tried many
different combinations and none of them seem to work.

 

Since the effect of changing the pivoting perturbation may not be obvious, I
also tried setting the number of OpenMP threads to use within Pardiso (icntl
= 65, without the environmental variable MKL_NUM_THREADS present) and my
setting was ignored (Pardiso defaulted to using the maximum number of cores
present).

 

Thanks in advance,

Matt Overholt



---
This email has been checked for viruses by Avast antivirus software.
https://www.avast.com/antivirus
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20160822/925d501c/attachment.html>