<div dir="ltr">Hi All,<div><br></div><div>I am using SLEPc to compute eigenvalues and MUMPS for factorization. </div><div><br></div><div>Please find attached:</div><div>1) A simple test program <font face="monospace">slepc_eps_mumps_test.cpp</font> that reads a given PETSc matrix and computes the smallest eigenvalues using MUMPS for factorization</div><div>2) <font face="arial, sans-serif">An example PETSc matrix </font><font face="monospace">MatA </font><font face="arial, sans-serif">of size 581343 rows (sending in .gz format via Google drive link, please extract it "</font><font face="monospace">gunzip MatA.gz</font><font face="arial, sans-serif">" before using). You should be able to reproduce this issue with any other matrix of a similar or bigger size.</font></div><div><br></div><div>I notice that when I run the attached test program in parallel with the attached test matrix the MUMPS statistics printed (using the -eps_view option in the command line) change with every run.</div><div><br></div><div>This is how I run the test:</div><div><font face="monospace"><b>$ </b>mpiexec -n 24 ./slepc_eps_mumps_test.o -nev 5 -f MatA -eps_view</font></div><div>and for example, the output of this includes the following MUMPS stats</div><div><font face="monospace">...</font></div><div><font face="monospace">PC Object: (st_) 24 MPI processes<br>    type: lu<br>      out-of-place factorization<br>      tolerance for zero pivot 2.22045e-14<br>      matrix ordering: external<br>      factor fill ratio given 0., needed 0.<br>        Factored matrix follows:<br>          Mat Object: 24 MPI processes<br>            type: mumps<br>            rows=581343, cols=581343<br>            package used to perform factorization: mumps<br>            <b>total: nonzeros=348236349, allocated nonzeros=348236349</b><br>              MUMPS run parameters:<br>                SYM (matrix type):                   0</font><br></div><div><font face="monospace">...</font></div><div><font face="monospace"><br></font></div><div><font face="arial, sans-serif">I ran this test 10 times as follows and got a different number of nonzeros (line highlighted above ) reported in each run. (If you save the full output and compare, you will notice many other differences, but I wouldn't have expected the nonzeros to change with every run.)</font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="monospace"><b>$</b> for i in `seq 1 10`; do echo "run $i :-----"; mpiexec -n 24 ./slepc_eps_mumps_test.o -nev 5 -f MatA -eps_view | grep -A 1 "factorization: mumps"; done</font></div><div><font face="monospace"><br></font></div><div><font face="monospace">run 1 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354789915, allocated nonzeros=354789915<br>run 2 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=359811101, allocated nonzeros=359811101<br>run 3 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354834871, allocated nonzeros=354834871<br>run 4 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354830397, allocated nonzeros=354830397<br>run 5 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=353942929, allocated nonzeros=353942929<br>run 6 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354147241, allocated nonzeros=354147241<br>run 7 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354980083, allocated nonzeros=354980083<br>run 8 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354980083, allocated nonzeros=354980083<br>run 9 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=354214219, allocated nonzeros=354214219<br>run 10 :-----<br>            package used to perform factorization: mumps<br>            total: nonzeros=355894047, allocated nonzeros=355894047<br></font></div><div><font face="monospace"><br></font></div><div><font face="arial, sans-serif">Can somebody please explain what causes these differences in MUMPS stats? </font></div><div><font face="arial, sans-serif"><br></font></div><div><font face="arial, sans-serif">Thanks,</font></div><div><font face="arial, sans-serif">Varun</font></div><div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:rgb(245,245,245);padding:5px;color:rgb(34,34,34);font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid rgb(221,221,221);line-height:1"><a href="https://drive.google.com/file/d/1gbi211awBK5nZtZ-dFiTmjBwStJQ_FKR/view?usp=drive_web" target="_blank" style="display:inline-block;max-width:366px;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0;border:none" aria-label="MatA.gz"><img style="vertical-align: bottom; border: none;" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="color:rgb(17,85,204);text-decoration:none;vertical-align:bottom">MatA.gz</span></a><img src="//ssl.gstatic.com/ui/v1/icons/common/x_8px.png" style="opacity: 0.55; cursor: pointer; float: right; position: relative; top: -1px; display: none;"></div><div><font face="monospace"><b><br></b></font></div></div>