Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Roman_F_
Beginner
31 Views

Pardisoinit

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
31 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