53 #ifndef MUELU_BRAESSSARAZINSMOOTHER_DECL_HPP_ 54 #define MUELU_BRAESSSARAZINSMOOTHER_DECL_HPP_ 57 #include <Teuchos_ParameterList.hpp> 60 #include <Xpetra_MapExtractor_fwd.hpp> 61 #include <Xpetra_MultiVectorFactory_fwd.hpp> 62 #include <Xpetra_CrsMatrixWrap_fwd.hpp> 63 #include <Xpetra_Matrix_fwd.hpp> 67 #include "MueLu_SmootherPrototype.hpp" 86 template <class Scalar = SmootherPrototype<>::scalar_type,
94 #undef MUELU_BRAESSSARAZINSMOOTHER_SHORT 134 void Apply(MultiVector& X,
const MultiVector& B,
bool InitialGuessIsZero =
false)
const;
137 RCP<SmootherPrototype>
Copy()
const;
179 #define MUELU_BRAESSSARAZINSMOOTHER_SHORT
RCP< Matrix > A00_
matrices
BraessSarazin smoother for 2x2 block matrices.
RCP< Vector > D_
Inverse to approximation to block (0,0). Here, D_ = omega*inv(diag(A(0,0)))
RCP< const MapExtractorClass > domainMapExtractor_
domain map extractor (from A_ generated by AFact)
std::string type_
smoother type
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level to an FancyOStream object.
RCP< Matrix > A_
block operator
Base class for smoother prototypes.
RCP< SmootherPrototype > Copy() const
Namespace for MueLu classes and methods.
LocalOrdinal local_ordinal_type
RCP< Matrix > S_
Schur complement.
RCP< const MapExtractorClass > rangeMapExtractor_
range map extractor (from A_ generated by AFact)
RCP< Matrix > A01_
Block (0,1) [typically, pressure gradient operator].
void DeclareInput(Level ¤tLevel) const
Input.
Class that holds all level-specific information.
virtual ~BraessSarazinSmoother()
Destructor.
RCP< Matrix > A11_
Block (1,1) [typically, pressure stabilization term or null block].
RCP< const ParameterList > GetValidParameterList() const
Input.
BraessSarazinSmoother()
Constructor.
RCP< const FactoryBase > AFact_
A Factory.
GlobalOrdinal global_ordinal_type
void AddFactoryManager(RCP< const FactoryManagerBase > FactManager, int pos=0)
Add a factory manager for BraessSarazin internal SchurComplement handling.
void Setup(Level ¤tLevel)
Setup routine.
RCP< Matrix > A10_
Block (1,0) [typically, divergence operator].
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the Braess Sarazin smoother.
Xpetra::MapExtractor< Scalar, LocalOrdinal, GlobalOrdinal, Node > MapExtractorClass
std::string description() const
Return a simple one-line description of this object.
Teuchos::RCP< SmootherBase > smoo_
Smoother for SchurComplement equation.
RCP< const FactoryManagerBase > FactManager_
Factory manager for creating the Schur Complement.