35 #ifndef _BLAZE_MATH_EXPRESSIONS_DVECSVECADDEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DVECSVECADDEXPR_H_
90 template<
typename VT1
129 template<
typename VT >
130 struct UseSMPAssign {
131 enum { value = ( !VT1::smpAssignable || !VT2::smpAssignable ) };
158 enum { vectorizable = 0 };
161 enum { smpAssignable = 0 };
197 inline ReturnType
at(
size_t index )
const {
198 if( index >=
lhs_.size() ) {
201 return (*
this)[index];
241 template<
typename T >
244 (
rhs_.canAlias( alias ) );
254 template<
typename T >
256 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
278 template<
typename VT >
286 addAssign( ~lhs, rhs.
rhs_ );
289 assign ( ~lhs, rhs.
lhs_ );
290 addAssign( ~lhs, rhs.
rhs_ );
308 template<
typename VT >
319 const ResultType tmp(
serial( rhs ) );
337 template<
typename VT >
338 friend inline void addAssign( DenseVector<VT,TF>& lhs,
const DVecSVecAddExpr& rhs )
344 addAssign( ~lhs, rhs.lhs_ );
345 addAssign( ~lhs, rhs.rhs_ );
366 template<
typename VT >
367 friend inline void subAssign( DenseVector<VT,TF>& lhs,
const DVecSVecAddExpr& rhs )
373 subAssign( ~lhs, rhs.lhs_ );
374 subAssign( ~lhs, rhs.rhs_ );
395 template<
typename VT >
396 friend inline void multAssign( DenseVector<VT,TF>& lhs,
const DVecSVecAddExpr& rhs )
406 const ResultType tmp(
serial( rhs ) );
407 multAssign( ~lhs, tmp );
430 template<
typename VT >
431 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
438 if( !IsComputation<VT1>::value &&
isSame( ~lhs, rhs.lhs_ ) ) {
463 template<
typename VT >
464 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
475 const ResultType tmp( rhs );
495 template<
typename VT >
496 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
527 template<
typename VT >
528 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
560 template<
typename VT >
561 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
572 const ResultType tmp( rhs );
629 template<
typename T1
632 inline const DVecSVecAddExpr<T1,T2,TF>
637 if( (~lhs).
size() != (~rhs).
size() ) {
672 template<
typename T1
675 inline const DVecSVecAddExpr<T2,T1,TF>
680 if( (~lhs).
size() != (~rhs).
size() ) {
710 template<
typename T1
714 inline const typename AddExprTrait< DVecSVecAddExpr<T1,T2,TF>, T3 >::Type
715 operator+(
const DVecSVecAddExpr<T1,T2,TF>& lhs,
const DenseVector<T3,TF>& rhs )
719 return ( lhs.leftOperand() + (~rhs) ) + lhs.rightOperand();
738 template<
typename T1
742 inline const typename SubExprTrait< DVecSVecAddExpr<T1,T2,TF>, T3 >::Type
743 operator-(
const DVecSVecAddExpr<T1,T2,TF>& lhs,
const DenseVector<T3,TF>& rhs )
747 return ( lhs.leftOperand() - (~rhs) ) + lhs.rightOperand();
763 template<
typename VT1,
typename VT2,
bool TF >
765 :
public Max< Size<VT1>, Size<VT2> >
781 template<
typename VT1,
typename VT2,
typename VT3 >
787 typedef typename SelectType< IsDenseVector<VT1>::value && IsColumnVector<VT1>::value &&
788 IsSparseVector<VT2>::value && IsColumnVector<VT2>::value &&
789 IsDenseVector<VT3>::value && IsColumnVector<VT3>::value
790 ,
typename DVecSVecAddExprTrait< typename DVecDVecAddExprTrait<VT1,VT3>::Type, VT2 >::Type
791 , INVALID_TYPE >::Type Type;
801 template<
typename VT1,
typename VT2,
typename VT3 >
807 typedef typename SelectType< IsDenseVector<VT1>::value && IsRowVector<VT1>::value &&
808 IsSparseVector<VT2>::value && IsRowVector<VT2>::value &&
809 IsDenseVector<VT3>::value && IsRowVector<VT3>::value
810 ,
typename TDVecTSVecAddExprTrait< typename TDVecTDVecAddExprTrait<VT1,VT3>::Type, VT2 >::Type
811 , INVALID_TYPE >::Type Type;
821 template<
typename VT1,
typename VT2,
typename VT3 >
827 typedef typename SelectType< IsDenseVector<VT1>::value && IsColumnVector<VT1>::value &&
828 IsSparseVector<VT2>::value && IsColumnVector<VT2>::value &&
829 IsDenseVector<VT3>::value && IsColumnVector<VT3>::value
830 ,
typename DVecSVecAddExprTrait< typename DVecDVecSubExprTrait<VT1,VT3>::Type, VT2 >::Type
831 , INVALID_TYPE >::Type Type;
841 template<
typename VT1,
typename VT2,
typename VT3 >
847 typedef typename SelectType< IsDenseVector<VT1>::value && IsRowVector<VT1>::value &&
848 IsSparseVector<VT2>::value && IsRowVector<VT2>::value &&
849 IsDenseVector<VT3>::value && IsRowVector<VT3>::value
850 ,
typename TDVecTSVecAddExprTrait< typename TDVecTDVecSubExprTrait<VT1,VT3>::Type, VT2 >::Type
851 , INVALID_TYPE >::Type Type;
861 template<
typename VT1,
typename VT2,
bool TF,
bool AF >
866 typedef typename AddExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
867 ,
typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exceptionThis macro encapsulates the default way of...
Definition: Exception.h:187
VT2::ResultType RT2
Result type of the right-hand side sparse vector expression.
Definition: DVecSVecAddExpr.h:100
Header file for the Max class template.
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DVecSVecAddExpr.h:255
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
RightOperand rhs_
Right-hand side sparse vector of the addition expression.
Definition: DVecSVecAddExpr.h:263
Evaluation of the return type of an addition expression.Via this type trait it is possible to evaluat...
Definition: AddExprTrait.h:104
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 basic type definitions.
VT2::ReturnType RN2
Return type of the right-hand side sparse vector expression.
Definition: DVecSVecAddExpr.h:102
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DVecSVecAddExpr.h:140
VT1::ResultType RT1
Result type of the left-hand side dense vector expression.
Definition: DVecSVecAddExpr.h:99
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:207
VT1::CompositeType CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecAddExpr.h:103
BLAZE_ALWAYS_INLINE 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:647
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2588
Header file for the IsRowVector type trait.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:259
Header file for the DenseVector base class.
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: DVecSVecAddExpr.h:197
Header file for the AddExprTrait class template.
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:721
Header file for the Computation base class.
Expression object for dense vector-sparse vector additions.The DVecSVecAddExpr class represents the c...
Definition: DVecSVecAddExpr.h:93
RightOperand rightOperand() const
Returns the right-hand side sparse vector operand.
Definition: DVecSVecAddExpr.h:230
Constraint on the data type.
VT1::TransposeType TT1
Transpose type of the left-hand side dense vector expression.
Definition: DVecSVecAddExpr.h:105
LeftOperand leftOperand() const
Returns the left-hand side dense vector operand.
Definition: DVecSVecAddExpr.h:220
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:261
VT2::TransposeType TT2
Transpose type of the right-hand side sparse vector expression.
Definition: DVecSVecAddExpr.h:106
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.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
Header file for the VecVecAddExpr base class.
DVecSVecAddExpr< VT1, VT2, TF > This
Type of this DVecSVecAddExpr instance.
Definition: DVecSVecAddExpr.h:138
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exceptionThis macro encapsulates the default way of Bla...
Definition: Exception.h:331
const DenseIterator< Type, AF > operator+(const DenseIterator< Type, AF > &it, ptrdiff_t inc)
Addition between a DenseIterator and an integral value.
Definition: DenseIterator.h:610
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_VECVECADDEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/vector ...
Definition: VecVecAddExpr.h:165
ResultType::ElementType ElementType
Resulting element type.
Definition: DVecSVecAddExpr.h:141
const ResultType CompositeType
Data type for composite expression templates.
Definition: DVecSVecAddExpr.h:147
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
const DenseIterator< Type, AF > operator-(const DenseIterator< Type, AF > &it, ptrdiff_t inc)
Subtraction between a DenseIterator and an integral value.
Definition: DenseIterator.h:642
Constraint on the data type.
LeftOperand lhs_
Left-hand side dense vector of the addition expression.
Definition: DVecSVecAddExpr.h:262
#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:2585
Constraint on the data type.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DVecSVecAddExpr.h:242
VT2::CompositeType CT2
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecAddExpr.h:104
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Base class for all vector/vector addition expression templates.The VecVecAddExpr class serves as a ta...
Definition: VecVecAddExpr.h:65
Header file for the EnableIf class template.
Header file for the serial shim.
Constraint on the data type.
EnableIf< IsDenseMatrix< MT1 > >::Type smpSubAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:160
Header file for the IsSparseVector type trait.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2587
AddExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DVecSVecAddExpr.h:119
Header file for run time assertion macros.
Base template for the AddTrait class.
Definition: AddTrait.h:138
EnableIf< IsDenseMatrix< MT1 > >::Type smpAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:98
Utility type for generic codes.
Header file for the addition trait.
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecAddExpr.h:150
DVecSVecAddExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the DVecSVecAddExpr class.
Definition: DVecSVecAddExpr.h:170
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DVecSVecAddExpr.h:144
size_t size() const
Returns the current size/dimension of the vector.
Definition: DVecSVecAddExpr.h:210
Header file for the IsDenseVector type trait.
#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
AddTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DVecSVecAddExpr.h:139
Header file for the IsComputation type trait class.
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
EnableIf< IsDenseMatrix< MT1 > >::Type smpAddAssign(Matrix< 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:129
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: DVecSVecAddExpr.h:184
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:118
#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:2583
Header file for the SubvectorExprTrait class template.
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecAddExpr.h:153
VT1::ReturnType RN1
Return type of the left-hand side dense vector expression.
Definition: DVecSVecAddExpr.h:101
Header file for exception macros.
Header file for the IsColumnVector type trait.
Header file for the SubExprTrait class template.
EnableIf< IsDenseVector< VT1 > >::Type smpMultAssign(Vector< 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:189
Header file for the Size type trait.
#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:81
#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
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.