47 #ifndef __Teko_BlockLowerTriInverseOp_hpp__
48 #define __Teko_BlockLowerTriInverseOp_hpp__
51 #include "Teko_BlockImplicitLinearOp.hpp"
99 virtual void implicitApply(
const BlockedMultiVector &x, BlockedMultiVector &y,
100 const double alpha = 1.0,
const double beta = 0.0)
const;
103 virtual void describe(Teuchos::FancyOStream &out_arg,
104 const Teuchos::EVerbosityLevel verbLevel)
const;
108 const BlockedLinearOp
L_;
111 Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> >
113 Teuchos::RCP<const Thyra::ProductVectorSpaceBase<double> >
117 mutable BlockedMultiVector srcScrap_;
118 mutable BlockedMultiVector dstScrap_;
126 inline LinearOp createBlockLowerTriInverseOp(BlockedLinearOp &U,
127 const std::vector<LinearOp> &invDiag) {
131 inline LinearOp createBlockLowerTriInverseOp(BlockedLinearOp &L,
132 const std::vector<LinearOp> &invDiag,
133 const std::string &str) {
134 Teuchos::RCP<Thyra::LinearOpBase<double> > result =
135 Teuchos::rcp(
new BlockLowerTriInverseOp(L, invDiag));
136 result->setObjectLabel(str);
virtual void implicitApply(const BlockedMultiVector &x, BlockedMultiVector &y, const double alpha=1.0, const double beta=0.0) const =0
Perform a matrix vector multiply with this implicitly defined blocked operator.
virtual VectorSpace domain() const
Domain space of this operator.
const BlockedLinearOp L_
operator
Teuchos::RCP< const Thyra::ProductVectorSpaceBase< double > > productDomain_
Domain vector space.
virtual void implicitApply(const BlockedMultiVector &x, BlockedMultiVector &y, const double alpha=1.0, const double beta=0.0) const
Perform a matrix vector multiply with this operator.
std::vector< LinearOp > invDiag_
(Approximate) Inverses of the diagonal operators
A virtual class that simplifies the construction of custom operators.
This linear operator computes the inverse of a lower triangular matrix.
Teuchos::RCP< const Thyra::ProductVectorSpaceBase< double > > productRange_
Range vector space.
virtual VectorSpace range() const
Range space of this operator.