47 #include "Teko_BlockImplicitLinearOp.hpp"
52 using Teuchos::rcp_const_cast;
53 using Teuchos::rcp_dynamic_cast;
54 using Teuchos::rcpFromRef;
56 using Thyra::ProductMultiVectorBase;
59 const BlockedMultiVector& x, BlockedMultiVector& y,
60 const double alpha,
const double beta)
const {
61 TEUCHOS_TEST_FOR_EXCEPTION(M_trans != Thyra::NOTRANS, std::runtime_error,
62 "Linear operators of inherited type BlockImplicitLinearOp "
63 "cannot handle conjugation (yet!)");
70 return (M_trans == Thyra::NOTRANS);
73 void BlockImplicitLinearOp::applyImpl(
const Thyra::EOpTransp M_trans,
74 const Thyra::MultiVectorBase<double>& x,
75 const Teuchos::Ptr<Thyra::MultiVectorBase<double> >& y,
76 const double alpha,
const double beta)
const {
78 RCP<const ProductMultiVectorBase<double> > src =
79 rcp_dynamic_cast<
const ProductMultiVectorBase<double> >(rcpFromRef(x));
80 BlockedMultiVector srcX = rcp_const_cast<ProductMultiVectorBase<double> >(src);
83 BlockedMultiVector destY = 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.