35 #ifndef _BLAZE_MATH_EXPRESSIONS_TSMATTDMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_TSMATTDMATMULTEXPR_H_
109 template<
typename MT1
111 class TSMatTDMatMultExpr :
public DenseMatrix< TSMatTDMatMultExpr<MT1,MT2>, true >
112 ,
private MatMatMultExpr
113 ,
private Computation
142 template<
typename T1,
typename T2,
typename T3 >
143 struct CanExploitSymmetry {
144 enum { value = ( IsSymmetric<T2>::value || IsSymmetric<T3>::value ) };
155 template<
typename T1,
typename T2,
typename T3 >
156 struct IsEvaluationRequired {
157 enum { value = ( evaluateLeft || evaluateRight ) &&
158 !CanExploitSymmetry<T1,T2,T3>::value };
169 template<
typename T1,
typename T2,
typename T3 >
170 struct UseOptimizedKernel {
171 enum { value = !IsResizable<typename T1::ElementType>::value };
181 template<
typename T1,
typename T2,
typename T3 >
182 struct UseDefaultKernel {
183 enum { value = !UseOptimizedKernel<T1,T2,T3>::value };
213 enum { vectorizable = 0 };
216 enum { smpAssignable = !evaluateLeft && MT1::smpAssignable &&
217 !evaluateRight && MT2::smpAssignable };
247 if(
lhs_.columns() != 0UL ) {
249 for(
size_t k=1UL; k<
lhs_.columns(); ++k ) {
277 return rhs_.columns();
307 template<
typename T >
309 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
319 template<
typename T >
321 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
331 return rhs_.isAligned();
364 template<
typename MT
384 TSMatTDMatMultExpr::selectAssignKernel( ~lhs, A, B );
404 template<
typename MT3
408 selectAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
414 const size_t block( 64UL );
416 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
417 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
418 for(
size_t i=0UL; i<A.columns(); ++i ) {
419 const ConstIterator
end( A.end(i) );
420 ConstIterator element( A.begin(i) );
421 for( ; element!=
end; ++element ) {
422 for(
size_t j=jj; j<jend; ++j ) {
424 C(element->index(),j) = element->value() * B(i,j);
426 C(element->index(),j) += element->value() * B(i,j);
450 template<
typename MT3
453 static inline typename EnableIf< UseOptimizedKernel<MT3,MT4,MT5> >::Type
454 selectAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
458 const size_t block( IsRowMajorMatrix<MT3>::value ? 128UL : 64UL );
460 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
461 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
462 for(
size_t i=0UL; i<A.rows(); ++i ) {
463 for(
size_t j=jj; j<jend; ++j ) {
467 for(
size_t i=0UL; i<A.columns(); ++i )
469 const ConstIterator
end( A.end(i) );
470 ConstIterator element( A.begin(i) );
472 const size_t nonzeros( A.nonZeros(i) );
474 const size_t kend( nonzeros &
size_t(-4) );
477 for(
size_t k=0UL; k<kend; k+=4UL ) {
478 const size_t i1( element->index() );
479 const ET1 v1( element->value() );
481 const size_t i2( element->index() );
482 const ET1 v2( element->value() );
484 const size_t i3( element->index() );
485 const ET1 v3( element->value() );
487 const size_t i4( element->index() );
488 const ET1 v4( element->value() );
491 for(
size_t j=jj; j<jend; ++j ) {
492 C(i1,j) += v1 * B(i,j);
493 C(i2,j) += v2 * B(i,j);
494 C(i3,j) += v3 * B(i,j);
495 C(i4,j) += v4 * B(i,j);
499 for( ; element!=
end; ++element ) {
500 for(
size_t j=jj; j<jend; ++j ) {
501 C(element->index(),j) += element->value() * B(i,j);
523 template<
typename MT
525 friend inline typename DisableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
530 typedef typename SelectType< SO, ResultType, OppositeType >::Type TmpType;
542 const TmpType tmp(
serial( rhs ) );
563 template<
typename MT >
564 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
574 if( IsSymmetric<MT1>::value && IsSymmetric<MT2>::value )
576 else if( IsSymmetric<MT1>::value )
599 template<
typename MT >
600 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
608 if( IsSymmetric<MT1>::value )
629 template<
typename MT
631 friend inline typename DisableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
649 TSMatTDMatMultExpr::selectAddAssignKernel( ~lhs, A, B );
668 template<
typename MT3
671 static inline void selectAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
675 const size_t block( IsRowMajorMatrix<MT3>::value ? 128UL : 64UL );
677 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
678 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
679 for(
size_t i=0UL; i<A.columns(); ++i )
681 const ConstIterator
end( A.end(i) );
682 ConstIterator element( A.begin(i) );
684 const size_t nonzeros( A.nonZeros(i) );
686 const size_t kend( nonzeros &
size_t(-4) );
689 for(
size_t k=0UL; k<kend; k+=4UL ) {
690 const size_t i1( element->index() );
691 const ET1 v1( element->value() );
693 const size_t i2( element->index() );
694 const ET1 v2( element->value() );
696 const size_t i3( element->index() );
697 const ET1 v3( element->value() );
699 const size_t i4( element->index() );
700 const ET1 v4( element->value() );
703 for(
size_t j=jj; j<jend; ++j ) {
704 C(i1,j) += v1 * B(i,j);
705 C(i2,j) += v2 * B(i,j);
706 C(i3,j) += v3 * B(i,j);
707 C(i4,j) += v4 * B(i,j);
711 for( ; element!=
end; ++element ) {
712 for(
size_t j=jj; j<jend; ++j ) {
713 C(element->index(),j) += element->value() * B(i,j);
737 template<
typename MT >
738 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
748 if( IsSymmetric<MT1>::value && IsSymmetric<MT2>::value )
750 else if( IsSymmetric<MT1>::value )
773 template<
typename MT >
774 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
782 if( IsSymmetric<MT1>::value )
807 template<
typename MT
809 friend inline typename DisableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
827 TSMatTDMatMultExpr::selectSubAssignKernel( ~lhs, A, B );
846 template<
typename MT3
849 static inline void selectSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
853 const size_t block( IsRowMajorMatrix<MT3>::value ? 128UL : 64UL );
855 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
856 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
857 for(
size_t i=0UL; i<A.columns(); ++i )
859 const ConstIterator
end( A.end(i) );
860 ConstIterator element( A.begin(i) );
862 const size_t nonzeros( A.nonZeros(i) );
864 const size_t kend( nonzeros &
size_t(-4) );
867 for(
size_t k=0UL; k<kend; k+=4UL ) {
868 const size_t i1( element->index() );
869 const ET1 v1( element->value() );
871 const size_t i2( element->index() );
872 const ET1 v2( element->value() );
874 const size_t i3( element->index() );
875 const ET1 v3( element->value() );
877 const size_t i4( element->index() );
878 const ET1 v4( element->value() );
881 for(
size_t j=jj; j<jend; ++j ) {
882 C(i1,j) -= v1 * B(i,j);
883 C(i2,j) -= v2 * B(i,j);
884 C(i3,j) -= v3 * B(i,j);
885 C(i4,j) -= v4 * B(i,j);
889 for( ; element!=
end; ++element ) {
890 for(
size_t j=jj; j<jend; ++j ) {
891 C(element->index(),j) -= element->value() * B(i,j);
916 template<
typename MT >
917 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
927 if( IsSymmetric<MT1>::value && IsSymmetric<MT2>::value )
929 else if( IsSymmetric<MT1>::value )
953 template<
typename MT >
954 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
962 if( IsSymmetric<MT1>::value )
998 template<
typename MT
1000 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1039 template<
typename MT
1041 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1046 typedef typename SelectType< SO, ResultType, OppositeType >::Type TmpType;
1058 const TmpType tmp( rhs );
1079 template<
typename MT >
1080 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1090 if( IsSymmetric<MT1>::value && IsSymmetric<MT2>::value )
1092 else if( IsSymmetric<MT1>::value )
1115 template<
typename MT >
1116 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1126 if( IsSymmetric<MT1>::value )
1150 template<
typename MT
1152 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1191 template<
typename MT >
1192 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1202 if( IsSymmetric<MT1>::value && IsSymmetric<MT2>::value )
1204 else if( IsSymmetric<MT1>::value )
1228 template<
typename MT >
1229 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1237 if( IsSymmetric<MT1>::value )
1265 template<
typename MT
1267 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1306 template<
typename MT >
1307 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1317 if( IsSymmetric<MT1>::value && IsSymmetric<MT2>::value )
1319 else if( IsSymmetric<MT1>::value )
1343 template<
typename MT >
1344 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1352 if( IsSymmetric<MT1>::value )
1422 template<
typename T1
1424 inline const TSMatTDMatMultExpr<T1,T2>
1430 throw std::invalid_argument(
"Matrix sizes do not match" );
1447 template<
typename MT1,
typename MT2 >
1465 template<
typename MT1,
typename MT2 >
1467 :
public Columns<MT2>
1483 template<
typename MT1,
typename MT2 >
1485 :
public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1501 template<
typename MT1,
typename MT2 >
1503 :
public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1519 template<
typename MT1,
typename MT2,
typename VT >
1524 typedef typename SelectType< IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1525 IsDenseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value &&
1526 IsDenseVector<VT>::value && IsColumnVector<VT>::value
1527 ,
typename TSMatDVecMultExprTrait< MT1, typename TDMatDVecMultExprTrait<MT2,VT>::Type >::Type
1528 , INVALID_TYPE >::Type Type;
1537 template<
typename MT1,
typename MT2,
typename VT >
1542 typedef typename SelectType< IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1543 IsDenseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value &&
1544 IsSparseVector<VT>::value && IsColumnVector<VT>::value
1545 ,
typename TSMatDVecMultExprTrait< MT1, typename TDMatSVecMultExprTrait<MT2,VT>::Type >::Type
1546 , INVALID_TYPE >::Type Type;
1555 template<
typename VT,
typename MT1,
typename MT2 >
1560 typedef typename SelectType< IsDenseVector<VT>::value && IsRowVector<VT>::value &&
1561 IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1562 IsDenseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value
1563 ,
typename TDVecTDMatMultExprTrait< typename TDVecTSMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1564 , INVALID_TYPE >::Type Type;
1573 template<
typename VT,
typename MT1,
typename MT2 >
1578 typedef typename SelectType< IsSparseVector<VT>::value && IsRowVector<VT>::value &&
1579 IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1580 IsDenseMatrix<MT2>::value && IsColumnMajorMatrix<MT2>::value
1581 ,
typename TSVecTDMatMultExprTrait< typename TSVecTSMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1582 , INVALID_TYPE >::Type Type;
1591 template<
typename MT1,
typename MT2,
bool AF >
1596 typedef typename MultExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1597 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1606 template<
typename MT1,
typename MT2 >
1611 typedef typename MultExprTrait< typename RowExprTrait<const MT1>::Type, MT2 >::Type Type;
1620 template<
typename MT1,
typename MT2 >
1625 typedef typename MultExprTrait< MT1, typename ColumnExprTrait<const MT2>::Type >::Type Type;
size_t columns() const
Returns the current number of columns of the matrix.
Definition: TSMatTDMatMultExpr.h:276
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
size_t rows() const
Returns the current number of rows of the matrix.
Definition: TSMatTDMatMultExpr.h:266
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
Expression object for transpose sparse matrix-transpose dense matrix multiplications.The TSMatTDMatMultExpr class represents the compile time expression for multiplications between a column-major sparse matrix and a column-major dense matrix.
Definition: Forward.h:143
TSMatTDMatMultExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the TSMatTDMatMultExpr class.
Definition: TSMatTDMatMultExpr.h:226
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
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.
LeftOperand leftOperand() const
Returns the left-hand side transpose sparse matrix operand.
Definition: TSMatTDMatMultExpr.h:286
Header file for the IsColumnMajorMatrix type trait.
Header file for the TSVecTSMatMultExprTrait class template.
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
TSMatTDMatMultExpr< MT1, MT2 > This
Type of this TSMatTDMatMultExpr instance.
Definition: TSMatTDMatMultExpr.h:190
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.
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: TSMatTDMatMultExpr.h:191
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
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: TSMatTDMatMultExpr.h:192
SelectType< evaluateRight, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side dense matrix operand.
Definition: TSMatTDMatMultExpr.h:208
const ElementType ReturnType
Return type for expression template evaluations.
Definition: TSMatTDMatMultExpr.h:195
MT1::ResultType RT1
Result type of the left-hand side sparse matrix expression.
Definition: TSMatTDMatMultExpr.h:117
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
RT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: TSMatTDMatMultExpr.h:119
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.
const size_t SMP_TSMATTDMATMULT_THRESHOLD
SMP column-major sparse matrix/column-major dense matrix multiplication threshold.This threshold specifies when a column-major sparse matrix/column-major dense matrix multiplication can be executed in parallel. In case the number of rows/columns of the target matrix is larger or equal to this threshold, the operation is executed in parallel. If the number of rows/columns is below this threshold the operation is executed single-threaded.
Definition: Thresholds.h:1087
Header file for the multiplication trait.
Header file for the IsSymmetric type trait.
#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
Header file for the TSVecTDMatMultExprTrait class template.
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
Header file for the TDMatSVecMultExprTrait class template.
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.
Header file for the IsLower type trait.
SelectType< evaluateLeft, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side sparse matrix operand.
Definition: TSMatTDMatMultExpr.h:205
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: TSMatTDMatMultExpr.h:308
#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
ResultType::ElementType ElementType
Resulting element type.
Definition: TSMatTDMatMultExpr.h:194
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: TSMatTDMatMultExpr.h:118
RightOperand rhs_
Right-hand side dense matrix of the multiplication expression.
Definition: TSMatTDMatMultExpr.h:348
Constraints on the storage order of matrix types.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
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.
MT1::CompositeType CT1
Composite type of the left-hand side sparse matrix expression.
Definition: TSMatTDMatMultExpr.h:121
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side sparse matrix expression.
Definition: TSMatTDMatMultExpr.h:199
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
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: TSMatTDMatMultExpr.h:296
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
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
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: TSMatTDMatMultExpr.h:122
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: TSMatTDMatMultExpr.h:241
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
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.
Header file for the isDefault shim.
BLAZE_ALWAYS_INLINE bool isDefault(const NonNumericProxy< MT > &proxy)
Returns whether the represented element is in default state.
Definition: NonNumericProxy.h:874
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: TSMatTDMatMultExpr.h:340
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.
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 isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: TSMatTDMatMultExpr.h:320
RT2::ElementType ET2
Element type of the right-hand side sparse matrix expression.
Definition: TSMatTDMatMultExpr.h:120
const ResultType CompositeType
Data type for composite expression templates.
Definition: TSMatTDMatMultExpr.h:196
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.
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
Header file for the TDMatDVecMultExprTrait 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:2473
Header file for the IsTrue value trait.
Header file for basic type definitions.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: TSMatTDMatMultExpr.h:193
Header file for the IsUpper type trait.
Header file for the IsColumnVector type trait.
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: TSMatTDMatMultExpr.h:202
Constraint on the data type.
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: TSMatTDMatMultExpr.h:330
Header file for the IsResizable 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
Header file for the TDVecTDMatMultExprTrait class template.
#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
Header file for the IsExpression type trait class.
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
LeftOperand lhs_
Left-hand side sparse matrix of the multiplication expression.
Definition: TSMatTDMatMultExpr.h:347