51 #ifndef IFPACK2_CHEBYSHEV_DECL_HPP 52 #define IFPACK2_CHEBYSHEV_DECL_HPP 59 #include "Ifpack2_Details_Chebyshev.hpp" 66 #include "Tpetra_CrsMatrix_decl.hpp" 68 #include <type_traits> 198 template<
class MatrixType>
201 typename MatrixType::local_ordinal_type,
202 typename MatrixType::global_ordinal_type,
203 typename MatrixType::node_type>,
205 typename MatrixType::local_ordinal_type,
206 typename MatrixType::global_ordinal_type,
207 typename MatrixType::node_type> >
229 typedef typename Teuchos::ScalarTraits<scalar_type>::magnitudeType
magnitude_type;
235 typedef Tpetra::RowMatrix<scalar_type,local_ordinal_type,global_ordinal_type,node_type>
row_matrix_type;
237 static_assert(std::is_same<MatrixType, row_matrix_type>::value,
238 "Ifpack2::Chebyshev: Please use MatrixType = Tpetra::RowMatrix.");
241 typedef Tpetra::Map<local_ordinal_type,global_ordinal_type,node_type>
map_type;
248 typedef Tpetra::Vector<scalar_type,local_ordinal_type,global_ordinal_type,node_type>
vector_type;
275 explicit Chebyshev (
const Teuchos::RCP<const row_matrix_type>& A);
433 return IsInitialized_;
510 setMatrix (
const Teuchos::RCP<const row_matrix_type>& A);
545 apply (
const Tpetra::MultiVector<scalar_type,local_ordinal_type,global_ordinal_type,node_type>& X,
546 Tpetra::MultiVector<scalar_type,local_ordinal_type,global_ordinal_type,node_type>& Y,
547 Teuchos::ETransp mode = Teuchos::NO_TRANS,
548 scalar_type alpha = Teuchos::ScalarTraits<scalar_type>::one(),
549 scalar_type beta = Teuchos::ScalarTraits<scalar_type>::zero())
const;
581 applyMat (
const Tpetra::MultiVector<scalar_type,local_ordinal_type,global_ordinal_type,node_type>& X,
582 Tpetra::MultiVector<scalar_type,local_ordinal_type,global_ordinal_type,node_type>& Y,
583 Teuchos::ETransp mode = Teuchos::NO_TRANS)
const;
590 Teuchos::RCP<const Teuchos::Comm<int> >
getComm()
const;
593 Teuchos::RCP<const row_matrix_type>
getMatrix()
const;
599 Teuchos::RCP<const Tpetra::CrsMatrix<scalar_type, local_ordinal_type, global_ordinal_type, node_type> >
637 void describe(Teuchos::FancyOStream &out,
const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default)
const;
646 template <
class NewMatrixType>
friend class Chebyshev;
667 template <
typename NewMatrixType>
668 Teuchos::RCP<Chebyshev<Tpetra::RowMatrix<typename NewMatrixType::scalar_type, typename NewMatrixType::local_ordinal_type, typename NewMatrixType::global_ordinal_type, typename NewMatrixType::node_type> > >
669 clone (
const Teuchos::RCP<const NewMatrixType>& A_newnode,
670 const Teuchos::ParameterList& params)
const;
677 typedef Teuchos::ScalarTraits<typename MatrixType::scalar_type> STS;
680 typedef Tpetra::MultiVector<scalar_type, local_ordinal_type, global_ordinal_type, node_type> MV;
696 applyImpl (
const MV& X,
698 Teuchos::ETransp mode,
700 scalar_type beta)
const;
726 mutable int NumApply_;
728 double InitializeTime_;
735 mutable double ApplyTime_;
737 double ComputeFlops_;
742 mutable double ApplyFlops_;
748 template <
typename MatrixType>
749 template <
typename NewMatrixType>
750 Teuchos::RCP<Chebyshev<Tpetra::RowMatrix<typename NewMatrixType::scalar_type, typename NewMatrixType::local_ordinal_type, typename NewMatrixType::global_ordinal_type, typename NewMatrixType::node_type> > >
752 clone (
const Teuchos::RCP<const NewMatrixType>& A_newnode,
753 const Teuchos::ParameterList& params)
const 756 typedef Tpetra::RowMatrix<
typename NewMatrixType::scalar_type,
757 typename NewMatrixType::local_ordinal_type,
758 typename NewMatrixType::global_ordinal_type,
759 typename NewMatrixType::node_type> new_row_matrix_type;
762 RCP<new_prec_type> prec (
new new_prec_type (A_newnode));
763 prec->setParameters (params);
771 #endif // IFPACK2_CHEBYSHEV_DECL_HPP Mix-in interface for preconditioners that can change their matrix after construction.
Definition: Ifpack2_Details_CanChangeMatrix.hpp:93
double getApplyFlops() const
The total number of floating-point operations taken by all calls to apply().
Definition: Ifpack2_Chebyshev_def.hpp:210
Teuchos::RCP< const Tpetra::CrsMatrix< scalar_type, local_ordinal_type, global_ordinal_type, node_type > > getCrsMatrix() const
Attempt to return the matrix A as a Tpetra::CrsMatrix.
Definition: Ifpack2_Chebyshev_def.hpp:126
MatrixType::global_ordinal_type global_ordinal_type
The type of global indices in the input MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:223
MatrixType::scalar_type scalar_type
The type of the entries of the input MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:217
double getApplyTime() const
The total time spent in all calls to apply().
Definition: Ifpack2_Chebyshev_def.hpp:198
virtual void setMatrix(const Teuchos::RCP< const row_matrix_type > &A)
Change the matrix to be preconditioned.
Definition: Ifpack2_Chebyshev_def.hpp:80
void setParameters(const Teuchos::ParameterList ¶ms)
Set (or reset) parameters.
Definition: Ifpack2_Chebyshev_def.hpp:92
double getComputeTime() const
The total time spent in all calls to compute().
Definition: Ifpack2_Chebyshev_def.hpp:192
int getNumApply() const
The total number of successful calls to apply().
Definition: Ifpack2_Chebyshev_def.hpp:180
MatrixType::node_type node_type
The Node type used by the input MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:226
Tpetra::Vector< scalar_type, local_ordinal_type, global_ordinal_type, node_type > vector_type
The Tpetra::Vector specialization matching MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:248
bool isInitialized() const
Definition: Ifpack2_Chebyshev_decl.hpp:432
void compute()
(Re)compute the left scaling, and (if applicable) estimate max and min eigenvalues of D_inv * A...
Definition: Ifpack2_Chebyshev_def.hpp:292
bool hasTransposeApply() const
Whether it's possible to apply the transpose of this operator.
Definition: Ifpack2_Chebyshev_def.hpp:162
Teuchos::RCP< Chebyshev< Tpetra::RowMatrix< typename NewMatrixType::scalar_type, typename NewMatrixType::local_ordinal_type, typename NewMatrixType::global_ordinal_type, typename NewMatrixType::node_type > > > clone(const Teuchos::RCP< const NewMatrixType > &A_newnode, const Teuchos::ParameterList ¶ms) const
Clone this object to one with a different Node type.
Definition: Ifpack2_Chebyshev_decl.hpp:752
double getComputeFlops() const
The total number of floating-point operations taken by all calls to compute().
Definition: Ifpack2_Chebyshev_def.hpp:204
Diagonally scaled Chebyshev iteration for Tpetra sparse matrices.
Definition: Ifpack2_Chebyshev_decl.hpp:199
void apply(const Tpetra::MultiVector< scalar_type, local_ordinal_type, global_ordinal_type, node_type > &X, Tpetra::MultiVector< scalar_type, local_ordinal_type, global_ordinal_type, node_type > &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS, scalar_type alpha=Teuchos::ScalarTraits< scalar_type >::one(), scalar_type beta=Teuchos::ScalarTraits< scalar_type >::zero()) const
Apply the preconditioner to X, returning the result in Y.
Definition: Ifpack2_Chebyshev_def.hpp:218
Interface for all Ifpack2 preconditioners.
Definition: Ifpack2_Preconditioner.hpp:107
bool isComputed() const
Definition: Ifpack2_Chebyshev_decl.hpp:479
std::string description() const
A simple one-line description of this object.
Definition: Ifpack2_Chebyshev_def.hpp:319
void initialize()
Initialize the preconditioner.
Definition: Ifpack2_Chebyshev_def.hpp:277
Declaration of interface for preconditioners that can change their matrix after construction.
Tpetra::Map< local_ordinal_type, global_ordinal_type, node_type > map_type
The Tpetra::Map specialization matching MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:238
Teuchos::RCP< const Teuchos::Comm< int > > getComm() const
The communicator over which the matrix is distributed.
Definition: Ifpack2_Chebyshev_def.hpp:101
MatrixType::local_ordinal_type local_ordinal_type
The type of local indices in the input MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:220
int getNumInitialize() const
The total number of successful calls to initialize().
Definition: Ifpack2_Chebyshev_def.hpp:168
Tpetra::RowMatrix< scalar_type, local_ordinal_type, global_ordinal_type, node_type > row_matrix_type
The Tpetra::RowMatrix specialization matching MatrixType.
Definition: Ifpack2_Chebyshev_decl.hpp:235
double getInitializeTime() const
The total time spent in all calls to initialize().
Definition: Ifpack2_Chebyshev_def.hpp:186
MatrixType::scalar_type getLambdaMaxForApply() const
The estimate of the maximum eigenvalue used in the apply().
Definition: Ifpack2_Chebyshev_def.hpp:511
Teuchos::ScalarTraits< scalar_type >::magnitudeType magnitude_type
The type of the magnitude (absolute value) of a matrix entry.
Definition: Ifpack2_Chebyshev_decl.hpp:229
void describe(Teuchos::FancyOStream &out, const Teuchos::EVerbosityLevel verbLevel=Teuchos::Describable::verbLevel_default) const
Print the object with some verbosity level to a Teuchos::FancyOStream.
Definition: Ifpack2_Chebyshev_def.hpp:348
Teuchos::RCP< const row_matrix_type > getMatrix() const
The matrix for which this is a preconditioner.
Definition: Ifpack2_Chebyshev_def.hpp:115
Teuchos::RCP< const map_type > getDomainMap() const
The Tpetra::Map representing the domain of this operator.
Definition: Ifpack2_Chebyshev_def.hpp:136
MatrixType matrix_type
The template parameter of this class.
Definition: Ifpack2_Chebyshev_decl.hpp:214
virtual ~Chebyshev()
Destructor.
Definition: Ifpack2_Chebyshev_def.hpp:75
Preconditioners and smoothers for Tpetra sparse matrices.
Definition: Ifpack2_AdditiveSchwarz_decl.hpp:72
Teuchos::RCP< const map_type > getRangeMap() const
The Tpetra::Map representing the range of this operator.
Definition: Ifpack2_Chebyshev_def.hpp:150
int getNumCompute() const
The total number of successful calls to compute().
Definition: Ifpack2_Chebyshev_def.hpp:174
void applyMat(const Tpetra::MultiVector< scalar_type, local_ordinal_type, global_ordinal_type, node_type > &X, Tpetra::MultiVector< scalar_type, local_ordinal_type, global_ordinal_type, node_type > &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS) const
Compute Y = Op(A)*X, where Op(A) is either A, , or .
Definition: Ifpack2_Chebyshev_def.hpp:258