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
MueLu::DefaultLocalOrdinal LocalOrdinal
MueLu::DefaultGlobalOrdinal GlobalOrdinal
BraessSarazin smoother for 2x2 block matrices.
RCP< Matrix > A01_
Block (0,1) [typically, pressure gradient operator].
RCP< Matrix > A10_
Block (1,0) [typically, divergence operator].
RCP< const ParameterList > GetValidParameterList() const
Input.
RCP< const MapExtractor > domainMapExtractor_
domain map extractor (from A_ generated by AFact)
RCP< Matrix > A_
block operator
RCP< const FactoryBase > AFact_
A Factory.
std::string type_
smoother type
RCP< Matrix > A11_
Block (1,1) [typically, pressure stabilization term or null block].
std::string description() const
Return a simple one-line description of this object.
void Apply(MultiVector &X, const MultiVector &B, bool InitialGuessIsZero=false) const
Apply the Braess Sarazin smoother.
RCP< Vector > D_
Inverse to approximation to block (0,0). Here, D_ = omega*inv(diag(A(0,0)))
RCP< const MapExtractor > rangeMapExtractor_
range map extractor (from A_ generated by AFact)
void DeclareInput(Level ¤tLevel) const
Input.
RCP< SmootherPrototype > Copy() const
RCP< Matrix > A00_
matrices
size_t getNodeSmootherComplexity() const
Get a rough estimate of cost per iteration.
void Setup(Level ¤tLevel)
Setup routine.
Teuchos::RCP< SmootherBase > smoo_
Smoother for SchurComplement equation.
RCP< const FactoryManagerBase > FactManager_
Factory manager for creating the Schur Complement.
RCP< Matrix > S_
Schur complement.
void print(Teuchos::FancyOStream &out, const VerbLevel verbLevel=Default) const
Print the object with some verbosity level to an FancyOStream object.
void AddFactoryManager(RCP< const FactoryManagerBase > FactManager, int pos=0)
Add a factory manager for BraessSarazin internal SchurComplement handling.
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.