47#ifndef MUELU_BLOCKEDGAUSSSEIDELSMOOTHER_DECL_HPP_
48#define MUELU_BLOCKEDGAUSSSEIDELSMOOTHER_DECL_HPP_
52#include <Teuchos_ParameterList.hpp>
54#include <Xpetra_Matrix_fwd.hpp>
55#include <Xpetra_BlockedCrsMatrix_fwd.hpp>
56#include <Xpetra_ReorderedBlockedCrsMatrix_fwd.hpp>
57#include <Xpetra_MultiVectorFactory_fwd.hpp>
58#include <Xpetra_MapExtractor_fwd.hpp>
62#include "MueLu_SmootherPrototype.hpp"
111 template <class Scalar = SmootherPrototype<>::scalar_type,
120#undef MUELU_BLOCKEDGAUSSSEIDELSMOOTHER_SHORT
167 void Apply(MultiVector& X,
const MultiVector& B,
bool InitialGuessIsZero =
false)
const;
170 RCP<SmootherPrototype>
Copy()
const;
199 std::vector<Teuchos::RCP<const SmootherBase> >
Inverse_;
219#define MUELU_BLOCKEDGAUSSSEIDELSMOOTHER_SHORT
MueLu::DefaultLocalOrdinal LocalOrdinal
MueLu::DefaultGlobalOrdinal GlobalOrdinal
block Gauss-Seidel method for blocked matrices
void DeclareInput(Level ¤tLevel) const
Input.
std::vector< Teuchos::RCP< const FactoryManagerBase > > FactManager_
vector of factory managers
std::vector< Teuchos::RCP< const SmootherBase > > Inverse_
vector of smoother/solver factories
RCP< const MapExtractorClass > domainMapExtractor_
domain map extractor (from A_ generated by AFact)
std::string description() const
Return a simple one-line description of this object.
virtual ~BlockedGaussSeidelSmoother()
Destructor.
size_t getNodeSmootherComplexity() const
Get a rough estimate of cost per iteration.
void AddFactoryManager(RCP< const FactoryManagerBase > FactManager, int pos)
Add a factory manager.
RCP< SmootherPrototype > Copy() const
RCP< const ParameterList > GetValidParameterList() const
Input.
RCP< const MapExtractorClass > rangeMapExtractor_
range map extractor (from A_ generated by AFact)
std::vector< bool > bIsBlockedOperator_
vector storing whether sub-block is a blocked operator (needed for nested blocked smoothers using Thy...
BlockedGaussSeidelSmoother()
Constructor.
Xpetra::MapExtractor< Scalar, LocalOrdinal, GlobalOrdinal, Node > MapExtractorClass
std::string type_
smoother type
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the direct solver. Solves the linear system AX=B using the constructed solver.
RCP< FactoryBase > AFact_
A Factory.
RCP< Matrix > A_
internal blocked operator "A" generated by AFact_
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level verbLevel to an FancyOStream object out.
void Setup(Level ¤tLevel)
Setup routine In the Setup method the Inverse_ vector is filled with the corresponding SmootherBase o...
Class that holds all level-specific information.
Base class for smoother prototypes.
GlobalOrdinal global_ordinal_type
LocalOrdinal local_ordinal_type
Namespace for MueLu classes and methods.