Hi everyone, I'm solving the eigenvalue problem of three bodies in the same program, it generates a three sets of matrices.
I installed PETSc as optimized version: Configure options --with-debugging=0 COPTFLAGS="-O2 -march=native -mtune=native" CXXOPTFLAGS="-O2 -march=native -mtune=native" FOPTFLAGS="-O2 -march=native -mtune=native" --download-mpich --download-superlu_dist --download-metis --download-parmetis --download-cmake --download-fblaslapack=1 --with-cxx-dialect=C++11 Then I installed SLEPc in the standard form and referring to PETSc optimized directory. I did NOT install SLEPc with --with-debugging=0, because I'm still testing my code. My matrices comes from DMCreateMatrix, the stiffness matrix and mass matrix. I use the function MatIsSymmetric to check if my matrices are symmetric or not, and always the matrices are symmetric (even when the program crash). For that reason I use: ierr = EPSSetProblemType(eps,EPS_GHEP); CHKERRQ(ierr); ierr = EPSSetType(eps,EPSLOBPCG); CHKERRQ(ierr); // because I need the smallest The problem is: sometimes the code works well for the three sets of matrices and I get the expected results, but sometimes it does not happen, it only works for the first sets of matrices, and then it crashes, and when that occurs the error message is: [0]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [0]PETSC ERROR: The inner product is not well defined: indefinite matrix [0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. [0]PETSC ERROR: Petsc Release Version 3.12.3, Jan, 03, 2020 [0]PETSC ERROR: ./comp on a linux-opt-02 named lnx by ayala Thu Feb 6 17:20:16 2020 [0]PETSC ERROR: Configure options --with-debugging=0 COPTFLAGS="-O2 -march=native -mtune=native" CXXOPTFLAGS="-O2 -march=native -mtune=native" FOPTFLAGS="-O2 -march=native -mtune=native" --download-mpich --download-superlu_dist --download-metis --download-parmetis --download-cmake --download-fblaslapack=1 --with-cxx-dialect=C++11 [0]PETSC ERROR: #1 BV_SafeSqrt() line 130 in /home/ayala/Documents/SLEPc/slepc-3.12.2/include/slepc/private/bvimpl.h [0]PETSC ERROR: #2 BVNorm_Private() line 473 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvglobal.c [0]PETSC ERROR: #3 BVNormColumn() line 718 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvglobal.c [0]PETSC ERROR: #4 BV_NormVecOrColumn() line 26 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c [0]PETSC ERROR: #5 BVOrthogonalizeCGS1() line 136 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c [0]PETSC ERROR: #6 BVOrthogonalizeGS() line 188 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c [0]PETSC ERROR: #7 BVOrthonormalizeColumn() line 416 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c [0]PETSC ERROR: #8 EPSSolve_LOBPCG() line 144 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/eps/impls/cg/lobpcg/lobpcg.c [0]PETSC ERROR: #9 EPSSolve() line 149 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/eps/interface/epssolve.c Number of iterations of the method: 0 Number of linear iterations of the method: 0 Solution method: lobpcg Number of requested eigenvalues: 6 Stopping condition: tol=1e-06, maxit=10000 Number of converged eigenpairs: 0 The problem appears even when I run the same compilation. Anyone have any suggestions to solve the problem? Kind regards.
