35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTDMATADDEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTDMATADDEXPR_H_
83 template<
typename MT1
135 enum { vectorizable = 0 };
138 enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
186 return lhs_.columns();
216 template<
typename T >
229 template<
typename T >
231 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
241 return lhs_.isAligned() &&
rhs_.isAligned();
251 return lhs_.canSMPAssign() ||
rhs_.canSMPAssign() ||
274 template<
typename MT
287 const size_t m( rhs.
rows() );
288 const size_t n( rhs.
columns() );
289 const size_t block( 16UL );
291 for(
size_t ii=0UL; ii<m; ii+=block ) {
292 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
293 for(
size_t jj=0UL; jj<n; jj+=block ) {
294 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
295 for(
size_t i=ii; i<iend; ++i ) {
296 for(
size_t j=jj; j<jend; ++j ) {
297 (~lhs)(i,j) = rhs.
lhs_(i,j) + rhs.
rhs_(i,j);
333 template<
typename MT
351 const TmpType tmp( rhs );
369 template<
typename MT
380 if( !RequiresEvaluation<MT1>::value && !RequiresEvaluation<MT2>::value )
382 const size_t m( rhs.rows() );
383 const size_t n( rhs.columns() );
384 const size_t block( 16UL );
386 for(
size_t ii=0UL; ii<m; ii+=block ) {
387 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
388 for(
size_t jj=0UL; jj<n; jj+=block ) {
389 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
390 for(
size_t i=ii; i<iend; ++i ) {
391 for(
size_t j=jj; j<jend; ++j ) {
392 (~lhs)(i,j) += rhs.lhs_(i,j) + rhs.rhs_(i,j);
426 template<
typename MT
437 if( !RequiresEvaluation<MT1>::value && !RequiresEvaluation<MT2>::value )
439 const size_t m( rhs.rows() );
440 const size_t n( rhs.columns() );
441 const size_t block( 16UL );
443 for(
size_t ii=0UL; ii<m; ii+=block ) {
444 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
445 for(
size_t jj=0UL; jj<n; jj+=block ) {
446 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
447 for(
size_t i=ii; i<iend; ++i ) {
448 for(
size_t j=jj; j<jend; ++j ) {
449 (~lhs)(i,j) -= rhs.lhs_(i,j) + rhs.rhs_(i,j);
528 template<
typename T1
530 inline const DMatTDMatAddExpr<T1,T2>
536 throw std::invalid_argument(
"Matrix sizes do not match" );
572 template<
typename T1
574 inline const DMatTDMatAddExpr<T2,T1>
580 throw std::invalid_argument(
"Matrix sizes do not match" );
597 template<
typename MT1,
typename MT2,
bool AF >
598 struct SubmatrixExprTrait< DMatTDMatAddExpr<MT1,MT2>, AF >
602 typedef typename AddExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
603 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
612 template<
typename MT1,
typename MT2 >
613 struct RowExprTrait< DMatTDMatAddExpr<MT1,MT2> >
617 typedef typename AddExprTrait< typename RowExprTrait<const MT1>::Type
618 ,
typename RowExprTrait<const MT2>::Type >::Type Type;
627 template<
typename MT1,
typename MT2 >
628 struct ColumnExprTrait< DMatTDMatAddExpr<MT1,MT2> >
632 typedef typename AddExprTrait< typename ColumnExprTrait<const MT1>::Type
633 ,
typename ColumnExprTrait<const MT2>::Type >::Type Type;
Evaluation of the return type of an addition expression.Via this type trait it is possible to evaluat...
Definition: AddExprTrait.h:103
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
const size_t SMP_DMATTDMATADD_THRESHOLD
SMP row-major dense matrix/column-major dense matrix addition threshold.This threshold represents the...
Definition: Thresholds.h:377
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
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.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTDMatAddExpr.h:116
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTDMatAddExpr.h:117
Header file for the sparse matrix SMP implementation.
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2384
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:249
AddTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTDMatAddExpr.h:115
Header file for the AddExprTrait class template.
Header file for the Computation base class.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTDMatAddExpr.h:217
Header file for the RequiresEvaluation type trait.
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:95
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.
DMatTDMatAddExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTDMatAddExpr class.
Definition: DMatTDMatAddExpr.h:147
Constraint on the data type.
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
MT1::ReturnType RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:93
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTDMatAddExpr.h:175
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 IsTemporary type trait class.
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTDMatAddExpr.h:124
#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 SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DMatTDMatAddExpr.h:121
Header file for the dense matrix SMP implementation.
Header file for the DenseMatrix base class.
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
Base class for all matrix/matrix addition expression templates.The MatMatAddExpr class serves as a ta...
Definition: MatMatAddExpr.h:65
const DenseIterator< Type > operator+(const DenseIterator< Type > &it, ptrdiff_t inc)
Addition between a DenseIterator and an integral value.
Definition: DenseIterator.h:556
AddExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatTDMatAddExpr.h:109
#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
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTDMatAddExpr.h:185
Constraints on the storage order of matrix types.
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatTDMatAddExpr.h:230
DMatTDMatAddExpr< MT1, MT2 > This
Type of this DMatTDMatAdd instance.
Definition: DMatTDMatAddExpr.h:114
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:92
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
Header file for the SubmatrixExprTrait class template.
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:130
#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 Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2383
LeftOperand lhs_
Left-hand side dense matrix of the addition expression.
Definition: DMatTDMatAddExpr.h:258
Header file for run time assertion macros.
Base template for the AddTrait class.
Definition: AddTrait.h:141
Header file for the addition trait.
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
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTDMatAddExpr.h:118
RightOperand rhs_
Right-hand side dense matrix of the addition expression.
Definition: DMatTDMatAddExpr.h:259
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: DMatTDMatAddExpr.h:205
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatTDMatAddExpr.h:195
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
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:94
Header file for the MatMatAddExpr base class.
#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
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:127
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:91
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTDMatAddExpr.h:163
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
#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
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatTDMatAddExpr.h:240
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
Header file for basic type definitions.
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatAddExpr.h:96
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
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTDMatAddExpr.h:250
Compile time check whether the given type is an expression template.This type trait class tests wheth...
Definition: IsExpression.h:87
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
Expression object for dense matrix-transpose dense matrix additions.The DMatTDMatAddExpr class repres...
Definition: DMatTDMatAddExpr.h:85