16 #ifndef dealii_trilinos_precondition_h
17 #define dealii_trilinos_precondition_h
22 #ifdef DEAL_II_WITH_TRILINOS
29 # include <Epetra_Map.h>
30 # include <Epetra_MpiComm.h>
31 # include <Epetra_MultiVector.h>
32 # include <Epetra_RowMatrix.h>
33 # include <Epetra_Vector.h>
34 # include <Teuchos_ParameterList.hpp>
40 class Ifpack_Preconditioner;
41 class Ifpack_Chebyshev;
44 class MultiLevelPreconditioner;
52 template <
typename number>
54 template <
typename number>
159 const ::Vector<double> &src)
const;
167 const ::LinearAlgebra::distributed::Vector<double> &src)
const;
175 const ::LinearAlgebra::distributed::Vector<double> &src)
const;
221 <<
"The sparse matrix the preconditioner is based on "
222 <<
"uses a map that is not compatible to the one in vector "
223 << arg1 <<
". Check preconditioner and matrix setup.");
371 const unsigned int overlap = 0,
465 const unsigned int overlap = 0,
551 const double omega = 1,
654 const double omega = 1,
656 const unsigned int overlap = 0,
765 const double omega = 1,
767 const unsigned int overlap = 0,
891 const unsigned int overlap = 0);
1010 const unsigned int overlap = 0);
1116 const unsigned int overlap = 0);
1377 std::vector<std::vector<bool>>(0),
1416 Teuchos::ParameterList & parameter_list,
1417 std::unique_ptr<Epetra_MultiVector> &distributed_constant_modes,
1418 const Epetra_RowMatrix &
matrix)
const;
1429 Teuchos::ParameterList & parameter_list,
1430 std::unique_ptr<Epetra_MultiVector> &distributed_constant_modes,
1440 Teuchos::ParameterList & parameter_list,
1441 std::unique_ptr<Epetra_MultiVector> &distributed_constant_modes,
1442 const Epetra_RowMatrix &
matrix)
const;
1451 Teuchos::ParameterList & parameter_list,
1452 std::unique_ptr<Epetra_MultiVector> &distributed_constant_modes,
1640 const Teuchos::ParameterList &ml_parameters);
1651 const Teuchos::ParameterList &ml_parameters);
1659 template <
typename number>
1661 initialize(const ::SparseMatrix<number> &deal_ii_sparse_matrix,
1663 const double drop_tolerance = 1
e-13,
1664 const ::SparsityPattern *use_this_sparsity =
nullptr);
1703 # if defined(DOXYGEN) || defined(DEAL_II_TRILINOS_WITH_MUELU)
1742 std::vector<std::vector<bool>>(0),
1901 Teuchos::ParameterList &muelu_parameters);
1910 Teuchos::ParameterList &muelu_parameters);
1918 template <
typename number>
1920 initialize(const ::SparseMatrix<number> &deal_ii_sparse_matrix,
1922 const double drop_tolerance = 1
e-13,
1923 const ::SparsityPattern *use_this_sparsity =
nullptr);
1992 const ::Vector<double> &src)
const override;
2000 const ::Vector<double> &src)
const override;
2008 const ::LinearAlgebra::distributed::Vector<double> &src)
2018 const ::LinearAlgebra::distributed::Vector<double> &src)
2056 Assert(dst.trilinos_partitioner().SameAs(
2059 Assert(src.trilinos_partitioner().SameAs(
2063 const int ierr =
preconditioner->ApplyInverse(src.trilinos_vector(),
2064 dst.trilinos_vector());
2071 Assert(dst.trilinos_partitioner().SameAs(
2074 Assert(src.trilinos_partitioner().SameAs(
2079 const int ierr =
preconditioner->ApplyInverse(src.trilinos_vector(),
2080 dst.trilinos_vector());
2097 const ::Vector<double> &src)
const
2103 Epetra_Vector tril_dst(
View,
2106 Epetra_Vector tril_src(
View,
2108 const_cast<double *
>(src.begin()));
2110 const int ierr =
preconditioner->ApplyInverse(tril_src, tril_dst);
2117 const ::Vector<double> &src)
const
2123 Epetra_Vector tril_dst(
View,
2126 Epetra_Vector tril_src(
View,
2128 const_cast<double *
>(src.begin()));
2131 const int ierr =
preconditioner->ApplyInverse(tril_src, tril_dst);
2147 Epetra_Vector tril_dst(
View,
2150 Epetra_Vector tril_src(
View,
2152 const_cast<double *
>(src.
begin()));
2154 const int ierr =
preconditioner->ApplyInverse(tril_src, tril_dst);
2167 Epetra_Vector tril_dst(
View,
2170 Epetra_Vector tril_src(
View,
2172 const_cast<double *
>(src.
begin()));
2175 const int ierr =
preconditioner->ApplyInverse(tril_src, tril_dst);
size_type locally_owned_size() const
std::shared_ptr< SparseMatrix > trilinos_matrix
size_type memory_consumption() const
virtual ~PreconditionAMGMueLu() override=default
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
~PreconditionAMG() override
std::shared_ptr< SparseMatrix > trilinos_matrix
size_type memory_consumption() const
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
IndexSet locally_owned_range_indices() const
Epetra_MpiComm communicator
virtual void Tvmult(::LinearAlgebra::distributed::Vector< double > &dst, const ::LinearAlgebra::distributed::Vector< double > &src) const
virtual void Tvmult(MPI::Vector &dst, const MPI::Vector &src) const
~PreconditionBase() override=default
virtual void vmult(MPI::Vector &dst, const MPI::Vector &src) const
Epetra_Operator & trilinos_operator() const
std::shared_ptr< Epetra_Map > vector_distributor
Teuchos::RCP< Epetra_Operator > preconditioner
virtual void Tvmult(::Vector< double > &dst, const ::Vector< double > &src) const
MPI_Comm get_mpi_communicator() const
virtual void vmult(::Vector< double > &dst, const ::Vector< double > &src) const
IndexSet locally_owned_domain_indices() const
virtual void vmult(::LinearAlgebra::distributed::Vector< double > &dst, const ::LinearAlgebra::distributed::Vector< double > &src) const
::types::global_dof_index size_type
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void vmult(LinearAlgebra::distributed::Vector< double > &dst, const ::LinearAlgebra::distributed::Vector< double > &src) const override
void vmult(MPI::Vector &dst, const MPI::Vector &src) const override
void Tvmult(MPI::Vector &dst, const MPI::Vector &src) const override
void Tvmult(LinearAlgebra::distributed::Vector< double > &dst, const ::LinearAlgebra::distributed::Vector< double > &src) const override
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
void initialize(const SparseMatrix &matrix, const AdditionalData &additional_data=AdditionalData())
#define DEAL_II_NAMESPACE_OPEN
#define DEAL_II_NAMESPACE_CLOSE
static ::ExceptionBase & ExcTrilinosError(int arg1)
static ::ExceptionBase & ExcNonMatchingMaps(std::string arg1)
#define Assert(cond, exc)
#define AssertDimension(dim1, dim2)
#define DeclException1(Exception1, type1, outsequence)
#define AssertThrow(cond, exc)
typename ::internal::FEInterfaceViews::ViewType< dim, spacedim, Extractor >::type View
@ matrix
Contents is actually a matrix.
SymmetricTensor< 2, dim, Number > e(const Tensor< 2, dim, Number > &F)
unsigned int global_dof_index
AdditionalData(const bool elliptic=true, const unsigned int n_cycles=1, const bool w_cycle=false, const double aggregation_threshold=1e-4, const std::vector< std::vector< bool >> &constant_modes=std::vector< std::vector< bool >>(0), const unsigned int smoother_sweeps=2, const unsigned int smoother_overlap=0, const bool output_details=false, const char *smoother_type="Chebyshev", const char *coarse_type="Amesos-KLU")
double aggregation_threshold
std::vector< std::vector< bool > > constant_modes
unsigned int smoother_sweeps
unsigned int smoother_overlap
const char * smoother_type
std::vector< std::vector< bool > > constant_modes
double aggregation_threshold
AdditionalData(const bool elliptic=true, const bool higher_order_elements=false, const unsigned int n_cycles=1, const bool w_cycle=false, const double aggregation_threshold=1e-4, const std::vector< std::vector< bool >> &constant_modes=std::vector< std::vector< bool >>(0), const unsigned int smoother_sweeps=2, const unsigned int smoother_overlap=0, const bool output_details=false, const char *smoother_type="Chebyshev", const char *coarse_type="Amesos-KLU")
unsigned int smoother_sweeps
bool higher_order_elements
unsigned int smoother_overlap
void set_operator_null_space(Teuchos::ParameterList ¶meter_list, std::unique_ptr< Epetra_MultiVector > &distributed_constant_modes, const Epetra_RowMatrix &matrix) const
const char * smoother_type
void set_parameters(Teuchos::ParameterList ¶meter_list, std::unique_ptr< Epetra_MultiVector > &distributed_constant_modes, const Epetra_RowMatrix &matrix) const
std::string block_creation_type
AdditionalData(const unsigned int block_size=1, const std::string &block_creation_type="linear", const double omega=1, const double min_diagonal=0, const unsigned int n_sweeps=1)
std::string block_creation_type
AdditionalData(const unsigned int block_size=1, const std::string &block_creation_type="linear", const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
AdditionalData(const unsigned int block_size=1, const std::string &block_creation_type="linear", const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
std::string block_creation_type
AdditionalData(const unsigned int overlap=0)
AdditionalData(const unsigned int degree=1, const double max_eigenvalue=10., const double eigenvalue_ratio=30., const double min_eigenvalue=1., const double min_diagonal=1e-12, const bool nonzero_starting=false)
AdditionalData(const unsigned int ic_fill=0, const double ic_atol=0., const double ic_rtol=1., const unsigned int overlap=0)
AdditionalData(const double ilut_drop=0., const unsigned int ilut_fill=0, const double ilut_atol=0., const double ilut_rtol=1., const unsigned int overlap=0)
AdditionalData(const unsigned int ilu_fill=0, const double ilu_atol=0., const double ilu_rtol=1., const unsigned int overlap=0)
AdditionalData(const double omega=1, const double min_diagonal=0, const unsigned int n_sweeps=1)
AdditionalData(const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)
AdditionalData(const double omega=1, const double min_diagonal=0, const unsigned int overlap=0, const unsigned int n_sweeps=1)