All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DMatDMatSubExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATDMATSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATDMATSUBEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <stdexcept>
50 #include <blaze/math/Intrinsics.h>
62 #include <blaze/util/Assert.h>
64 #include <blaze/util/EnableIf.h>
66 #include <blaze/util/SelectType.h>
67 #include <blaze/util/Types.h>
69 
70 
71 namespace blaze {
72 
73 //=================================================================================================
74 //
75 // CLASS DMATDMATSUBEXPR
76 //
77 //=================================================================================================
78 
79 //*************************************************************************************************
86 template< typename MT1 // Type of the left-hand side dense matrix
87  , typename MT2 // Type of the right-hand side dense matrix
88  , bool SO > // Storage order
89 class DMatDMatSubExpr : public DenseMatrix< DMatDMatSubExpr<MT1,MT2,SO>, SO >
90  , private MatMatSubExpr
91  , private Computation
92 {
93  private:
94  //**Type definitions****************************************************************************
95  typedef typename MT1::ResultType RT1;
96  typedef typename MT2::ResultType RT2;
97  typedef typename MT1::ReturnType RN1;
98  typedef typename MT2::ReturnType RN2;
99  typedef typename MT1::CompositeType CT1;
100  typedef typename MT2::CompositeType CT2;
101  typedef typename MT1::ElementType ET1;
102  typedef typename MT2::ElementType ET2;
103  //**********************************************************************************************
104 
105  //**Return type evaluation**********************************************************************
107 
112  enum { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
113 
116  //**********************************************************************************************
117 
118  //**Evaluation strategy*************************************************************************
120 
126  enum { useAssign = RequiresEvaluation<MT1>::value || RequiresEvaluation<MT2>::value || !returnExpr };
127 
129  template< typename VT >
131  struct UseAssign {
132  enum { value = useAssign };
133  };
135  //**********************************************************************************************
136 
137  public:
138  //**Type definitions****************************************************************************
145 
148 
151 
153  typedef typename SelectType< IsExpression<MT1>::value, const MT1, const MT1& >::Type LeftOperand;
154 
156  typedef typename SelectType< IsExpression<MT2>::value, const MT2, const MT2& >::Type RightOperand;
157  //**********************************************************************************************
158 
159  //**ConstIterator class definition**************************************************************
163  {
164  public:
165  //**Type definitions*************************************************************************
166  typedef std::random_access_iterator_tag IteratorCategory;
171 
172  // STL iterator requirements
178 
181 
184  //*******************************************************************************************
185 
186  //**Constructor******************************************************************************
192  explicit inline ConstIterator( LeftIteratorType left, RightIteratorType right )
193  : left_ ( left ) // Iterator to the current left-hand side element
194  , right_( right ) // Iterator to the current right-hand side element
195  {}
196  //*******************************************************************************************
197 
198  //**Addition assignment operator*************************************************************
204  inline ConstIterator& operator+=( size_t inc ) {
205  left_ += inc;
206  right_ += inc;
207  return *this;
208  }
209  //*******************************************************************************************
210 
211  //**Subtraction assignment operator**********************************************************
217  inline ConstIterator& operator-=( size_t dec ) {
218  left_ -= dec;
219  right_ -= dec;
220  return *this;
221  }
222  //*******************************************************************************************
223 
224  //**Prefix increment operator****************************************************************
230  ++left_;
231  ++right_;
232  return *this;
233  }
234  //*******************************************************************************************
235 
236  //**Postfix increment operator***************************************************************
241  inline const ConstIterator operator++( int ) {
242  return ConstIterator( left_++, right_++ );
243  }
244  //*******************************************************************************************
245 
246  //**Prefix decrement operator****************************************************************
252  --left_;
253  --right_;
254  return *this;
255  }
256  //*******************************************************************************************
257 
258  //**Postfix decrement operator***************************************************************
263  inline const ConstIterator operator--( int ) {
264  return ConstIterator( left_--, right_-- );
265  }
266  //*******************************************************************************************
267 
268  //**Element access operator******************************************************************
273  inline ReturnType operator*() const {
274  return (*left_) - (*right_);
275  }
276  //*******************************************************************************************
277 
278  //**Load function****************************************************************************
283  inline IntrinsicType load() const {
284  return left_.load() - right_.load();
285  }
286  //*******************************************************************************************
287 
288  //**Equality operator************************************************************************
294  inline bool operator==( const ConstIterator& rhs ) const {
295  return left_ == rhs.left_;
296  }
297  //*******************************************************************************************
298 
299  //**Inequality operator**********************************************************************
305  inline bool operator!=( const ConstIterator& rhs ) const {
306  return left_ != rhs.left_;
307  }
308  //*******************************************************************************************
309 
310  //**Less-than operator***********************************************************************
316  inline bool operator<( const ConstIterator& rhs ) const {
317  return left_ < rhs.left_;
318  }
319  //*******************************************************************************************
320 
321  //**Greater-than operator********************************************************************
327  inline bool operator>( const ConstIterator& rhs ) const {
328  return left_ > rhs.left_;
329  }
330  //*******************************************************************************************
331 
332  //**Less-or-equal-than operator**************************************************************
338  inline bool operator<=( const ConstIterator& rhs ) const {
339  return left_ <= rhs.left_;
340  }
341  //*******************************************************************************************
342 
343  //**Greater-or-equal-than operator***********************************************************
349  inline bool operator>=( const ConstIterator& rhs ) const {
350  return left_ >= rhs.left_;
351  }
352  //*******************************************************************************************
353 
354  //**Subtraction operator*********************************************************************
360  inline DifferenceType operator-( const ConstIterator& rhs ) const {
361  return left_ - rhs.left_;
362  }
363  //*******************************************************************************************
364 
365  //**Addition operator************************************************************************
372  friend inline const ConstIterator operator+( const ConstIterator& it, size_t inc ) {
373  return ConstIterator( it.left_ + inc, it.right_ + inc );
374  }
375  //*******************************************************************************************
376 
377  //**Addition operator************************************************************************
384  friend inline const ConstIterator operator+( size_t inc, const ConstIterator& it ) {
385  return ConstIterator( it.left_ + inc, it.right_ + inc );
386  }
387  //*******************************************************************************************
388 
389  //**Subtraction operator*********************************************************************
396  friend inline const ConstIterator operator-( const ConstIterator& it, size_t dec ) {
397  return ConstIterator( it.left_ - dec, it.right_ - dec );
398  }
399  //*******************************************************************************************
400 
401  private:
402  //**Member variables*************************************************************************
405  //*******************************************************************************************
406  };
407  //**********************************************************************************************
408 
409  //**Compilation flags***************************************************************************
411  enum { vectorizable = MT1::vectorizable && MT2::vectorizable &&
414 
416  enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
417  //**********************************************************************************************
418 
419  //**Constructor*********************************************************************************
425  explicit inline DMatDMatSubExpr( const MT1& lhs, const MT2& rhs )
426  : lhs_( lhs ) // Left-hand side dense matrix of the subtraction expression
427  , rhs_( rhs ) // Right-hand side dense matrix of the subtraction expression
428  {
429  BLAZE_INTERNAL_ASSERT( lhs.rows() == rhs.rows() , "Invalid number of rows" );
430  BLAZE_INTERNAL_ASSERT( lhs.columns() == rhs.columns(), "Invalid number of columns" );
431  }
432  //**********************************************************************************************
433 
434  //**Access operator*****************************************************************************
441  inline ReturnType operator()( size_t i, size_t j ) const {
442  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
443  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
444  return lhs_(i,j) - rhs_(i,j);
445  }
446  //**********************************************************************************************
447 
448  //**Load function*******************************************************************************
455  inline IntrinsicType load( size_t i, size_t j ) const {
456  typedef IntrinsicTrait<ElementType> IT;
457  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
458  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
459  BLAZE_INTERNAL_ASSERT( !SO || ( i % IT::size == 0UL ), "Invalid row access index" );
460  BLAZE_INTERNAL_ASSERT( SO || ( j % IT::size == 0UL ), "Invalid column access index" );
461  const IntrinsicType xmm1( lhs_.load(i,j) );
462  const IntrinsicType xmm2( rhs_.load(i,j) );
463  return xmm1 - xmm2;
464  }
465  //**********************************************************************************************
466 
467  //**Begin function******************************************************************************
473  inline ConstIterator begin( size_t i ) const {
474  return ConstIterator( lhs_.begin(i), rhs_.begin(i) );
475  }
476  //**********************************************************************************************
477 
478  //**End function********************************************************************************
484  inline ConstIterator end( size_t i ) const {
485  return ConstIterator( lhs_.end(i), rhs_.end(i) );
486  }
487  //**********************************************************************************************
488 
489  //**Rows function*******************************************************************************
494  inline size_t rows() const {
495  return lhs_.rows();
496  }
497  //**********************************************************************************************
498 
499  //**Columns function****************************************************************************
504  inline size_t columns() const {
505  return lhs_.columns();
506  }
507  //**********************************************************************************************
508 
509  //**Left operand access*************************************************************************
514  inline LeftOperand leftOperand() const {
515  return lhs_;
516  }
517  //**********************************************************************************************
518 
519  //**Right operand access************************************************************************
524  inline RightOperand rightOperand() const {
525  return rhs_;
526  }
527  //**********************************************************************************************
528 
529  //**********************************************************************************************
535  template< typename T >
536  inline bool canAlias( const T* alias ) const {
537  return ( IsExpression<MT1>::value && ( RequiresEvaluation<MT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
538  ( IsExpression<MT2>::value && ( RequiresEvaluation<MT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
539  }
540  //**********************************************************************************************
541 
542  //**********************************************************************************************
548  template< typename T >
549  inline bool isAliased( const T* alias ) const {
550  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
551  }
552  //**********************************************************************************************
553 
554  //**********************************************************************************************
559  inline bool isAligned() const {
560  return lhs_.isAligned() && rhs_.isAligned();
561  }
562  //**********************************************************************************************
563 
564  //**********************************************************************************************
569  inline bool canSMPAssign() const {
570  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
571  ( ( ( SO == rowMajor ) ? rows() : columns() ) > SMP_DMATDMATSUB_THRESHOLD );
572  }
573  //**********************************************************************************************
574 
575  private:
576  //**Member variables****************************************************************************
579  //**********************************************************************************************
580 
581  //**Assignment to dense matrices****************************************************************
595  template< typename MT // Type of the target dense matrix
596  , bool SO2 > // Storage order of the target dense matrix
597  friend inline typename EnableIf< UseAssign<MT> >::Type
598  assign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
599  {
601 
602  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
603  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
604 
605  if( !IsExpression<MT1>::value && (~lhs).isAliased( &rhs.lhs_ ) ) {
606  smpSubAssign( ~lhs, rhs.rhs_ );
607  }
608  else {
609  smpAssign ( ~lhs, rhs.lhs_ );
610  smpSubAssign( ~lhs, rhs.rhs_ );
611  }
612  }
614  //**********************************************************************************************
615 
616  //**Assignment to sparse matrices***************************************************************
630  template< typename MT // Type of the target sparse matrix
631  , bool SO2 > // Storage order of the target sparse matrix
632  friend inline typename EnableIf< UseAssign<MT> >::Type
633  assign( SparseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
634  {
636 
638 
645 
646  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
647  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
648 
649  const TmpType tmp( rhs );
650  smpAssign( ~lhs, tmp );
651  }
653  //**********************************************************************************************
654 
655  //**Addition assignment to dense matrices*******************************************************
669  template< typename MT // Type of the target dense matrix
670  , bool SO2 > // Storage order of the target dense matrix
671  friend inline typename EnableIf< UseAssign<MT> >::Type
672  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
673  {
675 
676  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
677  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
678 
679  smpAddAssign( ~lhs, rhs.lhs_ );
680  smpSubAssign( ~lhs, rhs.rhs_ );
681  }
683  //**********************************************************************************************
684 
685  //**Addition assignment to sparse matrices******************************************************
686  // No special implementation for the addition assignment to sparse matrices.
687  //**********************************************************************************************
688 
689  //**Subtraction assignment to dense matrices****************************************************
703  template< typename MT // Type of the target dense matrix
704  , bool SO2 > // Storage order of the target dense matrix
705  friend inline typename EnableIf< UseAssign<MT> >::Type
706  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
707  {
709 
710  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
711  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
712 
713  smpSubAssign( ~lhs, rhs.lhs_ );
714  smpAddAssign( ~lhs, rhs.rhs_ );
715  }
717  //**********************************************************************************************
718 
719  //**Subtraction assignment to sparse matrices***************************************************
720  // No special implementation for the subtraction assignment to sparse matrices.
721  //**********************************************************************************************
722 
723  //**Multiplication assignment to dense matrices*************************************************
724  // No special implementation for the multiplication assignment to dense matrices.
725  //**********************************************************************************************
726 
727  //**Multiplication assignment to sparse matrices************************************************
728  // No special implementation for the multiplication assignment to sparse matrices.
729  //**********************************************************************************************
730 
731  //**Compile time checks*************************************************************************
736  //**********************************************************************************************
737 };
738 //*************************************************************************************************
739 
740 
741 
742 
743 //=================================================================================================
744 //
745 // GLOBAL BINARY ARITHMETIC OPERATORS
746 //
747 //=================================================================================================
748 
749 //*************************************************************************************************
774 template< typename T1 // Type of the left-hand side dense matrix
775  , typename T2 // Type of the right-hand side dense matrix
776  , bool SO > // Storage order
777 inline const DMatDMatSubExpr<T1,T2,SO>
779 {
781 
782  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
783  throw std::invalid_argument( "Matrix sizes do not match" );
784 
785  return DMatDMatSubExpr<T1,T2,SO>( ~lhs, ~rhs );
786 }
787 //*************************************************************************************************
788 
789 
790 
791 
792 //=================================================================================================
793 //
794 // EXPRESSION TRAIT SPECIALIZATIONS
795 //
796 //=================================================================================================
797 
798 //*************************************************************************************************
800 template< typename MT1, typename MT2, bool SO, bool AF >
801 struct SubmatrixExprTrait< DMatDMatSubExpr<MT1,MT2,SO>, AF >
802 {
803  public:
804  //**********************************************************************************************
805  typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
806  , typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
807  //**********************************************************************************************
808 };
810 //*************************************************************************************************
811 
812 
813 //*************************************************************************************************
815 template< typename MT1, typename MT2, bool SO >
816 struct RowExprTrait< DMatDMatSubExpr<MT1,MT2,SO> >
817 {
818  public:
819  //**********************************************************************************************
820  typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
821  , typename RowExprTrait<const MT2>::Type >::Type Type;
822  //**********************************************************************************************
823 };
825 //*************************************************************************************************
826 
827 
828 //*************************************************************************************************
830 template< typename MT1, typename MT2, bool SO >
831 struct ColumnExprTrait< DMatDMatSubExpr<MT1,MT2,SO> >
832 {
833  public:
834  //**********************************************************************************************
835  typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
836  , typename ColumnExprTrait<const MT2>::Type >::Type Type;
837  //**********************************************************************************************
838 };
840 //*************************************************************************************************
841 
842 } // namespace blaze
843 
844 #endif
Pointer difference type of the Blaze library.
const ConstIterator operator++(int)
Post-increment operator.
Definition: DMatDMatSubExpr.h:241
ReturnType operator*() const
Direct access to the element at the current iterator position.
Definition: DMatDMatSubExpr.h:273
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
ConstIterator & operator++()
Pre-increment operator.
Definition: DMatDMatSubExpr.h:229
Header file for the subtraction trait.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatDMatSubExpr.h:536
bool operator>(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:327
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DMatDMatSubExpr.h:147
ElementType * PointerType
Pointer return type.
Definition: DMatDMatSubExpr.h:168
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
IntrinsicTrait< ElementType >::Type IntrinsicType
Resulting intrinsic element type.
Definition: DMatDMatSubExpr.h:144
#define BLAZE_CONSTRAINT_MUST_BE_MATRIX_WITH_STORAGE_ORDER(T, SO)
Constraint on the data type.In case the given data type T is not a dense or sparse matrix type and in...
Definition: StorageOrder.h:242
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:197
DMatDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatDMatSubExpr class.
Definition: DMatDMatSubExpr.h:425
#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
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatDMatSubExpr.h:441
Header file for the ColumnExprTrait class template.
Header file for the IsSame and IsStrictlySame type traits.
IntrinsicType load() const
Access to the intrinsic elements of the matrix.
Definition: DMatDMatSubExpr.h:283
Header file for the sparse matrix SMP implementation.
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2384
MT1::ReturnType RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:97
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:249
ConstIterator & operator--()
Pre-decrement operator.
Definition: DMatDMatSubExpr.h:251
PointerType pointer
Pointer return type.
Definition: DMatDMatSubExpr.h:175
Header file for the Computation base class.
Type relationship analysis.This class tests if the two data types A and B are equal. For this type comparison, the cv-qualifiers of both data types are ignored. If A and B are the same data type (ignoring the cv-qualifiers), then the value member enumeration is set to 1, the nested type definition Type is TrueType, and the class derives from TrueType. Otherwise value is set to 0, Type is FalseType, and the class derives from FalseType.
Definition: IsSame.h:158
MT2::ConstIterator RightIteratorType
ConstIterator type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:183
Header file for the RequiresEvaluation type trait.
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Iterator over the elements of the dense matrix.
Definition: DMatDMatSubExpr.h:162
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.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatDMatSubExpr.h:143
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:99
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
DifferenceType difference_type
Difference between two iterators.
Definition: DMatDMatSubExpr.h:177
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatDMatSubExpr.h:578
Expression object for dense matrix-dense matrix subtractions.The DMatDMatSubExpr class represents the...
Definition: DMatDMatSubExpr.h:89
Compile time type selection.The SelectType class template selects one of the two given types T1 and T...
Definition: SelectType.h:59
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatDMatSubExpr.h:494
Header file for the IsTemporary type trait class.
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatDMatSubExpr.h:140
Base class for all matrix/matrix subtraction expression templates.The MatMatSubExpr class serves as a...
Definition: MatMatSubExpr.h:65
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row i.
Definition: DMatDMatSubExpr.h:473
ElementType ValueType
Type of the underlying elements.
Definition: DMatDMatSubExpr.h:167
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2388
Header file for the dense matrix SMP implementation.
Header file for the DenseMatrix base class.
DMatDMatSubExpr< MT1, MT2, SO > This
Type of this DMatDMatSubExpr instance.
Definition: DMatDMatSubExpr.h:139
std::random_access_iterator_tag IteratorCategory
The iterator category.
Definition: DMatDMatSubExpr.h:166
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:153
RightIteratorType right_
Iterator to the current right-hand side element.
Definition: DMatDMatSubExpr.h:404
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 MatMatSubExpr base class.
ReferenceType reference
Reference return type.
Definition: DMatDMatSubExpr.h:176
const size_t SMP_DMATDMATSUB_THRESHOLD
SMP row-major dense matrix/row-major dense matrix subtraction threshold.This threshold represents the...
Definition: Thresholds.h:392
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatDMatSubExpr.h:569
ConstIterator(LeftIteratorType left, RightIteratorType right)
Constructor for the ConstIterator class.
Definition: DMatDMatSubExpr.h:192
#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 isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatDMatSubExpr.h:559
Constraints on the storage order of matrix types.
MT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:101
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row i.
Definition: DMatDMatSubExpr.h:484
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:103
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2382
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatDMatSubExpr.h:142
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:156
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.
bool operator<=(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:338
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatDMatSubExpr.h:141
ValueType value_type
Type of the underlying elements.
Definition: DMatDMatSubExpr.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
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:305
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two iterators.
Definition: DMatDMatSubExpr.h:360
Header file for the SubmatrixExprTrait class template.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2383
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:100
Intrinsic characteristics of data types.The IntrinsicTrait class template provides the intrinsic char...
Definition: IntrinsicTrait.h:748
Header file for run time assertion macros.
MT2::ElementType ET2
Element type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:102
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
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:294
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
friend const ConstIterator operator+(size_t inc, const ConstIterator &it)
Addition between an integral value and a ConstIterator.
Definition: DMatDMatSubExpr.h:384
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: DMatDMatSubExpr.h:98
ConstIterator & operator-=(size_t dec)
Subtraction assignment operator.
Definition: DMatDMatSubExpr.h:217
friend const ConstIterator operator-(const ConstIterator &it, size_t dec)
Subtraction between a ConstIterator and an integral value.
Definition: DMatDMatSubExpr.h:396
#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
bool operator>=(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:349
RightOperand rightOperand() const
Returns the right-hand side dense matrix operand.
Definition: DMatDMatSubExpr.h:524
Header file for all intrinsic functionality.
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DMatDMatSubExpr.h:170
IntrinsicType load(size_t i, size_t j) const
Access to the intrinsic elements of the matrix.
Definition: DMatDMatSubExpr.h:455
ElementType & ReferenceType
Reference return type.
Definition: DMatDMatSubExpr.h:169
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:96
const bool rowMajor
Storage order flag for row-major matrices.
Definition: StorageOrder.h:71
SelectType< useAssign, const ResultType, const DMatDMatSubExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: DMatDMatSubExpr.h:150
IteratorCategory iterator_category
The iterator category.
Definition: DMatDMatSubExpr.h:173
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatDMatSubExpr.h:514
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatDMatSubExpr.h:115
#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
LeftIteratorType left_
Iterator to the current left-hand side element.
Definition: DMatDMatSubExpr.h:403
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2379
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatDMatSubExpr.h:549
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.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:95
bool operator<(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:316
MT1::ConstIterator LeftIteratorType
ConstIterator type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:180
Base template for the SubTrait class.
Definition: SubTrait.h:141
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatDMatSubExpr.h:577
friend const ConstIterator operator+(const ConstIterator &it, size_t inc)
Addition between a ConstIterator and an integral value.
Definition: DMatDMatSubExpr.h:372
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatDMatSubExpr.h:504
const ConstIterator operator--(int)
Post-decrement operator.
Definition: DMatDMatSubExpr.h:263
Header file for the SubExprTrait class template.
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
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.
ConstIterator & operator+=(size_t inc)
Addition assignment operator.
Definition: DMatDMatSubExpr.h:204