<div dir="ltr"><div class="gmail_quote"><div dir="ltr">On Tue, Oct 30, 2018 at 1:18 PM Santiago Andres Triana via petsc-users <<a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Hi petsc-users,<div><br></div><div>I am solving a generalized eigenvalue problem using ex7 in $SLEPC_DIR/src/eps/examples/tutorials/. I provide the A and B matrices.</div><div>The program runs fine, with correct solutions on 12-core node and also on a mac laptop.</div><div><br></div><div>However, on a 16-core workstation running Debian testing (fresh install) and also a fresh install of petsc and slepc I get the following error:</div><div><br></div><div> $ mpiexec -n 2 ./ex7 -f1 A.petsc -f2 B.petsc -st_type sinvert -eps_nev 4 -eps_target -2e-3+1.01i</div><div><br></div><div>Generalized eigenproblem stored in file.</div><div><br></div><div> Reading COMPLEX matrices from binary files...</div><div>[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range</div><div>[1]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger</div><div>[1]PETSC ERROR: or see <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind</a></div><div>[1]PETSC ERROR: or try <a href="http://valgrind.org" target="_blank">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors</div><div>[1]PETSC ERROR: likely location of problem given in stack below</div><div>[1]PETSC ERROR: --------------------- Stack Frames ------------------------------------</div><div>[1]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,</div><div>[1]PETSC ERROR: INSTEAD the line number of the start of the function</div><div>[1]PETSC ERROR: is given.</div><div>[1]PETSC ERROR: [1] MatFactorNumeric_MUMPS line 1205 /home/spin2/petsc-3.10.2/src/mat/impls/aij/mpi/mumps/mumps.c</div><div>[1]PETSC ERROR: [1] MatLUFactorNumeric line 3054 /home/spin2/petsc-3.10.2/src/mat/interface/matrix.c</div><div>[1]PETSC ERROR: [1] PCSetUp_LU line 59 /home/spin2/petsc-3.10.2/src/ksp/pc/impls/factor/lu/lu.c</div><div>[1]PETSC ERROR: [1] PCSetUp line 894 /home/spin2/petsc-3.10.2/src/ksp/pc/interface/precon.c</div><div>[1]PETSC ERROR: [1] KSPSetUp line 304 /home/spin2/petsc-3.10.2/src/ksp/ksp/interface/itfunc.c</div><div>[1]PETSC ERROR: [1] STSetUp_Sinvert line 96 /home/spin2/slepc-3.10.1/src/sys/classes/st/impls/sinvert/sinvert.c</div><div>[1]PETSC ERROR: [1] STSetUp line 233 /home/spin2/slepc-3.10.1/src/sys/classes/st/interface/stsolve.c</div><div>[1]PETSC ERROR: [1] EPSSetUp line 104 /home/spin2/slepc-3.10.1/src/eps/interface/epssetup.c</div><div>[1]PETSC ERROR: [1] EPSSolve line 129 /home/spin2/slepc-3.10.1/src/eps/interface/epssolve.c</div><div>[1]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------</div><div>[1]PETSC ERROR: Signal received</div><div>[1]PETSC ERROR: See <a href="http://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">http://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.</div><div>[1]PETSC ERROR: Petsc Release Version 3.10.2, Oct, 09, 2018 </div><div>[1]PETSC ERROR: ./ex7 on a arch-linux2-c-opt named wobble-wkst-as by spin2 Tue Oct 30 17:40:51 2018</div><div>[1]PETSC ERROR: Configure options --download-mpich -with-scalar-type=complex --download-mumps --download-parmetis --download-metis --download-scalapack --download-fblaslapack --with-debugging=1 --download-superlu_dist --download-ptscotch</div><div>[1]PETSC ERROR: #1 User provided function() line 0 in unknown file</div><div>application called MPI_Abort(MPI_COMM_WORLD, 59) - process 1</div><div><br></div><div><br></div><div><br></div><div>the expected output is the following (on a compute node running petsc-3.9.2 and also on a mac laptop running petsc-3.10.2):</div><div><br></div><div><div>$ mpiexec -n 2 ./ex7 -f1 A.petsc -f2 B.petsc -st_type sinvert -eps_nev 4 -eps_target -2e-3+1.01i</div><div><br></div><div>Generalized eigenproblem stored in file.</div><div><br></div><div> Reading COMPLEX matrices from binary files...</div><div> Number of iterations of the method: 2</div><div> Number of linear iterations of the method: 27</div><div> Solution method: krylovschur</div><div><br></div><div> Number of requested eigenvalues: 4</div><div> Stopping condition: tol=1e-08, maxit=63157</div><div> Linear eigensolve converged (4 eigenpairs) due to CONVERGED_TOL; iterations 2</div><div> ---------------------- --------------------</div><div> k ||Ax-kBx||/||kx||</div><div> ---------------------- --------------------</div><div> -0.002806+1.009827i 2.00821e-19</div><div> -0.002980+1.008417i 8.08359e-17</div><div> -0.002676+1.011755i 9.49342e-18</div><div> -0.003201+1.007367i 1.50869e-16</div><div> ---------------------- --------------------</div></div><div><br></div><div><br></div><div>Just in case, the matrices can be downloaded from here if any one wants to give them a try</div><div><a href="https://www.dropbox.com/s/ejpa9owkv8tjnwi/A.petsc?dl=0" target="_blank">https://www.dropbox.com/s/ejpa9owkv8tjnwi/A.petsc?dl=0</a><br></div><div><a href="https://www.dropbox.com/s/urjtxaezl0cv3om/B.petsc?dl=0" target="_blank">https://www.dropbox.com/s/urjtxaezl0cv3om/B.petsc?dl=0</a></div></div></div></div></div></div></blockquote><div><br></div><div>Its not the matrices. Would you be willing to reconfigure and test with SuperLU_dist? Since it consistently fails</div><div>in the MUMPS factorization, it seems like it is either a bug in MUMPS or a bug in our interface to MUMPS. Valgrind</div><div>should show us which one (as Barry suggests), but running with SuperLU_dist should get you going in the mean time.</div><div><br></div><div> Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div>I tried different petsc/slepc versions to no avail, including an OS reinstall. So any help would be highly appreciated. Thanks in advance!</div><div><br></div><div>Santiago</div></div></div></div></div></div>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener</div><div><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>