17#ifdef DEAL_II_WITH_PETSC
27# include <petscconf.h>
60 PetscErrorCode ierr = PCDestroy(&
pc);
70 PetscErrorCode ierr = PCApply(
pc, src, dst);
79 PetscErrorCode ierr = PCApplyTranspose(
pc, src, dst);
88 PetscErrorCode ierr = PCSetUp(
pc);
95 return PetscObjectComm(
reinterpret_cast<PetscObject
>(
pc));
106 PetscErrorCode ierr = PetscObjectGetComm(
107 reinterpret_cast<PetscObject
>(
static_cast<const Mat &
>(matrix)), &
comm);
112 ierr = PCSetOperators(
pc, matrix, matrix);
120 PetscErrorCode ierr = PCCreate(
comm, &
pc);
145 PetscErrorCode ierr = PCCreate(
comm, &
pc);
167 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCJACOBI));
170 ierr = PCSetFromOptions(
pc);
202 PetscErrorCode ierr = PCCreate(
comm, &
pc);
223 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCBJACOBI));
226 ierr = PCSetFromOptions(
pc);
277 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCSOR));
284 ierr = PCSetFromOptions(
pc);
321 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCSOR));
329 ierr = PCSORSetSymmetric(
pc, SOR_SYMMETRIC_SWEEP);
332 ierr = PCSetFromOptions(
pc);
369 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCICC));
376 ierr = PCSetFromOptions(
pc);
413 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCILU));
420 ierr = PCSetFromOptions(
pc);
428 const bool symmetric_operator,
429 const double strong_threshold,
430 const double max_row_sum,
431 const unsigned int aggressive_coarsening_num_levels,
432 const bool output_details,
436 const unsigned int n_sweeps_coarse,
438 const unsigned int max_iter,
440 : symmetric_operator(symmetric_operator)
441 , strong_threshold(strong_threshold)
442 , max_row_sum(max_row_sum)
443 , aggressive_coarsening_num_levels(aggressive_coarsening_num_levels)
444 , output_details(output_details)
445 , relaxation_type_up(relaxation_type_up)
446 , relaxation_type_down(relaxation_type_down)
447 , relaxation_type_coarse(relaxation_type_coarse)
448 , n_sweeps_coarse(n_sweeps_coarse)
456# ifdef DEAL_II_PETSC_WITH_HYPRE
467 std::string string_type;
469 switch (relaxation_type)
472 string_type =
"Jacobi";
476 string_type =
"sequential-Gauss-Seidel";
480 string_type =
"seqboundary-Gauss-Seidel";
483 string_type =
"SOR/Jacobi";
487 string_type =
"backward-SOR/Jacobi";
491 string_type =
"symmetric-SOR/Jacobi";
495 string_type =
" l1scaled-SOR/Jacobi";
499 string_type =
"Gaussian-elimination";
503 string_type =
"l1-Gauss-Seidel";
507 string_type =
"backward-l1-Gauss-Seidel";
513 string_type =
"Chebyshev";
516 string_type =
"FCF-Jacobi";
520 string_type =
"l1scaled-Jacobi";
523 string_type =
"None";
548 PetscErrorCode ierr = PCCreate(
comm, &
pc);
551# ifdef DEAL_II_PETSC_WITH_HYPRE
555 ExcMessage(
"Your PETSc installation does not include a copy of "
556 "the hypre package necessary for this preconditioner."));
575# ifdef DEAL_II_PETSC_WITH_HYPRE
576 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCHYPRE));
579 ierr = PCHYPRESetType(
pc,
"boomeramg");
623 auto relaxation_type_is_symmetric =
640 ExcMessage(
"Use a symmetric smoother for relaxation_type_up"));
645 ExcMessage(
"Use a symmetric smoother for relaxation_type_down"));
650 ExcMessage(
"Use a symmetric smoother for relaxation_type_coarse"));
671 ierr = PCSetFromOptions(
pc);
675 ExcMessage(
"Your PETSc installation does not include a copy of "
676 "the hypre package necessary for this preconditioner."));
686# ifdef DEAL_II_PETSC_WITH_HYPRE
696 (void)additional_data_;
698 ExcMessage(
"Your PETSc installation does not include a copy of "
699 "the hypre package necessary for this preconditioner."));
707 const unsigned int symmetric,
708 const unsigned int n_levels,
709 const double threshold,
711 const bool output_details)
712 : symmetric(symmetric)
714 , threshold(threshold)
716 , output_details(output_details)
744# ifdef DEAL_II_PETSC_WITH_HYPRE
747 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCHYPRE));
750 ierr = PCHYPRESetType(
pc,
"parasails");
761 "ParaSails parameter symmetric can only be equal to 0, 1, 2!"));
763 std::stringstream ssStream;
769 ssStream <<
"nonsymmetric";
781 ssStream <<
"nonsymmetric,SPD";
789 "ParaSails parameter symmetric can only be equal to 0, 1, 2!"));
805 ierr = PCSetFromOptions(
pc);
811 ExcMessage(
"Your PETSc installation does not include a copy of "
812 "the hypre package necessary for this preconditioner."));
843 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCNONE));
846 ierr = PCSetFromOptions(
pc);
854 const double zero_pivot,
855 const double damping)
857 , zero_pivot(zero_pivot)
887 PetscErrorCode ierr = PCSetType(
pc,
const_cast<char *
>(PCLU));
900 ierr = PCSetFromOptions(
pc);
908 const bool use_vertices,
909 const bool use_edges,
910 const bool use_faces,
911 const bool symmetric,
913 : use_vertices(use_vertices)
914 , use_edges(use_edges)
915 , use_faces(use_faces)
916 , symmetric(symmetric)
937 PetscErrorCode ierr = PCCreate(
comm, &
pc);
959# if DEAL_II_PETSC_VERSION_GTE(3, 10, 0)
960 PetscErrorCode ierr = PCSetType(pc,
const_cast<char *
>(PCBDDC));
966 MatType current_type;
970 ierr = PCGetOperators(pc, &A, &P);
972 ierr = PCGetUseAmat(pc, &flg);
975 ierr = MatGetType(flg ? A : P, ¤t_type);
978 strcmp(current_type, MATIS) == 0,
980 "Matrix must be of IS type. For this, the variant of reinit that includes the active dofs must be used."));
984 std::stringstream ssStream;
986 if (additional_data.use_vertices)
990 if (additional_data.use_edges)
994 if (additional_data.use_faces)
998 if (additional_data.symmetric)
1002 if (additional_data.coords.size() > 0)
1006 std::vector<PetscReal> coords_petsc(additional_data.coords.size() *
1008 for (
unsigned int i = 0, j = 0; i < additional_data.coords.size(); ++i)
1010 for (j = 0; j < dim; ++j)
1011 coords_petsc[dim * i + j] = additional_data.coords[i][j];
1014 ierr = PCSetCoordinates(pc,
1016 additional_data.coords.size(),
1017 coords_petsc.data());
1023 ierr = PCSetCoordinates(pc, 0, 0,
nullptr);
1028 ierr = PCSetFromOptions(pc);
1032 false,
ExcMessage(
"BDDC preconditioner requires PETSc 3.10.0 or newer"));
1045 additional_data = additional_data_;
1047 create_pc_with_mat(matrix_);
1066 PetscErrorCode ierr;
1069 ierr = PCDestroy(&
pc);
1074 ierr = PCSetType(
pc, PCSHELL);
1076 ierr = PCShellSetContext(
pc,
static_cast<void *
>(
this));
1084 ierr = PCShellSetName(
pc,
"deal.II user solve");
1092 PetscErrorCode ierr;
1093 ierr = PCSetOperators(
pc, matrix, matrix);
1098# define PetscCall(code) \
1101 PetscErrorCode ierr = (code); \
1110 PetscFunctionBeginUser;
1115 PetscFunctionReturn(PETSC_SUCCESS);
1123 PetscFunctionBeginUser;
1124 PetscCall(PCShellGetContext(ppc, &ctx));
1129 PetscObjectComm((PetscObject)ppc),
1131 "Failure in ::PETScWrappers::PreconditionShell::pcapply. Missing std::function vmult");
1135 const int lineno = __LINE__;
1138 user->vmult(dst, src);
1147 PetscObjectComm((PetscObject)ppc),
1152 PETSC_ERROR_INITIAL,
1153 "Failure in pcapply from ::PETScWrappers::NonlinearSolver");
1156 PetscFunctionReturn(PETSC_SUCCESS);
1164 PetscFunctionBeginUser;
1165 PetscCall(PCShellGetContext(ppc, &ctx));
1170 PetscObjectComm((PetscObject)ppc),
1172 "Failure in ::PETScWrappers::PreconditionShell::pcapply_transpose. Missing std::function vmultT");
1176 const int lineno = __LINE__;
1179 user->vmultT(dst, src);
1188 PetscObjectComm((PetscObject)ppc),
1193 PETSC_ERROR_INITIAL,
1194 "Failure in pcapply_transpose from ::PETScWrappers::NonlinearSolver");
1197 PetscFunctionReturn(PETSC_SUCCESS);
AdditionalData additional_data
void Tvmult(VectorBase &dst, const VectorBase &src) const
MPI_Comm get_mpi_communicator() const
const PC & get_pc() const
void create_pc_with_comm(const MPI_Comm)
virtual ~PreconditionBase()
void vmult(VectorBase &dst, const VectorBase &src) const
void create_pc_with_mat(const MatrixBase &)
PreconditionBlockJacobi()
AdditionalData additional_data
AdditionalData additional_data
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
AdditionalData additional_data
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
AdditionalData additional_data
void initialize(const MatrixBase &matrix, const AdditionalData &additional_data=AdditionalData())
static PetscErrorCode pcapply_transpose(PC pc, Vec src, Vec dst)
static PetscErrorCode pcsetup(PC pc)
void initialize(const MPI_Comm comm)
PreconditionShell()=default
static PetscErrorCode pcapply(PC pc, Vec src, Vec dst)
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
#define DEAL_II_NOT_IMPLEMENTED()
#define Assert(cond, exc)
static ::ExceptionBase & ExcInvalidState()
static ::ExceptionBase & RecoverableUserCallbackError()
static ::ExceptionBase & ExcMessage(std::string arg1)
#define AssertThrow(cond, exc)
void set_option_value(const std::string &name, const std::string &value)
PetscErrorCode pc_set_failed_reason(PC pc, PCFailedReason reason)
void petsc_increment_state_counter(Vec v)
AdditionalData(const bool use_vertices=true, const bool use_edges=false, const bool use_faces=false, const bool symmetric=false, const std::vector< Point< dim > > coords={})
RelaxationType relaxation_type_up
unsigned int n_sweeps_coarse
AdditionalData(const bool symmetric_operator=false, const double strong_threshold=0.25, const double max_row_sum=0.9, const unsigned int aggressive_coarsening_num_levels=0, const bool output_details=false, const RelaxationType relaxation_type_up=RelaxationType::SORJacobi, const RelaxationType relaxation_type_down=RelaxationType::SORJacobi, const RelaxationType relaxation_type_coarse=RelaxationType::GaussianElimination, const unsigned int n_sweeps_coarse=1, const double tol=0.0, const unsigned int max_iter=1, const bool w_cycle=false)
unsigned int aggressive_coarsening_num_levels
RelaxationType relaxation_type_down
RelaxationType relaxation_type_coarse
AdditionalData(const unsigned int levels=0)
AdditionalData(const unsigned int levels=0)
AdditionalData(const double pivoting=1.e-6, const double zero_pivot=1.e-12, const double damping=0.0)
AdditionalData(const unsigned int symmetric=1, const unsigned int n_levels=1, const double threshold=0.1, const double filter=0.05, const bool output_details=false)
AdditionalData(const double omega=1)
AdditionalData(const double omega=1)