<div dir="ltr"><span class="gmail-tlid-translation gmail-translation" lang="en"><span title="" class="gmail-">Hi everyone,</span></span><br><br>I'm solving the eigenvalue problem of three bodies in the same program, it generates a three sets of matrices.<br><br>I installed PETSc as optimized version:<br><br>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<br><br>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.<br><br>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:<br><br>ierr = EPSSetProblemType(eps,EPS_GHEP); CHKERRQ(ierr);<br>ierr = EPSSetType(eps,EPSLOBPCG); CHKERRQ(ierr); // because I need the smallest<br><br>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:<br><br>[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>[0]PETSC ERROR: The inner product is not well defined: indefinite matrix<br>[0]PETSC ERROR: See <a href="https://www.mcs.anl.gov/petsc/documentation/faq.html" target="_blank">https://www.mcs.anl.gov/petsc/documentation/faq.html</a> for trouble shooting.<br>[0]PETSC ERROR: Petsc Release Version 3.12.3, Jan, 03, 2020 <br>[0]PETSC ERROR: ./comp on a linux-opt-02 named lnx by ayala Thu Feb  6 17:20:16 2020<br>[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<br>[0]PETSC ERROR: #1 BV_SafeSqrt() line 130 in /home/ayala/Documents/SLEPc/slepc-3.12.2/include/slepc/private/bvimpl.h<br>[0]PETSC ERROR: #2 BVNorm_Private() line 473 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvglobal.c<br>[0]PETSC ERROR: #3 BVNormColumn() line 718 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvglobal.c<br>[0]PETSC ERROR: #4 BV_NormVecOrColumn() line 26 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c<br>[0]PETSC ERROR: #5 BVOrthogonalizeCGS1() line 136 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c<br>[0]PETSC ERROR: #6 BVOrthogonalizeGS() line 188 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c<br>[0]PETSC ERROR: #7 BVOrthonormalizeColumn() line 416 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/sys/classes/bv/interface/bvorthog.c<br>[0]PETSC ERROR: #8 EPSSolve_LOBPCG() line 144 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/eps/impls/cg/lobpcg/lobpcg.c<br>[0]PETSC ERROR: #9 EPSSolve() line 149 in /home/ayala/Documents/SLEPc/slepc-3.12.2/src/eps/interface/epssolve.c<br><br>Number of iterations of the method: 0<br>Number of linear iterations of the method: 0<br>Solution method: lobpcg<br><br>Number of requested eigenvalues: 6<br>Stopping condition: tol=1e-06, maxit=10000<br>Number of converged eigenpairs: 0<br><br>The problem appears even when I run the same compilation.<br><br>Anyone have any suggestions to solve the problem?<br><br>Kind regards.</div>