46#ifndef MUELU_SHIFTEDLAPLACIANOPERATOR_DECL_HPP
47#define MUELU_SHIFTEDLAPLACIANOPERATOR_DECL_HPP
51#ifdef HAVE_MUELU_TPETRA
52#include <Tpetra_Operator.hpp>
53#include <Tpetra_MultiVector_decl.hpp>
56#include "MueLu_Utilities.hpp"
65 template <class Scalar = Tpetra::Operator<>::scalar_type,
66 class LocalOrdinal =
typename Tpetra::Operator<Scalar>::local_ordinal_type,
67 class GlobalOrdinal =
typename Tpetra::Operator<Scalar, LocalOrdinal>::global_ordinal_type,
68 class Node =
typename Tpetra::Operator<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>
70 :
public Tpetra::Operator<Scalar,LocalOrdinal,GlobalOrdinal,Node>
73 typedef Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>
Matrix;
74 typedef Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node>
CrsMatrix;
75 typedef Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>
MV;
76 typedef Tpetra::Operator<Scalar,LocalOrdinal,GlobalOrdinal,Node>
OP;
124 Teuchos::RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> >
getDomainMap()
const;
127 Teuchos::RCP<const Tpetra::Map<LocalOrdinal,GlobalOrdinal,Node> >
getRangeMap()
const;
135 void apply(
const Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& X,
136 Tpetra::MultiVector<Scalar,LocalOrdinal,GlobalOrdinal,Node>& Y,
137 Teuchos::ETransp mode = Teuchos::NO_TRANS,
138 Scalar alpha = Teuchos::ScalarTraits<Scalar>::one(),
139 Scalar beta = Teuchos::ScalarTraits<Scalar>::one())
const;
146 RCP<MueLu::Hierarchy<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
Hierarchy_;
147 RCP< Xpetra::Matrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
R_,
P_,
A_;
148 RCP< Tpetra::CrsMatrix<Scalar,LocalOrdinal,GlobalOrdinal,Node> >
Ac_;
MueLu::DefaultLocalOrdinal LocalOrdinal
MueLu::DefaultScalar Scalar
MueLu::DefaultGlobalOrdinal GlobalOrdinal
Provides methods to build a multigrid hierarchy and apply multigrid cycles.
Wraps an existing MueLu::Hierarchy as a Tpetra::Operator, with an optional two-level correction....
RCP< Teuchos::ParameterList > BelosList_
RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > P_
ShiftedLaplacianOperator(const RCP< MueLu::Hierarchy< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &H)
Constructor.
Teuchos::RCP< const Tpetra::Map< LocalOrdinal, GlobalOrdinal, Node > > getRangeMap() const
Returns the Tpetra::Map object associated with the range of this operator.
RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > A_
bool hasTransposeApply() const
Indicates whether this operator supports applying the adjoint operator.
Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > MV
RCP< Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > R_
ShiftedLaplacianOperator(const RCP< MueLu::Hierarchy< Scalar, LocalOrdinal, GlobalOrdinal, Node > > &H, const RCP< Matrix > A, int cycles, int iters, int option, double tol)
Auxiliary Constructor.
RCP< MueLu::Hierarchy< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Hierarchy_
MueLu::Utilities< Scalar, LocalOrdinal, GlobalOrdinal, Node > MUtils
Teuchos::RCP< const Tpetra::Map< LocalOrdinal, GlobalOrdinal, Node > > getDomainMap() const
Returns the Tpetra::Map object associated with the domain of this operator.
RCP< Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > > Ac_
Xpetra::Matrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > Matrix
Tpetra::CrsMatrix< Scalar, LocalOrdinal, GlobalOrdinal, Node > CrsMatrix
virtual ~ShiftedLaplacianOperator()
Destructor.
void apply(const Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &X, Tpetra::MultiVector< Scalar, LocalOrdinal, GlobalOrdinal, Node > &Y, Teuchos::ETransp mode=Teuchos::NO_TRANS, Scalar alpha=Teuchos::ScalarTraits< Scalar >::one(), Scalar beta=Teuchos::ScalarTraits< Scalar >::one()) const
Returns in Y the result of a Tpetra::Operator applied to a Tpetra::MultiVector X.
Tpetra::Operator< Scalar, LocalOrdinal, GlobalOrdinal, Node > OP
Namespace for MueLu classes and methods.