46 #ifndef MUELU_IFPACKSMOOTHER_HPP 47 #define MUELU_IFPACKSMOOTHER_HPP 50 #if defined(HAVE_MUELU_EPETRA) && defined(HAVE_MUELU_IFPACK) 52 #include <Teuchos_ParameterList.hpp> 53 #include <Xpetra_Matrix.hpp> 54 #include <Xpetra_MultiVectorFactory_fwd.hpp> 58 #include "MueLu_SmootherPrototype.hpp" 75 template <class Node = typename SmootherPrototype<double,int,int>::node_type>
80 #undef MUELU_IFPACKSMOOTHER_SHORT 122 IfpackSmoother(std::string
const & type, Teuchos::ParameterList
const & paramList = Teuchos::ParameterList(), LO
const &overlap=0);
156 void Apply(MultiVector& X,
const MultiVector& B,
bool InitialGuessIsZero =
false)
const;
163 RCP<SmootherPrototype>
Copy()
const;
184 void SetPrecParameters(
const Teuchos::ParameterList& list = Teuchos::ParameterList())
const;
205 template <
class Scalar,
class LocalOrdinal,
class GlobalOrdinal,
class Node>
206 RCP<MueLu::SmootherPrototype<Scalar, LocalOrdinal, GlobalOrdinal, Node> >
208 const Teuchos::ParameterList& paramList = Teuchos::ParameterList (),
209 const LocalOrdinal& overlap = 0)
211 TEUCHOS_TEST_FOR_EXCEPTION(
true,
Exceptions::RuntimeError,
"IfpackSmoother cannot be used with Scalar != double, LocalOrdinal != int, GlobalOrdinal != int");
212 TEUCHOS_UNREACHABLE_RETURN(Teuchos::null);
216 #if defined(HAVE_MUELU_EPETRA) 218 inline RCP<MueLu::SmootherPrototype<double, int, int, Xpetra::EpetraNode> >
219 GetIfpackSmoother<double, int, int, Xpetra::EpetraNode> (
const std::string& type,
const Teuchos::ParameterList& paramList,
const int& overlap) {
226 #define MUELU_IFPACKSMOOTHER_SHORT 227 #endif // HAVE_MUELU_EPETRA && HAVE_MUELU_IFPACK 228 #endif // MUELU_IFPACKSMOOTHER_HPP void Setup(Level ¤tLevel)
Set up the smoother.
void SetPrecParameters(const Teuchos::ParameterList &list=Teuchos::ParameterList()) const
virtual ~IfpackSmoother()
Destructor.
Base class for smoother prototypes.
RCP< SmootherPrototype > Copy() const
Namespace for MueLu classes and methods.
IfpackSmoother(std::string const &type, Teuchos::ParameterList const ¶mList=Teuchos::ParameterList(), LO const &overlap=0)
Constructor.
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level to an FancyOStream object.
RCP< Ifpack_Preconditioner > prec_
pointer to Ifpack solver object
std::string type_
ifpack-specific key phrase that denote smoother type
Class that holds all level-specific information.
LO overlap_
overlap when using the smoother in additive Schwarz mode
RCP< MueLu::SmootherPrototype< Scalar, LocalOrdinal, GlobalOrdinal, Node > > GetIfpackSmoother(const std::string &type="", const Teuchos::ParameterList ¶mList=Teuchos::ParameterList(), const LocalOrdinal &overlap=0)
RCP< Matrix > A_
Matrix. Not used directly, but held inside of prec_. So we have to keep an RCP pointer to it! ...
void DeclareInput(Level ¤tLevel) const
Input.
void SetParameterList(const Teuchos::ParameterList ¶mList)
size_t getNodeSmootherComplexity() const
Get a rough estimate of cost per iteration.
std::string description() const
Return a simple one-line description of this object.
Exception throws to report errors in the internal logical of the program.
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the preconditioner.
Class that encapsulates Ifpack smoothers.