<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=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 23 May 2021, at 6:17 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com" class="">knepley@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><meta charset="UTF-8" class=""><div dir="ltr" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><div dir="ltr" class="">On Sun, May 23, 2021 at 11:17 AM Karl Yang <<a href="mailto:y.juntao@hotmail.com" class="">y.juntao@hotmail.com</a>> wrote:<br class=""></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><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></blockquote><div class=""><br class=""></div><div class="">This says</div><div class=""><br class=""></div><div class="">On return from DMUMPS, INFOG(1)= -9<br class=""> On return from DMUMPS, INFOG(2)= 22<br class=""></div><div class=""><br class=""></div><div class="">that the internal work array for MUMPS is too small. I am not sure which option controls that.</div></div></div></div></blockquote><div><br class=""></div><div>-mat_mumps_icntl_14</div><div><br class=""></div><div>Juntao, you can usually troubleshoot MUMPS error codes by looking at <a href="http://mumps.enseeiht.fr/doc/userguide_5.4.0.pdf#page=92" class="">http://mumps.enseeiht.fr/doc/userguide_5.4.0.pdf#page=92</a></div><div><br class=""></div><div>Thanks,</div><div>Pierre</div><br class=""><blockquote type="cite" class=""><div class=""><div dir="ltr" style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;" class=""><div class="gmail_quote"><div class=""> Thanks,</div><div class=""><br class=""></div><div class=""> Matt</div><div class=""> </div><blockquote class="gmail_quote" style="margin: 0px 0px 0px 0.8ex; border-left-width: 1px; border-left-style: solid; border-left-color: rgb(204, 204, 204); padding-left: 1ex;"><div class="">regards</div><div class="">Juntao</div><br class=""><div class="">MUMPS options<span 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=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>N NNZ LIW INFO(1)</div><div class=""> <span class="Apple-converted-space"> </span>240 2448 5137 0</div><div class="">Matrix entries: IRN() ICN()</div><div class=""> <span class="Apple-converted-space"> </span>1 1 1 2 1 3</div><div class=""> <span class="Apple-converted-space"> </span>1 4 1 5 1 6</div><div class=""> <span class="Apple-converted-space"> </span>1 7 1 8 1 9</div><div class=""> <span class="Apple-converted-space"> </span>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<span class="Apple-converted-space"> </span></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 ...<span class="Apple-converted-space"> </span></div><div class=""> Estimations with standard Full-Rank (FR) factorization:</div><div class=""> <span class="Apple-converted-space"> </span>Total space in MBytes, IC factorization (INFOG(17)): 0</div><div class=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>ZERO PIVOT DETECTION ON, THRESHOLD = 2.8931920285365730E-020</div><div class=""> INFINITE FIXATION<span class="Apple-converted-space"> </span></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=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>--- (10) = 0</div><div class=""> <span class="Apple-converted-space"> </span>--- (11) = 0</div><div class=""> <span class="Apple-converted-space"> </span>--- (20) = 0</div><div class=""> <span class="Apple-converted-space"> </span>--- (21) = 0</div><div class=""> <span class="Apple-converted-space"> </span>--- (30) = 0</div><div class=""> <span class="Apple-converted-space"> </span>--- (35) = 0</div><br class=""><br class=""><div class=""> Vector solution for column 1</div><div class=""> RHS</div><div class=""> <span class="Apple-converted-space"> </span>-7.828363D-02 -3.255337D+00 1.054729D+00 1.379822D-01 -3.892113D-01</div><div class=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>.. Time in forward (fwd) step = 0.000053</div><div class=""> <span class="Apple-converted-space"> </span>.. 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=""> <span class="Apple-converted-space"> </span>type: preonly</div><div class=""> <span class="Apple-converted-space"> </span>maximum iterations=10000, initial guess is zero</div><div class=""> <span class="Apple-converted-space"> </span>tolerances: relative=1e-05, absolute=1e-50, divergence=10000.</div><div class=""> <span class="Apple-converted-space"> </span>left preconditioning</div><div class=""> <span class="Apple-converted-space"> </span>using NONE norm type for convergence test</div><div class="">PC Object: 1 MPI processes</div><div class=""> <span class="Apple-converted-space"> </span>type: cholesky</div><div class=""> <span class="Apple-converted-space"> </span>out-of-place factorization</div><div class=""> <span class="Apple-converted-space"> </span>tolerance for zero pivot 2.22045e-14</div><div class=""> <span class="Apple-converted-space"> </span>matrix ordering: natural</div><div class=""> <span class="Apple-converted-space"> </span>factor fill ratio given 0., needed 0.</div><div class=""> <span class="Apple-converted-space"> </span>Factored matrix follows:</div><div class=""> <span class="Apple-converted-space"> </span>Mat Object: 1 MPI processes</div><div class=""> <span class="Apple-converted-space"> </span>type: mumps</div><div class=""> <span class="Apple-converted-space"> </span>rows=240, cols=240</div><div class=""> <span class="Apple-converted-space"> </span>package used to perform factorization: mumps</div><div class=""> <span class="Apple-converted-space"> </span>total: nonzeros=3750, allocated nonzeros=3750</div><div class=""> <span class="Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls=0</div><div class=""> <span class="Apple-converted-space"> </span>MUMPS run parameters:</div><div class=""> <span class="Apple-converted-space"> </span>SYM (matrix type): 2<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>PAR (host participation): 1<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(1) (output for error): 1<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(2) (output of diagnostic msg): 1<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(3) (output for global info): 6<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(4) (level of printing): 3<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(5) (input mat struct): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(6) (matrix prescaling): 7<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(7) (sequential matrix ordering):2<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(8) (scaling strategy): 77<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(10) (max num of refinements): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(11) (error analysis): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(12) (efficiency control): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(13) (efficiency control): 1<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(14) (percentage of estimated workspace increase): 20<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(18) (input mat struct): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(19) (Schur complement info): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(20) (rhs sparse pattern): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(21) (solution struct): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(22) (in-core/out-of-core facility): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(23) (max size of memory can be allocated locally):0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(24) (detection of null pivot rows): 1<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(25) (computation of a null space basis): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(26) (Schur options for rhs or solution): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(27) (experimental parameter): -32<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(28) (use parallel or sequential ordering): 1<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(29) (parallel ordering): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(30) (user-specified set of entries in inv(A)): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(31) (factors is discarded in the solve phase): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(33) (compute determinant): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(35) (activate BLR based factorization): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(36) (choice of BLR factorization variant): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>ICNTL(38) (estimated compression rate of LU factors): 333<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>CNTL(1) (relative pivoting threshold): 0.01<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>CNTL(2) (stopping criterion of refinement): 1.49012e-08<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>CNTL(3) (absolute pivoting threshold): 0.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>CNTL(4) (value of static pivoting): -1.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>CNTL(5) (fixation for null pivots): 0.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>CNTL(7) (dropping parameter for BLR): 0.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>RINFO(1) (local estimated flops for the elimination after analysis):<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>[0] 71368.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>RINFO(2) (local estimated flops for the assembly after factorization):<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>[0] 5976.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>RINFO(3) (local estimated flops for the elimination after factorization):<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>[0] 119716.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFO(15) (estimated size of (in MB) MUMPS internal data for running numerical factorization):<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>[0] 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFO(16) (size of (in MB) MUMPS internal data used during numerical factorization):<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>[0] 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFO(23) (num of pivots eliminated on this processor after factorization):<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>[0] 240<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>RINFOG(1) (global estimated flops for the elimination after analysis): 71368.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>RINFOG(2) (global estimated flops for the assembly after factorization): 5976.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>RINFOG(3) (global estimated flops for the elimination after factorization): 119716.<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>(RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0.,0.)*(2^0)</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(3) (estimated real workspace for factors on all processors after analysis): 4641<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(4) (estimated integer workspace for factors on all processors after analysis): 2816<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(5) (estimated maximum front size in the complete tree): 38<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(6) (number of nodes in the complete tree): 56<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(7) (ordering option effectively use after analysis): 2<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(8) (structural symmetry in percent of the permuted matrix after analysis): 100<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(9) (total real/complex workspace to store the matrix factors after factorization): 6193<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(10) (total integer space store the matrix factors after factorization): 3036<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(11) (order of largest frontal matrix after factorization): 42<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(12) (number of off-diagonal pivots): 79<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(13) (number of delayed pivots after factorization): 110<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(14) (number of memory compress after factorization): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(15) (number of steps of iterative refinement after solution): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(16) (estimated size (in MB) of all MUMPS internal data for factorization after analysis: value on the most memory consuming processor): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(17) (estimated size of all MUMPS internal data for factorization after analysis: sum over all processors): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(18) (size of all MUMPS internal data allocated during factorization: value on the most memory consuming processor): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(19) (size of all MUMPS internal data allocated during factorization: sum over all processors): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(20) (estimated number of entries in the factors): 3750<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(21) (size in MB of memory effectively used during factorization - value on the most memory consuming processor): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(22) (size in MB of memory effectively used during factorization - sum over all processors): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(23) (after analysis: value of ICNTL(6) effectively used): 5<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(24) (after analysis: value of ICNTL(12) effectively used): 3<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(25) (after factorization: number of pivots modified by static pivoting): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(28) (after factorization: number of null pivots encountered): 0</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(29) (after factorization: effective number of entries in the factors (sum over all processors)): 4896</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(30, 31) (after solution: size in Mbytes of memory used during solution phase): 0, 0</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(32) (after analysis: type of analysis done): 1</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(33) (value used for ICNTL(8)): -2</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(34) (exponent of the determinant if determinant is requested): 0</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(35) (after factorization: number of entries taking into account BLR factor compression - sum over all processors): 4896</div><div class=""> <span class="Apple-converted-space"> </span>INFOG(36) (after analysis: estimated size of all MUMPS internal data for running BLR in-core - value on the most memory consuming processor): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(37) (after analysis: estimated size of all MUMPS internal data for running BLR in-core - sum over all processors): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>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<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>INFOG(39) (after analysis: estimated size of all MUMPS internal data for running BLR out-of-core - sum over all processors): 0<span class="Apple-converted-space"> </span></div><div class=""> <span class="Apple-converted-space"> </span>linear system matrix = precond matrix:</div><div class=""> <span class="Apple-converted-space"> </span>Mat Object: 1 MPI processes</div><div class=""> <span class="Apple-converted-space"> </span>type: seqaij</div><div class=""> <span class="Apple-converted-space"> </span>rows=240, cols=240</div><div class=""> <span class="Apple-converted-space"> </span>total: nonzeros=4656, allocated nonzeros=4656</div><div class=""> <span class="Apple-converted-space"> </span>total number of mallocs used during MatSetValues calls=0</div><div class=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>N NNZ LIW INFO(1)</div><div class=""> <span class="Apple-converted-space"> </span>960 9792 20545 0</div><div class="">Matrix entries: IRN() ICN()</div><div class=""> <span class="Apple-converted-space"> </span>1 1 1 2 1 3</div><div class=""> <span class="Apple-converted-space"> </span>1 4 1 5 1 6</div><div class=""> <span class="Apple-converted-space"> </span>1 7 1 8 1 9</div><div class=""> <span class="Apple-converted-space"> </span>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<span class="Apple-converted-space"> </span></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 ...<span class="Apple-converted-space"> </span></div><div class=""> Estimations with standard Full-Rank (FR) factorization:</div><div class=""> <span class="Apple-converted-space"> </span>Total space in MBytes, IC factorization (INFOG(17)): 1</div><div class=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>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=""> <span class="Apple-converted-space"> </span>ZERO PIVOT DETECTION ON, THRESHOLD = 2.9434468577175697E-020</div><div class=""> INFINITE FIXATION<span class="Apple-converted-space"> </span></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=""> <span class="Apple-converted-space"> </span>what(): <span class="Apple-converted-space"> </span></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="">*** <span class="Apple-converted-space"> </span><a href="mailto:fenics-support@googlegroups.com" target="_blank" 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="">***<span class="Apple-converted-space"> </span></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 alt="Sent from Mailspring" width="0" height="0" src="https://link.getmailspring.com/open/B6CC5A9F-3387-4A0A-BFBF-F5F0B322C09C@getmailspring.com?me=05842e25&recipient=cGV0c2MtdXNlcnNAbWNzLmFubC5nb3Y%3D" style="border: 0px; width: 0px; height: 0px;" class=""></blockquote></div><br clear="all" class=""><div class=""><br class=""></div>--<span class="Apple-converted-space"> </span><br class=""><div dir="ltr" class="gmail_signature"><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class=""><div dir="ltr" class=""><div class="">What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br class="">-- Norbert Wiener</div><div class=""><br class=""></div><div class=""><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank" class="">https://www.cse.buffalo.edu/~knepley/</a></div></div></div></div></div></div></div></div></div></blockquote></div><br class=""></body></html>