35 #ifndef _BLAZE_MATH_EXPRESSIONS_TDVECTSMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_TDVECTSMATMULTEXPR_H_
94 class TDVecTSMatMultExpr :
public DenseVector< TDVecTSMatMultExpr<VT,MT>, true >
95 ,
private TVecMatMultExpr
124 enum { useAssign = ( evaluateVector || evaluateMatrix ) };
129 template<
typename VT2 >
132 enum { value = useAssign };
143 template<
typename T1 >
144 struct UseSMPAssign {
145 enum { value = useAssign };
176 enum { vectorizable = 0 };
179 enum { smpAssignable = !evaluateVector && VT::smpAssignable &&
180 !evaluateMatrix && MT::smpAssignable };
212 const ConstIterator
end( A.end(index) );
213 ConstIterator element( A.begin(index) );
216 if( element !=
end ) {
217 res = x[element->index()] * element->value();
219 for( ; element!=
end; ++element )
220 res += x[element->index()] * element->value();
237 inline ReturnType
at(
size_t index )
const {
238 if( index >=
mat_.columns() ) {
241 return (*
this)[index];
251 return mat_.columns();
281 template<
typename T >
283 return (
vec_.isAliased( alias ) ||
mat_.isAliased( alias ) );
293 template<
typename T >
295 return (
vec_.isAliased( alias ) ||
mat_.isAliased( alias ) );
305 return vec_.isAligned();
315 return (
size() > SMP_TDVECTSMATMULT_THRESHOLD );
340 template<
typename VT2 >
350 if( rhs.mat_.rows() == 0UL ) {
355 LT x(
serial( rhs.vec_ ) );
356 RT A(
serial( rhs.mat_ ) );
382 template<
typename VT2 >
394 const ResultType tmp(
serial( rhs ) );
414 template<
typename VT2 >
424 if( rhs.mat_.rows() == 0UL ) {
428 LT x(
serial( rhs.vec_ ) );
429 RT A(
serial( rhs.mat_ ) );
459 template<
typename VT2 >
469 if( rhs.mat_.rows() == 0UL ) {
473 LT x(
serial( rhs.vec_ ) );
474 RT A(
serial( rhs.mat_ ) );
504 template<
typename VT2 >
516 const ResultType tmp(
serial( rhs ) );
539 template<
typename VT2 >
549 if( rhs.mat_.rows() == 0UL ) {
580 template<
typename VT2 >
592 const ResultType tmp( rhs );
611 template<
typename VT2 >
621 if( rhs.mat_.rows() == 0UL ) {
655 template<
typename VT2 >
665 if( rhs.mat_.rows() == 0UL ) {
699 template<
typename VT2 >
711 const ResultType tmp( rhs );
772 template<
typename T1
774 inline const typename DisableIf< IsMatMatMultExpr<T2>, TDVecTSMatMultExpr<T1,T2> >::Type
779 if( (~vec).
size() != (~mat).
rows() ) {
798 template<
typename VT,
typename MT >
815 template<
typename VT,
typename MT >
832 template<
typename VT,
typename MT,
bool AF >
837 typedef typename MultExprTrait< typename SubvectorExprTrait<const VT,AF>::Type
838 ,
typename SubmatrixExprTrait<const MT,AF>::Type >::Type Type;
SelectType< evaluateVector, const VRT, VCT >::Type LT
Composite type of the left-hand side dense vector expression.
Definition: TDVecTSMatMultExpr.h:168
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exceptionThis macro encapsulates the default way of...
Definition: Exception.h:187
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
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:7820
Header file for basic type definitions.
BLAZE_ALWAYS_INLINE MT::Iterator end(Matrix< MT, SO > &matrix, size_t i)
Returns an iterator just past the last element of row/column i.
Definition: Matrix.h:250
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:207
const ElementType ReturnType
Return type for expression template evaluations.
Definition: TDVecTSMatMultExpr.h:156
Expression object for transpose dense vector-transpose sparse matrix multiplications.The TDVecTSMatMultExpr class represents the compile time expression for multiplications between transpose dense vectors and column-major sparse matrices.
Definition: Forward.h:149
void reset(const DiagonalProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: DiagonalProxy.h:507
RightOperand mat_
Right-hand side sparse matrix of the multiplication expression.
Definition: TDVecTSMatMultExpr.h:322
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2588
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:259
BLAZE_ALWAYS_INLINE size_t rows(const Matrix< MT, SO > &matrix)
Returns the current number of rows of the matrix.
Definition: Matrix.h:308
Header file for the DenseVector base class.
TDVecTSMatMultExpr< VT, MT > This
Type of this TDVecTSMatMultExpr instance.
Definition: TDVecTSMatMultExpr.h:152
friend EnableIf< UseSMPAssign< VT2 > >::Type smpAssign(SparseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
SMP assignment of a transpose dense vector-transpose sparse matrix multiplication to a sparse vector ...
Definition: TDVecTSMatMultExpr.h:582
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:721
Header file for the Computation base class.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: TDVecTSMatMultExpr.h:282
Header file for the RequiresEvaluation type trait.
VT::ResultType VRT
Result type of the left-hand side dense vector expression.
Definition: TDVecTSMatMultExpr.h:100
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:117
Constraint on the data type.
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: TDVecTSMatMultExpr.h:200
MT::ResultType MRT
Result type of the right-hand side sparse matrix expression.
Definition: TDVecTSMatMultExpr.h:101
friend EnableIf< UseAssign< VT2 > >::Type addAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
Addition assignment of a transpose dense vector-transpose sparse matrix multiplication to a dense vec...
Definition: TDVecTSMatMultExpr.h:416
Constraint on the data type.
Header file for the MultExprTrait class template.
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:90
friend EnableIf< UseSMPAssign< VT2 > >::Type smpMultAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
SMP multiplication assignment of a transpose dense vector-transpose sparse matrix multiplication to a...
Definition: TDVecTSMatMultExpr.h:701
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 multiplication trait.
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: TDVecTSMatMultExpr.h:314
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
friend EnableIf< UseSMPAssign< VT2 > >::Type smpAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
SMP assignment of a transpose dense vector-transpose sparse matrix multiplication to a dense vector (...
Definition: TDVecTSMatMultExpr.h:541
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a column-major dense or sparse matri...
Definition: ColumnMajorMatrix.h:79
LeftOperand leftOperand() const
Returns the left-hand side dense vector operand.
Definition: TDVecTSMatMultExpr.h:260
size_t size() const
Returns the current size/dimension of the vector.
Definition: TDVecTSMatMultExpr.h:250
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2592
Header file for the IsMatMatMultExpr type trait class.
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exceptionThis macro encapsulates the default way of Bla...
Definition: Exception.h:331
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: TDVecTSMatMultExpr.h:294
Header file for the Columns type trait.
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: TDVecTSMatMultExpr.h:304
friend EnableIf< UseAssign< VT2 > >::Type assign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
Assignment of a transpose dense vector-transpose sparse matrix multiplication to a dense vector ( )...
Definition: TDVecTSMatMultExpr.h:342
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
Header file for the IsAligned type trait.
LeftOperand vec_
Left-hand side dense vector of the multiplication expression.
Definition: TDVecTSMatMultExpr.h:321
RightOperand rightOperand() const
Returns the right-hand side transpose sparse matrix operand.
Definition: TDVecTSMatMultExpr.h:270
#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
SelectType< IsExpression< VT >::value, const VT, const VT & >::Type LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: TDVecTSMatMultExpr.h:162
friend EnableIf< UseSMPAssign< VT2 > >::Type smpSubAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
SMP subtraction assignment of a transpose dense vector-transpose sparse matrix multiplication to a de...
Definition: TDVecTSMatMultExpr.h:657
Constraint on the data type.
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
ResultType::ElementType ElementType
Resulting element type.
Definition: TDVecTSMatMultExpr.h:155
Header file for the EnableIf class template.
Header file for the serial shim.
friend EnableIf< UseAssign< VT2 > >::Type assign(SparseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
Assignment of a transpose dense vector-transpose sparse matrix multiplication to a sparse vector ( )...
Definition: TDVecTSMatMultExpr.h:384
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: TDVecTSMatMultExpr.h:154
Header file for the SubmatrixExprTrait class template.
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
Header file for run time assertion macros.
Base template for the MultTrait class.
Definition: MultTrait.h:138
friend EnableIf< UseAssign< VT2 > >::Type multAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
Multiplication assignment of a transpose dense vector-transpose sparse matrix multiplication to a den...
Definition: TDVecTSMatMultExpr.h:506
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: TDVecTSMatMultExpr.h:237
Header file for the reset shim.
Constraint on the data type.
friend EnableIf< UseSMPAssign< VT2 > >::Type smpAddAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
SMP addition assignment of a transpose dense vector-transpose sparse matrix multiplication to a dense...
Definition: TDVecTSMatMultExpr.h:613
Header file for the TVecMatMultExpr base class.
Constraints on the storage order of matrix types.
MT::CompositeType MCT
Composite type of the right-hand side sparse matrix expression.
Definition: TDVecTSMatMultExpr.h:103
Header file for the RemoveReference type trait.
SelectType< IsExpression< MT >::value, const MT, const MT & >::Type RightOperand
Composite type of the right-hand side sparse matrix expression.
Definition: TDVecTSMatMultExpr.h:165
TDVecTSMatMultExpr(const VT &vec, const MT &mat)
Constructor for the TDVecTSMatMultExpr class.
Definition: TDVecTSMatMultExpr.h:186
#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
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_TVECMATMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/matrix ...
Definition: TVecMatMultExpr.h:166
MultTrait< VRT, MRT >::Type ResultType
Result type for expression template evaluations.
Definition: TDVecTSMatMultExpr.h:153
Header file for the IsComputation type trait class.
SelectType< useAssign, const ResultType, const TDVecTSMatMultExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: TDVecTSMatMultExpr.h:159
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:118
VT::CompositeType VCT
Composite type of the left-hand side dense vector expression.
Definition: TDVecTSMatMultExpr.h:102
#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:2583
#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: RowVector.h:79
Header file for the SubvectorExprTrait class template.
Header file for exception macros.
SelectType< evaluateMatrix, const MRT, MCT >::Type RT
Composite type of the right-hand side sparse matrix expression.
Definition: TDVecTSMatMultExpr.h:171
Header file for the Size type trait.
Header file for the thresholds for matrix/vector and matrix/matrix multiplications.
#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
#define BLAZE_CONSTRAINT_MUST_BE_SPARSE_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a sparse, N-dimensional matrix type...
Definition: SparseMatrix.h:79
Constraint on the transpose flag of vector types.
Header file for the IsExpression type trait class.
friend EnableIf< UseAssign< VT2 > >::Type subAssign(DenseVector< VT2, true > &lhs, const TDVecTSMatMultExpr &rhs)
Subtraction assignment of a transpose dense vector-transpose sparse matrix multiplication to a dense ...
Definition: TDVecTSMatMultExpr.h:461
Header file for the FunctionTrace class.