All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DMatDMatAddExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATDMATADDEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATDMATADDEXPR_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 DMATDMATADDEXPR
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 DMatDMatAddExpr : public DenseMatrix< DMatDMatAddExpr<MT1,MT2,SO>, SO >
97  , private MatMatAddExpr
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 DMatDMatAddExpr( const MT1& lhs, const MT2& rhs )
448  : lhs_( lhs ) // Left-hand side dense matrix of the addition expression
449  , rhs_( rhs ) // Right-hand side dense matrix of the addition 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  return lhs_.load(i,j) + rhs_.load(i,j);
484  }
485  //**********************************************************************************************
486 
487  //**Begin function******************************************************************************
493  inline ConstIterator begin( size_t i ) const {
494  return ConstIterator( lhs_.begin(i), rhs_.begin(i) );
495  }
496  //**********************************************************************************************
497 
498  //**End function********************************************************************************
504  inline ConstIterator end( size_t i ) const {
505  return ConstIterator( lhs_.end(i), rhs_.end(i) );
506  }
507  //**********************************************************************************************
508 
509  //**Rows function*******************************************************************************
514  inline size_t rows() const {
515  return lhs_.rows();
516  }
517  //**********************************************************************************************
518 
519  //**Columns function****************************************************************************
524  inline size_t columns() const {
525  return lhs_.columns();
526  }
527  //**********************************************************************************************
528 
529  //**Left operand access*************************************************************************
534  inline LeftOperand leftOperand() const {
535  return lhs_;
536  }
537  //**********************************************************************************************
538 
539  //**Right operand access************************************************************************
544  inline RightOperand rightOperand() const {
545  return rhs_;
546  }
547  //**********************************************************************************************
548 
549  //**********************************************************************************************
555  template< typename T >
556  inline bool canAlias( const T* alias ) const {
557  return ( IsExpression<MT1>::value && ( RequiresEvaluation<MT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
558  ( IsExpression<MT2>::value && ( RequiresEvaluation<MT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
559  }
560  //**********************************************************************************************
561 
562  //**********************************************************************************************
568  template< typename T >
569  inline bool isAliased( const T* alias ) const {
570  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
571  }
572  //**********************************************************************************************
573 
574  //**********************************************************************************************
579  inline bool isAligned() const {
580  return lhs_.isAligned() && rhs_.isAligned();
581  }
582  //**********************************************************************************************
583 
584  //**********************************************************************************************
589  inline bool canSMPAssign() const {
590  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
591  ( ( ( SO == rowMajor ) ? rows() : columns() ) > SMP_DMATDMATADD_THRESHOLD );
592  }
593  //**********************************************************************************************
594 
595  private:
596  //**Member variables****************************************************************************
599  //**********************************************************************************************
600 
601  //**Assignment to dense matrices****************************************************************
615  template< typename MT // Type of the target dense matrix
616  , bool SO2 > // Storage order of the target dense matrix
617  friend inline typename EnableIf< UseAssign<MT> >::Type
618  assign( DenseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
619  {
621 
622  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
623  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
624 
625  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
626  addAssign( ~lhs, rhs.rhs_ );
627  }
628  else if( !IsExpression<MT2>::value && isSame( ~lhs, rhs.rhs_ ) ) {
629  addAssign( ~lhs, rhs.lhs_ );
630  }
631  else {
632  assign ( ~lhs, rhs.lhs_ );
633  addAssign( ~lhs, rhs.rhs_ );
634  }
635  }
637  //**********************************************************************************************
638 
639  //**Assignment to sparse matrices***************************************************************
653  template< typename MT // Type of the target sparse matrix
654  , bool SO2 > // Storage order of the target sparse matrix
655  friend inline typename EnableIf< UseAssign<MT> >::Type
656  assign( SparseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
657  {
659 
661 
668 
669  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
670  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
671 
672  const TmpType tmp( serial( rhs ) );
673  assign( ~lhs, tmp );
674  }
676  //**********************************************************************************************
677 
678  //**Addition assignment to dense matrices*******************************************************
692  template< typename MT // Type of the target dense matrix
693  , bool SO2 > // Storage order of the target dense matrix
694  friend inline typename EnableIf< UseAssign<MT> >::Type
695  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
696  {
698 
699  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
700  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
701 
702  addAssign( ~lhs, rhs.lhs_ );
703  addAssign( ~lhs, rhs.rhs_ );
704  }
706  //**********************************************************************************************
707 
708  //**Addition assignment to sparse matrices******************************************************
709  // No special implementation for the addition assignment to sparse matrices.
710  //**********************************************************************************************
711 
712  //**Subtraction assignment to dense matrices****************************************************
726  template< typename MT // Type of the target dense matrix
727  , bool SO2 > // Storage order of the target dense matrix
728  friend inline typename EnableIf< UseAssign<MT> >::Type
729  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
730  {
732 
733  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
734  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
735 
736  subAssign( ~lhs, rhs.lhs_ );
737  subAssign( ~lhs, rhs.rhs_ );
738  }
740  //**********************************************************************************************
741 
742  //**Subtraction assignment to sparse matrices***************************************************
743  // No special implementation for the subtraction assignment to sparse matrices.
744  //**********************************************************************************************
745 
746  //**Multiplication assignment to dense matrices*************************************************
747  // No special implementation for the multiplication assignment to dense matrices.
748  //**********************************************************************************************
749 
750  //**Multiplication assignment to sparse matrices************************************************
751  // No special implementation for the multiplication assignment to sparse matrices.
752  //**********************************************************************************************
753 
754  //**SMP assignment to dense matrices************************************************************
768  template< typename MT // Type of the target dense matrix
769  , bool SO2 > // Storage order of the target dense matrix
770  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
771  smpAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
772  {
774 
775  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
776  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
777 
778  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
779  smpAddAssign( ~lhs, rhs.rhs_ );
780  }
781  else if( !IsExpression<MT2>::value && isSame( ~lhs, rhs.rhs_ ) ) {
782  smpAddAssign( ~lhs, rhs.lhs_ );
783  }
784  else {
785  smpAssign ( ~lhs, rhs.lhs_ );
786  smpAddAssign( ~lhs, rhs.rhs_ );
787  }
788  }
790  //**********************************************************************************************
791 
792  //**SMP assignment to sparse matrices***********************************************************
806  template< typename MT // Type of the target sparse matrix
807  , bool SO2 > // Storage order of the target sparse matrix
808  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
809  smpAssign( SparseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
810  {
812 
813  typedef typename SelectType< SO == SO2, ResultType, OppositeType >::Type TmpType;
814 
821 
822  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
823  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
824 
825  const TmpType tmp( rhs );
826  smpAssign( ~lhs, tmp );
827  }
829  //**********************************************************************************************
830 
831  //**SMP addition assignment to dense matrices***************************************************
845  template< typename MT // Type of the target dense matrix
846  , bool SO2 > // Storage order of the target dense matrix
847  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
848  smpAddAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
849  {
851 
852  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
853  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
854 
855  smpAddAssign( ~lhs, rhs.lhs_ );
856  smpAddAssign( ~lhs, rhs.rhs_ );
857  }
859  //**********************************************************************************************
860 
861  //**SMP addition assignment to sparse matrices**************************************************
862  // No special implementation for the SMP addition assignment to sparse matrices.
863  //**********************************************************************************************
864 
865  //**SMP subtraction assignment to dense matrices************************************************
879  template< typename MT // Type of the target dense matrix
880  , bool SO2 > // Storage order of the target dense matrix
881  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
882  smpSubAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatAddExpr& rhs )
883  {
885 
886  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
887  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
888 
889  smpSubAssign( ~lhs, rhs.lhs_ );
890  smpSubAssign( ~lhs, rhs.rhs_ );
891  }
893  //**********************************************************************************************
894 
895  //**SMP subtraction assignment to sparse matrices***********************************************
896  // No special implementation for the SMP subtraction assignment to sparse matrices.
897  //**********************************************************************************************
898 
899  //**SMP multiplication assignment to dense matrices*********************************************
900  // No special implementation for the SMP multiplication assignment to dense matrices.
901  //**********************************************************************************************
902 
903  //**SMP multiplication assignment to sparse matrices********************************************
904  // No special implementation for the SMP multiplication assignment to sparse matrices.
905  //**********************************************************************************************
906 
907  //**Compile time checks*************************************************************************
914  //**********************************************************************************************
915 };
916 //*************************************************************************************************
917 
918 
919 
920 
921 //=================================================================================================
922 //
923 // GLOBAL BINARY ARITHMETIC OPERATORS
924 //
925 //=================================================================================================
926 
927 //*************************************************************************************************
954 template< typename T1 // Type of the left-hand side dense matrix
955  , typename T2 // Type of the right-hand side dense matrix
956  , bool SO > // Storage order
957 inline const DMatDMatAddExpr<T1,T2,SO>
959 {
961 
962  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
963  throw std::invalid_argument( "Matrix sizes do not match" );
964 
965  return DMatDMatAddExpr<T1,T2,SO>( ~lhs, ~rhs );
966 }
967 //*************************************************************************************************
968 
969 
970 
971 
972 //=================================================================================================
973 //
974 // ROWS SPECIALIZATIONS
975 //
976 //=================================================================================================
977 
978 //*************************************************************************************************
980 template< typename MT1, typename MT2, bool SO >
981 struct Rows< DMatDMatAddExpr<MT1,MT2,SO> >
982  : public Max< Rows<MT1>, Rows<MT2> >::Type
983 {};
985 //*************************************************************************************************
986 
987 
988 
989 
990 //=================================================================================================
991 //
992 // COLUMNS SPECIALIZATIONS
993 //
994 //=================================================================================================
995 
996 //*************************************************************************************************
998 template< typename MT1, typename MT2, bool SO >
999 struct Columns< DMatDMatAddExpr<MT1,MT2,SO> >
1000  : public Max< Columns<MT1>, Columns<MT2> >::Type
1001 {};
1003 //*************************************************************************************************
1004 
1005 
1006 
1007 
1008 //=================================================================================================
1009 //
1010 // ISSYMMETRIC SPECIALIZATIONS
1011 //
1012 //=================================================================================================
1013 
1014 //*************************************************************************************************
1016 template< typename MT1, typename MT2, bool SO >
1017 struct IsSymmetric< DMatDMatAddExpr<MT1,MT2,SO> >
1018  : public IsTrue< IsSymmetric<MT1>::value && IsSymmetric<MT2>::value >
1019 {};
1021 //*************************************************************************************************
1022 
1023 
1024 
1025 
1026 //=================================================================================================
1027 //
1028 // ISLOWER SPECIALIZATIONS
1029 //
1030 //=================================================================================================
1031 
1032 //*************************************************************************************************
1034 template< typename MT1, typename MT2, bool SO >
1035 struct IsLower< DMatDMatAddExpr<MT1,MT2,SO> >
1036  : public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1037 {};
1039 //*************************************************************************************************
1040 
1041 
1042 
1043 
1044 //=================================================================================================
1045 //
1046 // ISUPPER SPECIALIZATIONS
1047 //
1048 //=================================================================================================
1049 
1050 //*************************************************************************************************
1052 template< typename MT1, typename MT2, bool SO >
1053 struct IsUpper< DMatDMatAddExpr<MT1,MT2,SO> >
1054  : public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1055 {};
1057 //*************************************************************************************************
1058 
1059 
1060 
1061 
1062 //=================================================================================================
1063 //
1064 // EXPRESSION TRAIT SPECIALIZATIONS
1065 //
1066 //=================================================================================================
1067 
1068 //*************************************************************************************************
1070 template< typename MT1, typename MT2, bool SO, bool AF >
1071 struct SubmatrixExprTrait< DMatDMatAddExpr<MT1,MT2,SO>, AF >
1072 {
1073  public:
1074  //**********************************************************************************************
1075  typedef typename AddExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1076  , typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1077  //**********************************************************************************************
1078 };
1080 //*************************************************************************************************
1081 
1082 
1083 //*************************************************************************************************
1085 template< typename MT1, typename MT2, bool SO >
1086 struct RowExprTrait< DMatDMatAddExpr<MT1,MT2,SO> >
1087 {
1088  public:
1089  //**********************************************************************************************
1090  typedef typename AddExprTrait< typename RowExprTrait<const MT1>::Type
1091  , typename RowExprTrait<const MT2>::Type >::Type Type;
1092  //**********************************************************************************************
1093 };
1095 //*************************************************************************************************
1096 
1097 
1098 //*************************************************************************************************
1100 template< typename MT1, typename MT2, bool SO >
1101 struct ColumnExprTrait< DMatDMatAddExpr<MT1,MT2,SO> >
1102 {
1103  public:
1104  //**********************************************************************************************
1105  typedef typename AddExprTrait< typename ColumnExprTrait<const MT1>::Type
1106  , typename ColumnExprTrait<const MT2>::Type >::Type Type;
1107  //**********************************************************************************************
1108 };
1110 //*************************************************************************************************
1111 
1112 } // namespace blaze
1113 
1114 #endif
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row i.
Definition: DMatDMatAddExpr.h:504
Pointer difference type of the Blaze library.
Header file for the Max class template.
bool operator>(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatAddExpr.h:349
ConstIterator & operator++()
Pre-increment operator.
Definition: DMatDMatAddExpr.h:251
ConstIterator & operator+=(size_t inc)
Addition assignment operator.
Definition: DMatDMatAddExpr.h:226
ConstIterator(LeftIteratorType left, RightIteratorType right)
Constructor for the ConstIterator class.
Definition: DMatDMatAddExpr.h:214
RightOperand rightOperand() const
Returns the right-hand side dense matrix operand.
Definition: DMatDMatAddExpr.h:544
Evaluation of the return type of an addition expression.Via this type trait it is possible to evaluat...
Definition: AddExprTrait.h:104
PointerType pointer
Pointer return type.
Definition: DMatDMatAddExpr.h:197
ConstIterator & operator-=(size_t dec)
Subtraction assignment operator.
Definition: DMatDMatAddExpr.h:239
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DMatDMatAddExpr.h:192
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.
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:258
bool operator<=(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatAddExpr.h:360
MT1::ElementType ET1
Element type of the left-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:108
#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
#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*() const
Direct access to the element at the current iterator position.
Definition: DMatDMatAddExpr.h:295
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(size_t i, size_t j) const
Access to the intrinsic elements of the matrix.
Definition: DMatDMatAddExpr.h:477
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DMatDMatAddExpr.h:327
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:103
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
MT2::ConstIterator RightIteratorType
ConstIterator type of the right-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:205
ElementType * PointerType
Pointer return type.
Definition: DMatDMatAddExpr.h:190
Header file for the AddExprTrait class template.
std::random_access_iterator_tag IteratorCategory
The iterator category.
Definition: DMatDMatAddExpr.h:188
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatDMatAddExpr.h:589
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
Header file for the RequiresEvaluation type trait.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatDMatAddExpr.h:165
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:105
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:107
IntrinsicType load() const
Access to the intrinsic elements of the matrix.
Definition: DMatDMatAddExpr.h:305
Constraint on the data type.
Iterator over the elements of the dense matrix.
Definition: DMatDMatAddExpr.h:184
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatDMatAddExpr.h:164
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
RightIteratorType right_
Iterator to the current right-hand side element.
Definition: DMatDMatAddExpr.h:426
Expression object for dense matrix-dense matrix additions.The DMatDMatAddExpr class represents the co...
Definition: DMatDMatAddExpr.h:96
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:102
Compile time type selection.The SelectType class template selects one of the two given types T1 and T...
Definition: SelectType.h:59
Header file for the IsTemporary type trait class.
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row i.
Definition: DMatDMatAddExpr.h:493
Header file for the IsSymmetric type trait.
const ConstIterator operator--(int)
Post-decrement operator.
Definition: DMatDMatAddExpr.h:285
AddExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatDMatAddExpr.h:122
AddTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatDMatAddExpr.h:162
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
SelectType< useAssign, const ResultType, const DMatDMatAddExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: DMatDMatAddExpr.h:172
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatDMatAddExpr.h:569
ConstIterator & operator--()
Pre-decrement operator.
Definition: DMatDMatAddExpr.h:273
IntrinsicTrait< ElementType >::Type IntrinsicType
Resulting intrinsic element type.
Definition: DMatDMatAddExpr.h:166
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.
friend const ConstIterator operator+(const ConstIterator &it, size_t inc)
Addition between a ConstIterator and an integral value.
Definition: DMatDMatAddExpr.h:394
Base class for all matrix/matrix addition expression templates.The MatMatAddExpr class serves as a ta...
Definition: MatMatAddExpr.h:65
ReferenceType reference
Reference return type.
Definition: DMatDMatAddExpr.h:198
const DenseIterator< Type > operator+(const DenseIterator< Type > &it, ptrdiff_t inc)
Addition between a DenseIterator and an integral value.
Definition: DenseIterator.h:556
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:107
Header file for the IsLower type trait.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatDMatAddExpr.h:163
bool operator<(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatAddExpr.h:338
DMatDMatAddExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatDMatAddExpr class.
Definition: DMatDMatAddExpr.h:447
#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.
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DMatDMatAddExpr.h:316
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_MATMATADDEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid matrix/matrix ...
Definition: MatMatAddExpr.h:165
MT1::ConstIterator LeftIteratorType
ConstIterator type of the left-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:202
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
RightOperand rhs_
Right-hand side dense matrix of the addition expression.
Definition: DMatDMatAddExpr.h:598
Header file for the EnableIf class template.
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two iterators.
Definition: DMatDMatAddExpr.h:382
Header file for the serial shim.
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatDMatAddExpr.h:524
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatDMatAddExpr.h:556
MT2::ElementType ET2
Element type of the right-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:109
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
Intrinsic characteristics of data types.The IntrinsicTrait class template provides the intrinsic char...
Definition: IntrinsicTrait.h:749
Header file for run time assertion macros.
Base template for the AddTrait class.
Definition: AddTrait.h:142
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
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatDMatAddExpr.h:514
Header file for the addition trait.
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
ElementType & ReferenceType
Reference return type.
Definition: DMatDMatAddExpr.h:191
Constraint on the data type.
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
ElementType ValueType
Type of the underlying elements.
Definition: DMatDMatAddExpr.h:189
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DMatDMatAddExpr.h:169
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatDMatAddExpr.h:579
Header file for the MatMatAddExpr base class.
const ConstIterator operator++(int)
Post-increment operator.
Definition: DMatDMatAddExpr.h:263
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:106
#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
friend const ConstIterator operator+(size_t inc, const ConstIterator &it)
Addition between an integral value and a ConstIterator.
Definition: DMatDMatAddExpr.h:406
Header file for all intrinsic functionality.
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:178
LeftIteratorType left_
Iterator to the current left-hand side element.
Definition: DMatDMatAddExpr.h:425
const bool rowMajor
Storage order flag for row-major matrices.
Definition: StorageOrder.h:71
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
IteratorCategory iterator_category
The iterator category.
Definition: DMatDMatAddExpr.h:195
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 lhs_
Left-hand side dense matrix of the addition expression.
Definition: DMatDMatAddExpr.h:597
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatDMatAddExpr.h:534
#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
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2473
Header file for the IsTrue value trait.
Header file for basic type definitions.
friend const ConstIterator operator-(const ConstIterator &it, size_t dec)
Subtraction between a ConstIterator and an integral value.
Definition: DMatDMatAddExpr.h:418
const size_t SMP_DMATDMATADD_THRESHOLD
SMP row-major dense matrix/row-major dense matrix addition threshold.This threshold specifies when a ...
Definition: Thresholds.h:714
bool operator>=(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatAddExpr.h:371
Header file for the IsUpper type trait.
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatDMatAddExpr.h:463
MT1::ReturnType RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:104
Header file for the thresholds for matrix/vector and matrix/matrix multiplications.
DifferenceType difference_type
Difference between two iterators.
Definition: DMatDMatAddExpr.h:199
#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
DMatDMatAddExpr< MT1, MT2, SO > This
Type of this DMatDMatAdd instance.
Definition: DMatDMatAddExpr.h:161
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatAddExpr.h:175
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
ValueType value_type
Type of the underlying elements.
Definition: DMatDMatAddExpr.h:196