47 #ifndef MUELU_BRAESSSARAZINSMOOTHER_DECL_HPP_ 48 #define MUELU_BRAESSSARAZINSMOOTHER_DECL_HPP_ 51 #include <Teuchos_ParameterList.hpp> 54 #include <Xpetra_MapExtractor_fwd.hpp> 55 #include <Xpetra_MultiVectorFactory_fwd.hpp> 56 #include <Xpetra_CrsMatrixWrap_fwd.hpp> 57 #include <Xpetra_Matrix_fwd.hpp> 61 #include "MueLu_SmootherPrototype.hpp" 80 template <class Scalar = SmootherPrototype<>::scalar_type,
86 #undef MUELU_BRAESSSARAZINSMOOTHER_SHORT 121 void Apply(MultiVector& X,
const MultiVector& B,
bool InitialGuessIsZero =
false)
const;
124 RCP<SmootherPrototype>
Copy()
const;
147 std::string
type_ =
"Braess Sarazin";
153 RCP<Matrix>
A_ = Teuchos::null;
172 #define MUELU_BRAESSSARAZINSMOOTHER_SHORT
RCP< SmootherPrototype > Copy() const
MueLu::DefaultLocalOrdinal LocalOrdinal
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)))
std::string type_
smoother type
RCP< const MapExtractor > domainMapExtractor_
domain map extractor (from A_ generated by AFact)
RCP< Matrix > A_
block operator
Base class for smoother prototypes.
Namespace for MueLu classes and methods.
LocalOrdinal local_ordinal_type
RCP< Matrix > S_
Schur complement.
std::string description() const
Return a simple one-line description of this object.
RCP< Matrix > A01_
Block (0,1) [typically, pressure gradient operator].
MueLu::DefaultGlobalOrdinal GlobalOrdinal
Class that holds all level-specific information.
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the Braess Sarazin smoother.
RCP< Matrix > A11_
Block (1,1) [typically, pressure stabilization term or null block].
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level to an FancyOStream object.
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 DeclareInput(Level ¤tLevel) const
Input.
void Setup(Level ¤tLevel)
Setup routine.
RCP< Matrix > A10_
Block (1,0) [typically, divergence operator].
size_t getNodeSmootherComplexity() const
Get a rough estimate of cost per iteration.
Teuchos::RCP< SmootherBase > smoo_
Smoother for SchurComplement equation.
RCP< const ParameterList > GetValidParameterList() const
Input.
RCP< const MapExtractor > rangeMapExtractor_
range map extractor (from A_ generated by AFact)
RCP< const FactoryManagerBase > FactManager_
Factory manager for creating the Schur Complement.