MKL 11, minor 3, update 4, version 110304, build date 20160823 2017-07-17 15:00:06.925 DETL| Loading matrix from C:/tmp/matrix_explosion.data... Size is 1563 x 1563 . Starting the clock... Re-ordering the matrix using Pardiso. 2017-07-17 15:00:06.961 DBUG| non-zero iparm values: 2017-07-17 15:00:06.961 DBUG| iparm[0] = 1 2017-07-17 15:00:06.961 DBUG| iparm[1] = 2 2017-07-17 15:00:06.961 DBUG| iparm[9] = 8 2017-07-17 15:00:06.961 DBUG| iparm[10] = 1 2017-07-17 15:00:06.961 DBUG| iparm[12] = 1 2017-07-17 15:00:06.961 DBUG| iparm[17] = -1 2017-07-17 15:00:06.961 DBUG| iparm[20] = 1 2017-07-17 15:00:06.961 DBUG| iparm[26] = 1 Factoring the matrix using Pardiso. Solving the system using Pardiso. Stopping the clock... 2017-07-17 15:00:07.068 DETL| Infinity residual norm ||Ax-b|| is 0.000629389. ||b|| is 0.99939. 2017-07-17 15:00:07.068 DETL| Relative infinity residual norm ||Ax-b||/||b|| is 0.000629773. Total solver time: 0.125855