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

PARDISO, How to escape the occurrence of zero pivot error

suzuki__atsushi
Beginner
565 Views

Hi everyone,

I  have zero pivot error -4, but I know the error occurs in the condftion of my calculation.

So I want to have a solution  using zero pivot replacement (the value i_parm(10)) and without stopping calculation.

Is it possible and how should I set parameters?

 

0 Kudos
1 Reply
Kirill_V_Intel
Employee
565 Views

Hello atsushi,

Could you tell us more about your setup? Intel MKL PARDISO allows the user to enable pivoting (see iparm[9], iparm[20] at https://software.intel.com/en-us/mkl-developer-reference-c-pardiso-iparm-parameter). What is the matrix type in your test case? If it is for positive definite matrices then your input matrix should not be singular. Are you using matching? So, please give more information about the parameters which you're using for the solver, namely, iparm values and the matrix type.

Best,
Kirill

 

0 Kudos
Reply