35 #ifndef _BLAZE_MATH_EXPRESSIONS_TSVECTDMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_TSVECTDMATMULTEXPR_H_
101 class TSVecTDMatMultExpr :
public DenseVector< TSVecTDMatMultExpr<VT,MT>, true >
102 ,
private TVecMatMultExpr
103 ,
private Computation
131 enum { useAssign = ( evaluateVector || evaluateMatrix ) };
136 template<
typename VT2 >
139 enum { value = useAssign };
150 template<
typename T1 >
151 struct UseSMPAssign {
152 enum { value = useAssign };
183 enum { vectorizable = 0 };
186 enum { smpAssignable = !evaluateVector && VT::smpAssignable &&
187 !evaluateMatrix && MT::smpAssignable };
228 if( element != end ) {
229 res = element->value() *
mat_( element->index(), index );
231 for( ; element!=
end; ++element ) {
232 res += element->value() *
mat_( element->index(), index );
246 return mat_.columns();
276 template<
typename T >
278 return vec_.isAliased( alias ) ||
mat_.isAliased( alias );
288 template<
typename T >
290 return vec_.isAliased( alias ) ||
mat_.isAliased( alias );
300 return mat_.isAligned();
336 template<
typename VT2 >
347 LT x(
serial( rhs.vec_ ) );
348 if( x.nonZeros() == 0UL ) {
354 RT A(
serial( rhs.mat_ ) );
384 template<
typename VT2 >
396 const ResultType tmp(
serial( rhs ) );
417 template<
typename VT2 >
418 friend inline typename EnableIf< UseAssign<VT2> >::Type
428 LT x(
serial( rhs.vec_ ) );
429 if( x.nonZeros() == 0UL )
return;
432 RT A(
serial( rhs.mat_ ) );
464 template<
typename VT2 >
475 LT x(
serial( rhs.vec_ ) );
476 if( x.nonZeros() == 0UL )
return;
479 RT A(
serial( rhs.mat_ ) );
511 template<
typename VT2 >
523 const ResultType tmp(
serial( rhs ) );
547 template<
typename VT2 >
559 if( x.nonZeros() == 0UL ) {
594 template<
typename VT2 >
595 friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
606 const ResultType tmp( rhs );
626 template<
typename VT2 >
627 friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
638 if( x.nonZeros() == 0UL )
return;
672 template<
typename VT2 >
684 if( x.nonZeros() == 0UL )
return;
719 template<
typename VT2 >
731 const ResultType tmp( rhs );
793 template<
typename T1,
typename T2 >
794 inline const typename DisableIf< Or< IsSymmetric<T2>, IsMatMatMultExpr<T2> >
795 , TSVecTDMatMultExpr<T1,T2> >::Type
800 if( (~vec).
size() != (~mat).
rows() )
801 throw std::invalid_argument(
"Vector and matrix sizes do not match" );
831 template<
typename T1
833 inline const typename EnableIf< IsSymmetric<T2>,
typename MultExprTrait<T1,T2>::Type >::Type
834 operator*(
const SparseVector<T1,true>& vec,
const DenseMatrix<T2,true>& mat )
840 if( (~vec).
size() != (~mat).
rows() )
841 throw std::invalid_argument(
"Vector and matrix sizes do not match" );
843 return (~vec) *
trans( ~mat );
859 template<
typename MT,
typename VT >
877 template<
typename VT,
typename MT,
bool AF >
882 typedef typename MultExprTrait< typename SubvectorExprTrait<const VT,AF>::Type
883 ,
typename SubmatrixExprTrait<const MT,AF>::Type >::Type Type;
LeftOperand vec_
Left-hand side sparse vector of the multiplication expression.
Definition: TSVecTDMatMultExpr.h:316
const ElementType ReturnType
Return type for expression template evaluations.
Definition: TSVecTDMatMultExpr.h:163
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:8247
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:258
SelectType< evaluateVector, const VRT, VCT >::Type LT
Type for the assignment of the left-hand side sparse vector operand.
Definition: TSVecTDMatMultExpr.h:175
size_t size() const
Returns the current size/dimension of the vector.
Definition: TSVecTDMatMultExpr.h:245
SelectType< evaluateMatrix, const MRT, MCT >::Type RT
Type for the assignment of the left-hand side dense matrix operand.
Definition: TSVecTDMatMultExpr.h:178
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:209
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a dense, N-dimensional matrix type...
Definition: DenseMatrix.h:79
void reset(const DiagonalProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: DiagonalProxy.h:821
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2507
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:261
BLAZE_ALWAYS_INLINE size_t rows(const Matrix< MT, SO > &matrix)
Returns the current number of rows of the matrix.
Definition: Matrix.h:316
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: TSVecTDMatMultExpr.h:265
Header file for the DenseVector base class.
Compile time check for lower triangular matrices.This type trait tests whether or not the given templ...
Definition: IsLower.h:90
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:699
Header file for the Computation base class.
MT::ResultType MRT
Result type of the right-hand side dense matrix expression.
Definition: TSVecTDMatMultExpr.h:108
Compile time check for upper triangular matrices.This type trait tests whether or not the given templ...
Definition: IsUpper.h:90
Header file for the RequiresEvaluation type trait.
ResultType::ElementType ElementType
Resulting element type.
Definition: TSVecTDMatMultExpr.h:162
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.
VT::CompositeType VCT
Composite type of the left-hand side sparse vector expression.
Definition: TSVecTDMatMultExpr.h:109
Constraint on the data type.
friend EnableIf< UseSMPAssign< VT2 > >::Type smpMultAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
SMP multiplication assignment of a transpose sparse vector-transpose dense matrix multiplication to a...
Definition: TSVecTDMatMultExpr.h:721
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
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:263
friend EnableIf< UseSMPAssign< VT2 > >::Type smpSubAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
SMP subtraction assignment of a transpose sparse vector-transpose dense matrix multiplication to a de...
Definition: TSVecTDMatMultExpr.h:674
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.
Header file for the IsStrictlyUpper type trait.
Header file for the IsSymmetric type trait.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
#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: StorageOrder.h:161
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: TSVecTDMatMultExpr.h:299
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2511
Header file for the Or class template.
friend EnableIf< UseAssign< VT2 > >::Type multAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
Multiplication assignment of a transpose sparse vector-transpose dense matrix multiplication to a den...
Definition: TSVecTDMatMultExpr.h:513
Header file for the IsMatMatMultExpr type trait class.
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
const size_t SMP_TSVECTDMATMULT_THRESHOLD
SMP sparse vector/column-major dense matrix multiplication threshold.This threshold specifies when a ...
Definition: Thresholds.h:483
Header file for the Columns type trait.
SelectType< useAssign, const ResultType, const TSVecTDMatMultExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: TSVecTDMatMultExpr.h:166
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 IsLower type trait.
#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.
#define BLAZE_CONSTRAINT_MUST_NOT_BE_MATMATMULTEXPR_TYPE(T)
Constraint on the data type.In case the given data type T is a matrix/matrix multiplication expressio...
Definition: MatMatMultExpr.h:126
#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
Constraints on the storage order of matrix types.
Constraint on the data type.
Compile time check for strictly upper triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyUpper.h:86
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Header file for the EnableIf class template.
Header file for the IsStrictlyLower type trait.
Header file for the serial shim.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: TSVecTDMatMultExpr.h:277
TSVecTDMatMultExpr< VT, MT > This
Type of this TSVecTDMatMultExpr instance.
Definition: TSVecTDMatMultExpr.h:159
#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.
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: TSVecTDMatMultExpr.h:161
Header file for run time assertion macros.
EnableIf< IsDenseMatrix< MT1 > >::Type smpAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:98
friend EnableIf< UseAssign< VT2 > >::Type subAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
Subtraction assignment of a transpose sparse vector-transpose dense matrix multiplication to a dense ...
Definition: TSVecTDMatMultExpr.h:466
MultTrait< VRT, MRT >::Type ResultType
Result type for expression template evaluations.
Definition: TSVecTDMatMultExpr.h:160
friend EnableIf< UseAssign< VT2 > >::Type addAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
Addition assignment of a transpose sparse vector-transpose dense matrix multiplication to a dense vec...
Definition: TSVecTDMatMultExpr.h:419
Base template for the MultTrait class.
Definition: MultTrait.h:150
MT::CompositeType MCT
Composite type of the right-hand side dense matrix expression.
Definition: TSVecTDMatMultExpr.h:110
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
Header file for the reset shim.
Constraint on the data type.
TSVecTDMatMultExpr(const VT &vec, const MT &mat)
Constructor for the TSVecTDMatMultExpr class.
Definition: TSVecTDMatMultExpr.h:196
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: TSVecTDMatMultExpr.h:289
Header file for the TVecMatMultExpr base class.
RightOperand mat_
Right-hand side dense matrix of the multiplication expression.
Definition: TSVecTDMatMultExpr.h:317
Header file for the RemoveReference type trait.
Expression object for transpose sparse vector-transpose dense matrix multiplications.The TSVecTDMatMultExpr class represents the compile time expression for multiplications between transpose sparse vectors and column-major dense matrices.
Definition: Forward.h:149
LeftOperand leftOperand() const
Returns the left-hand side sparse vector operand.
Definition: TSVecTDMatMultExpr.h:255
Compile time check for strictly lower triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyLower.h:86
#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
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: TSVecTDMatMultExpr.h:210
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:937
Header file for the IsComputation type trait class.
friend EnableIf< UseSMPAssign< VT2 > >::Type smpAddAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
SMP addition assignment of a transpose sparse vector-transpose dense matrix multiplication to a dense...
Definition: TSVecTDMatMultExpr.h:628
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:2502
#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.
Header file for the IsUpper type trait.
Constraint on the data type.
SelectType< IsExpression< MT >::value, const MT, const MT & >::Type RightOperand
Composite type of the right-hand side sparse matrix expression.
Definition: TSVecTDMatMultExpr.h:172
SelectType< IsExpression< VT >::value, const VT, const VT & >::Type LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: TSVecTDMatMultExpr.h:169
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
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: TSVecTDMatMultExpr.h:309
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
VT::ResultType VRT
Result type of the left-hand side sparse vector expression.
Definition: TSVecTDMatMultExpr.h:107