35 #ifndef _BLAZE_MATH_EXPRESSIONS_DVECSVECMULTEXPR_H_ 36 #define _BLAZE_MATH_EXPRESSIONS_DVECSVECMULTEXPR_H_ 92 template<
typename VT1
96 :
public VecVecMultExpr< SparseVector< DVecSVecMultExpr<VT1,VT2,TF>, TF > >
118 static constexpr
bool returnExpr = !IsTemporary_v<RN1> && !IsTemporary_v<RN2>;
132 static constexpr
bool useAssign = ( RequiresEvaluation_v<VT1> || RequiresEvaluation_v<VT2> );
135 template<
typename VT >
137 static constexpr
bool UseAssign_v =
useAssign;
332 if( index >=
lhs_.size() ) {
335 return (*
this)[index];
364 inline size_t size() const noexcept {
375 return rhs_.nonZeros();
441 template<
typename T >
442 inline bool canAlias(
const T* alias )
const noexcept {
443 return (
lhs_.canAlias( alias ) ||
rhs_.canAlias( alias ) );
453 template<
typename T >
454 inline bool isAliased(
const T* alias )
const noexcept {
455 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
479 template<
typename VT >
494 for(
auto element=y.begin(); element!=y.end(); ++element )
495 (~lhs)[element->index()] = x[element->index()] * element->value();
514 template<
typename VT >
530 (~lhs).reserve( y.nonZeros() );
533 for(
auto element=y.begin(); element!=y.end(); ++element )
534 (~lhs).append( element->index(), x[element->index()] * element->value() );
553 template<
typename VT >
554 friend inline auto addAssign( DenseVector<VT,TF>& lhs,
const DVecSVecMultExpr& rhs )
555 -> EnableIf_t< UseAssign_v<VT> >
568 for(
auto element=y.begin(); element!=y.end(); ++element )
569 (~lhs)[element->index()] += x[element->index()] * element->value();
592 template<
typename VT >
593 friend inline auto subAssign( DenseVector<VT,TF>& lhs,
const DVecSVecMultExpr& rhs )
594 -> EnableIf_t< UseAssign_v<VT> >
607 for(
auto element=y.begin(); element!=y.end(); ++element )
608 (~lhs)[element->index()] -= x[element->index()] * element->value();
631 template<
typename VT >
632 friend inline auto multAssign( DenseVector<VT,TF>& lhs,
const DVecSVecMultExpr& rhs )
633 -> EnableIf_t< UseAssign_v<VT> >
646 const auto end( y.end() );
647 auto begin( y.begin() );
652 for( ; i<index; ++i )
658 for( ; i<(~lhs).
size(); ++i )
703 template<
typename VT1
706 , DisableIf_t< IsZero_v<VT2> >* =
nullptr >
707 inline const DVecSVecMultExpr<VT1,VT2,TF>
708 dvecsvecmult(
const DenseVector<VT1,TF>& lhs,
const SparseVector<VT2,TF>& rhs )
714 return DVecSVecMultExpr<VT1,VT2,TF>( ~lhs, ~rhs );
733 template<
typename VT1
736 , EnableIf_t< IsZero_v<VT2> >* =
nullptr >
737 inline decltype(
auto)
738 dvecsvecmult( const DenseVector<VT1,TF>& lhs, const SparseVector<VT2,TF>& rhs )
746 using ReturnType =
const MultTrait_t< ResultType_t<VT1>, ResultType_t<VT2> >;
784 template<
typename VT1
787 inline decltype(
auto)
792 if( (~lhs).
size() != (~rhs).
size() ) {
796 return dvecsvecmult( ~lhs, ~rhs );
811 template<
typename VT1,
typename VT2,
bool TF >
812 struct Size< DVecSVecMultExpr<VT1,VT2,TF>, 0UL >
813 :
public Maximum< Size<VT1,0UL>, Size<VT2,0UL> >
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.This macro encapsulates the default way o...
Definition: Exception.h:235
Pointer difference type of the Blaze library.
LeftOperand vec_
Left-hand side dense vector expression.
Definition: DVecSVecMultExpr.h:292
Header file for auxiliary alias declarations.
const Element operator*() const
Direct access to the sparse vector element at the current iterator position.
Definition: DVecSVecMultExpr.h:222
ConstIterator(LeftOperand vec, IteratorType it)
Constructor for the ConstIterator class.
Definition: DVecSVecMultExpr.h:200
Header file for the UNUSED_PARAMETER function template.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: DVecSVecMultExpr.h:442
ValueType & ReferenceType
Reference return type.
Definition: DVecSVecMultExpr.h:183
Header file for basic type definitions.
Header file for the SparseVector base class.
typename If< Condition, T1, T2 >::Type If_t
Auxiliary alias declaration for the If class template.The If_t alias declaration provides a convenien...
Definition: If.h:109
If_t< IsExpression_v< VT2 >, const VT2, const VT2 &> RightOperand
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecMultExpr.h:159
ConstIterator_t< RemoveReference_t< RightOperand > > IteratorType
Iterator type of the sparse vector expression.
Definition: DVecSVecMultExpr.h:178
typename T::ResultType ResultType_t
Alias declaration for nested ResultType type definitions.The ResultType_t alias declaration provides ...
Definition: Aliases.h:390
Header file for the serial shim.
decltype(std::declval< RN1 >() *std::declval< RN2 >()) ExprReturnType
Expression return type for the subscript operator.
Definition: DVecSVecMultExpr.h:121
const If_t< returnExpr, ExprReturnType, ElementType > ReturnType
Return type for expression template evaluations.
Definition: DVecSVecMultExpr.h:150
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: DVecSVecMultExpr.h:454
const ConstIterator * operator->() const
Direct access to the sparse vector element at the current iterator position.
Definition: DVecSVecMultExpr.h:232
void reset(const DiagonalProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: DiagonalProxy.h:591
If_t< IsExpression_v< VT1 >, const VT1, const VT1 &> LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecMultExpr.h:156
ConstIterator begin() const
Returns an iterator to the first non-zero element of the sparse vector.
Definition: DVecSVecMultExpr.h:344
Constraint on the data type.
ConstIterator lowerBound(size_t index) const
Returns an iterator to the first index not less then the given index.
Definition: DVecSVecMultExpr.h:397
Expression object for dense vector-sparse vector multiplications.The DVecSVecMultExpr class represent...
Definition: DVecSVecMultExpr.h:95
Header file for the Computation base class.
Header file for the reset shim.
Header file for the RequiresEvaluation type trait.
constexpr void UNUSED_PARAMETER(const Args &...)
Suppression of unused parameter warnings.
Definition: Unused.h:81
ElementType_t< ResultType > ElementType
Resulting element type.
Definition: DVecSVecMultExpr.h:147
typename T::ReturnType ReturnType_t
Alias declaration for nested ReturnType type definitions.The ReturnType_t alias declaration provides ...
Definition: Aliases.h:410
size_t size() const noexcept
Returns the current size/dimension of the vector.
Definition: DVecSVecMultExpr.h:364
RightOperand rhs_
Right-hand side sparse vector of the multiplication expression.
Definition: DVecSVecMultExpr.h:462
typename T::ElementType ElementType_t
Alias declaration for nested ElementType type definitions.The ElementType_t alias declaration provide...
Definition: Aliases.h:170
LeftOperand leftOperand() const noexcept
Returns the left-hand side dense vector operand.
Definition: DVecSVecMultExpr.h:420
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two expression iterators.
Definition: DVecSVecMultExpr.h:285
Header file for the Maximum class template.
typename EnableIf< Condition, T >::Type EnableIf_t
Auxiliary type for the EnableIf class template.The EnableIf_t alias declaration provides a convenient...
Definition: EnableIf.h:138
ConstIterator find(size_t index) const
Searches for a specific vector element.
Definition: DVecSVecMultExpr.h:385
Header file for the ValueIndexPair class.
ValueType * PointerType
Pointer return type.
Definition: DVecSVecMultExpr.h:182
ConstIterator & operator++()
Pre-increment operator.
Definition: DVecSVecMultExpr.h:211
Header file for the DisableIf class template.
Header file for the IsTemporary type trait class.
Header file for the multiplication trait.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
Header file for the If class template.
If_t< useAssign, const ResultType, const DVecSVecMultExpr &> CompositeType
Data type for composite expression templates.
Definition: DVecSVecMultExpr.h:153
#define BLAZE_CONSTRAINT_MUST_BE_ZERO_TYPE(T)
Constraint on the data type.In case the given data type T is not a zero vector or matrix type...
Definition: Zero.h:61
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exception.This macro encapsulates the default way of Bl...
Definition: Exception.h:331
CompositeType_t< VT2 > CT2
Composite type of the right-hand side sparse vector expression.
Definition: DVecSVecMultExpr.h:106
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:3086
MultTrait_t< RT1, RT2 > ResultType
Result type for expression template evaluations.
Definition: DVecSVecMultExpr.h:145
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:76
#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:61
ConstIterator end() const
Returns an iterator just past the last non-zero element of the sparse vector.
Definition: DVecSVecMultExpr.h:354
CompositeType_t< VT1 > CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecSVecMultExpr.h:105
RightOperand rightOperand() const noexcept
Returns the right-hand side sparse vector operand.
Definition: DVecSVecMultExpr.h:430
Constraint on the data type.
Constraint on the data type.
Header file for the exception macros of the math module.
Header file for the VecVecMultExpr base class.
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DVecSVecMultExpr.h:274
Constraint on the data type.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Header file for the EnableIf class template.
DVecSVecMultExpr(const VT1 &lhs, const VT2 &rhs) noexcept
Constructor for the DVecSVecMultExpr class.
Definition: DVecSVecMultExpr.h:304
typename MultTrait< T1, T2 >::Type MultTrait_t
Auxiliary alias declaration for the MultTrait class template.The MultTrait_t alias declaration provid...
Definition: MultTrait.h:240
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DVecSVecMultExpr.h:263
IteratorCategory iterator_category
The iterator category.
Definition: DVecSVecMultExpr.h:187
typename T::TransposeType TransposeType_t
Alias declaration for nested TransposeType type definitions.The TransposeType_t alias declaration pro...
Definition: Aliases.h:470
Iterator over the elements of the dense vector-sparse vector multiplication expression.
Definition: DVecSVecMultExpr.h:170
Header file for run time assertion macros.
typename T::CompositeType CompositeType_t
Alias declaration for nested CompositeType type definitions.The CompositeType_t alias declaration pro...
Definition: Aliases.h:90
ValueIndexPair< ElementType > Element
Element type of the sparse vector expression.
Definition: DVecSVecMultExpr.h:175
TransposeType_t< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: DVecSVecMultExpr.h:146
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: DVecSVecMultExpr.h:331
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DVecSVecMultExpr.h:184
ConstIterator upperBound(size_t index) const
Returns an iterator to the first index greater then the given index.
Definition: DVecSVecMultExpr.h:409
Header file for the IsZero type trait.
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:94
constexpr size_t size(const Matrix< MT, SO > &matrix) noexcept
Returns the total number of elements of the matrix.
Definition: Matrix.h:530
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: DVecSVecMultExpr.h:318
decltype(auto) serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:808
#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:81
TransposeType_t< VT2 > TT2
Transpose type of the right-hand side sparse vector expression.
Definition: DVecSVecMultExpr.h:108
Header file for the RemoveReference type trait.
Element ValueType
Type of the underlying pointers.
Definition: DVecSVecMultExpr.h:181
ReturnType value() const
Access to the current value of the sparse element.
Definition: DVecSVecMultExpr.h:242
std::forward_iterator_tag IteratorCategory
The iterator category.
Definition: DVecSVecMultExpr.h:180
typename T::ConstIterator ConstIterator_t
Alias declaration for nested ConstIterator type definitions.The ConstIterator_t alias declaration pro...
Definition: Aliases.h:110
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3081
#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:61
TransposeType_t< VT1 > TT1
Transpose type of the left-hand side dense vector expression.
Definition: DVecSVecMultExpr.h:107
Index-value-pair for sparse vectors and matrices.The ValueIndexPair class represents a single index-v...
Definition: ValueIndexPair.h:73
Constraint on the data type.
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_VECVECMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/vector ...
Definition: VecVecMultExpr.h:104
static constexpr bool returnExpr
Compilation switch for the selection of the subscript operator return type.
Definition: DVecSVecMultExpr.h:118
LeftOperand lhs_
Left-hand side dense vector of the multiplication expression.
Definition: DVecSVecMultExpr.h:461
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:66
ResultType_t< VT2 > RT2
Result type of the right-hand side sparse vector expression.
Definition: DVecSVecMultExpr.h:102
ResultType_t< VT1 > RT1
Result type of the left-hand side dense vector expression.
Definition: DVecSVecMultExpr.h:101
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:138
size_t nonZeros() const
Returns the number of non-zero elements in the sparse vector.
Definition: DVecSVecMultExpr.h:374
ReturnType_t< VT1 > RN1
Return type of the left-hand side dense vector expression.
Definition: DVecSVecMultExpr.h:103
static constexpr bool useAssign
Compilation switch for the evaluation strategy of the multiplication expression.
Definition: DVecSVecMultExpr.h:132
IteratorType it_
Iterator over the elements of the right-hand side sparse vector expression.
Definition: DVecSVecMultExpr.h:293
Header file for the Size type trait.
#define BLAZE_CONSTRAINT_MUST_NOT_BE_ZERO_TYPE(T)
Constraint on the data type.In case the given data type T is a zero vector or matrix type...
Definition: Zero.h:81
#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:63
Base class for all vector/vector multiplication expression templates.The VecVecMultExpr class serves ...
Definition: VecVecMultExpr.h:66
#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
static constexpr bool smpAssignable
Compilation switch for the expression template assignment strategy.
Definition: DVecSVecMultExpr.h:164
ReturnType_t< VT2 > RN2
Return type of the right-hand side sparse vector expression.
Definition: DVecSVecMultExpr.h:104
If_t< Less< T1, T2 >::value, T2, T1 > Maximum
Compile time value evaluation.The Maximum alias declaration selects the larger of the two given templ...
Definition: Maximum.h:73
Header file for the IsExpression type trait class.
Header file for the function trace functionality.
size_t index() const
Access to the current index of the sparse element.
Definition: DVecSVecMultExpr.h:252