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().
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().