Dear Users,<div><br></div><div>I&#39;m having trouble figuring out why the MUMPS solver is failing on a specific range of one of my parameters. When using the PETSc direct solver on a single processor I have no issues. There error is:</div>
<div><br></div><div><div>[0]PETSC ERROR: --------------------- Error Message ------------------------------------</div><div>[0]PETSC ERROR: Error in external library!</div><div>[0]PETSC ERROR: Error reported by MUMPS in numerical factorization phase: INFO(1)=-9, INFO(2)=13</div>
<div>!</div><div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48 CDT 2011</div><div>[0]PETSC ERROR: See docs/changes/index.html for recent updates.</div>
<div>[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.</div><div>[0]PETSC ERROR: See docs/index.html for manual pages.</div><div>[0]PETSC ERROR: ------------------------------------------------------------------------</div>
<div>[0]PETSC ERROR: ./cntor on a complex-c named hpc-1-14.local by abyrd Wed Aug 31 10:53:42 2011</div><div>[0]PETSC ERROR: Libraries linked from /panfs/storage.local/scs/home/abyrd/petsc-3.1-p8/complex-cpp-mumps/lib</div>
<div>[0]PETSC ERROR: Configure run at Mon Jul 11 15:28:42 2011</div><div>[0]PETSC ERROR: Configure options PETSC_ARCH=complex-cpp-mumps --with-cc=mpicc --with-fc=mpif90 --with-blas-lapack-dir=/usr/lib64 --with-shared --with-clanguage=c++ --with-scalar-type=complex --download-mumps=1 --download-blacs=1 --download-scalapack=1 --download-parmetis=1 --with-cxx=mpicxx</div>
<div>[0]PETSC ERROR: ------------------------------------------------------------------------</div><div>[0]PETSC ERROR: MatFactorNumeric_MUMPS() line 517 in src/mat/impls/aij/mpi/mumps/mumps.c</div><div>[0]PETSC ERROR: MatLUFactorNumeric() line 2587 in src/mat/interface/matrix.c</div>
<div>[0]PETSC ERROR: PCSetUp_LU() line 158 in src/ksp/pc/impls/factor/lu/lu.c</div><div>[0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c</div><div>[0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c</div>
<div>[0]PETSC ERROR: InvertHamiltonian() line 102 in WDinvert.h</div></div><div><br></div><div><br></div><div>I suspect it has something to do with the preconditioning or setup of the matrix I am trying to invert. The matrix becomes singular at energy = 0 eV, and is nearly singular for values close to that, but the code is failing on energies relatively far from that point. The affected energy interval is [<span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">-0.03095, 0.03095].</span></div>
<div><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; "><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">Is anyone able to point in the right direction to figure out what I&#39;m not setting up properly?</span></div>
<div><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; "><br></span></div><div><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">Respectfully,</span></div>
<div><span class="Apple-style-span" style="border-collapse: collapse; font-family: arial, sans-serif; font-size: 13px; ">Adam Byrd</span></div>