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

Pardisoinit

Roman_F_
Beginner
112 Views

I guess this is another bug: the function pardisoinit resets iparm[34] to zero.

I must reset iparm[34] = 1 between calling pardisoinit() and pardiso().

0 Kudos
1 Reply
mecej4
Black Belt
112 Views

One of the purposes of the existence of pardisoinit() is stated in the MKL manual as  "...set default iparm values in accordance with the matrix type". You do not need to call pardisoinit() at all; if you do, do not set any iparm[] values before calling pardisoinit. As you found out, many of those values will be set in pardisoinit().

Reply