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>
51 #include <blaze/math/Intrinsics.h>
67 #include <blaze/util/Assert.h>
69 #include <blaze/util/EnableIf.h>
71 #include <blaze/util/mpl/Max.h>
72 #include <blaze/util/SelectType.h>
73 #include <blaze/util/Types.h>
76 
77 
78 namespace blaze {
79 
80 //=================================================================================================
81 //
82 // CLASS DMATDMATSUBEXPR
83 //
84 //=================================================================================================
85 
86 //*************************************************************************************************
93 template< typename MT1 // Type of the left-hand side dense matrix
94  , typename MT2 // Type of the right-hand side dense matrix
95  , bool SO > // Storage order
96 class DMatDMatSubExpr : public DenseMatrix< DMatDMatSubExpr<MT1,MT2,SO>, SO >
97  , private MatMatSubExpr
98  , private Computation
99 {
100  private:
101  //**Type definitions****************************************************************************
102  typedef typename MT1::ResultType RT1;
103  typedef typename MT2::ResultType RT2;
104  typedef typename MT1::ReturnType RN1;
105  typedef typename MT2::ReturnType RN2;
106  typedef typename MT1::CompositeType CT1;
107  typedef typename MT2::CompositeType CT2;
108  typedef typename MT1::ElementType ET1;
109  typedef typename MT2::ElementType ET2;
110  //**********************************************************************************************
111 
112  //**Return type evaluation**********************************************************************
114 
119  enum { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
120 
123  //**********************************************************************************************
124 
125  //**Serial evaluation strategy******************************************************************
127 
133  enum { useAssign = RequiresEvaluation<MT1>::value || RequiresEvaluation<MT2>::value || !returnExpr };
134 
136  template< typename MT >
138  struct UseAssign {
139  enum { value = useAssign };
140  };
142  //**********************************************************************************************
143 
144  //**Parallel evaluation strategy****************************************************************
146 
152  template< typename MT >
153  struct UseSMPAssign {
154  enum { value = ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign };
155  };
157  //**********************************************************************************************
158 
159  public:
160  //**Type definitions****************************************************************************
167 
170 
173 
175  typedef typename SelectType< IsExpression<MT1>::value, const MT1, const MT1& >::Type LeftOperand;
176 
178  typedef typename SelectType< IsExpression<MT2>::value, const MT2, const MT2& >::Type RightOperand;
179  //**********************************************************************************************
180 
181  //**ConstIterator class definition**************************************************************
185  {
186  public:
187  //**Type definitions*************************************************************************
188  typedef std::random_access_iterator_tag IteratorCategory;
193 
194  // STL iterator requirements
200 
203 
206  //*******************************************************************************************
207 
208  //**Constructor******************************************************************************
214  explicit inline ConstIterator( LeftIteratorType left, RightIteratorType right )
215  : left_ ( left ) // Iterator to the current left-hand side element
216  , right_( right ) // Iterator to the current right-hand side element
217  {}
218  //*******************************************************************************************
219 
220  //**Addition assignment operator*************************************************************
226  inline ConstIterator& operator+=( size_t inc ) {
227  left_ += inc;
228  right_ += inc;
229  return *this;
230  }
231  //*******************************************************************************************
232 
233  //**Subtraction assignment operator**********************************************************
239  inline ConstIterator& operator-=( size_t dec ) {
240  left_ -= dec;
241  right_ -= dec;
242  return *this;
243  }
244  //*******************************************************************************************
245 
246  //**Prefix increment operator****************************************************************
252  ++left_;
253  ++right_;
254  return *this;
255  }
256  //*******************************************************************************************
257 
258  //**Postfix increment operator***************************************************************
263  inline const ConstIterator operator++( int ) {
264  return ConstIterator( left_++, right_++ );
265  }
266  //*******************************************************************************************
267 
268  //**Prefix decrement operator****************************************************************
274  --left_;
275  --right_;
276  return *this;
277  }
278  //*******************************************************************************************
279 
280  //**Postfix decrement operator***************************************************************
285  inline const ConstIterator operator--( int ) {
286  return ConstIterator( left_--, right_-- );
287  }
288  //*******************************************************************************************
289 
290  //**Element access operator******************************************************************
295  inline ReturnType operator*() const {
296  return (*left_) - (*right_);
297  }
298  //*******************************************************************************************
299 
300  //**Load function****************************************************************************
305  inline IntrinsicType load() const {
306  return left_.load() - right_.load();
307  }
308  //*******************************************************************************************
309 
310  //**Equality operator************************************************************************
316  inline bool operator==( const ConstIterator& rhs ) const {
317  return left_ == rhs.left_;
318  }
319  //*******************************************************************************************
320 
321  //**Inequality operator**********************************************************************
327  inline bool operator!=( const ConstIterator& rhs ) const {
328  return left_ != rhs.left_;
329  }
330  //*******************************************************************************************
331 
332  //**Less-than operator***********************************************************************
338  inline bool operator<( const ConstIterator& rhs ) const {
339  return left_ < rhs.left_;
340  }
341  //*******************************************************************************************
342 
343  //**Greater-than operator********************************************************************
349  inline bool operator>( const ConstIterator& rhs ) const {
350  return left_ > rhs.left_;
351  }
352  //*******************************************************************************************
353 
354  //**Less-or-equal-than operator**************************************************************
360  inline bool operator<=( const ConstIterator& rhs ) const {
361  return left_ <= rhs.left_;
362  }
363  //*******************************************************************************************
364 
365  //**Greater-or-equal-than operator***********************************************************
371  inline bool operator>=( const ConstIterator& rhs ) const {
372  return left_ >= rhs.left_;
373  }
374  //*******************************************************************************************
375 
376  //**Subtraction operator*********************************************************************
382  inline DifferenceType operator-( const ConstIterator& rhs ) const {
383  return left_ - rhs.left_;
384  }
385  //*******************************************************************************************
386 
387  //**Addition operator************************************************************************
394  friend inline const ConstIterator operator+( const ConstIterator& it, size_t inc ) {
395  return ConstIterator( it.left_ + inc, it.right_ + inc );
396  }
397  //*******************************************************************************************
398 
399  //**Addition operator************************************************************************
406  friend inline const ConstIterator operator+( size_t inc, const ConstIterator& it ) {
407  return ConstIterator( it.left_ + inc, it.right_ + inc );
408  }
409  //*******************************************************************************************
410 
411  //**Subtraction operator*********************************************************************
418  friend inline const ConstIterator operator-( const ConstIterator& it, size_t dec ) {
419  return ConstIterator( it.left_ - dec, it.right_ - dec );
420  }
421  //*******************************************************************************************
422 
423  private:
424  //**Member variables*************************************************************************
427  //*******************************************************************************************
428  };
429  //**********************************************************************************************
430 
431  //**Compilation flags***************************************************************************
433  enum { vectorizable = MT1::vectorizable && MT2::vectorizable &&
436 
438  enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
439  //**********************************************************************************************
440 
441  //**Constructor*********************************************************************************
447  explicit inline DMatDMatSubExpr( const MT1& lhs, const MT2& rhs )
448  : lhs_( lhs ) // Left-hand side dense matrix of the subtraction expression
449  , rhs_( rhs ) // Right-hand side dense matrix of the subtraction expression
450  {
451  BLAZE_INTERNAL_ASSERT( lhs.rows() == rhs.rows() , "Invalid number of rows" );
452  BLAZE_INTERNAL_ASSERT( lhs.columns() == rhs.columns(), "Invalid number of columns" );
453  }
454  //**********************************************************************************************
455 
456  //**Access operator*****************************************************************************
463  inline ReturnType operator()( size_t i, size_t j ) const {
464  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
465  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
466  return lhs_(i,j) - rhs_(i,j);
467  }
468  //**********************************************************************************************
469 
470  //**Load function*******************************************************************************
477  inline IntrinsicType load( size_t i, size_t j ) const {
478  typedef IntrinsicTrait<ElementType> IT;
479  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
480  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
481  BLAZE_INTERNAL_ASSERT( !SO || ( i % IT::size == 0UL ), "Invalid row access index" );
482  BLAZE_INTERNAL_ASSERT( SO || ( j % IT::size == 0UL ), "Invalid column access index" );
483  const IntrinsicType xmm1( lhs_.load(i,j) );
484  const IntrinsicType xmm2( rhs_.load(i,j) );
485  return xmm1 - xmm2;
486  }
487  //**********************************************************************************************
488 
489  //**Begin function******************************************************************************
495  inline ConstIterator begin( size_t i ) const {
496  return ConstIterator( lhs_.begin(i), rhs_.begin(i) );
497  }
498  //**********************************************************************************************
499 
500  //**End function********************************************************************************
506  inline ConstIterator end( size_t i ) const {
507  return ConstIterator( lhs_.end(i), rhs_.end(i) );
508  }
509  //**********************************************************************************************
510 
511  //**Rows function*******************************************************************************
516  inline size_t rows() const {
517  return lhs_.rows();
518  }
519  //**********************************************************************************************
520 
521  //**Columns function****************************************************************************
526  inline size_t columns() const {
527  return lhs_.columns();
528  }
529  //**********************************************************************************************
530 
531  //**Left operand access*************************************************************************
536  inline LeftOperand leftOperand() const {
537  return lhs_;
538  }
539  //**********************************************************************************************
540 
541  //**Right operand access************************************************************************
546  inline RightOperand rightOperand() const {
547  return rhs_;
548  }
549  //**********************************************************************************************
550 
551  //**********************************************************************************************
557  template< typename T >
558  inline bool canAlias( const T* alias ) const {
559  return ( IsExpression<MT1>::value && ( RequiresEvaluation<MT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
560  ( IsExpression<MT2>::value && ( RequiresEvaluation<MT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
561  }
562  //**********************************************************************************************
563 
564  //**********************************************************************************************
570  template< typename T >
571  inline bool isAliased( const T* alias ) const {
572  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
573  }
574  //**********************************************************************************************
575 
576  //**********************************************************************************************
581  inline bool isAligned() const {
582  return lhs_.isAligned() && rhs_.isAligned();
583  }
584  //**********************************************************************************************
585 
586  //**********************************************************************************************
591  inline bool canSMPAssign() const {
592  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
593  ( ( ( SO == rowMajor ) ? rows() : columns() ) > SMP_DMATDMATSUB_THRESHOLD );
594  }
595  //**********************************************************************************************
596 
597  private:
598  //**Member variables****************************************************************************
601  //**********************************************************************************************
602 
603  //**Assignment to dense matrices****************************************************************
617  template< typename MT // Type of the target dense matrix
618  , bool SO2 > // Storage order of the target dense matrix
619  friend inline typename EnableIf< UseAssign<MT> >::Type
620  assign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
621  {
623 
624  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
625  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
626 
627  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
628  subAssign( ~lhs, rhs.rhs_ );
629  }
630  else {
631  assign ( ~lhs, rhs.lhs_ );
632  subAssign( ~lhs, rhs.rhs_ );
633  }
634  }
636  //**********************************************************************************************
637 
638  //**Assignment to sparse matrices***************************************************************
652  template< typename MT // Type of the target sparse matrix
653  , bool SO2 > // Storage order of the target sparse matrix
654  friend inline typename EnableIf< UseAssign<MT> >::Type
655  assign( SparseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
656  {
658 
660 
667 
668  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
669  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
670 
671  const TmpType tmp( serial( rhs ) );
672  assign( ~lhs, tmp );
673  }
675  //**********************************************************************************************
676 
677  //**Addition assignment to dense matrices*******************************************************
691  template< typename MT // Type of the target dense matrix
692  , bool SO2 > // Storage order of the target dense matrix
693  friend inline typename EnableIf< UseAssign<MT> >::Type
694  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
695  {
697 
698  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
699  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
700 
701  addAssign( ~lhs, rhs.lhs_ );
702  subAssign( ~lhs, rhs.rhs_ );
703  }
705  //**********************************************************************************************
706 
707  //**Addition assignment to sparse matrices******************************************************
708  // No special implementation for the addition assignment to sparse matrices.
709  //**********************************************************************************************
710 
711  //**Subtraction assignment to dense matrices****************************************************
725  template< typename MT // Type of the target dense matrix
726  , bool SO2 > // Storage order of the target dense matrix
727  friend inline typename EnableIf< UseAssign<MT> >::Type
728  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
729  {
731 
732  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
733  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
734 
735  subAssign( ~lhs, rhs.lhs_ );
736  addAssign( ~lhs, rhs.rhs_ );
737  }
739  //**********************************************************************************************
740 
741  //**Subtraction assignment to sparse matrices***************************************************
742  // No special implementation for the subtraction assignment to sparse matrices.
743  //**********************************************************************************************
744 
745  //**Multiplication assignment to dense matrices*************************************************
746  // No special implementation for the multiplication assignment to dense matrices.
747  //**********************************************************************************************
748 
749  //**Multiplication assignment to sparse matrices************************************************
750  // No special implementation for the multiplication assignment to sparse matrices.
751  //**********************************************************************************************
752 
753  //**SMP assignment to dense matrices************************************************************
767  template< typename MT // Type of the target dense matrix
768  , bool SO2 > // Storage order of the target dense matrix
769  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
770  smpAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
771  {
773 
774  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
775  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
776 
777  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
778  smpSubAssign( ~lhs, rhs.rhs_ );
779  }
780  else {
781  smpAssign ( ~lhs, rhs.lhs_ );
782  smpSubAssign( ~lhs, rhs.rhs_ );
783  }
784  }
786  //**********************************************************************************************
787 
788  //**SMP assignment to sparse matrices***********************************************************
802  template< typename MT // Type of the target sparse matrix
803  , bool SO2 > // Storage order of the target sparse matrix
804  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
805  smpAssign( SparseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
806  {
808 
809  typedef typename SelectType< SO == SO2, ResultType, OppositeType >::Type TmpType;
810 
817 
818  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
819  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
820 
821  const TmpType tmp( rhs );
822  smpAssign( ~lhs, tmp );
823  }
825  //**********************************************************************************************
826 
827  //**SMP addition assignment to dense matrices***************************************************
841  template< typename MT // Type of the target dense matrix
842  , bool SO2 > // Storage order of the target dense matrix
843  friend inline typename EnableIf< UseAssign<MT> >::Type
844  smpAddAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
845  {
847 
848  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
849  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
850 
851  smpAddAssign( ~lhs, rhs.lhs_ );
852  smpSubAssign( ~lhs, rhs.rhs_ );
853  }
855  //**********************************************************************************************
856 
857  //**SMP addition assignment to sparse matrices**************************************************
858  // No special implementation for the SMP addition assignment to sparse matrices.
859  //**********************************************************************************************
860 
861  //**SMP subtraction assignment to dense matrices************************************************
875  template< typename MT // Type of the target dense matrix
876  , bool SO2 > // Storage order of the target dense matrix
877  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
878  smpSubAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSubExpr& rhs )
879  {
881 
882  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
883  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
884 
885  smpSubAssign( ~lhs, rhs.lhs_ );
886  smpAddAssign( ~lhs, rhs.rhs_ );
887  }
889  //**********************************************************************************************
890 
891  //**SMP subtraction assignment to sparse matrices***********************************************
892  // No special implementation for the SMP subtraction assignment to sparse matrices.
893  //**********************************************************************************************
894 
895  //**SMP multiplication assignment to dense matrices*********************************************
896  // No special implementation for the SMP multiplication assignment to dense matrices.
897  //**********************************************************************************************
898 
899  //**SMP multiplication assignment to sparse matrices********************************************
900  // No special implementation for the SMP multiplication assignment to sparse matrices.
901  //**********************************************************************************************
902 
903  //**Compile time checks*************************************************************************
909  //**********************************************************************************************
910 };
911 //*************************************************************************************************
912 
913 
914 
915 
916 //=================================================================================================
917 //
918 // GLOBAL BINARY ARITHMETIC OPERATORS
919 //
920 //=================================================================================================
921 
922 //*************************************************************************************************
947 template< typename T1 // Type of the left-hand side dense matrix
948  , typename T2 // Type of the right-hand side dense matrix
949  , bool SO > // Storage order
950 inline const DMatDMatSubExpr<T1,T2,SO>
952 {
954 
955  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
956  throw std::invalid_argument( "Matrix sizes do not match" );
957 
958  return DMatDMatSubExpr<T1,T2,SO>( ~lhs, ~rhs );
959 }
960 //*************************************************************************************************
961 
962 
963 
964 
965 //=================================================================================================
966 //
967 // ROWS SPECIALIZATIONS
968 //
969 //=================================================================================================
970 
971 //*************************************************************************************************
973 template< typename MT1, typename MT2, bool SO >
974 struct Rows< DMatDMatSubExpr<MT1,MT2,SO> >
975  : public Max< Rows<MT1>, Rows<MT2> >::Type
976 {};
978 //*************************************************************************************************
979 
980 
981 
982 
983 //=================================================================================================
984 //
985 // COLUMNS SPECIALIZATIONS
986 //
987 //=================================================================================================
988 
989 //*************************************************************************************************
991 template< typename MT1, typename MT2, bool SO >
992 struct Columns< DMatDMatSubExpr<MT1,MT2,SO> >
993  : public Max< Columns<MT1>, Columns<MT2> >::Type
994 {};
996 //*************************************************************************************************
997 
998 
999 
1000 
1001 //=================================================================================================
1002 //
1003 // ISSYMMETRIC SPECIALIZATIONS
1004 //
1005 //=================================================================================================
1006 
1007 //*************************************************************************************************
1009 template< typename MT1, typename MT2, bool SO >
1010 struct IsSymmetric< DMatDMatSubExpr<MT1,MT2,SO> >
1011  : public IsTrue< IsSymmetric<MT1>::value && IsSymmetric<MT2>::value >
1012 {};
1014 //*************************************************************************************************
1015 
1016 
1017 
1018 
1019 //=================================================================================================
1020 //
1021 // ISLOWER SPECIALIZATIONS
1022 //
1023 //=================================================================================================
1024 
1025 //*************************************************************************************************
1027 template< typename MT1, typename MT2, bool SO >
1028 struct IsLower< DMatDMatSubExpr<MT1,MT2,SO> >
1029  : public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1030 {};
1032 //*************************************************************************************************
1033 
1034 
1035 
1036 
1037 //=================================================================================================
1038 //
1039 // ISUPPER SPECIALIZATIONS
1040 //
1041 //=================================================================================================
1042 
1043 //*************************************************************************************************
1045 template< typename MT1, typename MT2, bool SO >
1046 struct IsUpper< DMatDMatSubExpr<MT1,MT2,SO> >
1047  : public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1048 {};
1050 //*************************************************************************************************
1051 
1052 
1053 
1054 
1055 //=================================================================================================
1056 //
1057 // EXPRESSION TRAIT SPECIALIZATIONS
1058 //
1059 //=================================================================================================
1060 
1061 //*************************************************************************************************
1063 template< typename MT1, typename MT2, bool SO, bool AF >
1064 struct SubmatrixExprTrait< DMatDMatSubExpr<MT1,MT2,SO>, AF >
1065 {
1066  public:
1067  //**********************************************************************************************
1068  typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1069  , typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1070  //**********************************************************************************************
1071 };
1073 //*************************************************************************************************
1074 
1075 
1076 //*************************************************************************************************
1078 template< typename MT1, typename MT2, bool SO >
1079 struct RowExprTrait< DMatDMatSubExpr<MT1,MT2,SO> >
1080 {
1081  public:
1082  //**********************************************************************************************
1083  typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
1084  , typename RowExprTrait<const MT2>::Type >::Type Type;
1085  //**********************************************************************************************
1086 };
1088 //*************************************************************************************************
1089 
1090 
1091 //*************************************************************************************************
1093 template< typename MT1, typename MT2, bool SO >
1094 struct ColumnExprTrait< DMatDMatSubExpr<MT1,MT2,SO> >
1095 {
1096  public:
1097  //**********************************************************************************************
1098  typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
1099  , typename ColumnExprTrait<const MT2>::Type >::Type Type;
1100  //**********************************************************************************************
1101 };
1103 //*************************************************************************************************
1104 
1105 } // namespace blaze
1106 
1107 #endif
Constraint on the data type.
Pointer difference type of the Blaze library.
Header file for the Max class template.
const ConstIterator operator++(int)
Post-increment operator.
Definition: DMatDMatSubExpr.h:263
ReturnType operator*() const
Direct access to the element at the current iterator position.
Definition: DMatDMatSubExpr.h:295
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.
ConstIterator & operator++()
Pre-increment operator.
Definition: DMatDMatSubExpr.h:251
Header file for the subtraction trait.
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:258
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatDMatSubExpr.h:558
bool operator>(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:349
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DMatDMatSubExpr.h:169
ElementType * PointerType
Pointer return type.
Definition: DMatDMatSubExpr.h:190
IntrinsicTrait< ElementType >::Type IntrinsicType
Resulting intrinsic element type.
Definition: DMatDMatSubExpr.h:166
#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:205
DMatDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatDMatSubExpr class.
Definition: DMatDMatSubExpr.h:447
#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:463
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.
Header file for the IsSame and IsStrictlySame type traits.
IntrinsicType load() const
Access to the intrinsic elements of the matrix.
Definition: DMatDMatSubExpr.h:305
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
MT1::ReturnType RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:104
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
ConstIterator & operator--()
Pre-decrement operator.
Definition: DMatDMatSubExpr.h:273
PointerType pointer
Pointer return type.
Definition: DMatDMatSubExpr.h:197
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
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:205
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:184
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:107
Constraint on the data type.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatDMatSubExpr.h:165
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:106
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
DifferenceType difference_type
Difference between two iterators.
Definition: DMatDMatSubExpr.h:199
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatDMatSubExpr.h:600
Expression object for dense matrix-dense matrix subtractions.The DMatDMatSubExpr class represents the...
Definition: DMatDMatSubExpr.h:96
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:516
Header file for the IsTemporary type trait class.
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatDMatSubExpr.h:162
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
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row i.
Definition: DMatDMatSubExpr.h:495
ElementType ValueType
Type of the underlying elements.
Definition: DMatDMatSubExpr.h:189
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
const size_t SMP_DMATDMATSUB_THRESHOLD
SMP row-major dense matrix/row-major dense matrix subtraction threshold.This threshold specifies when...
Definition: Thresholds.h:763
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.
DMatDMatSubExpr< MT1, MT2, SO > This
Type of this DMatDMatSubExpr instance.
Definition: DMatDMatSubExpr.h:161
std::random_access_iterator_tag IteratorCategory
The iterator category.
Definition: DMatDMatSubExpr.h:188
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:175
RightIteratorType right_
Iterator to the current right-hand side element.
Definition: DMatDMatSubExpr.h:426
Header file for the MatMatSubExpr base class.
ReferenceType reference
Reference return type.
Definition: DMatDMatSubExpr.h:198
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatDMatSubExpr.h:591
Header file for the IsLower type trait.
ConstIterator(LeftIteratorType left, RightIteratorType right)
Constructor for the ConstIterator class.
Definition: DMatDMatSubExpr.h:214
#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:581
Constraints on the storage order of matrix types.
MT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:108
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row i.
Definition: DMatDMatSubExpr.h:506
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:104
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatDMatSubExpr.h:164
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:178
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.
bool operator<=(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:360
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatDMatSubExpr.h:163
ValueType value_type
Type of the underlying elements.
Definition: DMatDMatSubExpr.h:196
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:327
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two iterators.
Definition: DMatDMatSubExpr.h:382
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.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2477
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:107
Intrinsic characteristics of data types.The IntrinsicTrait class template provides the intrinsic char...
Definition: IntrinsicTrait.h:749
Header file for run time assertion macros.
MT2::ElementType ET2
Element type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:109
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
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:316
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:406
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:105
ConstIterator & operator-=(size_t dec)
Subtraction assignment operator.
Definition: DMatDMatSubExpr.h:239
friend const ConstIterator operator-(const ConstIterator &it, size_t dec)
Subtraction between a ConstIterator and an integral value.
Definition: DMatDMatSubExpr.h:418
#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:371
RightOperand rightOperand() const
Returns the right-hand side dense matrix operand.
Definition: DMatDMatSubExpr.h:546
Header file for all intrinsic functionality.
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DMatDMatSubExpr.h:192
IntrinsicType load(size_t i, size_t j) const
Access to the intrinsic elements of the matrix.
Definition: DMatDMatSubExpr.h:477
ElementType & ReferenceType
Reference return type.
Definition: DMatDMatSubExpr.h:191
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:103
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:172
IteratorCategory iterator_category
The iterator category.
Definition: DMatDMatSubExpr.h:195
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
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatDMatSubExpr.h:536
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatDMatSubExpr.h:122
#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:425
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2473
Header file for the IsTrue value trait.
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatDMatSubExpr.h:571
Header file for basic type definitions.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:102
bool operator<(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatSubExpr.h:338
MT1::ConstIterator LeftIteratorType
ConstIterator type of the left-hand side dense matrix expression.
Definition: DMatDMatSubExpr.h:202
#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
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatDMatSubExpr.h:599
Header file for the IsUpper type trait.
friend const ConstIterator operator+(const ConstIterator &it, size_t inc)
Addition between a ConstIterator and an integral value.
Definition: DMatDMatSubExpr.h:394
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatDMatSubExpr.h:526
const ConstIterator operator--(int)
Post-decrement operator.
Definition: DMatDMatSubExpr.h:285
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
ConstIterator & operator+=(size_t inc)
Addition assignment operator.
Definition: DMatDMatSubExpr.h:226