Gabriele: This time, you get INFO(1)=-1, "An error occurred on processor INFO(2)=39", according to mumps manual. Please ask mumps developer about INFO(1)=-1. My guess is that your matrix might be numerical singular, causing numerical factorization fail.
Hong > > I am getting an error with a different number, but the message is the same: > > 59]PETSC ERROR: > ------------------------------------------------------------------------ > [59]PETSC ERROR: MatFactorNumeric_MUMPS() line 517 in > src/mat/impls/aij/mpi/mumps/mumps.c > [10]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [10]PETSC ERROR: Error in external library! > [10]PETSC ERROR: Error reported by MUMPS in numerical factorization phase: > INFO(1)=-1, INFO(2)=39 > > I tried rerunning with the flags below: > ?-log_summary -ksp__monitor -ksp_type preonly -ksp_view -ksp_max_it 5000 > -pc_type lu -pc_factor_mat_solver_package mumps -mat_mumps_icntl_14 100 > but still get the same error. Did you find out anything new in the meantime? > > Thanks and greetings, > Gabriele > > > ________________________________________ > From: petsc-users-bounces at mcs.anl.gov [petsc-users-bounces at mcs.anl.gov] > on behalf of Max Rudolph [rudolph at berkeley.edu] > Sent: Wednesday, August 31, 2011 3:34 PM > To: petsc-users at mcs.anl.gov > Subject: Re: [petsc-users] Mumps Error > > I'm not sure if you figured out a solution yet, but I think that you might > want to run with > > -mat_mumps_icntl_14 100 > > Max > > >>> 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 >>> <PETScCntor.zip> >> >> >> >> >> _______________________________________________ >> petsc-users mailing list >> petsc-users at mcs.anl.gov >> https://lists.mcs.anl.gov/mailman/listinfo/petsc-users > >
