[petsc-users] Mumps Error

[Adam Byrd]

Dear Users,

I'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:

[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Error in external library!
[0]PETSC ERROR: Error reported by MUMPS in numerical factorization phase:
INFO(1)=-9, INFO(2)=13
!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Release Version 3.1.0, Patch 8, Thu Mar 17 13:37:48
CDT 2011
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: ./cntor on a complex-c named hpc-1-14.local by abyrd Wed Aug
31 10:53:42 2011
[0]PETSC ERROR: Libraries linked from
/panfs/storage.local/scs/home/abyrd/petsc-3.1-p8/complex-cpp-mumps/lib
[0]PETSC ERROR: Configure run at Mon Jul 11 15:28:42 2011
[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
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: MatFactorNumeric_MUMPS() line 517 in
src/mat/impls/aij/mpi/mumps/mumps.c
[0]PETSC ERROR: MatLUFactorNumeric() line 2587 in src/mat/interface/matrix.c
[0]PETSC ERROR: PCSetUp_LU() line 158 in src/ksp/pc/impls/factor/lu/lu.c
[0]PETSC ERROR: PCSetUp() line 795 in src/ksp/pc/interface/precon.c
[0]PETSC ERROR: KSPSetUp() line 237 in src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: InvertHamiltonian() line 102 in WDinvert.h


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 [-0.03095,
0.03095].

Is anyone able to point in the right direction to figure out what I'm not
setting up properly?

Respectfully,
Adam Byrd
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110831/c0eb5e74/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: PETScCntor.zip
Type: application/zip
Size: 80941 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20110831/c0eb5e74/attachment-0001.zip>