42 #ifndef THYRA_MULTI_VECTOR_ADAPTER_BASE_DEF_HPP 43 #define THYRA_MULTI_VECTOR_ADAPTER_BASE_DEF_HPP 45 #include "Thyra_MultiVectorAdapterBase_decl.hpp" 46 #include "Thyra_ScalarProdVectorSpaceBase.hpp" 47 #include "Thyra_ScalarProdBase.hpp" 56 template<
class Scalar>
57 RCP<const VectorSpaceBase<Scalar> >
60 return rangeScalarProdVecSpc();
64 template<
class Scalar>
65 RCP<const VectorSpaceBase<Scalar> >
68 return domainScalarProdVecSpc();
75 template<
class Scalar>
78 if (ScalarTraits<Scalar>::isComplex)
84 template<
class Scalar>
102 RCP<const ScalarProdVectorSpaceBase<Scalar> > scalarProdVecSpc =
104 ? domainScalarProdVecSpc()
105 : rangeScalarProdVecSpc() );
106 RCP<const ScalarProdBase<Scalar> > scalarProd = scalarProdVecSpc->getScalarProd();
107 if (scalarProd->isEuclidean()) {
109 this->euclideanApply(M_trans, X, Y, alpha, beta);
113 RCP<MultiVectorBase<Scalar> > T = createMembers(X.
range(), X.
domain());
114 ::Thyra::apply(*scalarProd->getLinearOp(),
NOTRANS, X, T.ptr());
116 this->euclideanApply(M_trans, *T, Y, alpha, beta);
124 #endif // THYRA_MULTI_VECTOR_ADAPTER_BASE_DEF_HPP virtual RCP< const VectorSpaceBase< Scalar > > range() const =0
Return a smart pointer for the range space for this operator.
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
Use the non-transposed operator.
EOpTransp real_trans(EOpTransp transp)
Return NOTRANS or TRANS for real scalar valued operators and this also is used for determining struct...
Use the transposed operator with complex-conjugate clements (same as TRANS for real scalar types)...
bool opSupportedImpl(EOpTransp M_trans) const
RCP< const VectorSpaceBase< Scalar > > range() const
Returns this->rangeScalarProdVecSpc()
Interface for a collection of column vectors called a multi-vector.
void applyImpl(const EOpTransp M_trans, const MultiVectorBase< Scalar > &X, const Ptr< MultiVectorBase< Scalar > > &Y, const Scalar alpha, const Scalar beta) const
virtual RCP< const VectorSpaceBase< Scalar > > domain() const =0
Return a smart pointer for the domain space for this operator.
RCP< const VectorSpaceBase< Scalar > > domain() const
Returns this->domainScalarProdVecSpc()