42 #ifndef THYRA_LINEAR_OP_SCALAR_PROD_DEF_HPP 43 #define THYRA_LINEAR_OP_SCALAR_PROD_DEF_HPP 45 #include "Thyra_LinearOpScalarProd_decl.hpp" 46 #include "Thyra_ScalarProdBase.hpp" 47 #include "Thyra_LinearOpBase.hpp" 48 #include "Thyra_MultiVectorStdOps.hpp" 49 #include "Thyra_MultiVectorBase.hpp" 50 #include "Thyra_VectorSpaceBase.hpp" 59 template<
class Scalar>
64 template<
class Scalar>
68 this->initialize(op_in);
72 template<
class Scalar>
78 TEUCHOS_TEST_FOR_EXCEPT(is_null(op_in));
84 template<
class Scalar>
89 if (!is_null(op_out)) *op_out = op_;
97 template<
class Scalar>
104 template<
class Scalar>
107 const ArrayView<Scalar> &scalarProds_out
110 Teuchos::RCP<MultiVectorBase<Scalar> >
112 Thyra::apply(*op_,
NOTRANS,Y, T.ptr());
113 dots(X, *T, scalarProds_out);
117 template<
class Scalar>
118 RCP<const LinearOpBase<Scalar> >
128 #endif // THYRA_LINEAR_OP_SCALAR_PROD_DEF_HPP virtual bool isEuclideanImpl() const
Returns false.
virtual RCP< const VectorSpaceBase< Scalar > > range() const =0
Return a smart pointer for the range space for this operator.
Use the non-transposed operator.
RCP< const LinearOpBase< Scalar > > getLinearOpImpl() const
Interface for a collection of column vectors called a multi-vector.
void initialize(const RCP< const LinearOpBase< Scalar > > &op)
Base class for all linear operators.
void uninitialize(const Ptr< RCP< const LinearOpBase< Scalar > > > &op=Teuchos::null)
virtual RCP< const VectorSpaceBase< Scalar > > domain() const =0
Return a smart pointer for the domain space for this operator.
void scalarProdsImpl(const MultiVectorBase< Scalar > &X, const MultiVectorBase< Scalar > &Y, const ArrayView< Scalar > &scalarProds_out) const