<div dir="ltr">Hi everyone,<div><br></div><div>I've run into a strange problem in my Fortran 90 code where it runs fine with 1 processor, but then throws a segmentation fault on KSPSolve() when I try to run it in parallel.  I'm using PETSc 3.3 with the SuperLU direct solver for the sequential case and SuperLU_dist for the parallel case.  I've called KSPView before and after KSPSolve.  I'll put the KSPView output for the sequential and parallel cases and the crash info for the parallel case below (with some details of my system redacted).  Any help would be appreciated.  If you need any other information, I'm happy to provide it.</div>
<div><br></div><div>Thank you,</div><div><br></div><div style>~Brendan</div><div>------------------------------</div><div><br></div><div>KSPView() before sequential solve:  </div><div><br></div><div>KSP Object: 1 MPI processes</div>
<div>  type: preonly</div><div>  maximum iterations=10000, initial guess is zero</div><div>  tolerances:  relative=1e-05, absolute=1e-50, divergence=10000</div>
<div>  left preconditioning</div><div>  using DEFAULT norm type for convergence test</div><div>PC Object: 1 MPI processes</div><div>  type: lu</div><div>    LU: out-of-place factorization</div><div>    tolerance for zero pivot 2.22045e-14</div>

<div>    matrix ordering: nd</div><div>  linear system matrix = precond matrix:</div><div>  Matrix Object:   1 MPI processes</div><div>    type: seqaij</div><div>    rows=11760, cols=11760</div><div>    total: nonzeros=506586, allocated nonzeros=509061</div>

<div>    total number of mallocs used during MatSetValues calls =0</div><div>      not using I-node routines</div><div><br></div><div>KSPView() after sequential solve:</div><div>  <br></div><div><div>KSP Object: 1 MPI processes</div>
<div>  type: preonly</div>
<div>  maximum iterations=10000, initial guess is zero</div><div>  tolerances:  relative=1e-05, absolute=1e-50, divergence=10000</div><div>  left preconditioning</div><div>  using NONE norm type for convergence test</div>

<div>PC Object: 1 MPI processes</div><div>  type: lu</div><div>    LU: out-of-place factorization</div><div>    tolerance for zero pivot 2.22045e-14</div><div>    matrix ordering: nd</div><div>    factor fill ratio given 0, needed 0</div>

<div>      Factored matrix follows:</div><div>        Matrix Object:         1 MPI processes</div><div>          type: seqaij</div><div>          rows=11760, cols=11760</div><div>          package used to perform factorization: superlu</div>

<div>          total: nonzeros=0, allocated nonzeros=0</div><div>          total number of mallocs used during MatSetValues calls =0</div><div>            SuperLU run parameters:</div><div>              Equil: NO</div><div>

              ColPerm: 3</div><div>              IterRefine: 0</div><div>              SymmetricMode: NO</div><div>              DiagPivotThresh: 1</div><div>              PivotGrowth: NO</div><div>              ConditionNumber: NO</div>

<div>              RowPerm: 0</div><div>              ReplaceTinyPivot: NO</div><div>              PrintStat: NO</div><div>              lwork: 0</div><div>  linear system matrix = precond matrix:</div><div>  Matrix Object:   1 MPI processes</div>

<div>    type: seqaij</div><div>    rows=11760, cols=11760</div><div>    total: nonzeros=506586, allocated nonzeros=509061</div><div>    total number of mallocs used during MatSetValues calls =0</div><div>      not using I-node routines</div>

<div><br></div><div><br></div><div>KSPView() before parallel solve:  <br></div><div><br></div><div><div>KSP Object: 2 MPI processes</div><div>  type: preonly</div><div>  maximum iterations=10000, initial guess is zero</div>
<div>  tolerances:  relative=1e-05, absolute=1e-50, divergence=10000</div><div>  left preconditioning</div><div>  using DEFAULT norm type for convergence test</div><div>PC Object: 2 MPI processes</div><div>  type: lu</div>
<div>    LU: out-of-place factorization</div><div>    tolerance for zero pivot 2.22045e-14</div><div>    matrix ordering: natural</div><div>  linear system matrix = precond matrix:</div><div>      Solving Electron Matrix Equation</div>
<div>  Matrix Object:   2 MPI processes</div><div>    type: mpiaij</div><div>    rows=11760, cols=11760</div><div>    total: nonzeros=506586, allocated nonzeros=520821</div><div>    total number of mallocs used during MatSetValues calls =0</div>
<div>      not using I-node (on process 0) routines</div></div><div><br></div><div style>Crash info for parallel solve:</div><div style><br></div><div style><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[1]PETSC">http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind[1]PETSC</a> ERROR: or try <a href="http://valgrind.org">http://valgrind.org</a> on GNU/linux and Apple Mac OS X to find memory corruption errors</div>
<div>[1]PETSC ERROR: configure using --with-debugging=yes, recompile, link, and run </div><div>[1]PETSC ERROR: to get more information on the crash.</div><div>[1]PETSC ERROR: --------------------- Error Message ------------------------------------</div>
<div>[1]PETSC ERROR: Signal received!</div><div>[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: Petsc Release Version 3.3.0, Patch 6, Mon Feb 11 12:26:34 CST 2013 </div>
<div>[1]PETSC ERROR: See docs/changes/index.html for recent updates.</div><div>[1]PETSC ERROR: See docs/faq.html for hints about trouble shooting.</div><div>[1]PETSC ERROR: See docs/index.html for manual pages.</div><div>
[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: <redacted> on a path-ompi named <redacted></div><div>[1]PETSC ERROR: Libraries linked from <redacted></div>
<div>[1]PETSC ERROR: Configure run at Thu Mar 21 14:19:42 2013</div><div>[1]PETSC ERROR: Configure options --PETSC_ARCH=path-ompi --PETSC_DIR=<redacted> --CFLAGS="-fPIC -O -mp" --CXXFLAGS="-fPIC -O -mp" --FFLAGS="-fPIC -O -mp" --with-debugging=0 --with-dynamic-loadin=no --with-mpi=1 --with-mpi-dir=<redacted> --with-superlu=1 --with-superlu-dir=<redacted> --with-blas-lapack-lib="<redacted>" --with-scalapack=1 --with-scalapack-dir=<redacted> --with-superlu_dist=1 --with-superlu_dist-dir=<redacted> --with-metis=1 --with-metis-dir=<redacted> --with-parmetis=1 --with-parmetis-dir=<redacted> --with-blacs-lib="<redacted>" --with-blacs-include=<redacted> --with-hypre=1 --download-hypre=1</div>
<div>[1]PETSC ERROR: ------------------------------------------------------------------------</div><div>[1]PETSC ERROR: User provided function() line 0 in unknown directory unknown file</div><div><br></div></div><div><br>
</div><div><br></div></div></div>