35 #ifndef _BLAZE_MATH_EXPRESSIONS_SVECSVECADDEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SVECSVECADDEXPR_H_
88 template<
typename VT1
91 class SVecSVecAddExpr :
public SparseVector< SVecSVecAddExpr<VT1,VT2,TF>, TF >
92 ,
private VecVecAddExpr
127 template<
typename VT >
128 struct UseSMPAssign {
129 enum { value = VT::smpAssignable };
156 enum { smpAssignable = 0 };
228 template<
typename T >
230 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
240 template<
typename T >
242 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
265 template<
typename VT >
276 CT1 x(
serial( rhs.lhs_ ) );
277 CT2 y(
serial( rhs.rhs_ ) );
283 const LeftIterator lend( x.end() );
284 const RightIterator rend( y.end() );
286 for( LeftIterator l=x.begin(); l!=lend; ++l ) {
287 (~lhs)[l->index()] = l->value();
290 for( RightIterator r=y.begin(); r!=rend; ++r ) {
292 (~lhs)[r->index()] = r->value();
294 (~lhs)[r->index()] += r->value();
313 template<
typename VT >
324 CT1 x(
serial( rhs.lhs_ ) );
325 CT2 y(
serial( rhs.rhs_ ) );
331 const LeftIterator lend( x.end() );
332 const RightIterator rend( y.end() );
334 for( LeftIterator l=x.begin(); l!=lend; ++l ) {
335 (~lhs)[l->index()] = l->value();
338 for( RightIterator r=y.begin(); r!=rend; ++r ) {
339 (~lhs)[r->index()] += r->value();
357 template<
typename VT >
367 CT1 x(
serial( rhs.lhs_ ) );
368 CT2 y(
serial( rhs.rhs_ ) );
374 const LeftIterator lend( x.end() );
375 const RightIterator rend( y.end() );
377 LeftIterator l( x.begin() );
378 RightIterator r( y.begin() );
380 while( l != lend && r != rend )
382 if( l->index() < r->index() ) {
383 (~lhs).append( l->index(), l->value() );
386 else if( l->index() > r->index() ) {
387 (~lhs).append( r->index(), r->value() );
391 (~lhs).append( l->index(), l->value() + r->value() );
398 (~lhs).append( l->index(), l->value() );
403 (~lhs).append( r->index(), r->value() );
422 template<
typename VT >
451 template<
typename VT >
480 template<
typename VT >
491 const ResultType tmp(
serial( rhs ) );
523 template<
typename VT >
524 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
555 template<
typename VT >
556 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
588 template<
typename VT >
589 friend inline typename EnableIf< UseSMPAssign<VT> >::Type
600 const ResultType tmp( rhs );
655 template<
typename T1
658 inline const SVecSVecAddExpr<T1,T2,TF>
663 if( (~lhs).
size() != (~rhs).
size() )
664 throw std::invalid_argument(
"Vector sizes do not match" );
681 template<
typename VT1,
typename VT2,
bool TF >
682 struct Size< SVecSVecAddExpr<VT1,VT2,TF> >
683 :
public Max< Size<VT1>, Size<VT2> >::Type
699 template<
typename VT1,
typename VT2,
bool TF,
bool AF >
700 struct SubvectorExprTrait< SVecSVecAddExpr<VT1,VT2,TF>, AF >
704 typedef typename AddExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
705 ,
typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
size_t nonZeros() const
Returns the number of non-zero elements in the sparse vector.
Definition: SVecSVecAddExpr.h:197
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: SVecSVecAddExpr.h:142
Header file for the Max class template.
Header file for mathematical functions.
BLAZE_ALWAYS_INLINE 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:879
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
AddTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: SVecSVecAddExpr.h:137
VT1::ResultType RT1
Result type of the left-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:97
Header file for basic type definitions.
Header file for the SparseVector base class.
VT2::ReturnType RN2
Return type of the right-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:100
LeftOperand lhs_
Left-hand side sparse vector of the addition expression.
Definition: SVecSVecAddExpr.h:248
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:264
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:209
ResultType::ElementType ElementType
Resulting element type.
Definition: SVecSVecAddExpr.h:139
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: SVecSVecAddExpr.h:229
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: SVecSVecAddExpr.h:176
VT2::CompositeType CT2
Composite type of the right-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:102
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2507
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:261
Header file for the AddExprTrait class template.
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:148
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:699
Header file for the Computation base class.
VT1::TransposeType TT1
Transpose type of the left-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:103
bool isDefault(const DiagonalProxy< MT > &proxy)
Returns whether the represented element is in default state.
Definition: DiagonalProxy.h:861
Constraint on the data type.
LeftOperand leftOperand() const
Returns the left-hand side sparse vector operand.
Definition: SVecSVecAddExpr.h:207
Expression object for sparse vector-sparse vector additions.The SVecSVecAddExpr class represents the ...
Definition: Forward.h:117
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:263
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 DisableIf class template.
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.
VT1::CompositeType CT1
Composite type of the left-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:101
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2511
#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
const MT::ElementType min(const DenseMatrix< MT, SO > &dm)
Returns the smallest element of the dense matrix.
Definition: DenseMatrix.h:1602
BLAZE_ALWAYS_INLINE 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:635
VT2::ResultType RT2
Result type of the right-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:98
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.
size_t size() const
Returns the current size/dimension of the vector.
Definition: SVecSVecAddExpr.h:187
#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:2504
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:151
Header file for the EnableIf class template.
RightOperand rightOperand() const
Returns the right-hand side sparse vector operand.
Definition: SVecSVecAddExpr.h:217
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
VT2::TransposeType TT2
Transpose type of the right-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:104
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2506
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
SVecSVecAddExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the SVecSVecAddExpr class.
Definition: SVecSVecAddExpr.h:162
Header file for run time assertion macros.
Base template for the AddTrait class.
Definition: AddTrait.h:150
Header file for the addition trait.
BLAZE_ALWAYS_INLINE 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:742
VT1::ReturnType RN1
Return type of the left-hand side sparse vector expression.
Definition: SVecSVecAddExpr.h:99
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
Header file for the isDefault shim.
AddExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: SVecSVecAddExpr.h:117
Header file for the RemoveReference type trait.
Substitution Failure Is Not An Error (SFINAE) class.The DisableIf class template is an auxiliary tool...
Definition: DisableIf.h:184
RightOperand rhs_
Right-hand side sparse vector of the addition expression.
Definition: SVecSVecAddExpr.h:249
Header file for the IsComputation type trait class.
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
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:108
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: SVecSVecAddExpr.h:241
#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:2502
Header file for the SubvectorExprTrait class template.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: SVecSVecAddExpr.h:138
Header file for the IsResizable type trait.
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:238
const ResultType CompositeType
Data type for composite expression templates.
Definition: SVecSVecAddExpr.h:145
#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
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
SVecSVecAddExpr< VT1, VT2, TF > This
Type of this SVecSVecAddExpr instance.
Definition: SVecSVecAddExpr.h:136
BLAZE_ALWAYS_INLINE 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:849