oneMKL PARDISO complex matrix pivoting perturbation

[AnttiR]

I have been struggling with getting correct results with oneMKL PARDISO mtype=13 (complex and nonsymmetric matrix). The algorithm the solution is a part of is known to handle this type of matrixes, and the results are consistent with Python between numpy.linalg.solve and scipy.linalg.spsolve.

I had one specific problem with a small 4x4 matrix with the solve producing totally incorrect results. I accidentally discovered that setting iparm(10) (Fortran) to a very large value or 0 solved the problem as the number of perturbed pivots changed from one to zero (with the default value being 13).

My actual question is that is iparm(10)=0 a valid parameter for disabling the pivoting perturbation? It is not evident from the documentation (https://www.intel.com/content/www/us/en/docs/onemkl/developer-reference-fortran/2023-2/pardiso-iparm-parameter.html) if this is a valid entry or not.

[IntelSupport]

Hi,

Thanks for posting in Intel Communities.

Could you please provide us a sample reproducer, MKL version and OS details, so that we can replicate and investigate more at our end. Please confirm, is it okay to use C as well as Fortran?

Regards,

Jilani

[JilaniS_Intel]

Hi,

A gentle reminder:

Could you please provide us a sample reproducer, MKL version and OS details, so that we can replicate and investigate more at our end. Please confirm, is it okay to use C as well as Fortran?

Regards,

Jilani

[JilaniS_Intel]

Hi,

A gentle reminder:

We have not heard back from you. Could you please provide us with an update?

Regards,

Jilani