35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
92 template<
typename MT1
132 template<
typename MT >
135 enum { value = useAssign };
148 template<
typename MT >
149 struct UseSMPAssign {
150 enum { value = ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign };
178 enum { vectorizable = 0 };
181 enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
229 return lhs_.columns();
259 template<
typename T >
272 template<
typename T >
274 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
284 return lhs_.isAligned() &&
rhs_.isAligned();
294 return lhs_.canSMPAssign() ||
rhs_.canSMPAssign() ||
319 template<
typename MT
329 const size_t m( rhs.
rows() );
330 const size_t n( rhs.
columns() );
331 const size_t block( 16UL );
333 for(
size_t ii=0UL; ii<m; ii+=block ) {
334 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
335 for(
size_t jj=0UL; jj<n; jj+=block ) {
336 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
337 for(
size_t i=ii; i<iend; ++i ) {
338 for(
size_t j=jj; j<jend; ++j ) {
339 (~lhs)(i,j) = rhs.
lhs_(i,j) - rhs.
rhs_(i,j);
362 template<
typename MT
395 template<
typename MT
401 typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
413 const TmpType tmp(
serial( rhs ) );
434 template<
typename MT
436 friend inline typename DisableIf< UseAssign<MT> >::Type
444 const size_t m( rhs.rows() );
445 const size_t n( rhs.columns() );
446 const size_t block( 16UL );
448 for(
size_t ii=0UL; ii<m; ii+=block ) {
449 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
450 for(
size_t jj=0UL; jj<n; jj+=block ) {
451 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
452 for(
size_t i=ii; i<iend; ++i ) {
453 for(
size_t j=jj; j<jend; ++j ) {
454 (~lhs)(i,j) += rhs.lhs_(i,j) - rhs.rhs_(i,j);
478 template<
typename MT
480 friend inline typename EnableIf< UseAssign<MT> >::Type
513 template<
typename MT
515 friend inline typename DisableIf< UseAssign<MT> >::Type
523 const size_t m( rhs.rows() );
524 const size_t n( rhs.columns() );
525 const size_t block( 16UL );
527 for(
size_t ii=0UL; ii<m; ii+=block ) {
528 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
529 for(
size_t jj=0UL; jj<n; jj+=block ) {
530 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
531 for(
size_t i=ii; i<iend; ++i ) {
532 for(
size_t j=jj; j<jend; ++j ) {
533 (~lhs)(i,j) -= rhs.lhs_(i,j) - rhs.rhs_(i,j);
557 template<
typename MT
559 friend inline typename EnableIf< UseAssign<MT> >::Type
599 template<
typename MT
601 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
609 if( !IsExpression<MT1>::value &&
isSame( ~lhs, rhs.lhs_ ) ) {
634 template<
typename MT
636 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
641 typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
653 const TmpType tmp( rhs );
674 template<
typename MT
676 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
709 template<
typename MT
711 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
786 template<
typename T1
788 inline typename EnableIf< And< Not< IsSymmetric<T1> >, Not< IsSymmetric<T2> > >
789 ,
const DMatTDMatSubExpr<T1,T2> >::Type
795 throw std::invalid_argument(
"Matrix sizes do not match" );
831 template<
typename T1
833 inline typename EnableIf< And< Not< IsSymmetric<T1> >, Not< IsSymmetric<T2> > >
834 ,
const DMatTDMatSubExpr<T1,T2> >::Type
840 throw std::invalid_argument(
"Matrix sizes do not match" );
869 template<
typename T1
871 inline typename EnableIf< And< IsSymmetric<T1>, Not< IsSymmetric<T2> > >
872 ,
const typename SubExprTrait<T1,T2>::Type >::Type
873 operator-(
const DenseMatrix<T1,false>& lhs,
const DenseMatrix<T2,true>& rhs )
878 throw std::invalid_argument(
"Matrix sizes do not match" );
880 return trans( ~lhs ) - ~rhs;
900 template<
typename T1
902 inline typename EnableIf< IsSymmetric<T2>
903 ,
const typename SubExprTrait<T1,T2>::Type >::Type
904 operator-(
const DenseMatrix<T1,false>& lhs,
const DenseMatrix<T2,true>& rhs )
909 throw std::invalid_argument(
"Matrix sizes do not match" );
911 return (~lhs) -
trans( ~rhs );
931 template<
typename T1
933 inline typename EnableIf< And< Not< IsSymmetric<T1> >, IsSymmetric<T2> >
934 ,
const typename SubExprTrait<T1,T2>::Type >::Type
935 operator-(
const DenseMatrix<T1,true>& lhs,
const DenseMatrix<T2,false>& rhs )
940 throw std::invalid_argument(
"Matrix sizes do not match" );
942 return (~lhs) -
trans( ~rhs );
962 template<
typename T1
964 inline typename EnableIf< IsSymmetric<T1>
965 ,
const typename SubExprTrait<T1,T2>::Type >::Type
966 operator-(
const DenseMatrix<T1,true>& lhs,
const DenseMatrix<T2,false>& rhs )
971 throw std::invalid_argument(
"Matrix sizes do not match" );
973 return trans( ~lhs ) - (~rhs);
989 template<
typename MT1,
typename MT2 >
991 :
public Max< Rows<MT1>, Rows<MT2> >::Type
1007 template<
typename MT1,
typename MT2 >
1009 :
public Max< Columns<MT1>, Columns<MT2> >::Type
1025 template<
typename MT1,
typename MT2 >
1027 :
public IsTrue< IsSymmetric<MT1>::value && IsSymmetric<MT2>::value >
1043 template<
typename MT1,
typename MT2 >
1045 :
public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1061 template<
typename MT1,
typename MT2 >
1063 :
public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1079 template<
typename MT1,
typename MT2,
bool AF >
1084 typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1085 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1094 template<
typename MT1,
typename MT2 >
1099 typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
1100 ,
typename RowExprTrait<const MT2>::Type >::Type Type;
1109 template<
typename MT1,
typename MT2 >
1114 typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
1115 ,
typename ColumnExprTrait<const MT2>::Type >::Type Type;
Constraint on the data type.
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:103
Header file for the Max class template.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTDMatSubExpr.h:260
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 Rows type trait.
Header file for the subtraction trait.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:160
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: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
BLAZE_ALWAYS_INLINE 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:946
Header file for the ColumnExprTrait class template.
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatTDMatSubExpr.h:118
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:301
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
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:104
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:102
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:164
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:101
Constraint on the data type.
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:94
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTDMatSubExpr.h:167
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 IsTemporary type trait class.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:100
Header file for the IsSymmetric type trait.
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:273
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.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTDMatSubExpr.h:161
Header file for the MatMatSubExpr base class.
DMatTDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTDMatSubExpr class.
Definition: DMatTDMatSubExpr.h:190
Header file for the IsLower type trait.
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:302
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:158
#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:104
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 EnableIf class template.
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:283
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:105
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 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:2477
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
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
DMatTDMatSubExpr< MT1, MT2 > This
Type of this DMatTDMatSubExpr instance.
Definition: DMatTDMatSubExpr.h:157
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
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTDMatSubExpr.h:206
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTDMatSubExpr.h:218
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:170
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:173
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:932
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
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
#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:238
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2473
Header file for the IsTrue value trait.
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: DMatTDMatSubExpr.h:248
Header file for basic type definitions.
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_MATMATSUBEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid matrix/matrix ...
Definition: MatMatSubExpr.h:165
Base template for the SubTrait class.
Definition: SubTrait.h:142
Header file for the IsUpper type trait.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTDMatSubExpr.h:159
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTDMatSubExpr.h:293
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTDMatSubExpr.h:228
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.
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