35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
98 template<
typename MT1
138 template<
typename MT >
141 enum { value = useAssign };
154 template<
typename MT >
155 struct UseSMPAssign {
156 enum { value = ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign };
184 enum { vectorizable = 0 };
187 enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
235 return lhs_.columns();
265 template<
typename T >
278 template<
typename T >
280 return (
lhs_.isAliased( alias ) ||
rhs_.isAliased( alias ) );
290 return lhs_.isAligned() &&
rhs_.isAligned();
300 return lhs_.canSMPAssign() ||
rhs_.canSMPAssign() ||
325 template<
typename MT
335 const size_t m( rhs.
rows() );
336 const size_t n( rhs.
columns() );
337 const size_t block( 16UL );
339 for(
size_t ii=0UL; ii<m; ii+=block ) {
340 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
341 for(
size_t jj=0UL; jj<n; jj+=block ) {
342 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
343 for(
size_t i=ii; i<iend; ++i ) {
344 for(
size_t j=jj; j<jend; ++j ) {
345 (~lhs)(i,j) = rhs.
lhs_(i,j) - rhs.
rhs_(i,j);
368 template<
typename MT
401 template<
typename MT
407 typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
419 const TmpType tmp(
serial( rhs ) );
440 template<
typename MT
442 friend inline typename DisableIf< UseAssign<MT> >::Type
450 const size_t m( rhs.rows() );
451 const size_t n( rhs.columns() );
452 const size_t block( 16UL );
454 for(
size_t ii=0UL; ii<m; ii+=block ) {
455 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
456 for(
size_t jj=0UL; jj<n; jj+=block ) {
457 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
458 for(
size_t i=ii; i<iend; ++i ) {
459 for(
size_t j=jj; j<jend; ++j ) {
460 (~lhs)(i,j) += rhs.lhs_(i,j) - rhs.rhs_(i,j);
484 template<
typename MT
486 friend inline typename EnableIf< UseAssign<MT> >::Type
519 template<
typename MT
521 friend inline typename DisableIf< UseAssign<MT> >::Type
529 const size_t m( rhs.rows() );
530 const size_t n( rhs.columns() );
531 const size_t block( 16UL );
533 for(
size_t ii=0UL; ii<m; ii+=block ) {
534 const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
535 for(
size_t jj=0UL; jj<n; jj+=block ) {
536 const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
537 for(
size_t i=ii; i<iend; ++i ) {
538 for(
size_t j=jj; j<jend; ++j ) {
539 (~lhs)(i,j) -= rhs.lhs_(i,j) - rhs.rhs_(i,j);
563 template<
typename MT
565 friend inline typename EnableIf< UseAssign<MT> >::Type
605 template<
typename MT
607 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
615 if( !IsExpression<MT1>::value &&
isSame( ~lhs, rhs.lhs_ ) ) {
640 template<
typename MT
642 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
647 typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
659 const TmpType tmp( rhs );
680 template<
typename MT
682 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
715 template<
typename MT
717 friend inline typename EnableIf< UseSMPAssign<MT> >::Type
792 template<
typename T1
794 inline typename EnableIf< And< Not< IsSymmetric<T1> >, Not< IsSymmetric<T2> > >
795 ,
const DMatTDMatSubExpr<T1,T2> >::Type
801 throw std::invalid_argument(
"Matrix sizes do not match" );
837 template<
typename T1
839 inline typename EnableIf< And< Not< IsSymmetric<T1> >, Not< IsSymmetric<T2> > >
840 ,
const DMatTDMatSubExpr<T1,T2> >::Type
846 throw std::invalid_argument(
"Matrix sizes do not match" );
875 template<
typename T1
877 inline typename EnableIf< And< IsSymmetric<T1>, Not< IsSymmetric<T2> > >
878 ,
const typename SubExprTrait<T1,T2>::Type >::Type
879 operator-(
const DenseMatrix<T1,false>& lhs,
const DenseMatrix<T2,true>& rhs )
884 throw std::invalid_argument(
"Matrix sizes do not match" );
886 return trans( ~lhs ) - ~rhs;
906 template<
typename T1
908 inline typename EnableIf< IsSymmetric<T2>
909 ,
const typename SubExprTrait<T1,T2>::Type >::Type
910 operator-(
const DenseMatrix<T1,false>& lhs,
const DenseMatrix<T2,true>& rhs )
915 throw std::invalid_argument(
"Matrix sizes do not match" );
917 return (~lhs) -
trans( ~rhs );
937 template<
typename T1
939 inline typename EnableIf< And< Not< IsSymmetric<T1> >, IsSymmetric<T2> >
940 ,
const typename SubExprTrait<T1,T2>::Type >::Type
941 operator-(
const DenseMatrix<T1,true>& lhs,
const DenseMatrix<T2,false>& rhs )
946 throw std::invalid_argument(
"Matrix sizes do not match" );
948 return (~lhs) -
trans( ~rhs );
968 template<
typename T1
970 inline typename EnableIf< IsSymmetric<T1>
971 ,
const typename SubExprTrait<T1,T2>::Type >::Type
972 operator-(
const DenseMatrix<T1,true>& lhs,
const DenseMatrix<T2,false>& rhs )
977 throw std::invalid_argument(
"Matrix sizes do not match" );
979 return trans( ~lhs ) - (~rhs);
995 template<
typename MT1,
typename MT2 >
997 :
public Max< Rows<MT1>, Rows<MT2> >::Type
1013 template<
typename MT1,
typename MT2 >
1015 :
public Max< Columns<MT1>, Columns<MT2> >::Type
1031 template<
typename MT1,
typename MT2 >
1033 :
public IsTrue< IsSymmetric<MT1>::value && IsSymmetric<MT2>::value >
1049 template<
typename MT1,
typename MT2 >
1051 :
public IsTrue< And< IsLower<MT1>, IsLower<MT2> >::value >
1067 template<
typename MT1,
typename MT2 >
1069 :
public IsTrue< And< IsUniLower<MT1>, IsStrictlyLower<MT2> >::value >
1085 template<
typename MT1,
typename MT2 >
1087 :
public IsTrue< And< IsStrictlyLower<MT1>, IsStrictlyLower<MT2> >::value >
1103 template<
typename MT1,
typename MT2 >
1105 :
public IsTrue< And< IsUpper<MT1>, IsUpper<MT2> >::value >
1121 template<
typename MT1,
typename MT2 >
1123 :
public IsTrue< And< IsUniUpper<MT1>, IsStrictlyUpper<MT2> >::value >
1139 template<
typename MT1,
typename MT2 >
1141 :
public IsTrue< And< IsStrictlyUpper<MT1>, IsStrictlyUpper<MT2> >::value >
1157 template<
typename MT1,
typename MT2,
bool AF >
1162 typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1163 ,
typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1172 template<
typename MT1,
typename MT2 >
1177 typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
1178 ,
typename RowExprTrait<const MT2>::Type >::Type Type;
1187 template<
typename MT1,
typename MT2 >
1192 typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
1193 ,
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:109
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:266
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 IsUniUpper type trait.
Header file for the subtraction trait.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:166
Header file for basic type definitions.
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:209
#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:124
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2507
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:261
Header file for the And class template.
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:307
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:699
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:110
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:108
Header file for the IsUniLower type trait.
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:170
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:107
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:100
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTDMatSubExpr.h:173
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:106
Header file for the IsStrictlyUpper type trait.
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
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
#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:279
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 Not class template.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTDMatSubExpr.h:167
Header file for the MatMatSubExpr base class.
DMatTDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTDMatSubExpr class.
Definition: DMatTDMatSubExpr.h:196
Header file for the IsLower type trait.
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:308
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:164
#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 IsStrictlyLower type trait.
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:289
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:111
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:2506
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:163
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:212
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTDMatSubExpr.h:224
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:176
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:179
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:937
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:244
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2502
Header file for the IsTrue value trait.
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: DMatTDMatSubExpr.h:254
#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:150
Header file for the IsUpper type trait.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTDMatSubExpr.h:165
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTDMatSubExpr.h:299
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTDMatSubExpr.h:234
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