SLEPcでmatrix-free、LOBPCGまたはRQCGを用いることができないというエラー

[t-sakashita]

routine=lobpcg
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: No support for this operation for this object type
[0]PETSC ERROR: Not coded for matrix type shell
[0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.15.1, Jun 17, 2021 
[0]PETSC ERROR: Unknown Name on a  named sakashitatatsuyanoMacBook-Pro.local by sakashitatatsuya Thu Jun 24 11:26:31 2021
[0]PETSC ERROR: Configure options --prefix=/Users/sakashitatatsuya/rokko/petsc/petsc-3.15.1-1/Debug --with-cc=mpicc --COPTFLAGS="-g -O0" --with-cxx=mpicxx --CXXOPTFLAGS="-g -O0" --with-fc=mpif77 --FOPTFLAGS="-g -O0" --with-mpiexec=mpiexec --with-blaslapack-lib="-framework Accelerate" --with-c++-support=1 --with-debugging=1
[0]PETSC ERROR: #1 MatGetDiagonalBlock() at /Users/sakashitatatsuya/build/petsc-3.15.1-build-Debug/src/mat/interface/matrix.c:337
[0]PETSC ERROR: #2 PCSetUp_BJacobi() at /Users/sakashitatatsuya/build/petsc-3.15.1-build-Debug/src/ksp/pc/impls/bjacobi/bjacobi.c:121
[0]PETSC ERROR: #3 PCSetUp() at /Users/sakashitatatsuya/build/petsc-3.15.1-build-Debug/src/ksp/pc/interface/precon.c:1015
[0]PETSC ERROR: #4 KSPSetUp() at /Users/sakashitatatsuya/build/petsc-3.15.1-build-Debug/src/ksp/ksp/interface/itfunc.c:406
[0]PETSC ERROR: #5 KSPMatSolve() at /Users/sakashitatatsuya/build/petsc-3.15.1-build-Debug/src/ksp/ksp/interface/itfunc.c:1223
[0]PETSC ERROR: #6 STMatMatSolve() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/sys/classes/st/interface/stsles.c:229
[0]PETSC ERROR: #7 STApplyMat_Generic() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/sys/classes/st/interface/stsolve.c:82
[0]PETSC ERROR: #8 STApplyMat() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/sys/classes/st/interface/stsolve.c:117
[0]PETSC ERROR: #9 EPSSolve_LOBPCG() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/eps/impls/cg/lobpcg/lobpcg.c:313
[0]PETSC ERROR: #10 EPSSolve() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/eps/interface/epssolve.c:149
[0]PETSC ERROR: #11 EPSGetConverged() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/eps/interface/epssolve.c:264
[0]PETSC ERROR: #12 EPSGetEigenvalue() at /Users/sakashitatatsuya/build/slepc-3.15.1-build-Debug/src/eps/interface/epssolve.c:453
/Users/sakashitatatsuya/development/rokko/test/sparse/laplacian_mfree_variable_local_rows_mpi.cpp:66: Failure
The difference between eigval and th_eigval is 3.9765608475606973, which exceeds eigval*eps, where
eigval evaluates to 6.9287851074886541e-310,
th_eigval evaluates to 3.9765608475606973, and
eigval*eps evaluates to 6.9287864985905584e-318.

1プロセスでは、正常終了した。 複数プロセスでエラーが起こる。

[t-sakashita]

https://lists.mcs.anl.gov/pipermail/petsc-users/2016-January/028060.html

[t-sakashita]

https://www.mcs.anl.gov/petsc/petsc-current/docs/manualpages/Mat/MatGetDiagonalBlock.html

[t-sakashita]

3ce7725eb2b5af1c55eb190699087a9f6533e22e