<div dir="ltr">Ok, thank you very much for your help and links to articles.</div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">вс, 15 мая 2022 г. в 20:02, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es">jroman@dsic.upv.es</a>>:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Preconditioned eigensolvers avoid an explicit inversion. You can try Davidson methods <a href="https://doi.org/10.1145/2543696" rel="noreferrer" target="_blank">https://doi.org/10.1145/2543696</a> or LOBPCG if the problem is symmetric-definite. However, for reasonable convergence these methods usually require a good preconditioner, so if your matrix is implicit you still need to explicitly build some type of approximation that can be used to build the preconditioner.<br>
<br>
Jose<br>
<br>
<br>
> El 15 may 2022, a las 18:54, Миша Сальников <<a href="mailto:rolekst095@gmail.com" target="_blank">rolekst095@gmail.com</a>> escribió:<br>
> <br>
> Thanks a lot for your answer. I found a solution in the documentation. So if I understand correctly SLEPc during eigenvalue iterations of Ax = k Bx is trying to calculate the inverse of B or y = Bx. So the question is don't you know maybe there exist methods or ways for generalized eigenproblems without inverting of B? Again thanks a lot for the answer     <br>
> <br>
> вс, 15 мая 2022 г. в 17:52, Jose E. Roman <<a href="mailto:jroman@dsic.upv.es" target="_blank">jroman@dsic.upv.es</a>>:<br>
> You cannot compute the LU factorization of a shell matrix. You will have to use a preconditioned iterative linear solver, see section 3.4.1 of the SLEPc users guide.<br>
> <br>
> Jose<br>
> <br>
> <br>
> > El 15 may 2022, a las 14:39, Миша Сальников <<a href="mailto:rolekst095@gmail.com" target="_blank">rolekst095@gmail.com</a>> escribió:<br>
> > <br>
> > Hi, I am developing a program for solving generalized eigenvalue problem Ax = k Bx when both operators A and B are matrix-free. I found some examples on SLEPc documentation but and there was and example where problem Ax = k x are solving matrix free. So I tried to change this example a little bit and added operator B as identity matrix-free operator. But for now I am getting an error. My source code is attached. And the error message is below<br>
> > Can you please explain to me what I am doing wrong? <br>
> > <br>
> > [0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------<br>
> > [0]PETSC ERROR: See <a href="https://petsc.org/release/overview/linear_solve_table/" rel="noreferrer" target="_blank">https://petsc.org/release/overview/linear_solve_table/</a> for possible LU and Cholesky solvers<br>
> > [0]PETSC ERROR: Could not locate a solver type for factorization type LU and matrix type shell.<br>
> > [0]PETSC ERROR: See <a href="https://petsc.org/release/faq/" rel="noreferrer" target="_blank">https://petsc.org/release/faq/</a> for trouble shooting.<br>
> > [0]PETSC ERROR: Petsc Release Version 3.15.5, Sep 29, 2021 <br>
> > [0]PETSC ERROR: ./exe_eps_impl_matr on a  named LAPTOP-7DN2DH7N by misha Sun May 15 15:21:59 2022<br>
> > [0]PETSC ERROR: Configure options --build=x86_64-linux-gnu --prefix=/usr --includedir=${prefix}/include --mandir=${prefix}/share/man --infodir=${prefix}/share/info --sysconfdir=/etc --localstatedir=/var --with-option-checking=0 --with-silent-rules=0 --libdir=${prefix}/lib/x86_64-linux-gnu --runstatedir=/run --with-maintainer-mode=0 --with-dependency-tracking=0 --with-64-bit-indices --with-debugging=0 --shared-library-extension=64_real --with-shared-libraries --with-pic=1 --with-cc=mpicc --with-cxx=mpicxx --with-fc=mpif90 --with-cxx-dialect=C++11 --with-opencl=1 --with-blas-lib=-lblas --with-lapack-lib=-llapack --with-scalapack=1 --with-scalapack-lib=-lscalapack-openmpi --with-ptscotch=1 --with-ptscotch-include=/usr/include/scotch --with-ptscotch-lib="-lptesmumps -lptscotch -lptscotcherr" --with-fftw=1 --with-fftw-include="[]" --with-fftw-lib="-lfftw3 -lfftw3_mpi" --with-yaml=1 --with-hdf5-include=/usr/include/hdf5/openmpi --with-hdf5-lib="-L/usr/lib/x86_64-linux-gnu/hdf5/openmpi -L/usr/lib/x86_64-linux-gnu/openmpi/lib -lhdf5 -lmpi" --CXX_LINKER_FLAGS=-Wl,--no-as-needed --with-hypre=1 --with-hypre-include=/usr/include/hypre64m --with-hypre-lib=-lHYPRE64m --with-mumps=1 --with-mumps-include="[]" --with-mumps-lib="-ldmumps_64 -lzmumps_64 -lsmumps_64 -lcmumps_64 -lmumps_common_64 -lpord_64" --with-suitesparse=1 --with-suitesparse-include=/usr/include/suitesparse --with-suitesparse-lib="-lumfpack -lamd -lcholmod -lklu" --with-zoltan=1 --with-zoltan-include=/usr/include/trilinos --with-zoltan-lib=-ltrilinos_zoltan --prefix=/usr/lib/petscdir/petsc64-3.15/x86_64-linux-gnu-real --PETSC_ARCH=x86_64-linux-gnu-real-64 CFLAGS="-g -O2 -ffile-prefix-map=/build/petsc-tFq6Xk/petsc-3.15.5+dfsg1=. -flto=auto -ffat-lto-objects -flto=auto -ffat-lto-objects -fstack-protector-strong -Wformat -Werror=format-security -fPIC" CXXFLAGS="-g -O2 -ffile-prefix-map=/build/petsc-tFq6Xk/petsc-3.15.5+dfsg1=. -flto=auto -ffat-lto-objects -flto=auto -ffat-lto-objects -fstack-protector-strong -Wformat -Werror=format-security -fPIC" FCFLAGS="-g -O2 -ffile-prefix-map=/build/petsc-tFq6Xk/petsc-3.15.5+dfsg1=. -flto=auto -ffat-lto-objects -flto=auto -ffat-lto-objects -fstack-protector-strong -fPIC -ffree-line-length-0" FFLAGS="-g -O2 -ffile-prefix-map=/build/petsc-tFq6Xk/petsc-3.15.5+dfsg1=. -flto=auto -ffat-lto-objects -flto=auto -ffat-lto-objects -fstack-protector-strong -fPIC -ffree-line-length-0" CPPFLAGS="-Wdate-time -D_FORTIFY_SOURCE=2" LDFLAGS="-Wl,-Bsymbolic-functions -flto=auto -ffat-lto-objects -flto=auto -Wl,-z,relro -fPIC" MAKEFLAGS=w<br>
> > [0]PETSC ERROR: #1 MatGetFactor() at ./src/mat/interface/matrix.c:4756<br>
> > [0]PETSC ERROR: #2 PCSetUp_LU() at ./src/ksp/pc/impls/factor/lu/lu.c:82<br>
> > [0]PETSC ERROR: #3 PCSetUp() at ./src/ksp/pc/interface/precon.c:1017<br>
> > [0]PETSC ERROR: #4 KSPSetUp() at ./src/ksp/ksp/interface/itfunc.c:406<br>
> > [0]PETSC ERROR: #5 STSetUp_Shift() at ./src/sys/classes/st/impls/shift/shift.c:107<br>
> > [0]PETSC ERROR: #6 STSetUp() at ./src/sys/classes/st/interface/stsolve.c:582<br>
> > [0]PETSC ERROR: #7 EPSSetUp() at ./src/eps/interface/epssetup.c:350<br>
> > [0]PETSC ERROR: #8 EPSSolve() at ./src/eps/interface/epssolve.c:136<br>
> > [0]PETSC ERROR: #9 main() at ex_implicit_matr.cpp:95<br>
> > [0]PETSC ERROR: No PETSc Option Table entries<br>
> > [0]PETSC ERROR: ----------------End of Error Message -------send entire error message to petsc-maint@mcs.anl.gov----------<br>
> > --------------------------------------------------------------------------<br>
> > MPI_ABORT was invoked on rank 0 in communicator MPI_COMM_WORLD<br>
> > with errorcode 92.<br>
> > <br>
> > NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.<br>
> > You may or may not see output from other processes, depending on<br>
> > exactly when Open MPI kills them.<br>
> > <ex_implicit_matr.cpp><br>
> <br>
<br>
</blockquote></div>