35 #ifndef _BLAZE_MATH_EXPRESSIONS_DVECTSVECMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DVECTSVECMULTEXPR_H_
90 template<
typename VT1
134 template<
typename MT >
137 enum { value = useAssign };
147 template<
typename T1,
typename T2,
typename T3 >
148 struct UseVectorizedKernel {
149 enum { value = T1::vectorizable && T2::vectorizable &&
162 template<
typename T1,
typename T2,
typename T3 >
163 struct UseDefaultKernel {
164 enum { value = !UseVectorizedKernel<T1,T2,T3>::value };
198 enum { smpAssignable = 0 };
277 return v_ *
it_->value();
407 return lhs_.size() *
rhs_.nonZeros();
419 return rhs_.nonZeros();
449 template<
typename T >
451 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
461 template<
typename T >
463 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
488 template<
typename MT >
507 const ConstIterator
begin( y.begin() );
508 const ConstIterator
end ( y.end() );
510 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
511 for( ConstIterator element=
begin; element!=
end; ++element ) {
512 (~lhs)(i,element->index()) = x[i] * element->
value();
532 template<
typename MT >
548 DVecTSVecMultExpr::selectAssignKernel( ~lhs, x, y );
567 template<
typename MT
570 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
571 selectAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
575 const ConstIterator
begin( y.begin() );
576 const ConstIterator
end ( y.end() );
578 for( ConstIterator element=
begin; element!=
end; ++element ) {
579 for(
size_t i=0UL; i<(~A).
rows(); ++i ) {
580 (~A)(i,element->index()) = x[i] * element->
value();
601 template<
typename MT
604 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
605 selectAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
609 typedef IntrinsicTrait<ElementType> IT;
610 typedef typename IT::Type IntrinsicType;
612 const ConstIterator
begin( y.begin() );
613 const ConstIterator
end ( y.end() );
615 for( ConstIterator element=
begin; element!=
end; ++element )
617 const IntrinsicType y1(
set( element->value() ) );
619 for(
size_t i=0UL; i<(~A).
rows(); i+=IT::size ) {
620 (~A).
store( i, element->index(), x.load(i) * y1 );
642 template<
typename MT >
643 friend inline typename EnableIf< UseAssign<MT> >::Type
661 const ConstIterator
begin( y.begin() );
662 const ConstIterator
end ( y.end() );
664 for(
size_t i=0UL; i<x.size(); ++i ) {
666 (~lhs).reserve( i, y.nonZeros() );
667 for( ConstIterator element=
begin; element!=
end; ++element ) {
668 (~lhs).append( i, element->index(), x[i] * element->value() );
689 template<
typename MT >
707 const ConstIterator
begin( y.begin() );
708 const ConstIterator
end ( y.end() );
710 for( ConstIterator element=
begin; element!=
end; ++element ) {
712 (~lhs).reserve( element->index(), x.size() );
713 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
714 (~lhs).append( i, element->index(), x[i] * element->value() );
738 template<
typename MT >
739 friend inline typename EnableIf< UseAssign<MT> >::Type
757 const ConstIterator
begin( y.begin() );
758 const ConstIterator
end ( y.end() );
760 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
761 for( ConstIterator element=
begin; element!=
end; ++element ) {
762 (~lhs)(i,element->index()) += x[i] * element->
value();
782 template<
typename MT >
798 DVecTSVecMultExpr::selectAddAssignKernel( ~lhs, x, y );
817 template<
typename MT
820 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
821 selectAddAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
825 const ConstIterator
begin( y.begin() );
826 const ConstIterator
end ( y.end() );
828 for( ConstIterator element=
begin; element!=
end; ++element ) {
829 for(
size_t i=0UL; i<(~A).
rows(); ++i ) {
830 (~A)(i,element->index()) += x[i] * element->
value();
851 template<
typename MT
854 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
855 selectAddAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
859 typedef IntrinsicTrait<ElementType> IT;
860 typedef typename IT::Type IntrinsicType;
862 const ConstIterator
begin( y.begin() );
863 const ConstIterator
end ( y.end() );
865 for( ConstIterator element=
begin; element!=
end; ++element )
867 const IntrinsicType y1(
set( element->value() ) );
869 for(
size_t i=0UL; i<(~A).
rows(); i+=IT::size ) {
870 (~A).
store( i, element->index(), (~A).
load(i,element->index()) + x.load(i) * y1 );
897 template<
typename MT >
898 friend inline typename EnableIf< UseAssign<MT> >::Type
916 const ConstIterator
begin( y.begin() );
917 const ConstIterator
end ( y.end() );
919 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
920 for( ConstIterator element=
begin; element!=
end; ++element ) {
921 (~lhs)(i,element->index()) -= x[i] * element->
value();
941 template<
typename MT >
957 DVecTSVecMultExpr::selectSubAssignKernel( ~lhs, x, y );
976 template<
typename MT
979 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
980 selectSubAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
984 const ConstIterator
begin( y.begin() );
985 const ConstIterator
end ( y.end() );
987 for( ConstIterator element=
begin; element!=
end; ++element ) {
988 for(
size_t i=0UL; i<(~A).
rows(); ++i ) {
989 (~A)(i,element->index()) -= x[i] * element->
value();
1010 template<
typename MT
1013 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
1014 selectSubAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
1018 typedef IntrinsicTrait<ElementType> IT;
1019 typedef typename IT::Type IntrinsicType;
1021 const ConstIterator
begin( y.begin() );
1022 const ConstIterator
end ( y.end() );
1024 for( ConstIterator element=
begin; element!=
end; ++element )
1026 const IntrinsicType y1(
set( element->value() ) );
1028 for(
size_t i=0UL; i<(~A).
rows(); i+=IT::size ) {
1029 (~A).
store( i, element->index(), (~A).
load(i,element->index()) - x.load(i) * y1 );
1097 template<
typename T1
1099 inline const DVecTSVecMultExpr<T1,T2>
1119 template<
typename VT1,
typename VT2,
bool AF >
1124 typedef typename MultExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
1125 ,
typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
1134 template<
typename VT1,
typename VT2 >
1139 typedef typename MultExprTrait< typename VT1::ReturnType, VT2 >::Type Type;
1148 template<
typename VT1,
typename VT2 >
1153 typedef typename MultExprTrait< VT1, typename VT2::ReturnType >::Type Type;
ValueType * PointerType
Pointer return type.
Definition: DVecTSVecMultExpr.h:219
Pointer difference type of the Blaze library.
VT1::ReturnType RN1
Return type of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:100
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
Compile time check for numeric types.This type trait tests whether or not the given template paramete...
Definition: IsNumeric.h:98
PointerType pointer
Pointer return type.
Definition: DVecTSVecMultExpr.h:226
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
EnableIf< IsIntegral< T >, Load< T, sizeof(T)> >::Type::Type load(const T *address)
Loads a vector of integral values.
Definition: Load.h:222
VT2::ElementType ET2
Element type of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:105
ValueType & ReferenceType
Reference return type.
Definition: DVecTSVecMultExpr.h:220
Header file for the UNUSED_PARAMETER function template.
const DMatDMatMultExpr< T1, T2 > operator*(const DenseMatrix< T1, false > &lhs, const DenseMatrix< T2, false > &rhs)
Multiplication operator for the multiplication of two row-major dense matrices ( ).
Definition: DMatDMatMultExpr.h:4075
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DVecTSVecMultExpr.h:462
SelectType< IsComputation< VT1 >::value, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side dense vector operand.
Definition: DVecTSVecMultExpr.h:190
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DVecTSVecMultExpr.h:221
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DVecTSVecMultExpr.h:386
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a column dense or sparse vector type...
Definition: TransposeFlag.h:159
bool isDefault(const DynamicMatrix< Type, SO > &m)
Returns whether the given dense matrix is in default state.
Definition: DynamicMatrix.h:4622
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:197
Header file for the ColumnExprTrait class template.
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DVecTSVecMultExpr.h:308
const ConstIterator * operator->() const
Direct access to the sparse matrix element at the current iterator position.
Definition: DVecTSVecMultExpr.h:266
VT1::ElementType ET1
Element type of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:104
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2384
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:249
IteratorType it_
Iterator over the elements of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:327
void UNUSED_PARAMETER(const T1 &)
Suppression of unused parameter warnings.
Definition: Unused.h:84
size_t index() const
Access to the current index of the sparse element.
Definition: DVecTSVecMultExpr.h:286
std::forward_iterator_tag IteratorCategory
The iterator category.
Definition: DVecTSVecMultExpr.h:217
Header file for the Computation base class.
Type relationship analysis.This class tests if the two data types A and B are equal. For this type comparison, the cv-qualifiers of both data types are ignored. If A and B are the same data type (ignoring the cv-qualifiers), then the value member enumeration is set to 1, the nested type definition Type is TrueType, and the class derives from TrueType. Otherwise value is set to 0, Type is FalseType, and the class derives from FalseType.
Definition: IsSame.h:158
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:187
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DVecTSVecMultExpr.h:351
Element ValueType
Type of the underlying pointers.
Definition: DVecTSVecMultExpr.h:218
const Element operator*() const
Direct access to the sparse matrix element at the current iterator position.
Definition: DVecTSVecMultExpr.h:256
DVecTSVecMultExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the DVecTSVecMultExpr class.
Definition: DVecTSVecMultExpr.h:338
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:104
Constraint on the data type.
MultExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DVecTSVecMultExpr.h:118
Header file for the SparseMatrix base class.
Constraint on the data type.
Header file for the MultExprTrait class template.
ValueType value_type
Type of the underlying pointers.
Definition: DVecTSVecMultExpr.h:225
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row i.
Definition: DVecTSVecMultExpr.h:365
Header file for the ValueIndexPair class.
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.
Header file for the multiplication trait.
ConstIterator & operator++()
Pre-increment operator.
Definition: DVecTSVecMultExpr.h:245
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2388
LeftOperand lhs_
Left-hand side dense vector of the multiplication expression.
Definition: DVecTSVecMultExpr.h:469
ResultType::ElementType ElementType
Resulting element type.
Definition: DVecTSVecMultExpr.h:175
Expression object for dense vector-sparse vector outer products.The DVecTSVecMultExpr class represent...
Definition: DVecTSVecMultExpr.h:92
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DVecTSVecMultExpr.h:450
ConstIterator(LeftElement v, IteratorType it)
Constructor for the ConstIterator class.
Definition: DVecTSVecMultExpr.h:234
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DVecTSVecMultExpr.h:297
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:179
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
IteratorCategory iterator_category
The iterator category.
Definition: DVecTSVecMultExpr.h:224
#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.
Constraint on the data type.
DifferenceType difference_type
Difference between two iterators.
Definition: DVecTSVecMultExpr.h:228
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2382
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Header file for the VecTVecMultExpr base class.
VT1::CompositeType CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:102
Header file for the EnableIf class template.
RightOperand rightOperand() const
Returns the right-hand side sparse vector operand.
Definition: DVecTSVecMultExpr.h:438
Header file for the IsNumeric type trait.
LeftOperand leftOperand() const
Returns the left-hand side dense vector operand.
Definition: DVecTSVecMultExpr.h:428
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row i.
Definition: DVecTSVecMultExpr.h:376
Iterator over the elements of the dense vector-sparse vector outer product expression.
Definition: DVecTSVecMultExpr.h:204
Header file for the SubmatrixExprTrait class template.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2383
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
VT1::ResultType RT1
Result type of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:98
Intrinsic characteristics of data types.The IntrinsicTrait class template provides the intrinsic char...
Definition: IntrinsicTrait.h:748
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DVecTSVecMultExpr.h:396
Header file for run time assertion macros.
Base template for the MultTrait class.
Definition: MultTrait.h:141
RightOperand rhs_
Right-hand side sparse vector of the multiplication expression.
Definition: DVecTSVecMultExpr.h:470
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:209
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DVecTSVecMultExpr.h:174
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
DVecTSVecMultExpr< VT1, VT2 > This
Type of this DVecTSVecMultExpr instance.
Definition: DVecTSVecMultExpr.h:171
ValueIndexPair< ElementType > Element
Element type of the sparse matrix expression.
Definition: DVecTSVecMultExpr.h:209
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row.
Definition: DVecTSVecMultExpr.h:417
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:239
Header file for the isDefault shim.
ReferenceType reference
Reference return type.
Definition: DVecTSVecMultExpr.h:227
ReturnType value() const
Access to the current value of the sparse element.
Definition: DVecTSVecMultExpr.h:276
Header file for the RemoveReference type trait.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DVecTSVecMultExpr.h:173
Base class for all outer product expression templates.The VecTVecMultExpr class serves as a tag for a...
Definition: VecTVecMultExpr.h:66
ET1 LeftElement
Element type of the dense vector expression.
Definition: DVecTSVecMultExpr.h:212
Header file for all intrinsic functionality.
#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
VT2::ReturnType RN2
Return type of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:101
Index-value-pair for sparse vectors and matrices.The ValueIndexPair class represents a single index-v...
Definition: ValueIndexPair.h:69
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DVecTSVecMultExpr.h:172
Header file for the IsComputation type trait class.
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DVecTSVecMultExpr.h:178
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:105
#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:2379
LeftElement v_
Element of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:326
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two expression iterators.
Definition: DVecTSVecMultExpr.h:319
SelectType< useAssign, const ResultType, const DVecTSVecMultExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: DVecTSVecMultExpr.h:181
Header file for basic type definitions.
#define BLAZE_CONSTRAINT_MUST_BE_ROW_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a row dense or sparse vector type (i...
Definition: TransposeFlag.h:81
RemoveReference< RightOperand >::Type::ConstIterator IteratorType
Iterator type of the sparse vector expression.
Definition: DVecTSVecMultExpr.h:215
Header file for the SubvectorExprTrait class template.
VT2::CompositeType CT2
Composite type of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:103
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: DVecTSVecMultExpr.h:406
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:184
VT2::ResultType RT2
Result type of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:99
Evaluation of the resulting expression type of a multiplication.Via this type trait it is possible to...
Definition: MultExprTrait.h:137
#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
EnableIf< IsIntegral< T >, Set< T, sizeof(T)> >::Type::Type set(T value)
Sets all values in the vector to the given integral value.
Definition: Set.h:209
void store(float *address, const sse_float_t &value)
Aligned store of a vector of 'float' values.
Definition: Store.h:242
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
SelectType< IsComputation< VT2 >::value, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side dense vector operand.
Definition: DVecTSVecMultExpr.h:193