47 #include "Teko_BlockImplicitLinearOp.hpp" 51 using Teuchos::rcpFromRef;
52 using Teuchos::rcp_dynamic_cast;
53 using Teuchos::rcp_const_cast;
56 using Thyra::ProductMultiVectorBase;
60 return (M_trans == Thyra::NOTRANS);
63 void BlockImplicitLinearOp::applyImpl(
64 const Thyra::EOpTransp M_trans,
65 const Thyra::MultiVectorBase<double> & x,
66 const Teuchos::Ptr<Thyra::MultiVectorBase<double> > & y,
71 TEUCHOS_TEST_FOR_EXCEPTION(M_trans!=Thyra::NOTRANS, std::runtime_error,
72 "Linear operators of inherited type BlockImplicitLinearOp " 73 "cannot handle conjugation (yet!)");
76 RCP<const ProductMultiVectorBase<double> > src =
77 rcp_dynamic_cast<
const ProductMultiVectorBase<double> >(rcpFromRef(x));
78 BlockedMultiVector srcX = rcp_const_cast<ProductMultiVectorBase<double> >(src);
81 BlockedMultiVector destY =
82 rcp_dynamic_cast<ProductMultiVectorBase<double> >(rcpFromPtr(y));
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 bool opSupportedImpl(const Thyra::EOpTransp M_trans) const
Functions required by Thyra::LinearOpBase.