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 };
272 return v_ *
it_->value();
402 return lhs_.size() *
rhs_.nonZeros();
414 return rhs_.nonZeros();
444 template<
typename T >
446 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
456 template<
typename T >
458 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
483 template<
typename MT >
502 const ConstIterator
begin( y.begin() );
503 const ConstIterator
end ( y.end() );
505 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
506 for( ConstIterator element=
begin; element!=
end; ++element ) {
507 (~lhs)(i,element->index()) = x[i] * element->
value();
527 template<
typename MT >
543 DVecTSVecMultExpr::selectAssignKernel( ~lhs, x, y );
562 template<
typename MT
565 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
566 selectAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
570 const ConstIterator
begin( y.begin() );
571 const ConstIterator
end ( y.end() );
573 for( ConstIterator element=
begin; element!=
end; ++element ) {
574 for(
size_t i=0UL; i<(~A).
rows(); ++i ) {
575 (~A)(i,element->index()) = x[i] * element->
value();
596 template<
typename MT
599 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
600 selectAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
604 typedef IntrinsicTrait<ElementType> IT;
605 typedef typename IT::Type IntrinsicType;
607 const ConstIterator
begin( y.begin() );
608 const ConstIterator
end ( y.end() );
610 for( ConstIterator element=
begin; element!=
end; ++element )
612 const IntrinsicType y1(
set( element->value() ) );
614 for(
size_t i=0UL; i<(~A).
rows(); i+=IT::size ) {
615 (~A).
store( i, element->index(), x.load(i) * y1 );
637 template<
typename MT >
638 friend inline typename EnableIf< UseAssign<MT> >::Type
656 const ConstIterator
begin( y.begin() );
657 const ConstIterator
end ( y.end() );
659 for(
size_t i=0UL; i<x.size(); ++i ) {
661 (~lhs).reserve( i, y.nonZeros() );
662 for( ConstIterator element=
begin; element!=
end; ++element ) {
663 (~lhs).append( i, element->index(), x[i] * element->value() );
684 template<
typename MT >
702 const ConstIterator
begin( y.begin() );
703 const ConstIterator
end ( y.end() );
705 for( ConstIterator element=
begin; element!=
end; ++element ) {
707 (~lhs).reserve( element->index(), x.size() );
708 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
709 (~lhs).append( i, element->index(), x[i] * element->value() );
733 template<
typename MT >
734 friend inline typename EnableIf< UseAssign<MT> >::Type
752 const ConstIterator
begin( y.begin() );
753 const ConstIterator
end ( y.end() );
755 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
756 for( ConstIterator element=
begin; element!=
end; ++element ) {
757 (~lhs)(i,element->index()) += x[i] * element->
value();
777 template<
typename MT >
793 DVecTSVecMultExpr::selectAddAssignKernel( ~lhs, x, y );
812 template<
typename MT
815 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
816 selectAddAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
820 const ConstIterator
begin( y.begin() );
821 const ConstIterator
end ( y.end() );
823 for( ConstIterator element=
begin; element!=
end; ++element ) {
824 for(
size_t i=0UL; i<(~A).
rows(); ++i ) {
825 (~A)(i,element->index()) += x[i] * element->
value();
846 template<
typename MT
849 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
850 selectAddAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
854 typedef IntrinsicTrait<ElementType> IT;
855 typedef typename IT::Type IntrinsicType;
857 const ConstIterator
begin( y.begin() );
858 const ConstIterator
end ( y.end() );
860 for( ConstIterator element=
begin; element!=
end; ++element )
862 const IntrinsicType y1(
set( element->value() ) );
864 for(
size_t i=0UL; i<(~A).
rows(); i+=IT::size ) {
865 (~A).
store( i, element->index(), (~A).
load(i,element->index()) + x.load(i) * y1 );
892 template<
typename MT >
893 friend inline typename EnableIf< UseAssign<MT> >::Type
911 const ConstIterator
begin( y.begin() );
912 const ConstIterator
end ( y.end() );
914 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
915 for( ConstIterator element=
begin; element!=
end; ++element ) {
916 (~lhs)(i,element->index()) -= x[i] * element->
value();
936 template<
typename MT >
952 DVecTSVecMultExpr::selectSubAssignKernel( ~lhs, x, y );
971 template<
typename MT
974 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
975 selectSubAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
979 const ConstIterator
begin( y.begin() );
980 const ConstIterator
end ( y.end() );
982 for( ConstIterator element=
begin; element!=
end; ++element ) {
983 for(
size_t i=0UL; i<(~A).
rows(); ++i ) {
984 (~A)(i,element->index()) -= x[i] * element->
value();
1005 template<
typename MT
1008 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
1009 selectSubAssignKernel( DenseMatrix<MT,true>& A,
const VT3& x,
const VT4& y )
1013 typedef IntrinsicTrait<ElementType> IT;
1014 typedef typename IT::Type IntrinsicType;
1016 const ConstIterator
begin( y.begin() );
1017 const ConstIterator
end ( y.end() );
1019 for( ConstIterator element=
begin; element!=
end; ++element )
1021 const IntrinsicType y1(
set( element->value() ) );
1023 for(
size_t i=0UL; i<(~A).
rows(); i+=IT::size ) {
1024 (~A).
store( i, element->index(), (~A).
load(i,element->index()) - x.load(i) * y1 );
1092 template<
typename T1
1094 inline const DVecTSVecMultExpr<T1,T2>
1114 template<
typename VT1,
typename VT2 >
1119 typedef typename MultExprTrait< typename SubvectorExprTrait<const VT1>::Type
1120 ,
typename SubvectorExprTrait<const VT2>::Type >::Type Type;
1129 template<
typename VT1,
typename VT2 >
1134 typedef typename MultExprTrait< typename VT1::ReturnType, VT2 >::Type Type;
1143 template<
typename VT1,
typename VT2 >
1148 typedef typename MultExprTrait< VT1, typename VT2::ReturnType >::Type Type;
ValueType * PointerType
Pointer return type.
Definition: DVecTSVecMultExpr.h:214
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:221
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:215
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:3703
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DVecTSVecMultExpr.h:457
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:216
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DVecTSVecMultExpr.h:381
#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:4555
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:196
Header file for the ColumnExprTrait class template.
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DVecTSVecMultExpr.h:303
const ConstIterator * operator->() const
Direct access to the sparse matrix element at the current iterator position.
Definition: DVecTSVecMultExpr.h:261
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:2375
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:248
IteratorType it_
Iterator over the elements of the right-hand side sparse vector expression.
Definition: DVecTSVecMultExpr.h:322
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:281
std::forward_iterator_tag IteratorCategory
The iterator category.
Definition: DVecTSVecMultExpr.h:212
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:346
Element ValueType
Type of the underlying pointers.
Definition: DVecTSVecMultExpr.h:213
const Element operator*() const
Direct access to the sparse matrix element at the current iterator position.
Definition: DVecTSVecMultExpr.h:251
DVecTSVecMultExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the DVecTSVecMultExpr class.
Definition: DVecTSVecMultExpr.h:333
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:220
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row i.
Definition: DVecTSVecMultExpr.h:360
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:240
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2379
LeftOperand lhs_
Left-hand side dense vector of the multiplication expression.
Definition: DVecTSVecMultExpr.h:464
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:445
ConstIterator(LeftElement v, IteratorType it)
Constructor for the ConstIterator class.
Definition: DVecTSVecMultExpr.h:229
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DVecTSVecMultExpr.h:292
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:219
#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:223
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2373
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:433
Header file for the IsNumeric type trait.
LeftOperand leftOperand() const
Returns the left-hand side dense vector operand.
Definition: DVecTSVecMultExpr.h:423
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row i.
Definition: DVecTSVecMultExpr.h:371
Iterator over the elements of the dense vector-sparse vector outer product expression.
Definition: DVecTSVecMultExpr.h:199
Header file for the SubmatrixExprTrait class template.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2374
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:648
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DVecTSVecMultExpr.h:391
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:465
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:204
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row.
Definition: DVecTSVecMultExpr.h:412
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:222
ReturnType value() const
Access to the current value of the sparse element.
Definition: DVecTSVecMultExpr.h:271
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:207
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:2370
LeftElement v_
Element of the left-hand side dense vector expression.
Definition: DVecTSVecMultExpr.h:321
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two expression iterators.
Definition: DVecTSVecMultExpr.h:314
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:210
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:401
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