47 #include "Teko_ImplicitLinearOp.hpp"
52 using Teuchos::rcp_const_cast;
53 using Teuchos::rcp_dynamic_cast;
54 using Teuchos::rcpFromRef;
56 using Thyra::MultiVectorBase;
59 return (M_trans == Thyra::NOTRANS);
62 void ImplicitLinearOp::applyImpl(
const Thyra::EOpTransp M_trans,
63 const Thyra::MultiVectorBase<double>& x,
64 const Teuchos::Ptr<Thyra::MultiVectorBase<double> >& y,
65 const double alpha,
const double beta)
const {
66 TEUCHOS_TEST_FOR_EXCEPTION(M_trans != Thyra::NOTRANS, std::runtime_error,
67 "Linear operators of inherited type Teko::ImplicitLinearOp "
68 "cannot handle conjugation (yet!)");
70 MultiVector srcX = rcp_const_cast<MultiVectorBase<double> >(rcpFromRef(x));
71 MultiVector destY = rcp_dynamic_cast<MultiVectorBase<double> >(rcpFromRef(*y));
virtual bool opSupportedImpl(const Thyra::EOpTransp M_trans) const
Functions required by Thyra::LinearOpBase.
virtual void implicitApply(const MultiVector &x, MultiVector &y, const double alpha=1.0, const double beta=0.0) const =0
Perform a matrix vector multiply with this implicitly defined blocked operator.