35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATSMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATSMATMULTEXPR_H_
103 template<
typename MT1
134 template<
typename T1,
typename T2,
typename T3 >
135 struct UseSMPAssignKernel {
136 enum { value = evaluateLeft || evaluateRight };
147 template<
typename T1,
typename T2,
typename T3 >
148 struct UseOptimizedKernel {
149 enum { value = !UseSMPAssignKernel<T1,T2,T3>::value &&
150 !IsResizable<typename T1::ElementType>::value };
160 template<
typename T1,
typename T2,
typename T3 >
161 struct UseDefaultKernel {
162 enum { value = !UseSMPAssignKernel<T1,T2,T3>::value &&
163 !UseOptimizedKernel<T1,T2,T3>::value };
193 enum { vectorizable = 0 };
196 enum { smpAssignable = !evaluateLeft && !evaluateRight };
226 if(
lhs_.columns() != 0UL ) {
228 for(
size_t k=1UL; k<
lhs_.columns(); ++k ) {
256 return rhs_.columns();
286 template<
typename T >
288 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
298 template<
typename T >
300 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
310 return lhs_.isAligned();
343 template<
typename MT >
361 DMatSMatMultExpr::selectRowMajorAssignKernel( ~lhs, A, B );
380 template<
typename MT3
384 selectRowMajorAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
388 for(
size_t i=0UL; i<A.rows(); ++i ) {
389 for(
size_t j=0UL; j<C.columns(); ++j ) {
392 for(
size_t j=0UL; j<B.rows(); ++j ) {
393 ConstIterator element( B.begin(j) );
394 const ConstIterator end( B.end(j) );
395 for( ; element!=end; ++element ) {
397 C(i,element->index()) = A(i,j) * element->value();
399 C(i,element->index()) += A(i,j) * element->value();
421 template<
typename MT3
424 static inline typename EnableIf< UseOptimizedKernel<MT3,MT4,MT5> >::Type
425 selectRowMajorAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
429 BLAZE_INTERNAL_ASSERT( ( A.rows() - ( A.rows() % 4UL ) ) == ( A.rows() & size_t(-4) ),
"Invalid end calculation" );
430 const size_t last( A.rows() & size_t(-4) );
432 for(
size_t i=0UL; i<last; i+=4UL ) {
433 for(
size_t j=0UL; j<C.columns(); ++j ) {
439 for(
size_t j=0UL; j<B.rows(); ++j ) {
440 ConstIterator element( B.begin(j) );
441 const ConstIterator end( B.end(j) );
442 for( ; element!=end; ++element ) {
443 C(i ,element->index()) += A(i ,j) * element->value();
444 C(i+1UL,element->index()) += A(i+1UL,j) * element->value();
445 C(i+2UL,element->index()) += A(i+2UL,j) * element->value();
446 C(i+3UL,element->index()) += A(i+3UL,j) * element->value();
451 for(
size_t i=last; i<A.rows(); ++i ) {
452 for(
size_t j=0UL; j<C.columns(); ++j ) {
455 for(
size_t j=0UL; j<B.rows(); ++j ) {
456 ConstIterator element( B.begin(j) );
457 const ConstIterator end( B.end(j) );
458 for( ; element!=end; ++element ) {
459 C(i,element->index()) += A(i,j) * element->value();
480 template<
typename MT3
483 static inline typename EnableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
484 selectRowMajorAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
504 template<
typename MT >
522 DMatSMatMultExpr::selectColumnMajorAssignKernel( ~lhs, A, B );
541 template<
typename MT3
544 static inline typename EnableIf< UseDefaultKernel<MT3,MT4,MT5> >::Type
545 selectColumnMajorAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
551 for(
size_t i=0UL; i<B.rows(); ++i ) {
552 for(
size_t jj=0UL; jj<A.rows(); jj+=8UL ) {
553 const size_t jend( ( jj+8UL > A.rows() )?( A.rows() ):( jj+8UL ) );
554 ConstIterator element( B.begin(i) );
555 const ConstIterator end( B.end(i) );
556 for( ; element!=end; ++element ) {
557 for(
size_t j=jj; j<jend; ++j ) {
559 C(j,element->index()) = A(j,i) * element->value();
561 C(j,element->index()) += A(j,i) * element->value();
584 template<
typename MT3
587 static inline typename EnableIf< UseOptimizedKernel<MT3,MT4,MT5> >::Type
588 selectColumnMajorAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
594 for(
size_t i=0UL; i<B.rows(); ++i ) {
595 for(
size_t jj=0UL; jj<A.rows(); jj+=8UL ) {
596 const size_t jend( ( jj+8UL > A.rows() )?( A.rows() ):( jj+8UL ) );
597 ConstIterator element( B.begin(i) );
598 const ConstIterator end( B.end(i) );
599 for( ; element!=end; ++element ) {
600 for(
size_t j=jj; j<jend; ++j ) {
601 C(j,element->index()) += A(j,i) * element->value();
623 template<
typename MT3
626 static inline typename EnableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
627 selectColumnMajorAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
647 template<
typename MT
653 typedef typename SelectType< SO, OppositeType, ResultType >::Type TmpType;
665 const TmpType tmp( rhs );
684 template<
typename MT >
702 DMatSMatMultExpr::selectRowMajorAddAssignKernel( ~lhs, A, B );
721 template<
typename MT3
724 static inline typename DisableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
725 selectRowMajorAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
729 const size_t last( A.rows() & size_t(-4) );
731 for(
size_t i=0UL; i<last; i+=4UL ) {
732 for(
size_t j=0UL; j<B.rows(); ++j ) {
733 ConstIterator element( B.begin(j) );
734 const ConstIterator end( B.end(j) );
735 for( ; element!=end; ++element ) {
736 C(i ,element->index()) += A(i ,j) * element->value();
737 C(i+1UL,element->index()) += A(i+1UL,j) * element->value();
738 C(i+2UL,element->index()) += A(i+2UL,j) * element->value();
739 C(i+3UL,element->index()) += A(i+3UL,j) * element->value();
744 for(
size_t i=last; i<A.rows(); ++i ) {
745 for(
size_t j=0UL; j<B.rows(); ++j ) {
746 ConstIterator element( B.begin(j) );
747 const ConstIterator end( B.end(j) );
748 for( ; element!=end; ++element ) {
749 C(i,element->index()) += A(i,j) * element->value();
771 template<
typename MT3
774 static inline typename EnableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
775 selectRowMajorAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
795 template<
typename MT >
813 DMatSMatMultExpr::selectColumnMajorAddAssignKernel( ~lhs, A, B );
832 template<
typename MT3
835 static inline typename DisableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
836 selectColumnMajorAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
840 for(
size_t i=0UL; i<B.rows(); ++i ) {
841 for(
size_t jj=0UL; jj<A.rows(); jj+=8UL ) {
842 const size_t jend( ( jj+8UL > A.rows() )?( A.rows() ):( jj+8UL ) );
843 ConstIterator element( B.begin(i) );
844 const ConstIterator end( B.end(i) );
845 for( ; element!=end; ++element ) {
846 for(
size_t j=jj; j<jend; ++j ) {
847 C(j,element->index()) += A(j,i) * element->value();
870 template<
typename MT3
873 static inline typename EnableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
874 selectColumnMajorAddAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
898 template<
typename MT >
916 DMatSMatMultExpr::selectRowMajorSubAssignKernel( ~lhs, A, B );
935 template<
typename MT3
938 static inline typename DisableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
939 selectRowMajorSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
943 const size_t last( A.rows() & size_t(-4) );
945 for(
size_t i=0UL; i<last; i+=4UL ) {
946 for(
size_t j=0UL; j<B.rows(); ++j ) {
947 ConstIterator element( B.begin(j) );
948 const ConstIterator end( B.end(j) );
949 for( ; element!=end; ++element ) {
950 C(i ,element->index()) -= A(i ,j) * element->value();
951 C(i+1UL,element->index()) -= A(i+1UL,j) * element->value();
952 C(i+2UL,element->index()) -= A(i+2UL,j) * element->value();
953 C(i+3UL,element->index()) -= A(i+3UL,j) * element->value();
958 for(
size_t i=last; i<A.rows(); ++i ) {
959 for(
size_t j=0UL; j<B.rows(); ++j ) {
960 ConstIterator element( B.begin(j) );
961 const ConstIterator end( B.end(j) );
962 for( ; element!=end; ++element ) {
963 C(i,element->index()) -= A(i,j) * element->value();
985 template<
typename MT3
988 static inline typename EnableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
989 selectRowMajorSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
1009 template<
typename MT >
1027 DMatSMatMultExpr::selectColumnMajorSubAssignKernel( ~lhs, A, B );
1046 template<
typename MT3
1049 static inline typename DisableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
1050 selectColumnMajorSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
1054 for(
size_t i=0UL; i<B.rows(); ++i ) {
1055 for(
size_t jj=0UL; jj<A.rows(); jj+=8UL ) {
1056 const size_t jend( ( jj+8UL > A.rows() )?( A.rows() ):( jj+8UL ) );
1057 ConstIterator element( B.begin(i) );
1058 const ConstIterator end( B.end(i) );
1059 for( ; element!=end; ++element ) {
1060 for(
size_t j=jj; j<jend; ++j ) {
1061 C(j,element->index()) -= A(j,i) * element->value();
1084 template<
typename MT3
1087 static inline typename EnableIf< UseSMPAssignKernel<MT3,MT4,MT5> >::Type
1088 selectColumnMajorSubAssignKernel( MT3& C,
const MT4& A,
const MT5& B )
1156 template<
typename T1
1158 inline const DMatSMatMultExpr<T1,T2>
1164 throw std::invalid_argument(
"Matrix sizes do not match" );
1181 template<
typename MT1,
typename MT2,
typename VT >
1182 struct DMatDVecMultExprTrait< DMatSMatMultExpr<MT1,MT2>, VT >
1186 typedef typename SelectType< IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1187 IsSparseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value &&
1188 IsDenseVector<VT>::value && IsColumnVector<VT>::value
1189 ,
typename DMatDVecMultExprTrait< MT1, typename SMatDVecMultExprTrait<MT2,VT>::Type >::Type
1190 , INVALID_TYPE >::Type Type;
1199 template<
typename MT1,
typename MT2,
typename VT >
1200 struct DMatSVecMultExprTrait< DMatSMatMultExpr<MT1,MT2>, VT >
1204 typedef typename SelectType< IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1205 IsSparseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value &&
1206 IsSparseVector<VT>::value && IsColumnVector<VT>::value
1207 ,
typename DMatSVecMultExprTrait< MT1, typename SMatSVecMultExprTrait<MT2,VT>::Type >::Type
1208 , INVALID_TYPE >::Type Type;
1217 template<
typename VT,
typename MT1,
typename MT2 >
1218 struct TDVecDMatMultExprTrait< VT, DMatSMatMultExpr<MT1,MT2> >
1222 typedef typename SelectType< IsDenseVector<VT>::value && IsRowVector<VT>::value &&
1223 IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1224 IsSparseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value
1225 ,
typename TDVecSMatMultExprTrait< typename TDVecDMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1226 , INVALID_TYPE >::Type Type;
1235 template<
typename VT,
typename MT1,
typename MT2 >
1236 struct TSVecDMatMultExprTrait< VT, DMatSMatMultExpr<MT1,MT2> >
1240 typedef typename SelectType< IsSparseVector<VT>::value && IsRowVector<VT>::value &&
1241 IsDenseMatrix<MT1>::value && IsRowMajorMatrix<MT1>::value &&
1242 IsSparseMatrix<MT2>::value && IsRowMajorMatrix<MT2>::value
1243 ,
typename TDVecSMatMultExprTrait< typename TDVecDMatMultExprTrait<VT,MT1>::Type, MT2 >::Type
1244 , INVALID_TYPE >::Type Type;
1253 template<
typename MT1,
typename MT2,
bool AF >
1254 struct SubmatrixExprTrait< DMatSMatMultExpr<MT1,MT2>, AF >
1258 typedef typename MultExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1259 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1268 template<
typename MT1,
typename MT2 >
1269 struct RowExprTrait< DMatSMatMultExpr<MT1,MT2> >
1273 typedef typename MultExprTrait< typename RowExprTrait<const MT1>::Type, MT2 >::Type Type;
1282 template<
typename MT1,
typename MT2 >
1283 struct ColumnExprTrait< DMatSMatMultExpr<MT1,MT2> >
1287 typedef typename MultExprTrait< MT1, typename ColumnExprTrait<const MT2>::Type >::Type Type;
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatSMatMultExpr.h:173
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
const size_t SMP_DMATSMATMULT_THRESHOLD
SMP row-major dense matrix/row-major sparse matrix multiplication threshold.This threshold represents...
Definition: Thresholds.h:485
RightOperand rhs_
Right-hand side sparse matrix of the multiplication expression.
Definition: DMatSMatMultExpr.h:327
Header file for the SMatDVecMultExprTrait class template.
void reset(DynamicMatrix< Type, SO > &m)
Resetting the given dense matrix.
Definition: DynamicMatrix.h:4579
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:4075
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatSMatMultExpr.h:176
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatSMatMultExpr.h:220
const ElementType ReturnType
Return type for expression template evaluations.
Definition: DMatSMatMultExpr.h:175
SelectType< evaluateRight, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side sparse matrix operand.
Definition: DMatSMatMultExpr.h:188
void smpSubAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:151
Header file for the IsSparseMatrix type trait.
bool isDefault(const DynamicMatrix< Type, SO > &m)
Returns whether the given dense matrix is in default state.
Definition: DynamicMatrix.h:4622
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:197
#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 sparse matrix SMP implementation.
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2384
Header file for the IsRowVector type trait.
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:249
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatSMatMultExpr.h:255
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatSMatMultExpr.h:309
MT2::ResultType RT2
Result type of the right-hand side sparse matrix expression.
Definition: DMatSMatMultExpr.h:112
Header file for the TDVecSMatMultExprTrait class template.
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatSMatMultExpr.h:171
Header file for the Computation base class.
Header file for the MatMatMultExpr base class.
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.
Constraint on the data type.
Constraint on the data type.
Header file for the MultExprTrait class template.
void smpAddAssign(DenseMatrix< 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:121
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:90
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side sparse matrix expression.
Definition: DMatSMatMultExpr.h:182
LeftOperand lhs_
Left-hand side dense matrix of the multiplication expression.
Definition: DMatSMatMultExpr.h:326
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatSMatMultExpr.h:245
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.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatSMatMultExpr.h:111
#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
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2388
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatSMatMultExpr.h:265
Header file for the dense matrix SMP implementation.
Header file for the DenseMatrix base class.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatSMatMultExpr.h:172
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.
RT2::ElementType ET2
Element type of the right-hand side sparse matrix expression.
Definition: DMatSMatMultExpr.h:114
Header file for the SMatSVecMultExprTrait class template.
#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
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatSMatMultExpr.h:319
Constraints on the storage order of matrix types.
SelectType< evaluateLeft, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side dense matrix operand.
Definition: DMatSMatMultExpr.h:185
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2382
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.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatSMatMultExpr.h:174
void smpAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:91
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatSMatMultExpr.h:115
RT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: DMatSMatMultExpr.h:113
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
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatSMatMultExpr.h:287
Header file for run time assertion macros.
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatSMatMultExpr.h:299
Utility type for generic codes.
Base template for the MultTrait class.
Definition: MultTrait.h:141
Expression object for dense matrix-sparse matrix multiplications.The DMatSMatMultExpr class represent...
Definition: DMatSMatMultExpr.h:105
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
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.
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
Header file for the isDefault shim.
#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
MT2::CompositeType CT2
Composite type of the right-hand side sparse matrix expression.
Definition: DMatSMatMultExpr.h:116
Header file for the IsDenseVector type trait.
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatSMatMultExpr.h:179
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:2379
size_t columns(const Matrix< MT, SO > &m)
Returns the current number of columns of the matrix.
Definition: Matrix.h:154
DMatSMatMultExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatSMatMultExpr class.
Definition: DMatSMatMultExpr.h:205
Header file for basic type definitions.
Header file for the TSVecDMatMultExprTrait class template.
DMatSMatMultExpr< MT1, MT2 > This
Type of this DMatSMatMultExpr instance.
Definition: DMatSMatMultExpr.h:170
Header file for the DMatSVecMultExprTrait class template.
Header file for the IsColumnVector type trait.
Header file for the IsResizable type trait.
Header file for the thresholds for matrix/vector and matrix/matrix multiplications.
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
#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.
RightOperand rightOperand() const
Returns the right-hand side sparse matrix operand.
Definition: DMatSMatMultExpr.h:275