=== PARDISO: solving a symmetric indefinite system === Matrix checker is turned ON 1-based array indexing is turned ON PARDISO double precision computation is turned ON METIS algorithm at reorder step is turned ON Scaling is turned ON Matching is turned ON Summary: ( reordering phase ) ================ Times: ====== Time spent in calculations of symmetric matrix portrait (fulladj): 0.000174 s Time spent in reordering of the initial matrix (reorder) : 0.000003 s Time spent in symbolic factorization (symbfct) : 0.002724 s Time spent in data preparations for factorization (parlist) : 0.000195 s Time spent in allocation of internal data structures (malloc) : 0.001730 s Time spent in matching/scaling : 0.043543 s Time spent in additional calculations : 0.004403 s Total time spent : 0.052773 s Statistics: =========== Parallel Direct Factorization is running on 1 OpenMP < Linear system Ax = b > number of equations: 3429 number of non-zeros in A: 90125 number of non-zeros in A (%): 0.766497 number of right-hand sides: 1 < Factors L and U > number of columns for each panel: 80 number of independent subgraphs: 0 < Preprocessing with state of the art partitioning metis> number of supernodes: 1398 size of largest supernode: 97 number of non-zeros in L: 398526 number of non-zeros in U: 1 number of non-zeros in L+U: 398527 === PARDISO is running in In-Core mode, because iparam(60)=0 === Percentage of computed non-zeros for LL^T factorization 1 % 2 % 3 % 4 % 5 % 6 % 7 % 8 % 9 % 10 % 11 % 12 % 13 % 14 % 15 % 16 % 17 % 19 % 20 % 21 % 22 % 23 % 24 % 25 % 26 % 29 % 30 % 32 % 33 % 34 % 35 % 36 % 37 % 38 % 39 % 40 % 41 % 42 % 43 % 44 % 45 % 47 % 50 % 54 % 55 % 56 % 57 % 58 % 59 % 60 % 61 % 62 % 63 % 64 % 65 % 66 % 67 % 68 % 72 % 73 % 74 % 75 % 76 % 77 % 78 % 79 % 80 % 81 % 82 % 83 % 84 % 85 % 86 % 87 % 88 % 89 % 93 % 94 % 95 % 97 % 99 % 100 % === PARDISO: solving a symmetric indefinite system === Single-level factorization algorithm is turned ON Summary: ( factorization phase ) ================ Times: ====== Time spent in copying matrix to internal data structure (A to LU): 0.000000 s Time spent in factorization step (numfct) : 0.032004 s Time spent in allocation of internal data structures (malloc) : 0.000830 s Time spent in additional calculations : 0.000002 s Total time spent : 0.032835 s Statistics: =========== Parallel Direct Factorization is running on 1 OpenMP < Linear system Ax = b > number of equations: 3429 number of non-zeros in A: 90125 number of non-zeros in A (%): 0.766497 number of right-hand sides: 1 < Factors L and U > number of columns for each panel: 80 number of independent subgraphs: 0 < Preprocessing with state of the art partitioning metis> number of supernodes: 1398 size of largest supernode: 97 number of non-zeros in L: 398526 number of non-zeros in U: 1 number of non-zeros in L+U: 398527 gflop for the numerical factorization: 0.064288 gflop/s for the numerical factorization: 2.008787 Number of negative eigenvalues according to Pardiso is 1311. However, the correct number for this matrix should be 1309. Eigenvalues computed by Eigen (dense): [-243.214,..., -1.23553e-05, 0.00443743,..., 223.752] number of negative eigenvalues is 1309