35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTSMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTSMATMULTEXPR_H_
99 template<
typename MT1
138 enum { vectorizable = 0 };
171 if(
lhs_.columns() == 0UL )
179 const ConstIterator end( B.end(j) );
180 ConstIterator element( B.begin(j) );
187 tmp =
lhs_(i,element->index()) * element->value();
189 for( ; element!=end; ++element )
190 tmp +=
lhs_(i,element->index()) * element->value();
196 for(
size_t k=1UL; k<
lhs_.columns(); ++k ) {
221 return rhs_.columns();
251 template<
typename T >
253 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
263 template<
typename T >
265 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
287 template<
typename MT
305 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
306 for(
size_t j=0UL; j<(~lhs).
columns(); ++j )
308 ConstIterator element( B.begin(j) );
309 const ConstIterator end( B.end(j) );
311 if( element == end ) {
312 reset( (~lhs)(i,j) );
316 (~lhs)(i,j) = A(i,element->index()) * element->value();
318 for( ; element!=end; ++element )
319 (~lhs)(i,j) += A(i,element->index()) * element->value();
339 template<
typename MT
357 const TmpType tmp( rhs );
376 template<
typename MT
394 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
395 for(
size_t j=0UL; j<(~lhs).
columns(); ++j )
397 ConstIterator element( B.begin(j) );
398 const ConstIterator end( B.end(j) );
400 for( ; element!=end; ++element )
401 (~lhs)(i,j) += A(i,element->index()) * element->value();
425 template<
typename MT
443 for(
size_t i=0UL; i<(~lhs).
rows(); ++i ) {
444 for(
size_t j=0UL; j<(~lhs).
columns(); ++j )
446 ConstIterator element( B.begin(j) );
447 const ConstIterator end( B.end(j) );
449 for( ; element!=end; ++element )
450 (~lhs)(i,j) -= A(i,element->index()) * element->value();
519 template<
typename T1
521 inline const DMatTSMatMultExpr<T1,T2>
527 throw std::invalid_argument(
"Matrix sizes do not match" );
544 template<
typename MT1,
typename MT2,
typename VT >
545 struct DMatDVecMultExprTrait< DMatTSMatMultExpr<MT1,MT2>, VT >
549 typedef typename SelectType< IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
550 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value &&
551 IsDenseVector<VT>::value && IsColumnVector<VT>::value
552 ,
typename DMatDVecMultExprTrait< MT1, typename TSMatDVecMultExprTrait<MT2,VT>::Type >::Type
553 , INVALID_TYPE >::Type Type;
562 template<
typename MT1,
typename MT2,
typename VT >
563 struct DMatSVecMultExprTrait< DMatTSMatMultExpr<MT1,MT2>, VT >
567 typedef typename SelectType< IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
568 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value &&
569 IsSparseVector<VT>::value && IsColumnVector<VT>::value
570 ,
typename DMatSVecMultExprTrait< MT1, typename TSMatSVecMultExprTrait<MT2,VT>::Type >::Type
571 , INVALID_TYPE >::Type Type;
580 template<
typename VT,
typename MT1,
typename MT2 >
581 struct TDVecDMatMultExprTrait< VT, DMatTSMatMultExpr<MT1,MT2> >
585 typedef typename SelectType< IsDenseVector<VT>::value && IsRowVector<VT>::value &&
586 IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
587 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value
588 ,
typename TDVecTSMatMultExprTrait< typename TDVecDMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
589 , INVALID_TYPE >::Type Type;
598 template<
typename VT,
typename MT1,
typename MT2 >
599 struct TSVecDMatMultExprTrait< VT, DMatTSMatMultExpr<MT1,MT2> >
603 typedef typename SelectType< IsSparseVector<VT>::value && IsRowVector<VT>::value &&
604 IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
605 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value
606 ,
typename TDVecTSMatMultExprTrait< typename TSVecDMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
607 , INVALID_TYPE >::Type Type;
616 template<
typename MT1,
typename MT2 >
617 struct SubmatrixExprTrait< DMatTSMatMultExpr<MT1,MT2> >
621 typedef typename MultExprTrait< typename SubmatrixExprTrait<const MT1>::Type
622 ,
typename SubmatrixExprTrait<const MT2>::Type >::Type Type;
631 template<
typename MT1,
typename MT2 >
632 struct RowExprTrait< DMatTSMatMultExpr<MT1,MT2> >
636 typedef typename MultExprTrait< typename RowExprTrait<const MT1>::Type, MT2 >::Type Type;
645 template<
typename MT1,
typename MT2 >
646 struct ColumnExprTrait< DMatTSMatMultExpr<MT1,MT2> >
650 typedef typename MultExprTrait< MT1, typename ColumnExprTrait<const MT2>::Type >::Type Type;
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTSMatMultExpr.h:119
void reset(DynamicMatrix< Type, SO > &m)
Resetting the given dense matrix.
Definition: DynamicMatrix.h:4512
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTSMatMultExpr.h:162
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
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTSMatMultExpr.h:124
MT2::CompositeType CT2
Composite type of the right-hand side sparse matrix expression.
Definition: DMatTSMatMultExpr.h:110
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTSMatMultExpr.h:252
Header file for the IsSparseMatrix type trait.
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:196
#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
Header file for the ColumnExprTrait class template.
Header file for the IsColumnMajorMatrix type trait.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTSMatMultExpr.h:118
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2375
Header file for the IsRowVector type trait.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:248
Header file for the Computation base class.
Header file for the MatMatMultExpr base class.
Header file for the RequiresEvaluation type trait.
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatTSMatMultExpr.h:230
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.
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTSMatMultExpr.h:121
Constraint on the data type.
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
DMatTSMatMultExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTSMatMultExpr class.
Definition: DMatTSMatMultExpr.h:147
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 multiplication trait.
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTSMatMultExpr.h:210
SelectType< IsComputation< MT2 >::value, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side sparse matrix operand.
Definition: DMatTSMatMultExpr.h:133
#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
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTSMatMultExpr.h:109
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2379
Header file for the TDVecTSMatMultExprTrait class template.
SelectType< IsComputation< MT1 >::value, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side dense matrix operand.
Definition: DMatTSMatMultExpr.h:130
Header file for the DenseMatrix base class.
Header file for the TSMatDVecMultExprTrait class template.
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
Header file for the DMatDVecMultExprTrait class template.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTSMatMultExpr.h:107
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTSMatMultExpr.h:117
#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.
RightOperand rhs_
Right-hand side sparse matrix of the multiplication expression.
Definition: DMatTSMatMultExpr.h:272
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
Header file for the IsDenseMatrix type trait.
Base class for all matrix/matrix multiplication expression templates.The MatMatMultExpr class serves ...
Definition: MatMatMultExpr.h:65
Header file for the IsSparseVector type trait.
Header file for the SubmatrixExprTrait class template.
#define BLAZE_CONSTRAINT_MUST_BE_ROW_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a row-major dense or sparse matrix t...
Definition: StorageOrder.h:81
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.
Utility type for generic codes.
Base template for the MultTrait class.
Definition: MultTrait.h:141
const ElementType ReturnType
Return type for expression template evaluations.
Definition: DMatTSMatMultExpr.h:120
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTSMatMultExpr.h:116
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
RightOperand rightOperand() const
Returns the right-hand side transpose sparse matrix operand.
Definition: DMatTSMatMultExpr.h:240
Header file for the reset shim.
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
Expression object for dense matrix-transpose sparse matrix multiplications.The DMatTSMatMultExpr clas...
Definition: DMatTSMatMultExpr.h:101
Header file for the RemoveReference type trait.
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatTSMatMultExpr.h:264
#define BLAZE_CONSTRAINT_MATRICES_MUST_HAVE_SAME_STORAGE_ORDER(T1, T2)
Constraint on the data type.In case either of the two given data types T1 or T2 is not a matrix type ...
Definition: StorageOrder.h:283
Header file for the IsDenseVector type trait.
LeftOperand lhs_
Left-hand side dense matrix of the multiplication expression.
Definition: DMatTSMatMultExpr.h:271
Header file for the IsRowMajorMatrix type trait.
Header file for the IsComputation type trait class.
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
Header file for the TDVecDMatMultExprTrait class template.
#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
size_t columns(const Matrix< MT, SO > &m)
Returns the current number of columns of the matrix.
Definition: Matrix.h:154
Header file for basic type definitions.
Header file for the TSVecDMatMultExprTrait class template.
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTSMatMultExpr.h:220
MT2::ResultType RT2
Result type of the right-hand side sparse matrix expression.
Definition: DMatTSMatMultExpr.h:108
Header file for the DMatSVecMultExprTrait class template.
Header file for the IsColumnVector type trait.
size_t rows(const Matrix< MT, SO > &m)
Returns the current number of rows of the matrix.
Definition: Matrix.h:138
#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
DMatTSMatMultExpr< MT1, MT2 > This
Type of this DMatTSMatMultExpr instance.
Definition: DMatTSMatMultExpr.h:115
#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
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side sparse matrix expression.
Definition: DMatTSMatMultExpr.h:127
Header file for the IsExpression type trait class.
Header file for the TSMatSVecMultExprTrait class template.
Header file for the FunctionTrace class.