42#ifndef THYRA_DIAGONAL_LINEAR_OP_DEF_HPP
43#define THYRA_DIAGONAL_LINEAR_OP_DEF_HPP
46#include "Thyra_DefaultDiagonalLinearOp_decl.hpp"
47#include "Thyra_MultiVectorStdOps.hpp"
48#include "Thyra_VectorStdOps.hpp"
49#include "Thyra_VectorBase.hpp"
50#include "Thyra_AssertOp.hpp"
99 initialize(createMember(space));
103template<
class Scalar>
108 diag_.initialize(diag);
112template<
class Scalar>
117 diag_.initialize(diag);
121template<
class Scalar>
124 diag_.uninitialize();
131template<
class Scalar>
134 return diag_.isConst();
138template<
class Scalar>
142 return diag_.getNonconstObj();
146template<
class Scalar>
150 return diag_.getConstObj();
157template<
class Scalar>
161 return diag_.getConstObj()->space();
165template<
class Scalar>
169 return diag_.getConstObj()->space();
173template<
class Scalar>
187template<
class Scalar>
194template<
class Scalar>
207 "DefaultDiagonalLinearOp<Scalar>::apply(...)",*
this, M_trans, X, &*Y
213 if( beta != ST::one() ) scale<Scalar>(beta, Y);
219 for (
Ordinal col_j = 0; col_j < m; ++col_j) {
224 ele_wise_prod( alpha, *diag_.getConstObj(), *x, y.
ptr() );
227 ele_wise_conj_prod( alpha, *diag_.getConstObj(), *x, y.
ptr() );
231 ele_wise_prod( alpha, *diag_.getConstObj(), *x, y.
ptr() );
Default concrete LinearOpBase subclass for diagonal linear operators.
void uninitialize()
Uninitialize.
RCP< const VectorSpaceBase< Scalar > > range() const
Returns this->getDiag()->space().
void initialize(const RCP< const VectorSpaceBase< Scalar > > &space)
Initialize given a vector space which allocates a vector internally.
RCP< const VectorBase< Scalar > > getDiag() const
RCP< VectorBase< Scalar > > getNonconstDiag()
RCP< const VectorSpaceBase< Scalar > > domain() const
Returns this->getDiag()->space().
DefaultDiagonalLinearOp()
Constructs to uninitialized.
RCP< const LinearOpBase< Scalar > > clone() const
bool opSupportedImpl(EOpTransp M_trans) const
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.
Interface for a collection of column vectors called a multi-vector.
RCP< const VectorBase< Scalar > > col(Ordinal j) const
Calls colImpl().
Abstract interface for finite-dimensional dense vectors.
Abstract interface for objects that represent a space for vectors.
#define TEUCHOS_TEST_FOR_EXCEPT(throw_exception_test)
#define THYRA_ASSERT_LINEAR_OP_MULTIVEC_APPLY_SPACES(FUNC_NAME, M, M_T, X, Y)
This is a very useful macro that should be used to validate that the spaces for the multi-vector vers...
EOpTransp
Enumeration for determining how a linear operator is applied. `*.
Teuchos::Ordinal Ordinal
Type for the dimension of a vector space. `*.
@ TRANS
Use the transposed operator.
@ NOTRANS
Use the non-transposed operator.
TEUCHOS_DEPRECATED RCP< T > rcp(T *p, Dealloc_T dealloc, bool owns_mem)