47 #ifndef MUELU_BLOCKEDJACOBISMOOTHER_DECL_HPP_
48 #define MUELU_BLOCKEDJACOBISMOOTHER_DECL_HPP_
62 #include "MueLu_SmootherPrototype.hpp"
113 template <class Scalar = SmootherPrototype<>::scalar_type,
114 class LocalOrdinal =
typename SmootherPrototype<Scalar>::local_ordinal_type,
115 class GlobalOrdinal =
typename SmootherPrototype<Scalar, LocalOrdinal>::global_ordinal_type,
116 class Node =
typename SmootherPrototype<Scalar, LocalOrdinal, GlobalOrdinal>::node_type>
122 #undef MUELU_BLOCKEDJACOBISMOOTHER_SHORT
170 void Apply(MultiVector& X,
const MultiVector& B,
bool InitialGuessIsZero =
false)
const;
201 std::vector<Teuchos::RCP<const SmootherBase> >
Inverse_;
221 #define MUELU_BLOCKEDJACOBISMOOTHER_SHORT
RCP< const ParameterList > GetValidParameterList() const
Input.
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level to an FancyOStream object.
std::string description() const
Return a simple one-line description of this object.
Base class for smoother prototypes.
std::vector< Teuchos::RCP< const SmootherBase > > Inverse_
vector of smoother/solver factories
void Setup(Level ¤tLevel)
Setup routine In the Setup method the Inverse_ vector is filled with the corresponding SmootherBase o...
RCP< const MapExtractorClass > domainMapExtractor_
RCP< FactoryBase > AFact_
A Factory.
std::vector< bool > bIsBlockedOperator_
vector storing whether sub-block is a blocked operator (needed for nested blocked smoothers using Thy...
Class that holds all level-specific information.
void DeclareInput(Level ¤tLevel) const
Input.
std::string type_
smoother type
Xpetra::MapExtractor< Scalar, LocalOrdinal, GlobalOrdinal, Node > MapExtractorClass
std::vector< Teuchos::RCP< const FactoryManagerBase > > FactManager_
vector of factory managers
size_t getNodeSmootherComplexity() const
Get a rough estimate of cost per iteration.
block Jacobi method for blocked matrices
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the direct solver.
BlockedJacobiSmoother()
Constructor.
RCP< const MapExtractorClass > rangeMapExtractor_
domain map extractor (from A_ generated by AFact)
RCP< Matrix > A_
internal blocked operator "A" generated by AFact_
virtual ~BlockedJacobiSmoother()
Destructor.
RCP< SmootherPrototype > Copy() const
void AddFactoryManager(RCP< const FactoryManagerBase > FactManager, int pos)
Add a factory manager.