35 #ifndef _BLAZE_MATH_EXPRESSIONS_TSMATDMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_TSMATDMATMULTEXPR_H_
114 template<
typename MT1
116 class TSMatDMatMultExpr :
public DenseMatrix< TSMatDMatMultExpr<MT1,MT2>, true >
117 ,
private MatMatMultExpr
118 ,
private Computation
148 template<
typename T1,
typename T2,
typename T3 >
149 struct CanExploitSymmetry {
150 enum { value = IsRowMajorMatrix<T1>::value && IsSymmetric<T2>::value };
161 template<
typename T1,
typename T2,
typename T3 >
162 struct IsEvaluationRequired {
163 enum { value = ( evaluateLeft || evaluateRight ) &&
164 !CanExploitSymmetry<T1,T2,T3>::value };
175 template<
typename T1,
typename T2,
typename T3 >
176 struct UseOptimizedKernel {
177 enum { value = !IsResizable<typename T1::ElementType>::value };
187 template<
typename T1,
typename T2,
typename T3 >
188 struct UseDefaultKernel {
189 enum { value = !UseOptimizedKernel<T1,T2,T3>::value };
219 enum { vectorizable = 0 };
222 enum { smpAssignable = !evaluateLeft && MT1::smpAssignable &&
223 !evaluateRight && MT2::smpAssignable };
253 if(
lhs_.columns() != 0UL ) {
255 for(
size_t k=1UL; k<
lhs_.columns(); ++k ) {
283 return rhs_.columns();
313 template<
typename T >
315 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
325 template<
typename T >
327 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
337 return rhs_.isAligned();
370 template<
typename MT
390 TSMatDMatMultExpr::selectAssignKernel( ~lhs, A, B );
410 template<
typename MT3
414 selectAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
422 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
423 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
424 for(
size_t i=0UL; i<A.columns(); ++i ) {
425 const ConstIterator
end( A.end(i) );
426 ConstIterator element( A.begin(i) );
427 for( ; element!=
end; ++element ) {
428 for(
size_t j=jj; j<jend; ++j ) {
430 C(element->index(),j) = element->value() * B(i,j);
432 C(element->index(),j) += element->value() * B(i,j);
456 template<
typename MT3
459 static inline typename EnableIf< UseOptimizedKernel<MT3,MT4,MT5> >::Type
460 selectAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
464 const size_t block( IsRowMajorMatrix<MT3>::value ? 256UL : 8UL );
466 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
467 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
468 for(
size_t i=0UL; i<A.rows(); ++i ) {
469 for(
size_t j=jj; j<jend; ++j ) {
473 for(
size_t i=0UL; i<A.columns(); ++i )
475 const ConstIterator
end( A.end(i) );
476 ConstIterator element( A.begin(i) );
478 const size_t nonzeros( A.nonZeros(i) );
480 const size_t kend( nonzeros &
size_t(-4) );
483 for(
size_t k=0UL; k<kend; k+=4UL ) {
484 const size_t i1( element->index() );
485 const ET1 v1( element->value() );
487 const size_t i2( element->index() );
488 const ET1 v2( element->value() );
490 const size_t i3( element->index() );
491 const ET1 v3( element->value() );
493 const size_t i4( element->index() );
494 const ET1 v4( element->value() );
497 for(
size_t j=jj; j<jend; ++j ) {
498 C(i1,j) += v1 * B(i,j);
499 C(i2,j) += v2 * B(i,j);
500 C(i3,j) += v3 * B(i,j);
501 C(i4,j) += v4 * B(i,j);
505 for( ; element!=
end; ++element ) {
506 for(
size_t j=jj; j<jend; ++j ) {
507 C(element->index(),j) += element->value() * B(i,j);
529 template<
typename MT
531 friend inline typename DisableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
536 typedef typename SelectType< SO, ResultType, OppositeType >::Type TmpType;
548 const TmpType tmp(
serial( rhs ) );
569 template<
typename MT >
570 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
598 template<
typename MT
617 TSMatDMatMultExpr::selectAddAssignKernel( ~lhs, A, B );
636 template<
typename MT3
639 static inline void selectAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
643 const size_t block( IsRowMajorMatrix<MT3>::value ? 256UL : 8UL );
645 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
646 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
647 for(
size_t i=0UL; i<A.columns(); ++i )
649 const ConstIterator
end( A.end(i) );
650 ConstIterator element( A.begin(i) );
652 const size_t nonzeros( A.nonZeros(i) );
654 const size_t kend( nonzeros &
size_t(-4) );
657 for(
size_t k=0UL; k<kend; k+=4UL ) {
658 const size_t i1( element->index() );
659 const ET1 v1( element->value() );
661 const size_t i2( element->index() );
662 const ET1 v2( element->value() );
664 const size_t i3( element->index() );
665 const ET1 v3( element->value() );
667 const size_t i4( element->index() );
668 const ET1 v4( element->value() );
671 for(
size_t j=jj; j<jend; ++j ) {
672 C(i1,j) += v1 * B(i,j);
673 C(i2,j) += v2 * B(i,j);
674 C(i3,j) += v3 * B(i,j);
675 C(i4,j) += v4 * B(i,j);
679 for( ; element!=
end; ++element ) {
680 for(
size_t j=jj; j<jend; ++j ) {
681 C(element->index(),j) += element->value() * B(i,j);
705 template<
typename MT >
706 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
738 template<
typename MT
757 TSMatDMatMultExpr::selectSubAssignKernel( ~lhs, A, B );
776 template<
typename MT3
779 static inline void selectSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
783 const size_t block( IsRowMajorMatrix<MT3>::value ? 256UL : 8UL );
785 for(
size_t jj=0UL; jj<B.columns(); jj+=block ) {
786 const size_t jend( ( jj+block > B.columns() )?( B.columns() ):( jj+block ) );
787 for(
size_t i=0UL; i<A.columns(); ++i )
789 const ConstIterator
end( A.end(i) );
790 ConstIterator element( A.begin(i) );
792 const size_t nonzeros( A.nonZeros(i) );
794 const size_t kend( nonzeros &
size_t(-4) );
797 for(
size_t k=0UL; k<kend; k+=4UL ) {
798 const size_t i1( element->index() );
799 const ET1 v1( element->value() );
801 const size_t i2( element->index() );
802 const ET1 v2( element->value() );
804 const size_t i3( element->index() );
805 const ET1 v3( element->value() );
807 const size_t i4( element->index() );
808 const ET1 v4( element->value() );
811 for(
size_t j=jj; j<jend; ++j ) {
812 C(i1,j) -= v1 * B(i,j);
813 C(i2,j) -= v2 * B(i,j);
814 C(i3,j) -= v3 * B(i,j);
815 C(i4,j) -= v4 * B(i,j);
819 for( ; element!=
end; ++element ) {
820 for(
size_t j=jj; j<jend; ++j ) {
821 C(element->index(),j) -= element->value() * B(i,j);
845 template<
typename MT >
846 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
889 template<
typename MT
891 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
930 template<
typename MT
932 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
937 typedef typename SelectType< SO, ResultType, OppositeType >::Type TmpType;
949 const TmpType tmp( rhs );
970 template<
typename MT >
971 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1000 template<
typename MT
1002 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1040 template<
typename MT >
1041 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1074 template<
typename MT
1076 friend inline typename EnableIf< IsEvaluationRequired<MT,MT1,MT2> >::Type
1114 template<
typename MT >
1115 friend inline typename EnableIf< CanExploitSymmetry<MT,MT1,MT2> >::Type
1192 template<
typename T1
1194 inline const TSMatDMatMultExpr<T1,T2>
1200 throw std::invalid_argument(
"Matrix sizes do not match" );
1217 template<
typename MT1,
typename MT2 >
1235 template<
typename MT1,
typename MT2 >
1237 :
public Columns<MT2>
1253 template<
typename MT1,
typename MT2 >
1255 :
public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1271 template<
typename MT1,
typename MT2 >
1273 :
public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1289 template<
typename MT1,
typename MT2,
typename VT >
1294 typedef typename SelectType< IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1295 IsDenseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value &&
1296 IsDenseVector<VT>::value && IsColumnVector<VT>::value
1297 ,
typename TSMatDVecMultExprTrait< MT1, typename DMatDVecMultExprTrait<MT2,VT>::Type >::Type
1298 , INVALID_TYPE >::Type Type;
1307 template<
typename MT1,
typename MT2,
typename VT >
1312 typedef typename SelectType< IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1313 IsDenseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value &&
1314 IsSparseVector<VT>::value && IsColumnVector<VT>::value
1315 ,
typename TSMatDVecMultExprTrait< MT1, typename DMatSVecMultExprTrait<MT2,VT>::Type >::Type
1316 , INVALID_TYPE >::Type Type;
1325 template<
typename VT,
typename MT1,
typename MT2 >
1330 typedef typename SelectType< IsDenseVector<VT>::value && IsRowVector<VT>::value &&
1331 IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1332 IsDenseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value
1333 ,
typename TDVecDMatMultExprTrait< typename TDVecTSMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1334 , INVALID_TYPE >::Type Type;
1343 template<
typename VT,
typename MT1,
typename MT2 >
1348 typedef typename SelectType< IsSparseVector<VT>::value && IsRowVector<VT>::value &&
1349 IsSparseMatrix<MT1>::value && IsColumnMajorMatrix<MT1>::value &&
1350 IsDenseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value
1351 ,
typename TDVecDMatMultExprTrait< typename TDVecTSMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1352 , INVALID_TYPE >::Type Type;
1361 template<
typename MT1,
typename MT2,
bool AF >
1366 typedef typename MultExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1367 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1376 template<
typename MT1,
typename MT2 >
1381 typedef typename MultExprTrait< typename RowExprTrait<const MT1>::Type, MT2 >::Type Type;
1390 template<
typename MT1,
typename MT2 >
1395 typedef typename MultExprTrait< MT1, typename ColumnExprTrait<const MT2>::Type >::Type Type;
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: TSMatDMatMultExpr.h:198
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
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: TSMatDMatMultExpr.h:197
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.
size_t columns() const
Returns the current number of columns of the matrix.
Definition: TSMatDMatMultExpr.h:282
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: TSMatDMatMultExpr.h:314
TSMatDMatMultExpr< MT1, MT2 > This
Type of this TSMatDMatMultExpr instance.
Definition: TSMatDMatMultExpr.h:196
Header file for the IsColumnMajorMatrix type trait.
Header file for the TSVecTSMatMultExprTrait class template.
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: TSMatDMatMultExpr.h:336
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
LeftOperand leftOperand() const
Returns the left-hand side transpose sparse matrix operand.
Definition: TSMatDMatMultExpr.h:292
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: TSMatDMatMultExpr.h:199
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.
MT1::CompositeType CT1
Composite type of the left-hand side sparse matrix expression.
Definition: TSMatDMatMultExpr.h:126
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
LeftOperand lhs_
Left-hand side sparse matrix of the multiplication expression.
Definition: TSMatDMatMultExpr.h:353
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: TSMatDMatMultExpr.h:326
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
SelectType< evaluateLeft, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side sparse matrix operand.
Definition: TSMatDMatMultExpr.h:211
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.
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: TSMatDMatMultExpr.h:346
RT2::ElementType ET2
Element type of the right-hand side sparse matrix expression.
Definition: TSMatDMatMultExpr.h:125
Header file for the multiplication trait.
RightOperand rhs_
Right-hand side dense matrix of the multiplication expression.
Definition: TSMatDMatMultExpr.h:354
Header file for the IsSymmetric type trait.
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: TSMatDMatMultExpr.h:127
Compile time check for row-major matrix types.This type trait tests whether or not the given template...
Definition: IsRowMajorMatrix.h:104
#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
RightOperand rightOperand() const
Returns the right-hand side dense matrix operand.
Definition: TSMatDMatMultExpr.h:302
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.
MT1::ResultType RT1
Result type of the left-hand side sparse matrix expression.
Definition: TSMatDMatMultExpr.h:122
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.
TSMatDMatMultExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the TSMatDMatMultExpr class.
Definition: TSMatDMatMultExpr.h:232
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: TSMatDMatMultExpr.h:247
SelectType< evaluateRight, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side dense matrix operand.
Definition: TSMatDMatMultExpr.h:214
const size_t SMP_TSMATDMATMULT_THRESHOLD
SMP column-major sparse matrix/row-major dense matrix multiplication threshold.This threshold specifi...
Definition: Thresholds.h:1064
Header file for the DMatDVecMultExprTrait class template.
Header file for the IsLower type trait.
#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
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: TSMatDMatMultExpr.h:123
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: TSMatDMatMultExpr.h:208
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.
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side sparse matrix expression.
Definition: TSMatDMatMultExpr.h:205
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
RT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: TSMatDMatMultExpr.h:124
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
const ResultType CompositeType
Data type for composite expression templates.
Definition: TSMatDMatMultExpr.h:202
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
const ElementType ReturnType
Return type for expression template evaluations.
Definition: TSMatDMatMultExpr.h:201
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
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
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.
size_t rows() const
Returns the current number of rows of the matrix.
Definition: TSMatDMatMultExpr.h:272
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.
Expression object for transpose sparse matrix-dense matrix multiplications.The TSMatDMatMultExpr clas...
Definition: Forward.h:137
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
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.
ResultType::ElementType ElementType
Resulting element type.
Definition: TSMatDMatMultExpr.h:200
Header file for basic type definitions.
Header file for the TSVecDMatMultExprTrait class template.
Header file for the IsUpper type trait.
Header file for the DMatSVecMultExprTrait class template.
Header file for the IsColumnVector type trait.
Constraint on the data type.
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