16 #ifndef dealii_petsc_solver_h
17 #define dealii_petsc_solver_h
22 #ifdef DEAL_II_WITH_PETSC
27 # include <petscksp.h>
31 # ifdef DEAL_II_WITH_SLEPC
39 # ifdef DEAL_II_WITH_SLEPC
43 class TransformationBase;
54 class PreconditionBase;
175 static PetscErrorCode
177 const PetscInt iteration,
178 const PetscReal residual_norm,
179 KSPConvergedReason *reason,
218 # ifdef DEAL_II_WITH_SLEPC
268 DEAL_II_DEPRECATED_EARLY
270 const MPI_Comm & mpi_communicator,
320 DEAL_II_DEPRECATED_EARLY
322 const MPI_Comm & mpi_communicator,
370 DEAL_II_DEPRECATED_EARLY
372 const MPI_Comm & mpi_communicator,
421 DEAL_II_DEPRECATED_EARLY
423 const MPI_Comm & mpi_communicator,
489 DEAL_II_DEPRECATED_EARLY
491 const MPI_Comm & mpi_communicator,
541 DEAL_II_DEPRECATED_EARLY
543 const MPI_Comm & mpi_communicator,
592 DEAL_II_DEPRECATED_EARLY
594 const MPI_Comm & mpi_communicator,
643 DEAL_II_DEPRECATED_EARLY
645 const MPI_Comm & mpi_communicator,
699 DEAL_II_DEPRECATED_EARLY
701 const MPI_Comm & mpi_communicator,
749 DEAL_II_DEPRECATED_EARLY
751 const MPI_Comm & mpi_communicator,
801 DEAL_II_DEPRECATED_EARLY
803 const MPI_Comm & mpi_communicator,
857 DEAL_II_DEPRECATED_EARLY
859 const MPI_Comm & mpi_communicator,
920 DEAL_II_DEPRECATED_EARLY
922 const MPI_Comm & mpi_communicator,
955 static PetscErrorCode
957 const PetscInt iteration,
958 const PetscReal residual_norm,
959 KSPConvergedReason *reason,
virtual void set_solver_type(KSP &ksp) const =0
SolverControl & control() const
SolverControl & solver_control
void initialize(const PreconditionBase &preconditioner)
SolverBase(SolverControl &cn)
static PetscErrorCode convergence_test(KSP ksp, const PetscInt iteration, const PetscReal residual_norm, KSPConvergedReason *reason, void *solver_control)
virtual ~SolverBase()=default
void set_prefix(const std::string &prefix)
void solve(const MatrixBase &A, VectorBase &x, const VectorBase &b, const PreconditionBase &preconditioner)
std::unique_ptr< SolverData > solver_data
SolverBiCG(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
virtual void set_solver_type(KSP &ksp) const override
const AdditionalData additional_data
SolverBicgstab(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
SolverCGS(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
SolverCG(SolverControl &cn, const AdditionalData &data=AdditionalData())
virtual void set_solver_type(KSP &ksp) const override
virtual void set_solver_type(KSP &ksp) const override
const AdditionalData additional_data
SolverCR(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
SolverChebychev(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
SolverGMRES(SolverControl &cn, const AdditionalData &data=AdditionalData())
virtual void set_solver_type(KSP &ksp) const override
SolverLSQR(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
SolverPreOnly(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
virtual void set_solver_type(KSP &ksp) const override
const AdditionalData additional_data
SolverRichardson(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
SolverTCQMR(SolverControl &cn, const AdditionalData &data=AdditionalData())
const AdditionalData additional_data
virtual void set_solver_type(KSP &ksp) const override
SolverTFQMR(SolverControl &cn, const AdditionalData &data=AdditionalData())
void solve(const MatrixBase &A, VectorBase &x, const VectorBase &b)
virtual void set_solver_type(KSP &ksp) const override
const AdditionalData additional_data
void set_symmetric_mode(const bool flag)
static PetscErrorCode convergence_test(KSP ksp, const PetscInt iteration, const PetscReal residual_norm, KSPConvergedReason *reason, void *solver_control)
SparseDirectMUMPS(SolverControl &cn, const AdditionalData &data=AdditionalData())
std::unique_ptr< SolverDataMUMPS > solver_data
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)
AdditionalData(const unsigned int restart_parameter=30, const bool right_preconditioning=false)
bool right_preconditioning
unsigned int restart_parameter
AdditionalData(const double omega=1)