47 #include "Teko_MultPreconditionerFactory.hpp"
54 Teko::BlockedMultiVector &y,
const double ,
55 const double )
const {
62 Teko::MultiVector MOne_r = Teko::deepcopy(r);
63 Teko::MultiVector t = Teko::deepcopy(r);
64 Teko::MultiVector w = Teko::toMultiVector(y);
66 Teko::applyOp(M1_, r, MOne_r);
67 Teko::applyOp(A_, MOne_r, t);
68 Teko::update(1., r, -1., t);
69 Teko::applyOp(M2_, t, w);
70 Teko::update(1., MOne_r, 1., w);
74 MultPreconditionerFactory ::MultPreconditionerFactory(
75 const RCP<const Teko::BlockPreconditionerFactory> &FirstFactory,
76 const RCP<const Teko::BlockPreconditionerFactory> &SecondFactory)
77 : FirstFactory_(FirstFactory), SecondFactory_(SecondFactory) {}
79 MultPreconditionerFactory::MultPreconditionerFactory() {}
85 FirstFactory_->buildPreconditionerState());
87 SecondFactory_->buildPreconditionerState());
96 TEUCHOS_ASSERT(MyState != 0);
98 Teko::LinearOp M1 = FirstFactory_->buildPreconditionerOperator(blockOp, *MyState->StateOne_);
99 Teko::LinearOp M2 = SecondFactory_->buildPreconditionerOperator(blockOp, *MyState->StateTwo_);
129 std::string aStr =
"", bStr =
"";
132 aStr = pl.get<std::string>(
"Preconditioner A");
133 bStr = pl.get<std::string>(
"Preconditioner B");
135 RCP<const Teuchos::ParameterList> aSettings = invLib->getParameterList(aStr);
136 RCP<const Teuchos::ParameterList> bSettings = invLib->getParameterList(bStr);
139 std::string aType = aSettings->get<std::string>(
"Preconditioner Type");
140 RCP<Teko::PreconditionerFactory> precA =
142 aType, aSettings->sublist(
"Preconditioner Settings"), invLib);
145 std::string bType = bSettings->get<std::string>(
"Preconditioner Type");
146 RCP<Teko::PreconditionerFactory> precB =
148 bType, bSettings->sublist(
"Preconditioner Settings"), invLib);
virtual Teuchos::RCP< Teko::PreconditionerState > buildPreconditionerState() const
Build the MultPrecondState object.
virtual void initializeFromParameterList(const Teuchos::ParameterList &pl)
Initialize from a parameter list.
Abstract class which block preconditioner factories in Teko should be based on.
An implementation of a state object for block preconditioners.
Teuchos::RCP< const InverseLibrary > getInverseLibrary() const
Get the inverse library used by this preconditioner factory.
Teko::LinearOp buildPreconditionerOperator(Teko::BlockedLinearOp &blo, Teko::BlockPreconditionerState &state) const
Function inherited from Teko::BlockPreconditionerFactory.
static Teuchos::RCP< PreconditionerFactory > buildPreconditionerFactory(const std::string &name, const Teuchos::ParameterList &settings, const Teuchos::RCP< const InverseLibrary > &invLib=Teuchos::null)
Builder function for creating preconditioner factories (yes this is a factory factory).
virtual void implicitApply(const Teko::BlockedMultiVector &r, Teko::BlockedMultiVector &y, const double alpha=1.0, const double beta=0.0) const
Perform a matrix vector multiply with this implicitly defined blocked operator.