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