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