<div dir="ltr"><div dir="ltr">On Sun, May 23, 2021 at 11:17 AM Karl Yang <<a href="mailto:y.juntao@hotmail.com">y.juntao@hotmail.com</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>Hello,</div><br><div>I am using MUMPS direct solver for my project. I used the following options for solving my problem and it works in most cases. But for some cases I encounter a divergence error. But I think it is actually error due to MUMPS?</div><br><div>I'm not sure how to debug the error. It is appreciated if anyone familiar with MUMPS solver to offer me some guidance.</div></blockquote><div><br></div><div>This says</div><div><br></div><div>On return from DMUMPS, INFOG(1)=              -9<br> On return from DMUMPS, INFOG(2)=              22<br></div><div><br></div><div>that the internal work array for MUMPS is too small. I am not sure which option controls that.</div><div><br></div><div>  Thanks,</div><div><br></div><div>     Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div>regards</div><div>Juntao</div><br><div>MUMPS options<span>:</span></div><code><pre style="background-color:rgba(0,0,0,0.05);padding:0.2em 1em"><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-ksp_type"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"preonly"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-pc_type"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"cholesky"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-pc_factor_mat_solver_type"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"mumps"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_1"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"1"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_2"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"1"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_3"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"1"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_4"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"3"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_28"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"1"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_7"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"2"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div><div><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(220,220,170)">PetscOptionsSetValue</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">(</span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(86,156,214)">NULL</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"-mat_mumps_icntl_24"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">, </span></font><span style="color:rgb(212,212,212)"><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(206,145,120)">"1"</span></font></span><font style="font-family:"Droid Sans Mono",monospace,monospace,"Droid Sans Fallback""><span style="color:rgb(212,212,212)">);</span></font></div></pre></code><div> </div><br><code><pre style="background-color:rgba(0,0,0,0.05);padding:0.2em 1em">Log output from MUMPS and error message from PETSC at the bottom</pre></code><div>Entering DMUMPS 5.2.1 from C interface with JOB, N, NNZ =   1         240           2448</div><div>      executing #MPI =      1, without OMP</div><br><div> =================================================</div><div> MUMPS compiled with option -Dmetis</div><div> MUMPS compiled with option -Dptscotch</div><div> MUMPS compiled with option -Dscotch</div><div> This MUMPS version includes code for SAVE_RESTORE</div><div> =================================================</div><div>L D L^T Solver for general symmetric matrices</div><div>Type of parallelism: Working host</div><br><div> ****** ANALYSIS STEP ********</div><br><div> Scaling will be computed during analysis</div><div>Compute maximum matching (Maximum Transversal):  5</div><div> ... JOB =  5: MAXIMIZE PRODUCT DIAGONAL AND SCALE</div><br><div>Entering analysis phase with ...</div><div>                N        NNZ         LIW       INFO(1)</div><div>                240       2448        5137             0</div><div>Matrix entries:    IRN()   ICN()</div><div>           1      1           1      2           1      3</div><div>           1      4           1      5           1      6</div><div>           1      7           1      8           1      9</div><div>           1     10</div><div> Average density of rows/columns =   18</div><div> Average density of rows/columns =   18</div><div> Ordering based on AMF </div><div> Constrained Ordering based on AMF</div><div> Average density of rows/columns =   18</div><div> Average density of rows/columns =   18</div><div>NFSIZ(.)  =     0    38    14     0    33    33     0     0     0     0</div><br><div>FILS (.)  =     0   148     4   -96   224   163    20   -43     8     1</div><br><div>FRERE(.)  =   241    -5    -6   241     0    -2   241   241   241   241</div><br><br><div>Leaving analysis phase with  ...</div><div> INFOG(1)                                       =               0</div><div> INFOG(2)                                       =               0</div><div> -- (20) Number of entries in factors (estim.)  =            3750</div><div> --  (3) Real space for factors    (estimated)  =            4641</div><div> --  (4) Integer space for factors (estimated)  =            2816</div><div> --  (5) Maximum frontal size      (estimated)  =              38</div><div> --  (6) Number of nodes in the tree            =              56</div><div> -- (32) Type of analysis effectively used      =               1</div><div> --  (7) Ordering option effectively used       =               2</div><div> ICNTL(6) Maximum transversal option            =               0</div><div> ICNTL(7) Pivot order option                    =               2</div><div> ICNTL(14) Percentage of memory relaxation      =              20</div><div> Number of level 2 nodes                        =               0</div><div> Number of split nodes                          =               0</div><div> RINFOG(1) Operations during elimination (estim)= 7.137D+04</div><div> Ordering compressed/constrained (ICNTL(12))    =               3</div><br><div> MEMORY ESTIMATIONS ... </div><div> Estimations with standard Full-Rank (FR) factorization:</div><div>    Total space in MBytes, IC factorization      (INFOG(17)):           0</div><div>    Total space in MBytes,  OOC factorization    (INFOG(27)):           0</div><br><div> Elapsed time in analysis driver=       0.0016</div><br><div>Entering DMUMPS 5.2.1 from C interface with JOB, N, NNZ =   2         240           2448</div><div>      executing #MPI =      1, without OMP</div><br><br><br><div>****** FACTORIZATION STEP ********</div><br><div> GLOBAL STATISTICS PRIOR NUMERICAL FACTORIZATION ...</div><div> Number of working processes                =               1</div><div> ICNTL(22) Out-of-core option               =               0</div><div> ICNTL(35) BLR activation (eff. choice)     =               0</div><div> ICNTL(14) Memory relaxation                =              20</div><div> INFOG(3) Real space for factors (estimated)=            4641</div><div> INFOG(4) Integer space for factors (estim.)=            2816</div><div> Maximum frontal size (estimated)           =              38</div><div> Number of nodes in the tree                =              56</div><div> Memory allowed (MB -- 0: N/A )             =               0</div><div> Memory provided by user, sum of LWK_USER   =               0</div><div> Relative threshold for pivoting, CNTL(1)   =      0.1000D-01</div><div>  ZERO PIVOT DETECTION ON, THRESHOLD          =   2.8931920285365730E-020</div><div> INFINITE FIXATION </div><div> Effective size of S     (based on INFO(39))=                 7981</div><div> Elapsed time to reformat/distribute matrix =      0.0001</div><div> ** Memory allocated, total in Mbytes           (INFOG(19)):           0</div><div> ** Memory effectively used, total in Mbytes    (INFOG(22)):           0</div><div> ** Memory dynamically allocated for CB, total in Mbytes   :           0</div><br><div> Elapsed time for factorization             =      0.0006</div><br><div>Leaving factorization with ...</div><div> RINFOG(2)  Operations in node assembly     = 5.976D+03</div><div> ------(3)  Operations in node elimination  = 1.197D+05</div><div> INFOG (9)  Real space for factors          =            6193</div><div> INFOG(10)  Integer space for factors       =            3036</div><div> INFOG(11)  Maximum front size              =              42</div><div> INFOG(29)  Number of entries in factors    =            4896</div><div> INFOG(12)  Number of negative pivots       =              79</div><div> INFOG(13)  Number of delayed pivots        =             110</div><div> Number of 2x2 pivots in type 1 nodes       =               1</div><div> Number of 2X2 pivots in type 2 nodes       =               0</div><div> Nb of null pivots detected by ICNTL(24)    =               0</div><div> INFOG(28)  Estimated deficiency            =               0</div><div> INFOG(14)  Number of memory compress       =               0</div><br><div> Elapsed time in factorization driver=       0.0009</div><br><div>Entering DMUMPS 5.2.1 from C interface with JOB, N, NNZ =   3         240           2448</div><div>      executing #MPI =      1, without OMP</div><br><br><br><div> ****** SOLVE & CHECK STEP ********</div><br><div> GLOBAL STATISTICS PRIOR SOLVE PHASE ...........</div><div> Number of right-hand-sides                    =           1</div><div> Blocking factor for multiple rhs              =           1</div><div> ICNTL (9)                                     =           1</div><div>  --- (10)                                     =           0</div><div>  --- (11)                                     =           0</div><div>  --- (20)                                     =           0</div><div>  --- (21)                                     =           0</div><div>  --- (30)                                     =           0</div><div>  --- (35)                                     =           0</div><br><br><div> Vector solution for column            1</div><div> RHS</div><div>  -7.828363D-02 -3.255337D+00  1.054729D+00  1.379822D-01 -3.892113D-01</div><div>   1.433990D-01  1.089250D+00  2.252611D+00  3.215399D+00 -6.788806D-02</div><div> ** Space in MBYTES used for solve                        :         0</div><br><div> Leaving solve with ...</div><div> Time to build/scatter RHS        =       0.000003</div><div> Time in solution step (fwd/bwd)  =       0.000167</div><div>  .. Time in forward (fwd) step   =          0.000053</div><div>  .. Time in backward (bwd) step  =          0.000093</div><div> Time to gather solution(cent.sol)=       0.000000</div><div> Time to copy/scale dist. solution=       0.000000</div><br><div> Elapsed time in solve driver=       0.0004</div><div>*** Warning: Verbose output for PETScKrylovSolver not implemented, calling PETSc KSPView directly.</div><div>KSP Object: 1 MPI processes</div><div>  type: preonly</div><div>  maximum iterations=10000, initial guess is zero</div><div>  tolerances:  relative=1e-05, absolute=1e-50, divergence=10000.</div><div>  left preconditioning</div><div>  using NONE norm type for convergence test</div><div>PC Object: 1 MPI processes</div><div>  type: cholesky</div><div>    out-of-place factorization</div><div>    tolerance for zero pivot 2.22045e-14</div><div>    matrix ordering: natural</div><div>    factor fill ratio given 0., needed 0.</div><div>      Factored matrix follows:</div><div>        Mat Object: 1 MPI processes</div><div>          type: mumps</div><div>          rows=240, cols=240</div><div>          package used to perform factorization: mumps</div><div>          total: nonzeros=3750, allocated nonzeros=3750</div><div>          total number of mallocs used during MatSetValues calls=0</div><div>            MUMPS run parameters:</div><div>              SYM (matrix type):                   2 </div><div>              PAR (host participation):            1 </div><div>              ICNTL(1) (output for error):         1 </div><div>              ICNTL(2) (output of diagnostic msg): 1 </div><div>              ICNTL(3) (output for global info):   6 </div><div>              ICNTL(4) (level of printing):        3 </div><div>              ICNTL(5) (input mat struct):         0 </div><div>              ICNTL(6) (matrix prescaling):        7 </div><div>              ICNTL(7) (sequential matrix ordering):2 </div><div>              ICNTL(8) (scaling strategy):        77 </div><div>              ICNTL(10) (max num of refinements):  0 </div><div>              ICNTL(11) (error analysis):          0 </div><div>              ICNTL(12) (efficiency control):                         0 </div><div>              ICNTL(13) (efficiency control):                         1 </div><div>              ICNTL(14) (percentage of estimated workspace increase): 20 </div><div>              ICNTL(18) (input mat struct):                           0 </div><div>              ICNTL(19) (Schur complement info):                      0 </div><div>              ICNTL(20) (rhs sparse pattern):                         0 </div><div>              ICNTL(21) (solution struct):                            0 </div><div>              ICNTL(22) (in-core/out-of-core facility):               0 </div><div>              ICNTL(23) (max size of memory can be allocated locally):0 </div><div>              ICNTL(24) (detection of null pivot rows):               1 </div><div>              ICNTL(25) (computation of a null space basis):          0 </div><div>              ICNTL(26) (Schur options for rhs or solution):          0 </div><div>              ICNTL(27) (experimental parameter):                     -32 </div><div>              ICNTL(28) (use parallel or sequential ordering):        1 </div><div>              ICNTL(29) (parallel ordering):                          0 </div><div>              ICNTL(30) (user-specified set of entries in inv(A)):    0 </div><div>              ICNTL(31) (factors is discarded in the solve phase):    0 </div><div>              ICNTL(33) (compute determinant):                        0 </div><div>              ICNTL(35) (activate BLR based factorization):           0 </div><div>              ICNTL(36) (choice of BLR factorization variant):        0 </div><div>              ICNTL(38) (estimated compression rate of LU factors):   333 </div><div>              CNTL(1) (relative pivoting threshold):      0.01 </div><div>              CNTL(2) (stopping criterion of refinement): 1.49012e-08 </div><div>              CNTL(3) (absolute pivoting threshold):      0. </div><div>              CNTL(4) (value of static pivoting):         -1. </div><div>              CNTL(5) (fixation for null pivots):         0. </div><div>              CNTL(7) (dropping parameter for BLR):       0. </div><div>              RINFO(1) (local estimated flops for the elimination after analysis): </div><div>                [0] 71368. </div><div>              RINFO(2) (local estimated flops for the assembly after factorization): </div><div>                [0]  5976. </div><div>              RINFO(3) (local estimated flops for the elimination after factorization): </div><div>                [0]  119716. </div><div>              INFO(15) (estimated size of (in MB) MUMPS internal data for running numerical factorization): </div><div>              [0] 0 </div><div>              INFO(16) (size of (in MB) MUMPS internal data used during numerical factorization): </div><div>                [0] 0 </div><div>              INFO(23) (num of pivots eliminated on this processor after factorization): </div><div>                [0] 240 </div><div>              RINFOG(1) (global estimated flops for the elimination after analysis): 71368. </div><div>              RINFOG(2) (global estimated flops for the assembly after factorization): 5976. </div><div>              RINFOG(3) (global estimated flops for the elimination after factorization): 119716. </div><div>              (RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0.,0.)*(2^0)</div><div>              INFOG(3) (estimated real workspace for factors on all processors after analysis): 4641 </div><div>              INFOG(4) (estimated integer workspace for factors on all processors after analysis): 2816 </div><div>              INFOG(5) (estimated maximum front size in the complete tree): 38 </div><div>              INFOG(6) (number of nodes in the complete tree): 56 </div><div>              INFOG(7) (ordering option effectively use after analysis): 2 </div><div>              INFOG(8) (structural symmetry in percent of the permuted matrix after analysis): 100 </div><div>              INFOG(9) (total real/complex workspace to store the matrix factors after factorization): 6193 </div><div>              INFOG(10) (total integer space store the matrix factors after factorization): 3036 </div><div>              INFOG(11) (order of largest frontal matrix after factorization): 42 </div><div>              INFOG(12) (number of off-diagonal pivots): 79 </div><div>              INFOG(13) (number of delayed pivots after factorization): 110 </div><div>              INFOG(14) (number of memory compress after factorization): 0 </div><div>              INFOG(15) (number of steps of iterative refinement after solution): 0 </div><div>              INFOG(16) (estimated size (in MB) of all MUMPS internal data for factorization after analysis: value on the most memory consuming processor): 0 </div><div>              INFOG(17) (estimated size of all MUMPS internal data for factorization after analysis: sum over all processors): 0 </div><div>              INFOG(18) (size of all MUMPS internal data allocated during factorization: value on the most memory consuming processor): 0 </div><div>              INFOG(19) (size of all MUMPS internal data allocated during factorization: sum over all processors): 0 </div><div>              INFOG(20) (estimated number of entries in the factors): 3750 </div><div>              INFOG(21) (size in MB of memory effectively used during factorization - value on the most memory consuming processor): 0 </div><div>              INFOG(22) (size in MB of memory effectively used during factorization - sum over all processors): 0 </div><div>              INFOG(23) (after analysis: value of ICNTL(6) effectively used): 5 </div><div>              INFOG(24) (after analysis: value of ICNTL(12) effectively used): 3 </div><div>              INFOG(25) (after factorization: number of pivots modified by static pivoting): 0 </div><div>              INFOG(28) (after factorization: number of null pivots encountered): 0</div><div>              INFOG(29) (after factorization: effective number of entries in the factors (sum over all processors)): 4896</div><div>              INFOG(30, 31) (after solution: size in Mbytes of memory used during solution phase): 0, 0</div><div>              INFOG(32) (after analysis: type of analysis done): 1</div><div>              INFOG(33) (value used for ICNTL(8)): -2</div><div>              INFOG(34) (exponent of the determinant if determinant is requested): 0</div><div>              INFOG(35) (after factorization: number of entries taking into account BLR factor compression - sum over all processors): 4896</div><div>              INFOG(36) (after analysis: estimated size of all MUMPS internal data for running BLR in-core - value on the most memory consuming processor): 0 </div><div>              INFOG(37) (after analysis: estimated size of all MUMPS internal data for running BLR in-core - sum over all processors): 0 </div><div>              INFOG(38) (after analysis: estimated size of all MUMPS internal data for running BLR out-of-core - value on the most memory consuming processor): 0 </div><div>              INFOG(39) (after analysis: estimated size of all MUMPS internal data for running BLR out-of-core - sum over all processors): 0 </div><div>  linear system matrix = precond matrix:</div><div>  Mat Object: 1 MPI processes</div><div>    type: seqaij</div><div>    rows=240, cols=240</div><div>    total: nonzeros=4656, allocated nonzeros=4656</div><div>    total number of mallocs used during MatSetValues calls=0</div><div>      using I-node routines: found 167 nodes, limit used is 5</div><br><div>Entering DMUMPS 5.2.1 from C interface with JOB =  -2</div><div>      executing #MPI =      1, without OMP</div><div>rank: 0 coefficient: 0.132368</div><br><div>Entering DMUMPS 5.2.1 from C interface with JOB, N, NNZ =   1         960           9792</div><div>      executing #MPI =      1, without OMP</div><br><div> =================================================</div><div> MUMPS compiled with option -Dmetis</div><div> MUMPS compiled with option -Dptscotch</div><div> MUMPS compiled with option -Dscotch</div><div> This MUMPS version includes code for SAVE_RESTORE</div><div> =================================================</div><div>L D L^T Solver for general symmetric matrices</div><div>Type of parallelism: Working host</div><br><div> ****** ANALYSIS STEP ********</div><br><div> Scaling will be computed during analysis</div><div>Compute maximum matching (Maximum Transversal):  5</div><div> ... JOB =  5: MAXIMIZE PRODUCT DIAGONAL AND SCALE</div><br><div>Entering analysis phase with ...</div><div>                N        NNZ         LIW       INFO(1)</div><div>                960       9792       20545             0</div><div>Matrix entries:    IRN()   ICN()</div><div>           1      1           1      2           1      3</div><div>           1      4           1      5           1      6</div><div>           1      7           1      8           1      9</div><div>           1     10</div><div> Average density of rows/columns =   18</div><div> Average density of rows/columns =   18</div><div> Ordering based on AMF </div><div> Constrained Ordering based on AMF</div><div> Average density of rows/columns =   18</div><div> Average density of rows/columns =   18</div><div>NFSIZ(.)  =     0     0     0    58     0     0     0    73    14     0</div><br><div>FILS (.)  =     0  -747   -80   922   146     5     6   669     3     1</div><br><div>FRERE(.)  =   961   961   961     0   961   961   961    -4   -69   961</div><br><br><div>Leaving analysis phase with  ...</div><div> INFOG(1)                                       =               0</div><div> INFOG(2)                                       =               0</div><div> -- (20) Number of entries in factors (estim.)  =           20336</div><div> --  (3) Real space for factors    (estimated)  =           24094</div><div> --  (4) Integer space for factors (estimated)  =           12143</div><div> --  (5) Maximum frontal size      (estimated)  =              80</div><div> --  (6) Number of nodes in the tree            =             227</div><div> -- (32) Type of analysis effectively used      =               1</div><div> --  (7) Ordering option effectively used       =               2</div><div> ICNTL(6) Maximum transversal option            =               0</div><div> ICNTL(7) Pivot order option                    =               2</div><div> ICNTL(14) Percentage of memory relaxation      =              20</div><div> Number of level 2 nodes                        =               0</div><div> Number of split nodes                          =               0</div><div> RINFOG(1) Operations during elimination (estim)= 6.966D+05</div><div> Ordering compressed/constrained (ICNTL(12))    =               3</div><br><div> MEMORY ESTIMATIONS ... </div><div> Estimations with standard Full-Rank (FR) factorization:</div><div>    Total space in MBytes, IC factorization      (INFOG(17)):           1</div><div>    Total space in MBytes,  OOC factorization    (INFOG(27)):           1</div><br><div> Elapsed time in analysis driver=       0.0066</div><br><div>Entering DMUMPS 5.2.1 from C interface with JOB, N, NNZ =   2         960           9792</div><div>      executing #MPI =      1, without OMP</div><br><br><br><div>****** FACTORIZATION STEP ********</div><br><div> GLOBAL STATISTICS PRIOR NUMERICAL FACTORIZATION ...</div><div> Number of working processes                =               1</div><div> ICNTL(22) Out-of-core option               =               0</div><div> ICNTL(35) BLR activation (eff. choice)     =               0</div><div> ICNTL(14) Memory relaxation                =              20</div><div> INFOG(3) Real space for factors (estimated)=           24094</div><div> INFOG(4) Integer space for factors (estim.)=           12143</div><div> Maximum frontal size (estimated)           =              80</div><div> Number of nodes in the tree                =             227</div><div> Memory allowed (MB -- 0: N/A )             =               0</div><div> Memory provided by user, sum of LWK_USER   =               0</div><div> Relative threshold for pivoting, CNTL(1)   =      0.1000D-01</div><div>  ZERO PIVOT DETECTION ON, THRESHOLD          =   2.9434468577175697E-020</div><div> INFINITE FIXATION </div><div> Effective size of S     (based on INFO(39))=                31314</div><div> Elapsed time to reformat/distribute matrix =      0.0006</div><div> ** Memory allocated, total in Mbytes           (INFOG(19)):           1</div><div> ** Memory effectively used, total in Mbytes    (INFOG(22)):           1</div><div> ** Memory dynamically allocated for CB, total in Mbytes   :           0</div><br><div> Elapsed time for (failed) factorization    =      0.0032</div><br><div>Leaving factorization with ...</div><div> RINFOG(2)  Operations in node assembly     = 3.366D+04</div><div> ------(3)  Operations in node elimination  = 9.346D+05</div><div> INFOG (9)  Real space for factors          =           26980</div><div> INFOG(10)  Integer space for factors       =           13047</div><div> INFOG(11)  Maximum front size              =              84</div><div> INFOG(29)  Number of entries in factors    =           24047</div><div> INFOG(12)  Number of negative pivots       =             294</div><div> INFOG(13)  Number of delayed pivots        =             452</div><div> Number of 2x2 pivots in type 1 nodes       =               0</div><div> Number of 2X2 pivots in type 2 nodes       =               0</div><div> Nb of null pivots detected by ICNTL(24)    =               0</div><div> INFOG(28)  Estimated deficiency            =               0</div><div> INFOG(14)  Number of memory compress       =               1</div><br><div> Elapsed time in factorization driver=       0.0042</div><div> On return from DMUMPS, INFOG(1)=              -9</div><div> On return from DMUMPS, INFOG(2)=              22</div><div>terminate called after throwing an instance of 'std::runtime_error'</div><div>  what():  </div><br><div>*** -------------------------------------------------------------------------</div><div>*** DOLFIN encountered an error. If you are not able to resolve this issue</div><div>*** using the information listed below, you can ask for help at</div><div>***</div><div>***     <a href="mailto:fenics-support@googlegroups.com" target="_blank">fenics-support@googlegroups.com</a></div><div>***</div><div>*** Remember to include the error message listed below and, if possible,</div><div>*** include a *minimal* running example to reproduce the error.</div><div>***</div><div>*** -------------------------------------------------------------------------</div><div>*** Error:   Unable to solve linear system using PETSc Krylov solver.</div><div>*** Reason:  Solution failed to converge in 0 iterations (PETSc reason DIVERGED_PC_FAILED, residual norm ||r|| = 0.000000e+00).</div><div>*** Where:   This error was encountered inside PETScKrylovSolver.cpp.</div><div>*** Process: 0</div><div>*** </div><div>*** DOLFIN version: 2019.1.0</div><div>*** Git changeset:  74d7efe1e84d65e9433fd96c50f1d278fa3e3f3f</div><div>*** -------------------------------------------------------------------------</div><br><div>Aborted (core dumped)</div><br><img alt="Sent from Mailspring" width="0" height="0" style="border: 0px; width: 0px; height: 0px;" src="https://link.getmailspring.com/open/B6CC5A9F-3387-4A0A-BFBF-F5F0B322C09C@getmailspring.com?me=05842e25&recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D"></blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>