35 #ifndef _BLAZE_MATH_EXPRESSIONS_SVECTDVECMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SVECTDVECMULTEXPR_H_
97 template<
typename VT1
99 class SVecTDVecMultExpr :
public SparseMatrix< SVecTDVecMultExpr<VT1,VT2>, true >
100 ,
private VecTVecMultExpr
101 ,
private Computation
140 template<
typename MT >
143 enum { value = useAssign };
153 template<
typename T1,
typename T2,
typename T3 >
154 struct UseVectorizedKernel {
155 enum { value = T1::vectorizable && T3::vectorizable &&
168 template<
typename T1,
typename T2,
typename T3 >
169 struct UseDefaultKernel {
170 enum { value = !UseVectorizedKernel<T1,T2,T3>::value };
204 enum { smpAssignable = 0 };
283 return it_->value() *
v_;
413 return lhs_.nonZeros() *
rhs_.size();
425 return lhs_.nonZeros();
494 template<
typename T >
496 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
506 template<
typename T >
508 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
530 template<
typename MT >
540 LT x(
serial( rhs.lhs_ ) );
541 RT y(
serial( rhs.rhs_ ) );
548 SVecTDVecMultExpr::selectAssignKernel( ~lhs, x, y );
567 template<
typename MT
575 const ConstIterator
end( x.end() );
577 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
579 for(
size_t j=0UL; j<y.size(); ++j ) {
580 (~A)(element->index(),j) = element->
value() * y[j];
602 template<
typename MT
605 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
606 selectAssignKernel( DenseMatrix<MT,false>& A,
const VT3& x,
const VT4& y )
610 typedef IntrinsicTrait<ElementType> IT;
611 typedef typename IT::Type IntrinsicType;
613 const ConstIterator
begin( x.begin() );
614 const ConstIterator
end ( x.end() );
616 for( ConstIterator element=
begin; element!=
end; ++element )
618 const IntrinsicType x1(
set( element->value() ) );
621 (~A).
store( element->index(), j, x1 * y.load(j) );
644 template<
typename MT >
645 friend inline typename EnableIf< UseAssign<MT> >::Type
655 LT x(
serial( rhs.lhs_ ) );
656 RT y(
serial( rhs.rhs_ ) );
663 const ConstIterator
end( x.end() );
665 for(
size_t i=0UL; i<y.size(); ++i ) {
667 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
668 (~lhs)(element->index(),i) = element->
value() * y[i];
688 template<
typename MT >
700 LT x(
serial( rhs.lhs_ ) );
701 RT y(
serial( rhs.rhs_ ) );
708 const ConstIterator
end( x.end() );
710 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
712 (~lhs).reserve( element->index(), y.size() );
713 for(
size_t i=0UL; i<y.size(); ++i ) {
714 (~lhs).append( element->index(), i, element->value() * y[i] );
738 template<
typename MT >
739 friend inline typename EnableIf< UseAssign<MT> >::Type
749 LT x(
serial( rhs.lhs_ ) );
750 RT y(
serial( rhs.rhs_ ) );
757 const ConstIterator
end( x.end() );
759 for(
size_t i=0UL; i<y.size(); ++i ) {
761 (~lhs).reserve( i, x.nonZeros() );
762 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
763 (~lhs).append( element->index(), i, element->value() * y[i] );
784 template<
typename MT >
794 LT x(
serial( rhs.lhs_ ) );
795 RT y(
serial( rhs.rhs_ ) );
802 SVecTDVecMultExpr::selectAddAssignKernel( ~lhs, x, y );
821 template<
typename MT
824 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
825 selectAddAssignKernel( DenseMatrix<MT,false>& A,
const VT3& x,
const VT4& y )
829 const ConstIterator
end( x.end() );
831 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
833 for(
size_t i=0UL; i<y.size(); ++i ) {
834 (~A)(element->index(),i) += element->
value() * y[i];
856 template<
typename MT
859 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
860 selectAddAssignKernel( DenseMatrix<MT,false>& A,
const VT3& x,
const VT4& y )
864 typedef IntrinsicTrait<ElementType> IT;
865 typedef typename IT::Type IntrinsicType;
867 const ConstIterator
begin( x.begin() );
868 const ConstIterator
end ( x.end() );
870 for( ConstIterator element=
begin; element!=
end; ++element )
872 const IntrinsicType x1(
set( element->value() ) );
875 (~A).
store( element->index(), j, (~A).
load(element->index(),j) + x1 * y.load(j) );
898 template<
typename MT >
899 friend inline typename EnableIf< UseAssign<MT> >::Type
909 LT x(
serial( rhs.lhs_ ) );
910 RT y(
serial( rhs.rhs_ ) );
917 const ConstIterator
end( x.end() );
919 for(
size_t i=0UL; i<y.size(); ++i ) {
921 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
922 (~lhs)(element->index(),i) += element->
value() * y[i];
947 template<
typename MT >
957 LT x(
serial( rhs.lhs_ ) );
958 RT y(
serial( rhs.rhs_ ) );
965 SVecTDVecMultExpr::selectSubAssignKernel( ~lhs, x, y );
984 template<
typename MT
987 static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
988 selectSubAssignKernel( DenseMatrix<MT,false>& A,
const VT3& x,
const VT4& y )
992 const ConstIterator
end( x.end() );
994 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
996 for(
size_t i=0UL; i<y.size(); ++i ) {
997 (~A)(element->index(),i) -= element->
value() * y[i];
1019 template<
typename MT
1022 static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
1023 selectSubAssignKernel( DenseMatrix<MT,false>& A,
const VT3& x,
const VT4& y )
1027 typedef IntrinsicTrait<ElementType> IT;
1028 typedef typename IT::Type IntrinsicType;
1030 const ConstIterator
begin( x.begin() );
1031 const ConstIterator
end ( x.end() );
1033 for( ConstIterator element=
begin; element!=
end; ++element )
1035 const IntrinsicType x1(
set( element->value() ) );
1038 (~A).
store( element->index(), j, (~A).
load(element->index(),j) - x1 * y.load(j) );
1061 template<
typename MT >
1062 friend inline typename EnableIf< UseAssign<MT> >::Type
1072 LT x(
serial( rhs.lhs_ ) );
1073 RT y(
serial( rhs.rhs_ ) );
1080 const ConstIterator
end( x.end() );
1082 for(
size_t i=0UL; i<y.size(); ++i ) {
1084 for( ConstIterator element=x.begin(); element!=
end; ++element ) {
1085 (~lhs)(element->index(),i) -= element->
value() * y[i];
1155 template<
typename T1
1157 inline const SVecTDVecMultExpr<T1,T2>
1177 template<
typename VT1,
typename VT2 >
1194 template<
typename VT1,
typename VT2 >
1211 template<
typename VT1,
typename VT2,
bool AF >
1216 typedef typename MultExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
1217 ,
typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
1226 template<
typename VT1,
typename VT2 >
1231 typedef typename MultExprTrait< typename VT1::ReturnType, VT2 >::Type Type;
1240 template<
typename VT1,
typename VT2 >
1245 typedef typename MultExprTrait< VT1, typename VT2::ReturnType >::Type Type;
std::forward_iterator_tag IteratorCategory
The iterator category.
Definition: SVecTDVecMultExpr.h:223
Pointer difference type of the Blaze library.
ValueType value_type
Type of the underlying pointers.
Definition: SVecTDVecMultExpr.h:231
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
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 Rows type trait.
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:4838
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SVecTDVecMultExpr.h:179
Iterator over the elements of the sparse vector-dense vector outer product expression.
Definition: SVecTDVecMultExpr.h:210
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: SVecTDVecMultExpr.h:495
RightOperand rightOperand() const
Returns the right-hand side dense vector operand.
Definition: SVecTDVecMultExpr.h:483
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:258
#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
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:205
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:193
Header file for the ColumnExprTrait class template.
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of column i.
Definition: SVecTDVecMultExpr.h:371
ConstIterator lowerBound(size_t i, size_t j) const
Returns an iterator to the first index not less then the given index.
Definition: SVecTDVecMultExpr.h:449
Expression object for sparse vector-dense vector outer products.The SVecTDVecMultExpr class represent...
Definition: Forward.h:121
VT1::ResultType RT1
Result type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:105
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
size_t index() const
Access to the current index of the sparse element.
Definition: SVecTDVecMultExpr.h:292
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
BLAZE_ALWAYS_INLINE EnableIf< IsIntegral< T >, Load< T, sizeof(T)> >::Type::Type load(const T *address)
Loads a vector of integral values.
Definition: Load.h:224
void UNUSED_PARAMETER(const T1 &)
Suppression of unused parameter warnings.
Definition: Unused.h:81
SelectType< IsComputation< VT1 >::value, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side dense vector operand.
Definition: SVecTDVecMultExpr.h:196
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of column i.
Definition: SVecTDVecMultExpr.h:382
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:695
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
ConstIterator(IteratorType it, RightElement v)
Constructor for the ConstIterator class.
Definition: SVecTDVecMultExpr.h:240
const Element operator*() const
Direct access to the sparse matrix element at the current iterator position.
Definition: SVecTDVecMultExpr.h:262
RightOperand rhs_
Right-hand side dense vector of the multiplication expression.
Definition: SVecTDVecMultExpr.h:515
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Constraint on the data type.
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: SVecTDVecMultExpr.h:227
DifferenceType difference_type
Difference between two iterators.
Definition: SVecTDVecMultExpr.h:234
Header file for the SparseMatrix base class.
Constraint on the data type.
Header file for the MultExprTrait class template.
IteratorType it_
Iterator over the elements of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:332
SVecTDVecMultExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the SVecTDVecMultExpr class.
Definition: SVecTDVecMultExpr.h:344
Header file for the ValueIndexPair class.
SelectType< IsComputation< VT2 >::value, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side dense vector operand.
Definition: SVecTDVecMultExpr.h:199
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.
RemoveReference< LeftOperand >::Type::ConstIterator IteratorType
Iterator type of the sparse vector expression.
Definition: SVecTDVecMultExpr.h:218
size_t columns() const
Returns the current number of columns of the matrix.
Definition: SVecTDVecMultExpr.h:402
VT2::CompositeType CT2
Composite type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:110
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: SVecTDVecMultExpr.h:314
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_VECTVECMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/vector ...
Definition: VecTVecMultExpr.h:161
MultExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: SVecTDVecMultExpr.h:125
Element ValueType
Type of the underlying pointers.
Definition: SVecTDVecMultExpr.h:224
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
Header file for the Columns type trait.
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
ET2 RightElement
Element type of the dense vector expression.
Definition: SVecTDVecMultExpr.h:221
Constraint on the data type.
VT2::ResultType RT2
Result type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:106
size_t rows() const
Returns the current number of rows of the matrix.
Definition: SVecTDVecMultExpr.h:392
RightElement v_
Element of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:333
Constraint on the data type.
SelectType< useAssign, const ResultType, const SVecTDVecMultExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: SVecTDVecMultExpr.h:187
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
Constraint on the data type.
SVecTDVecMultExpr< VT1, VT2 > This
Type of this SVecTDVecMultExpr instance.
Definition: SVecTDVecMultExpr.h:177
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 sparse vector expression.
Definition: SVecTDVecMultExpr.h:109
Header file for the EnableIf class template.
Header file for the serial shim.
ReferenceType reference
Reference return type.
Definition: SVecTDVecMultExpr.h:233
ConstIterator find(size_t i, size_t j) const
Searches for a specific matrix element.
Definition: SVecTDVecMultExpr.h:436
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: SVecTDVecMultExpr.h:303
Header file for the IsNumeric type trait.
VT2::ReturnType RN2
Return type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:108
BLAZE_ALWAYS_INLINE 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:211
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: SVecTDVecMultExpr.h:184
#define BLAZE_CONSTRAINT_MUST_NOT_BE_SYMMETRIC_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is a symmetric matrix type, a compilation error is created.
Definition: Symmetric.h:116
Header file for the SubmatrixExprTrait class template.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2477
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
IteratorCategory iterator_category
The iterator category.
Definition: SVecTDVecMultExpr.h:230
Intrinsic characteristics of data types.The IntrinsicTrait class template provides the intrinsic char...
Definition: IntrinsicTrait.h:749
Header file for run time assertion macros.
Base template for the MultTrait class.
Definition: MultTrait.h:142
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
ReturnType value() const
Access to the current value of the sparse element.
Definition: SVecTDVecMultExpr.h:282
LeftOperand leftOperand() const
Returns the left-hand side sparse vector operand.
Definition: SVecTDVecMultExpr.h:473
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
VT2::ElementType ET2
Element type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:112
ResultType::ElementType ElementType
Resulting element type.
Definition: SVecTDVecMultExpr.h:181
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: SVecTDVecMultExpr.h:507
Header file for the isDefault shim.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: SVecTDVecMultExpr.h:180
BLAZE_ALWAYS_INLINE bool isDefault(const NonNumericProxy< MT > &proxy)
Returns whether the represented element is in default state.
Definition: NonNumericProxy.h:874
Constraint on the data type.
ValueType * PointerType
Pointer return type.
Definition: SVecTDVecMultExpr.h:225
ConstIterator upperBound(size_t i, size_t j) const
Returns an iterator to the first index greater then the given index.
Definition: SVecTDVecMultExpr.h:462
VT1::ElementType ET1
Element type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:111
ValueIndexPair< ElementType > Element
Element type of the sparse matrix expression.
Definition: SVecTDVecMultExpr.h:215
#define BLAZE_CONSTRAINT_MUST_NOT_REQUIRE_EVALUATION(T)
Constraint on the data type.In case the given data type T requires an intermediate evaluation within ...
Definition: RequiresEvaluation.h:118
Constraint on the data type.
Header file for the RemoveReference type trait.
LeftOperand lhs_
Left-hand side sparse vector of the multiplication expression.
Definition: SVecTDVecMultExpr.h:514
Header file for all intrinsic functionality.
PointerType pointer
Pointer return type.
Definition: SVecTDVecMultExpr.h:232
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two expression iterators.
Definition: SVecTDVecMultExpr.h:325
#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
Index-value-pair for sparse vectors and matrices.The ValueIndexPair class represents a single index-v...
Definition: ValueIndexPair.h:70
Header file for the IsComputation type trait class.
const ConstIterator * operator->() const
Direct access to the sparse matrix element at the current iterator position.
Definition: SVecTDVecMultExpr.h:272
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:2473
Header file for basic type definitions.
ConstIterator & operator++()
Pre-increment operator.
Definition: SVecTDVecMultExpr.h:251
#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
Header file for the SubvectorExprTrait class template.
ValueType & ReferenceType
Reference return type.
Definition: SVecTDVecMultExpr.h:226
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: SVecTDVecMultExpr.h:178
VT1::ReturnType RN1
Return type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:107
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified column.
Definition: SVecTDVecMultExpr.h:423
BLAZE_ALWAYS_INLINE EnableIf< IsIntegral< T > >::Type store(T *address, const typename Store< T, sizeof(T)>::Type &value)
Aligned store of a vector of integral values.
Definition: Store.h:225
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SVecTDVecMultExpr.h:357
Header file for the Size type trait.
Evaluation of the resulting expression type of a multiplication.Via this type trait it is possible to...
Definition: MultExprTrait.h:137
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:190
#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.
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SVecTDVecMultExpr.h:412
Header file for the FunctionTrace class.
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