[petsc-users] SLEPc Eigenvalue problem - Zero Pivot

Fri May 15 07:44:01 CDT 2015

You have to choose MUMPS as the solver. ONCE AGAIN, read section 3.4.1 of the users manual.
Jose

El 15/05/2015, a las 13:48, venkatesh g escribió:

> Ok, I followed this process and compiled ex7 successfully. 
> 
> Now I give the command: './ex7 -f1 A11 -f2 B11 -eps_nev 4 -eps_target 0.5 -st_type sinvert'
> 
> The error is the zero pivot problem:
> 
>  Reading COMPLEX matrices from binary files...
> [0]PETSC ERROR: --------------------- Error Message ------------------------------------
> [0]PETSC ERROR: Detected zero pivot in LU factorization:
> see http://www.mcs.anl.gov/petsc/documentation/faq.html#ZeroPivot!
> [0]PETSC ERROR: Zero pivot row 420 value 0 tolerance 2.22045e-14!
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: Petsc Release Version 3.3.0, Patch 7, Sat May 11 22:15:24 CDT 2013 
> [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: ./ex7 on a linux-gnu named earth.ceas.iisc.ernet.in by venkatesh Fri May 15 17:07:02 2015
> [0]PETSC ERROR: Libraries linked from /cluster/share/venkatesh/petsc-3.3-p7/linux-gnu-intel-cmplx-batch/lib
> [0]PETSC ERROR: Configure run at Fri May 15 15:10:57 2015
> [0]PETSC ERROR: Configure options --with-cc=gcc --with-fc=ifort --with-blas-lapack-dir=/opt/intel/composer_xe_2013.5.192/mkl/lib/intel64 --download-mpich=/cluster/share/apps/mpich2-1.4.1p1.tar.gz --download-mumps=/cluster/share/apps/MUMPS_4.10.0-p3.tar.gz --download-scalapack=/cluster/share/apps/SCALAPACK-1.7.tar.gz --download-blacs=/cluster/share/apps/blacs-dev.tar.gz --download-parmetis=/cluster/share/apps/parmetis-4.0.2-p3.tar.gz --download-metis=/cluster/share/apps/metis-5.0.2-p3.tar.gz --download-cmake=/cluster/share/apps/cmake-2.8.7.tar.gz --with-scalar-type=complex --with-fortran-kernels=generic --with-large-file-io=1
> [0]PETSC ERROR: ------------------------------------------------------------------------
> [0]PETSC ERROR: MatPivotCheck_none() line 583 in /cluster/share/venkatesh/petsc-3.3-p7/include/petsc-private/matimpl.h
> [0]PETSC ERROR: MatPivotCheck() line 602 in /cluster/share/venkatesh/petsc-3.3-p7/include/petsc-private/matimpl.h
> [0]PETSC ERROR: MatLUFactorNumeric_SeqAIJ() line 570 in /cluster/share/venkatesh/petsc-3.3-p7/src/mat/impls/aij/seq/aijfact.c
> [0]PETSC ERROR: MatLUFactorNumeric() line 2803 in /cluster/share/venkatesh/petsc-3.3-p7/src/mat/interface/matrix.c
> [0]PETSC ERROR: PCSetUp_LU() line 160 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/impls/factor/lu/lu.c
> [0]PETSC ERROR: PCSetUp() line 832 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: KSPSetUp() line 278 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: PCSetUp_Redundant() line 176 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/impls/redundant/redundant.c
> [0]PETSC ERROR: PCSetUp() line 832 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/pc/interface/precon.c
> [0]PETSC ERROR: KSPSetUp() line 278 in /cluster/share/venkatesh/petsc-3.3-p7/src/ksp/ksp/interface/itfunc.c
> [0]PETSC ERROR: STSetUp_Sinvert() line 153 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/st/impls/sinvert/sinvert.c
> [0]PETSC ERROR: STSetUp() line 280 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/st/interface/stsolve.c
> [0]PETSC ERROR: EPSSetUp() line 204 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/eps/interface/setup.c
> [0]PETSC ERROR: EPSSolve() line 109 in /cluster/share/venkatesh/petsc-3.3-p7/slepc-3.3-p4/src/eps/interface/solve.c
> [0]PETSC ERROR: main() line 148 in ex7.c
> application called MPI_Abort(MPI_COMM_WORLD, 71) - process 0
> [unset]: aborting job:
> application called MPI_Abort(MPI_COMM_WORLD, 71) - process 0
> 
> Venkatesh
> 
> On Fri, May 15, 2015 at 4:22 PM, Jose E. Roman <jroman at dsic.upv.es> wrote:
> 
> El 15/05/2015, a las 12:49, venkatesh g escribió:
> 
> > Hi,
> > 'make test' was successful on the $PETSC_DIR.
> > I configured MUMPS, metis, parmetis, scalapack through '--download-mumps=/path/MUMPS_4.10.0-p3.tar.gz' and etc..
> >
> > I will anyway download the new version and start installing it now., Is it due to some flag issue ? I have attached my makefile.
> >
> > Venkatesh
> >
> 
> You must link with all libraries, otherwise you will get undefined symbols. The easiest way is to include PETSc and SLEPc makefile definitions. See for instance section 1.4.2 of SLEPc's users manual.
> 
> Jose
> 
> 