35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
82 template<
typename MT1
122 template<
typename MT >
125 enum { value = useAssign };
138 template<
typename MT >
139 struct UseSMPAssign {
140 enum { value = ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign };
168 enum { vectorizable = 0 };
171 enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
219 return lhs_.columns();
249 template<
typename T >
262 template<
typename T >
264 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
274 return lhs_.isAligned() &&
rhs_.isAligned();
284 return lhs_.canSMPAssign() ||
rhs_.canSMPAssign() ||
309 template<
typename MT
319 const size_t m( rhs.
rows() );
320 const size_t n( rhs.
columns() );
321 const size_t block( 16UL );
323 for(
size_t ii=0UL; ii<m; ii+=block ) {
324 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
325 for(
size_t jj=0UL; jj<n; jj+=block ) {
326 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
327 for(
size_t i=ii; i<iend; ++i ) {
328 for(
size_t j=jj; j<jend; ++j ) {
329 (~lhs)(i,j) = rhs.
lhs_(i,j) - rhs.
rhs_(i,j);
352 template<
typename MT
385 template<
typename MT
391 typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
403 const TmpType tmp(
serial( rhs ) );
424 template<
typename MT
426 friend inline typename DisableIf< UseAssign<MT> >::Type
434 const size_t m( rhs.rows() );
435 const size_t n( rhs.columns() );
436 const size_t block( 16UL );
438 for(
size_t ii=0UL; ii<m; ii+=block ) {
439 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
440 for(
size_t jj=0UL; jj<n; jj+=block ) {
441 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
442 for(
size_t i=ii; i<iend; ++i ) {
443 for(
size_t j=jj; j<jend; ++j ) {
444 (~lhs)(i,j) += rhs.lhs_(i,j) - rhs.rhs_(i,j);
468 template<
typename MT
470 friend inline typename EnableIf< UseAssign<MT> >::Type
503 template<
typename MT
505 friend inline typename DisableIf< UseAssign<MT> >::Type
513 const size_t m( rhs.rows() );
514 const size_t n( rhs.columns() );
515 const size_t block( 16UL );
517 for(
size_t ii=0UL; ii<m; ii+=block ) {
518 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
519 for(
size_t jj=0UL; jj<n; jj+=block ) {
520 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
521 for(
size_t i=ii; i<iend; ++i ) {
522 for(
size_t j=jj; j<jend; ++j ) {
523 (~lhs)(i,j) -= rhs.lhs_(i,j) - rhs.rhs_(i,j);
547 template<
typename MT
549 friend inline typename EnableIf< UseAssign<MT> >::Type
589 template<
typename MT
591 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
599 if( !IsExpression<MT1>::value &&
isSame( ~lhs, rhs.lhs_ ) ) {
624 template<
typename MT
626 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
631 typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
643 const TmpType tmp( rhs );
664 template<
typename MT
666 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
699 template<
typename MT
701 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
775 template<
typename T1
777 inline const DMatTDMatSubExpr<T1,T2>
783 throw std::invalid_argument(
"Matrix sizes do not match" );
819 template<
typename T1
821 inline const DMatTDMatSubExpr<T1,T2>
827 throw std::invalid_argument(
"Matrix sizes do not match" );
844 template<
typename MT1,
typename MT2,
bool AF >
849 typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
850 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
859 template<
typename MT1,
typename MT2 >
864 typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
865 ,
typename RowExprTrait<const MT2>::Type >::Type Type;
874 template<
typename MT1,
typename MT2 >
879 typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
880 ,
typename ColumnExprTrait<const MT2>::Type >::Type Type;
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:93
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTDMatSubExpr.h:250
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
Header file for the subtraction trait.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:150
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:152
const size_t SMP_DMATTDMATSUB_THRESHOLD
SMP row-major dense matrix/column-major dense matrix subtraction threshold.This threshold specifies w...
Definition: Thresholds.h:788
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:199
#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.
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatTDMatSubExpr.h:108
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2408
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:251
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:291
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:690
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:94
Header file for the Computation base class.
Header file for the RequiresEvaluation type trait.
MT1::ReturnType RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:92
bool isSame(const Matrix< MT1, SO1 > &a, const Matrix< MT2, SO2 > &b)
Returns whether the two given matrices represent the same observable state.
Definition: Matrix.h:409
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Constraint on the data type.
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:154
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:91
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:122
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:90
Expression object for dense matrix-transpose dense matrix subtractions.The DMatTDMatSubExpr class rep...
Definition: DMatTDMatSubExpr.h:84
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTDMatSubExpr.h:157
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.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:90
Base class for all matrix/matrix subtraction expression templates.The MatMatSubExpr class serves as a...
Definition: MatMatSubExpr.h:65
#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
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatTDMatSubExpr.h:263
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:271
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTDMatSubExpr.h:151
Header file for the MatMatSubExpr base class.
DMatTDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTDMatSubExpr class.
Definition: DMatTDMatSubExpr.h:180
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:292
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:148
#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.
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:103
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 serial shim.
#define BLAZE_CONSTRAINT_MATRICES_MUST_HAVE_DIFFERENT_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:325
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatTDMatSubExpr.h:273
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:95
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:92
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 Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2407
Header file for run time assertion macros.
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:301
DMatTDMatSubExpr< MT1, MT2 > This
Type of this DMatTDMatSubExpr instance.
Definition: DMatTDMatSubExpr.h:147
const DenseIterator< Type > operator-(const DenseIterator< Type > &it, ptrdiff_t inc)
Subtraction between a DenseIterator and an integral value.
Definition: DenseIterator.h:585
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
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:331
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTDMatSubExpr.h:196
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTDMatSubExpr.h:208
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:160
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
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.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
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatTDMatSubExpr.h:228
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2403
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: DMatTDMatSubExpr.h:238
Header file for basic type definitions.
Base template for the SubTrait class.
Definition: SubTrait.h:141
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTDMatSubExpr.h:149
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTDMatSubExpr.h:283
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTDMatSubExpr.h:218
Header file for the SubExprTrait class template.
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
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.