35 #ifndef _BLAZE_MATH_EXPRESSIONS_DVECSVECSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DVECSVECSUBEXPR_H_
85 template<
typename VT1
124 template<
typename VT >
125 struct UseSMPAssign {
126 enum { value = ( !VT1::smpAssignable || !VT2::smpAssignable ) };
153 enum { vectorizable = 0 };
156 enum { smpAssignable = 0 };
221 template<
typename T >
224 (
rhs_.canAlias( alias ) );
234 template<
typename T >
236 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
258 template<
typename VT >
288 template<
typename VT >
317 template<
typename VT >
346 template<
typename VT >
375 template<
typename VT >
410 template<
typename VT >
411 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
418 if( !IsComputation<VT1>::value &&
isSame( ~lhs, rhs.lhs_ ) ) {
443 template<
typename VT >
444 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
475 template<
typename VT >
476 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
508 template<
typename VT >
509 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
541 template<
typename VT >
542 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
609 template<
typename T1
612 inline const DVecSVecSubExpr<T1,T2,TF>
617 if( (~lhs).
size() != (~rhs).
size() )
618 throw std::invalid_argument(
"Vector sizes do not match" );
646 template<
typename T1
650 inline const typename AddExprTrait< DVecSVecSubExpr<T1,T2,TF>, T3 >::Type
651 operator+(
const DVecSVecSubExpr<T1,T2,TF>& lhs,
const DenseVector<T3,TF>& rhs )
655 return ( lhs.leftOperand() + (~rhs) ) - lhs.rightOperand();
674 template<
typename T1
678 inline const typename SubExprTrait< DVecSVecSubExpr<T1,T2,TF>, T3 >::Type
679 operator-(
const DVecSVecSubExpr<T1,T2,TF>& lhs,
const DenseVector<T3,TF>& rhs )
683 return ( lhs.leftOperand() - (~rhs) ) - lhs.rightOperand();
699 template<
typename VT1,
typename VT2,
typename VT3 >
705 typedef typename SelectType< IsDenseVector<VT1>::value && IsColumnVector<VT1>::value &&
706 IsSparseVector<VT2>::value && IsColumnVector<VT2>::value &&
707 IsDenseVector<VT3>::value && IsColumnVector<VT3>::value
708 ,
typename DVecSVecSubExprTrait< typename DVecDVecAddExprTrait<VT1,VT3>::Type, VT2 >::Type
709 , INVALID_TYPE >::Type Type;
719 template<
typename VT1,
typename VT2,
typename VT3 >
725 typedef typename SelectType< IsDenseVector<VT1>::value && IsRowVector<VT1>::value &&
726 IsSparseVector<VT2>::value && IsRowVector<VT2>::value &&
727 IsDenseVector<VT3>::value && IsRowVector<VT3>::value
728 ,
typename TDVecTSVecSubExprTrait< typename TDVecTDVecAddExprTrait<VT1,VT3>::Type, VT2 >::Type
729 , INVALID_TYPE >::Type Type;
739 template<
typename VT1,
typename VT2,
typename VT3 >
745 typedef typename SelectType< IsDenseVector<VT1>::value && IsColumnVector<VT1>::value &&
746 IsSparseVector<VT2>::value && IsColumnVector<VT2>::value &&
747 IsDenseVector<VT3>::value && IsColumnVector<VT3>::value
748 ,
typename DVecSVecSubExprTrait< typename DVecDVecSubExprTrait<VT1,VT3>::Type, VT2 >::Type
749 , INVALID_TYPE >::Type Type;
759 template<
typename VT1,
typename VT2,
typename VT3 >
765 typedef typename SelectType< IsDenseVector<VT1>::value && IsRowVector<VT1>::value &&
766 IsSparseVector<VT2>::value && IsRowVector<VT2>::value &&
767 IsDenseVector<VT3>::value && IsRowVector<VT3>::value
768 ,
typename TDVecTSVecSubExprTrait< typename TDVecTDVecSubExprTrait<VT1,VT3>::Type, VT2 >::Type
769 , INVALID_TYPE >::Type Type;
779 template<
typename VT1,
typename VT2,
bool TF,
bool AF >
784 typedef typename SubExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
785 ,
typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
ResultType::ElementType ElementType
Resulting element type.
Definition: DVecSVecSubExpr.h:136
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
VT2::ReturnType RN2
Return type of the right-hand side sparse vector expression.
Definition: DVecSVecSubExpr.h:97
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
Header file for the subtraction trait.
VT1::ResultType RT1
Result type of the left-hand side dense vector expression.
Definition: DVecSVecSubExpr.h:94
RightOperand rightOperand() const
Returns the right-hand side sparse vector operand.
Definition: DVecSVecSubExpr.h:210
void smpSubAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:152
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:199
void smpMultAssign(DenseVector< VT1, TF1 > &lhs, const Vector< VT2, TF2 > &rhs)
Default implementation of the SMP multiplication assignment of a vector to a dense vector...
Definition: DenseVector.h:179
RightOperand rhs_
Right-hand side sparse vector of the subtraction expression.
Definition: DVecSVecSubExpr.h:243
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DVecSVecSubExpr.h:139
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2408
Header file for the IsRowVector type trait.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:251
Header file for the DenseVector base class.
Header file for the AddExprTrait class template.
DVecSVecSubExpr< VT1, VT2, TF > This
Type of this DVecSVecSubExpr instance.
Definition: DVecSVecSubExpr.h:133
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DVecSVecSubExpr.h:134
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:690
Header file for the Computation base class.
bool isSame(const Matrix< MT1, SO1 > &a, const Matrix< MT2, SO2 > &b)
Returns whether the two given matrices represent the same observable state.
Definition: Matrix.h:409
Constraint on the data type.
Expression object for dense vector-sparse vector subtractions.The DVecSVecSubExpr class represents th...
Definition: DVecSVecSubExpr.h:88
void smpAddAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP addition assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:122
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DVecSVecSubExpr.h:222
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:253
Compile time type selection.The SelectType class template selects one of the two given types T1 and T...
Definition: SelectType.h:59
Header file for the IsTemporary type trait class.
const ResultType CompositeType
Data type for composite expression templates.
Definition: DVecSVecSubExpr.h:142
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DVecSVecSubExpr.h:135
VT2::TransposeType TT2
Transpose type of the right-hand side sparse vector expression.
Definition: DVecSVecSubExpr.h:101
void assign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the assignment of a matrix to a matrix.
Definition: Matrix.h:271
const DenseIterator< Type > operator+(const DenseIterator< Type > &it, ptrdiff_t inc)
Addition between a DenseIterator and an integral value.
Definition: DenseIterator.h:556
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
#define BLAZE_CONSTRAINT_MUST_BE_SPARSE_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a sparse, N-dimensional vector type...
Definition: SparseVector.h:79
Constraint on the data type.
#define BLAZE_CONSTRAINT_MUST_BE_REFERENCE_TYPE(T)
Constraint on the data type.In case the given data type T is not a reference type, a compilation error is created.
Definition: Reference.h:78
CompressedMatrix< Type, false > TransposeType
Transpose type for expression template evaluations.
Definition: CompressedMatrix.h:2405
Constraint on the data type.
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:103
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DVecSVecSubExpr.h:114
VT2::CompositeType CT2
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecSubExpr.h:99
void multAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the multiplication assignment of a matrix to a matrix.
Definition: Matrix.h:361
Header file for the VecVecSubExpr base class.
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
LeftOperand leftOperand() const
Returns the left-hand side dense vector operand.
Definition: DVecSVecSubExpr.h:200
VT1::TransposeType TT1
Transpose type of the left-hand side dense vector expression.
Definition: DVecSVecSubExpr.h:100
Header file for the serial shim.
void smpAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:92
DVecSVecSubExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the DVecSVecSubExpr class.
Definition: DVecSVecSubExpr.h:165
Header file for the IsSparseVector type trait.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2407
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecSubExpr.h:148
Header file for run time assertion macros.
void addAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the addition assignment of a matrix to a matrix.
Definition: Matrix.h:301
const DenseIterator< Type > operator-(const DenseIterator< Type > &it, ptrdiff_t inc)
Subtraction between a DenseIterator and an integral value.
Definition: DenseIterator.h:585
void subAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the subtraction assignment of a matrix to matrix.
Definition: Matrix.h:331
VT1::ReturnType RN1
Return type of the left-hand side dense vector expression.
Definition: DVecSVecSubExpr.h:96
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecSubExpr.h:145
Header file for the IsDenseVector type trait.
LeftOperand lhs_
Left-hand side dense vector of the subtraction expression.
Definition: DVecSVecSubExpr.h:242
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a dense, N-dimensional vector type...
Definition: DenseVector.h:79
size_t size() const
Returns the current size/dimension of the vector.
Definition: DVecSVecSubExpr.h:190
Header file for the IsComputation type trait class.
VT2::ResultType RT2
Result type of the right-hand side sparse vector expression.
Definition: DVecSVecSubExpr.h:95
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:108
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:157
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2403
Header file for basic type definitions.
Header file for the SubvectorExprTrait class template.
Base template for the SubTrait class.
Definition: SubTrait.h:141
Header file for the IsColumnVector type trait.
Base class for all vector/vector subtraction expression templates.The VecVecSubExpr class serves as a...
Definition: VecVecSubExpr.h:65
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: DVecSVecSubExpr.h:179
Header file for the SubExprTrait class template.
#define BLAZE_CONSTRAINT_MUST_BE_VECTOR_WITH_TRANSPOSE_FLAG(T, TF)
Constraint on the data type.In case the given data type T is not a dense or sparse vector type and in...
Definition: TransposeFlag.h:238
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.In case of an invalid run time expression, the program execution is terminated. The BLAZE_INTERNAL_ASSERT macro can be disabled by setting the BLAZE_USER_ASSERTION flag to zero or by defining NDEBUG during the compilation.
Definition: Assert.h:101
Compile time check whether the given type is an expression template.This type trait class tests wheth...
Definition: IsExpression.h:87
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DVecSVecSubExpr.h:235
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
VT1::CompositeType CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecSubExpr.h:98