35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTSMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTSMATMULTEXPR_H_
110 template<
typename MT1
145 template<
typename T1,
typename T2,
typename T3 >
146 struct CanExploitSymmetry {
147 enum { value = ( IsRowMajorMatrix<T1>::value && IsSymmetric<T3>::value ) ||
148 ( IsColumnMajorMatrix<T1>::value && IsSymmetric<T2>::value ) };
159 template<
typename T1,
typename T2,
typename T3 >
160 struct IsEvaluationRequired {
161 enum { value = ( evaluateLeft || evaluateRight ) &&
162 !CanExploitSymmetry<T1,T2,T3>::value };
192 enum { vectorizable = 0 };
195 enum { smpAssignable = !evaluateLeft && MT1::smpAssignable &&
196 !evaluateRight && MT2::smpAssignable };
229 if(
lhs_.columns() == 0UL )
237 const ConstIterator
end( B.end(j) );
238 ConstIterator element( B.begin(j) );
245 tmp =
lhs_(i,element->index()) * element->value();
247 for( ; element!=
end; ++element )
248 tmp +=
lhs_(i,element->index()) * element->value();
254 for(
size_t k=1UL; k<
lhs_.columns(); ++k ) {
279 return rhs_.columns();
309 template<
typename T >
311 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
321 template<
typename T >
323 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
333 return lhs_.isAligned();
366 template<
typename MT
383 DMatTSMatMultExpr::selectAssignKernel( ~lhs, A, B );
402 template<
typename MT3
405 static inline void selectAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
409 for(
size_t i=0UL; i<C.rows(); ++i ) {
410 for(
size_t j=0UL; j<C.columns(); ++j )
412 ConstIterator element( B.begin(j) );
413 const ConstIterator
end( B.end(j) );
415 if( element ==
end ) {
420 C(i,j) = A(i,element->index()) * element->value();
422 for( ; element!=
end; ++element )
423 C(i,j) += A(i,element->index()) * element->value();
443 template<
typename MT
445 friend inline typename DisableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
450 typedef typename SelectType< SO, OppositeType, ResultType >::Type TmpType;
462 const TmpType tmp(
serial( rhs ) );
483 template<
typename MT >
484 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
512 template<
typename MT >
513 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
541 template<
typename MT
543 friend inline typename DisableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
558 DMatTSMatMultExpr::selectAddAssignKernel( ~lhs, A, B );
577 template<
typename MT3
580 static inline void selectAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
584 for(
size_t i=0UL; i<C.rows(); ++i ) {
585 for(
size_t j=0UL; j<C.columns(); ++j )
587 ConstIterator element( B.begin(j) );
588 const ConstIterator
end( B.end(j) );
590 for( ; element!=
end; ++element )
591 C(i,j) += A(i,element->index()) * element->value();
613 template<
typename MT >
614 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
642 template<
typename MT >
643 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
675 template<
typename MT
691 DMatTSMatMultExpr::selectSubAssignKernel( ~lhs, A, B );
710 template<
typename MT3
713 static inline void selectSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
717 for(
size_t i=0UL; i<C.rows(); ++i ) {
718 for(
size_t j=0UL; j<C.columns(); ++j )
720 ConstIterator element( B.begin(j) );
721 const ConstIterator
end( B.end(j) );
723 for( ; element!=
end; ++element )
724 C(i,j) -= A(i,element->index()) * element->value();
746 template<
typename MT >
747 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
775 template<
typename MT >
776 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
818 template<
typename MT
820 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
855 template<
typename MT
857 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
862 typedef typename SelectType< SO, OppositeType, ResultType >::Type TmpType;
874 const TmpType tmp( rhs );
895 template<
typename MT >
896 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
924 template<
typename MT >
925 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
956 template<
typename MT
958 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
993 template<
typename MT >
994 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1022 template<
typename MT >
1023 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1058 template<
typename MT
1060 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1095 template<
typename MT >
1096 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1124 template<
typename MT >
1125 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1203 template<
typename T1
1205 inline const DMatTSMatMultExpr<T1,T2>
1211 throw std::invalid_argument(
"Matrix sizes do not match" );
1228 template<
typename MT1,
typename MT2 >
1246 template<
typename MT1,
typename MT2 >
1248 :
public Columns<MT2>
1264 template<
typename MT1,
typename MT2 >
1266 :
public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1282 template<
typename MT1,
typename MT2 >
1284 :
public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1300 template<
typename MT1,
typename MT2,
typename VT >
1305 typedef typename SelectType< IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1306 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value &&
1307 IsDenseVector<VT>::value && IsColumnVector<VT>::value
1308 ,
typename DMatDVecMultExprTrait< MT1, typename TSMatDVecMultExprTrait<MT2,VT>::Type >::Type
1309 , INVALID_TYPE >::Type Type;
1318 template<
typename MT1,
typename MT2,
typename VT >
1323 typedef typename SelectType< IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1324 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value &&
1325 IsSparseVector<VT>::value && IsColumnVector<VT>::value
1326 ,
typename DMatSVecMultExprTrait< MT1, typename TSMatSVecMultExprTrait<MT2,VT>::Type >::Type
1327 , INVALID_TYPE >::Type Type;
1336 template<
typename VT,
typename MT1,
typename MT2 >
1341 typedef typename SelectType< IsDenseVector<VT>::value && IsRowVector<VT>::value &&
1342 IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1343 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value
1344 ,
typename TDVecTSMatMultExprTrait< typename TDVecDMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1345 , INVALID_TYPE >::Type Type;
1354 template<
typename VT,
typename MT1,
typename MT2 >
1359 typedef typename SelectType< IsSparseVector<VT>::value && IsRowVector<VT>::value &&
1360 IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1361 IsSparseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value
1362 ,
typename TDVecTSMatMultExprTrait< typename TSVecDMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1363 , INVALID_TYPE >::Type Type;
1372 template<
typename MT1,
typename MT2,
bool AF >
1377 typedef typename MultExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1378 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1387 template<
typename MT1,
typename MT2 >
1392 typedef typename MultExprTrait< typename RowExprTrait<const MT1>::Type, MT2 >::Type Type;
1401 template<
typename MT1,
typename MT2 >
1406 typedef typename MultExprTrait< MT1, typename ColumnExprTrait<const MT2>::Type >::Type Type;
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTSMatMultExpr.h:173
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTSMatMultExpr.h:220
Header file for the Rows type trait.
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
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTSMatMultExpr.h:178
MT2::CompositeType CT2
Composite type of the right-hand side sparse matrix expression.
Definition: DMatTSMatMultExpr.h:123
const size_t SMP_DMATTSMATMULT_THRESHOLD
SMP row-major dense matrix/column-major sparse matrix multiplication threshold.This threshold specifi...
Definition: Thresholds.h:949
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
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTSMatMultExpr.h:310
Header file for the IsSparseMatrix type trait.
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:205
#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:172
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
Header file for the IsRowVector type trait.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
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.
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:288
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:107
Constraint on the data type.
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTSMatMultExpr.h:175
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
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatTSMatMultExpr.h:332
DMatTSMatMultExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTSMatMultExpr class.
Definition: DMatTSMatMultExpr.h:205
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 IsSymmetric type trait.
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTSMatMultExpr.h:268
#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:122
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
Header file for the TDVecTSMatMultExprTrait class template.
Header file for the DenseMatrix base 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
Header file for the Columns type trait.
Header file for the TSMatDVecMultExprTrait class template.
RT2::ElementType ET2
Element type of the right-hand side sparse matrix expression.
Definition: DMatTSMatMultExpr.h:121
SelectType< evaluateLeft, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side dense matrix operand.
Definition: DMatTSMatMultExpr.h:184
Header file for the DMatDVecMultExprTrait class template.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTSMatMultExpr.h:118
Header file for the IsLower type trait.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTSMatMultExpr.h:171
#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.
RT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: DMatTSMatMultExpr.h:120
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
RightOperand rhs_
Right-hand side sparse matrix of the multiplication expression.
Definition: DMatTSMatMultExpr.h:350
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.
Header file for the EnableIf class template.
Header file for the serial shim.
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_MATMATMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid matrix/matrix ...
Definition: MatMatMultExpr.h:165
Base class for all matrix/matrix multiplication expression templates.The MatMatMultExpr class serves ...
Definition: MatMatMultExpr.h:65
EnableIf< IsDenseMatrix< MT1 > >::Type smpSubAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:160
Header file for the IsSparseVector type trait.
#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.
#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.
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
Utility type for generic codes.
Base template for the MultTrait class.
Definition: MultTrait.h:142
const ElementType ReturnType
Return type for expression template evaluations.
Definition: DMatTSMatMultExpr.h:174
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTSMatMultExpr.h:170
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
RightOperand rightOperand() const
Returns the right-hand side transpose sparse matrix operand.
Definition: DMatTSMatMultExpr.h:298
Header file for the reset shim.
SelectType< evaluateRight, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side sparse matrix operand.
Definition: DMatTSMatMultExpr.h:187
Expression object for dense matrix-transpose sparse matrix multiplications.The DMatTSMatMultExpr clas...
Definition: DMatTSMatMultExpr.h:112
BLAZE_ALWAYS_INLINE void reset(const NonNumericProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: NonNumericProxy.h:833
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:322
Substitution Failure Is Not An Error (SFINAE) class.The DisableIf class template is an auxiliary tool...
Definition: DisableIf.h:184
#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.
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTSMatMultExpr.h:342
LeftOperand lhs_
Left-hand side dense matrix of the multiplication expression.
Definition: DMatTSMatMultExpr.h:349
Header file for the IsRowMajorMatrix type trait.
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:932
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.
EnableIf< IsDenseMatrix< MT1 > >::Type smpAddAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP addition assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:129
#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 the IsTrue value trait.
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:278
Header file for the IsUpper type trait.
MT2::ResultType RT2
Result type of the right-hand side sparse matrix expression.
Definition: DMatTSMatMultExpr.h:119
Header file for the DMatSVecMultExprTrait class template.
Header file for the IsColumnVector type trait.
Constraint on the data type.
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
DMatTSMatMultExpr< MT1, MT2 > This
Type of this DMatTSMatMultExpr instance.
Definition: DMatTSMatMultExpr.h:169
#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:181
Header file for the IsExpression type trait class.
Header file for the TSMatSVecMultExprTrait class template.
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