Hi, > [0]PETSC ERROR: See > https://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html for > possible LU and Cholesky solvers > [0]PETSC ERROR: Could not locate a solver package for factorization type ILU > and matrix type mpiaij.
This means that ILU is not available for parallel matrices (mpiaij). You need to pick a different preconditioner. You can try PreconditionBlockJacobi which is a blocked ILU and probably what you wanted to do. On Tue, Mar 15, 2022 at 12:39 PM Hermes Sampedro <[email protected]> wrote: > > Dear all, > > I am getting a running error that I can not really understand when calling > PETScWrappers::PreconditionILU preconditioner(system_matrix). However, it > seems to work with PETScWrappers::PreconditionBlockJacobi > preconditioner(system_matrix). > > May I ask for help to understand/solve this issue? The solver function ad > error that I get is the following: > > void LaplaceProblem<dim>::solve() > > { > > PETScWrappers::MPI::Vector > completely_distributed_solution(locally_owned_dofs,mpi_communicator); > > SolverControl cn(completely_distributed_solution.size(), 1e-8 * > system_rhs.l2_norm()); > > PETScWrappers::SolverGMRES solver(cn, mpi_communicator); > > PETScWrappers::PreconditionILU preconditioner(system_matrix); > > solver.solve(system_matrix, completely_distributed_solution, > system_rhs, preconditioner); > > constraints.distribute(completely_distributed_solution); > > locally_relevant_solution = completely_distributed_solution; > > } > > > [0]PETSC ERROR: --------------------- Error Message > -------------------------------------------------------------- > > [0]PETSC ERROR: See > https://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html for > possible LU and Cholesky solvers > > [0]PETSC ERROR: Could not locate a solver package for factorization type ILU > and matrix type mpiaij. > > [0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for > trouble shooting. > > [0]PETSC ERROR: Petsc Release Version 3.13.1, May 02, 2020 > > [0]PETSC ERROR: ./waveLaplaceSolver on a named gbarlogin1 by hsllo Fri Mar > 11 11:05:23 2022 > > [0]PETSC ERROR: Configure options > --prefix=/zhome/32/9/115503/dealii-candi/petsc-3.13.1 --with-debugging=0 > --with-shared-libraries=1 --with-mpi=1 --with-x=0 --with-64-bit-indices=0 > --download-hypre=1 CC=mpicc CXX=mpicxx FC=mpif90 > --with-blaslapack-dir=/appl/OpenBLAS/0.3.17/XeonE5-2660v3/gcc-11.2.0/lib > --with-parmetis-dir=/zhome/32/9/115503/dealii-candi/parmetis-4.0.3 > --with-metis-dir=/zhome/32/9/115503/dealii-candi/parmetis-4.0.3 > --download-scalapack=1 --download-mumps=1 > > [0]PETSC ERROR: #1 MatGetFactor() line 4492 in > /zhome/32/9/115503/dealii-candi/tmp/build/petsc-lite-3.13.1/src/mat/interface/matrix.c > > [0]PETSC ERROR: #2 PCSetUp_ILU() line 133 in > /zhome/32/9/115503/dealii-candi/tmp/build/petsc-lite-3.13.1/src/ksp/pc/impls/factor/ilu/ilu.c > > [0]PETSC ERROR: #3 PCSetUp() line 894 in > /zhome/32/9/115503/dealii-candi/tmp/build/petsc-lite-3.13.1/src/ksp/pc/interface/precon.c > > ---------------------------------------------------- > > Exception on processing: > > -------------------------------------------------------- > > An error occurred in line <431> of file > </zhome/32/9/115503/dealii-candi/tmp/unpack/deal.II-v9.3.1/source/lac/petsc_precondition.cc> > in function > > void dealii::PETScWrappers::PreconditionILU::initialize(const > dealii::PETScWrappers::MatrixBase&, const > dealii::PETScWrappers::PreconditionILU::AdditionalData&) > > The violated condition was: > > ierr == 0 > > Additional information: > > deal.II encountered an error while calling a PETSc function. > > The description of the error provided by PETSc is "See > > https://www.mcs.anl.gov/petsc/documentation/linearsolvertable.html for > > possible LU and Cholesky solvers". > > The numerical value of the original error code is 92. > > > > > Thank you very much > > Regards, > > H > > -- > The deal.II project is located at http://www.dealii.org/ > For mailing list/forum options, see > https://groups.google.com/d/forum/dealii?hl=en > --- > You received this message because you are subscribed to the Google Groups > "deal.II User Group" group. > To unsubscribe from this group and stop receiving emails from it, send an > email to [email protected]. > To view this discussion on the web visit > https://groups.google.com/d/msgid/dealii/c4109bcf-cd4e-453e-9c81-9ed1230a235cn%40googlegroups.com. -- Timo Heister http://www.math.clemson.edu/~heister/ -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/CAMRj59HC9OzxGkNWSniENzVEhZFcKgWqo1eUH1ZykQm-atiSrQ%40mail.gmail.com.
