All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DVecDVecSubExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_DVECDVECSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DVECDVECSUBEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <stdexcept>
51 #include <blaze/math/Intrinsics.h>
62 #include <blaze/util/Assert.h>
64 #include <blaze/util/EnableIf.h>
66 #include <blaze/util/mpl/Max.h>
67 #include <blaze/util/SelectType.h>
68 #include <blaze/util/Types.h>
70 
71 
72 namespace blaze {
73 
74 //=================================================================================================
75 //
76 // CLASS DVECDVECSUBEXPR
77 //
78 //=================================================================================================
79 
80 //*************************************************************************************************
87 template< typename VT1 // Type of the left-hand side dense vector
88  , typename VT2 // Type of the right-hand side dense vector
89  , bool TF > // Transpose flag
90 class DVecDVecSubExpr : public DenseVector< DVecDVecSubExpr<VT1,VT2,TF>, TF >
91  , private VecVecSubExpr
92  , private Computation
93 {
94  private:
95  //**Type definitions****************************************************************************
96  typedef typename VT1::ResultType RT1;
97  typedef typename VT2::ResultType RT2;
98  typedef typename VT1::ReturnType RN1;
99  typedef typename VT2::ReturnType RN2;
100  typedef typename VT1::CompositeType CT1;
101  typedef typename VT2::CompositeType CT2;
102  typedef typename VT1::ElementType ET1;
103  typedef typename VT2::ElementType ET2;
104  //**********************************************************************************************
105 
106  //**Return type evaluation**********************************************************************
108 
113  enum { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
114 
117  //**********************************************************************************************
118 
119  //**Serial evaluation strategy******************************************************************
121 
127  enum { useAssign = ( RequiresEvaluation<VT1>::value || RequiresEvaluation<VT2>::value || !returnExpr ) };
128 
130  template< typename VT >
132  struct UseAssign {
133  enum { value = useAssign };
134  };
136  //**********************************************************************************************
137 
138  //**Parallel evaluation strategy****************************************************************
140 
146  template< typename VT >
147  struct UseSMPAssign {
148  enum { value = ( !VT1::smpAssignable || !VT2::smpAssignable ) && useAssign };
149  };
151  //**********************************************************************************************
152 
153  public:
154  //**Type definitions****************************************************************************
160 
163 
166 
168  typedef typename SelectType< IsExpression<VT1>::value, const VT1, const VT1& >::Type LeftOperand;
169 
171  typedef typename SelectType< IsExpression<VT2>::value, const VT2, const VT2& >::Type RightOperand;
172  //**********************************************************************************************
173 
174  //**ConstIterator class definition**************************************************************
178  {
179  public:
180  //**Type definitions*************************************************************************
181  typedef std::random_access_iterator_tag IteratorCategory;
186 
187  // STL iterator requirements
193 
196 
199  //*******************************************************************************************
200 
201  //**Constructor******************************************************************************
207  explicit inline ConstIterator( LeftIteratorType left, RightIteratorType right )
208  : left_ ( left ) // Iterator to the current left-hand side element
209  , right_( right ) // Iterator to the current right-hand side element
210  {}
211  //*******************************************************************************************
212 
213  //**Addition assignment operator*************************************************************
219  inline ConstIterator& operator+=( size_t inc ) {
220  left_ += inc;
221  right_ += inc;
222  return *this;
223  }
224  //*******************************************************************************************
225 
226  //**Subtraction assignment operator**********************************************************
232  inline ConstIterator& operator-=( size_t dec ) {
233  left_ -= dec;
234  right_ -= dec;
235  return *this;
236  }
237  //*******************************************************************************************
238 
239  //**Prefix increment operator****************************************************************
245  ++left_;
246  ++right_;
247  return *this;
248  }
249  //*******************************************************************************************
250 
251  //**Postfix increment operator***************************************************************
256  inline const ConstIterator operator++( int ) {
257  return ConstIterator( left_++, right_++ );
258  }
259  //*******************************************************************************************
260 
261  //**Prefix decrement operator****************************************************************
267  --left_;
268  --right_;
269  return *this;
270  }
271  //*******************************************************************************************
272 
273  //**Postfix decrement operator***************************************************************
278  inline const ConstIterator operator--( int ) {
279  return ConstIterator( left_--, right_-- );
280  }
281  //*******************************************************************************************
282 
283  //**Element access operator******************************************************************
288  inline ReturnType operator*() const {
289  return (*left_) - (*right_);
290  }
291  //*******************************************************************************************
292 
293  //**Load function****************************************************************************
298  inline IntrinsicType load() const {
299  return left_.load() - right_.load();
300  }
301  //*******************************************************************************************
302 
303  //**Equality operator************************************************************************
309  inline bool operator==( const ConstIterator& rhs ) const {
310  return left_ == rhs.left_;
311  }
312  //*******************************************************************************************
313 
314  //**Inequality operator**********************************************************************
320  inline bool operator!=( const ConstIterator& rhs ) const {
321  return left_ != rhs.left_;
322  }
323  //*******************************************************************************************
324 
325  //**Less-than operator***********************************************************************
331  inline bool operator<( const ConstIterator& rhs ) const {
332  return left_ < rhs.left_;
333  }
334  //*******************************************************************************************
335 
336  //**Greater-than operator********************************************************************
342  inline bool operator>( const ConstIterator& rhs ) const {
343  return left_ > rhs.left_;
344  }
345  //*******************************************************************************************
346 
347  //**Less-or-equal-than operator**************************************************************
353  inline bool operator<=( const ConstIterator& rhs ) const {
354  return left_ <= rhs.left_;
355  }
356  //*******************************************************************************************
357 
358  //**Greater-or-equal-than operator***********************************************************
364  inline bool operator>=( const ConstIterator& rhs ) const {
365  return left_ >= rhs.left_;
366  }
367  //*******************************************************************************************
368 
369  //**Subtraction operator*********************************************************************
375  inline DifferenceType operator-( const ConstIterator& rhs ) const {
376  return left_ - rhs.left_;
377  }
378  //*******************************************************************************************
379 
380  //**Addition operator************************************************************************
387  friend inline const ConstIterator operator+( const ConstIterator& it, size_t inc ) {
388  return ConstIterator( it.left_ + inc, it.right_ + inc );
389  }
390  //*******************************************************************************************
391 
392  //**Addition operator************************************************************************
399  friend inline const ConstIterator operator+( size_t inc, const ConstIterator& it ) {
400  return ConstIterator( it.left_ + inc, it.right_ + inc );
401  }
402  //*******************************************************************************************
403 
404  //**Subtraction operator*********************************************************************
411  friend inline const ConstIterator operator-( const ConstIterator& it, size_t dec ) {
412  return ConstIterator( it.left_ - dec, it.right_ - dec );
413  }
414  //*******************************************************************************************
415 
416  private:
417  //**Member variables*************************************************************************
420  //*******************************************************************************************
421  };
422  //**********************************************************************************************
423 
424  //**Compilation flags***************************************************************************
426  enum { vectorizable = VT1::vectorizable && VT2::vectorizable &&
429 
431  enum { smpAssignable = VT1::smpAssignable && VT2::smpAssignable };
432  //**********************************************************************************************
433 
434  //**Constructor*********************************************************************************
440  explicit inline DVecDVecSubExpr( const VT1& lhs, const VT2& rhs )
441  : lhs_( lhs ) // Left-hand side dense vector of the subtraction expression
442  , rhs_( rhs ) // Right-hand side dense vector of the subtraction expression
443  {
444  BLAZE_INTERNAL_ASSERT( lhs.size() == rhs.size(), "Invalid vector sizes" );
445  }
446  //**********************************************************************************************
447 
448  //**Subscript operator**************************************************************************
454  inline ReturnType operator[]( size_t index ) const {
455  BLAZE_INTERNAL_ASSERT( index < lhs_.size(), "Invalid vector access index" );
456  return lhs_[index] - rhs_[index];
457  }
458  //**********************************************************************************************
459 
460  //**Load function*******************************************************************************
466  inline IntrinsicType load( size_t index ) const {
467  typedef IntrinsicTrait<ElementType> IT;
468  BLAZE_INTERNAL_ASSERT( index < lhs_.size() , "Invalid vector access index" );
469  BLAZE_INTERNAL_ASSERT( index % IT::size == 0UL, "Invalid vector access index" );
470  const IntrinsicType xmm1( lhs_.load( index ) );
471  const IntrinsicType xmm2( rhs_.load( index ) );
472  return xmm1 - xmm2;
473  }
474  //**********************************************************************************************
475 
476  //**Begin function******************************************************************************
481  inline ConstIterator begin() const {
482  return ConstIterator( lhs_.begin(), rhs_.begin() );
483  }
484  //**********************************************************************************************
485 
486  //**End function********************************************************************************
491  inline ConstIterator end() const {
492  return ConstIterator( lhs_.end(), rhs_.end() );
493  }
494  //**********************************************************************************************
495 
496  //**Size function*******************************************************************************
501  inline size_t size() const {
502  return lhs_.size();
503  }
504  //**********************************************************************************************
505 
506  //**Left operand access*************************************************************************
511  inline LeftOperand leftOperand() const {
512  return lhs_;
513  }
514  //**********************************************************************************************
515 
516  //**Right operand access************************************************************************
521  inline RightOperand rightOperand() const {
522  return rhs_;
523  }
524  //**********************************************************************************************
525 
526  //**********************************************************************************************
532  template< typename T >
533  inline bool canAlias( const T* alias ) const {
534  return ( IsExpression<VT1>::value && ( RequiresEvaluation<VT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
535  ( IsExpression<VT2>::value && ( RequiresEvaluation<VT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
536  }
537  //**********************************************************************************************
538 
539  //**********************************************************************************************
545  template< typename T >
546  inline bool isAliased( const T* alias ) const {
547  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
548  }
549  //**********************************************************************************************
550 
551  //**********************************************************************************************
556  inline bool isAligned() const {
557  return lhs_.isAligned() && rhs_.isAligned();
558  }
559  //**********************************************************************************************
560 
561  //**********************************************************************************************
566  inline bool canSMPAssign() const {
567  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
569  }
570  //**********************************************************************************************
571 
572  private:
573  //**Member variables****************************************************************************
576  //**********************************************************************************************
577 
578  //**Assignment to dense vectors*****************************************************************
592  template< typename VT > // Type of the target dense vector
593  friend inline typename EnableIf< UseAssign<VT> >::Type
594  assign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
595  {
597 
598  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
599 
600  if( !IsComputation<VT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
601  subAssign( ~lhs, rhs.rhs_ );
602  }
603  else {
604  assign ( ~lhs, rhs.lhs_ );
605  subAssign( ~lhs, rhs.rhs_ );
606  }
607  }
609  //**********************************************************************************************
610 
611  //**Assignment to sparse vectors****************************************************************
625  template< typename VT > // Type of the target sparse vector
626  friend inline typename EnableIf< UseAssign<VT> >::Type
627  assign( SparseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
628  {
630 
634 
635  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
636 
637  const ResultType tmp( serial( rhs ) );
638  assign( ~lhs, tmp );
639  }
641  //**********************************************************************************************
642 
643  //**Addition assignment to dense vectors********************************************************
657  template< typename VT > // Type of the target dense vector
658  friend inline typename EnableIf< UseAssign<VT> >::Type
659  addAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
660  {
662 
663  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
664 
665  addAssign( ~lhs, rhs.lhs_ );
666  subAssign( ~lhs, rhs.rhs_ );
667  }
669  //**********************************************************************************************
670 
671  //**Addition assignment to sparse vectors*******************************************************
672  // No special implementation for the addition assignment to sparse vectors.
673  //**********************************************************************************************
674 
675  //**Subtraction assignment to dense vectors*****************************************************
689  template< typename VT > // Type of the target dense vector
690  friend inline typename EnableIf< UseAssign<VT> >::Type
691  subAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
692  {
694 
695  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
696 
697  subAssign( ~lhs, rhs.lhs_ );
698  addAssign( ~lhs, rhs.rhs_ );
699  }
701  //**********************************************************************************************
702 
703  //**Subtraction assignment to sparse vectors****************************************************
704  // No special implementation for the subtraction assignment to sparse vectors.
705  //**********************************************************************************************
706 
707  //**Multiplication assignment to dense vectors**************************************************
721  template< typename VT > // Type of the target dense vector
722  friend inline typename EnableIf< UseAssign<VT> >::Type
723  multAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
724  {
726 
730 
731  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
732 
733  const ResultType tmp( serial( rhs ) );
734  multAssign( ~lhs, tmp );
735  }
737  //**********************************************************************************************
738 
739  //**Multiplication assignment to sparse vectors*************************************************
740  // No special implementation for the multiplication assignment to sparse vectors.
741  //**********************************************************************************************
742 
743  //**SMP assignment to dense vectors*************************************************************
757  template< typename VT > // Type of the target dense vector
758  friend inline typename EnableIf< UseSMPAssign<VT> >::Type
759  smpAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
760  {
762 
763  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
764 
765  if( !IsComputation<VT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
766  smpSubAssign( ~lhs, rhs.rhs_ );
767  }
768  else {
769  smpAssign ( ~lhs, rhs.lhs_ );
770  smpSubAssign( ~lhs, rhs.rhs_ );
771  }
772  }
774  //**********************************************************************************************
775 
776  //**SMP assignment to sparse vectors************************************************************
790  template< typename VT > // Type of the target sparse vector
791  friend inline typename EnableIf< UseSMPAssign<VT> >::Type
792  smpAssign( SparseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
793  {
795 
799 
800  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
801 
802  const ResultType tmp( rhs );
803  smpAssign( ~lhs, tmp );
804  }
806  //**********************************************************************************************
807 
808  //**SMP addition assignment to dense vectors****************************************************
822  template< typename VT > // Type of the target dense vector
823  friend inline typename EnableIf< UseSMPAssign<VT> >::Type
824  smpAddAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
825  {
827 
828  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
829 
830  smpAddAssign( ~lhs, rhs.lhs_ );
831  smpSubAssign( ~lhs, rhs.rhs_ );
832  }
834  //**********************************************************************************************
835 
836  //**SMP addition assignment to sparse vectors***************************************************
837  // No special implementation for the SMP addition assignment to sparse vectors.
838  //**********************************************************************************************
839 
840  //**SMP subtraction assignment to dense vectors*************************************************
854  template< typename VT > // Type of the target dense vector
855  friend inline typename EnableIf< UseSMPAssign<VT> >::Type
856  smpSubAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
857  {
859 
860  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
861 
862  smpSubAssign( ~lhs, rhs.lhs_ );
863  smpAddAssign( ~lhs, rhs.rhs_ );
864  }
866  //**********************************************************************************************
867 
868  //**SMP subtraction assignment to sparse vectors************************************************
869  // No special implementation for the SMP subtraction assignment to sparse vectors.
870  //**********************************************************************************************
871 
872  //**SMP multiplication assignment to dense vectors**********************************************
886  template< typename VT > // Type of the target dense vector
887  friend inline typename EnableIf< UseSMPAssign<VT> >::Type
888  smpMultAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
889  {
891 
895 
896  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
897 
898  const ResultType tmp( rhs );
899  smpMultAssign( ~lhs, tmp );
900  }
902  //**********************************************************************************************
903 
904  //**SMP multiplication assignment to sparse vectors*********************************************
905  // No special implementation for the SMP multiplication assignment to sparse vectors.
906  //**********************************************************************************************
907 
908  //**Compile time checks*************************************************************************
916  //**********************************************************************************************
917 };
918 //*************************************************************************************************
919 
920 
921 
922 
923 //=================================================================================================
924 //
925 // GLOBAL BINARY ARITHMETIC OPERATORS
926 //
927 //=================================================================================================
928 
929 //*************************************************************************************************
953 template< typename T1 // Type of the left-hand side dense vector
954  , typename T2 // Type of the right-hand side dense vector
955  , bool TF > // Transpose flag
956 inline const DVecDVecSubExpr<T1,T2,TF>
958 {
960 
961  if( (~lhs).size() != (~rhs).size() )
962  throw std::invalid_argument( "Vector sizes do not match" );
963 
964  return DVecDVecSubExpr<T1,T2,TF>( ~lhs, ~rhs );
965 }
966 //*************************************************************************************************
967 
968 
969 
970 
971 //=================================================================================================
972 //
973 // SIZE SPECIALIZATIONS
974 //
975 //=================================================================================================
976 
977 //*************************************************************************************************
979 template< typename VT1, typename VT2, bool TF >
980 struct Size< DVecDVecSubExpr<VT1,VT2,TF> >
981  : public Max< Size<VT1>, Size<VT2> >::Type
982 {};
984 //*************************************************************************************************
985 
986 
987 
988 
989 //=================================================================================================
990 //
991 // EXPRESSION TRAIT SPECIALIZATIONS
992 //
993 //=================================================================================================
994 
995 //*************************************************************************************************
997 template< typename VT1, typename VT2, bool TF, bool AF >
998 struct SubvectorExprTrait< DVecDVecSubExpr<VT1,VT2,TF>, AF >
999 {
1000  public:
1001  //**********************************************************************************************
1002  typedef typename SubExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
1003  , typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
1004  //**********************************************************************************************
1005 };
1007 //*************************************************************************************************
1008 
1009 } // namespace blaze
1010 
1011 #endif
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DVecDVecSubExpr.h:185
RightIteratorType right_
Iterator to the current right-hand side element.
Definition: DVecDVecSubExpr.h:419
PointerType pointer
Pointer return type.
Definition: DVecDVecSubExpr.h:190
Pointer difference type of the Blaze library.
Header file for the Max class template.
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
BLAZE_ALWAYS_INLINE void multAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the multiplication assignment of a matrix to a matrix.
Definition: Matrix.h:879
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
ValueType value_type
Type of the underlying elements.
Definition: DVecDVecSubExpr.h:189
Header file for the subtraction trait.
ReferenceType reference
Reference return type.
Definition: DVecDVecSubExpr.h:191
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:258
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:171
const ConstIterator operator++(int)
Post-increment operator.
Definition: DVecDVecSubExpr.h:256
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: DVecDVecSubExpr.h:454
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:205
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:168
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 IsSame and IsStrictlySame type traits.
DVecDVecSubExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the DVecDVecSubExpr class.
Definition: DVecDVecSubExpr.h:440
VT1::ReturnType RN1
Return type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:98
friend const ConstIterator operator+(const ConstIterator &it, size_t inc)
Addition between a ConstIterator and an integral value.
Definition: DVecDVecSubExpr.h:387
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
VT2::ReturnType RN2
Return type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:99
Header file for the DenseVector base class.
VT1::ElementType ET1
Element type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:102
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
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DVecDVecSubExpr.h:556
Header file for the RequiresEvaluation type trait.
ConstIterator & operator-=(size_t dec)
Subtraction assignment operator.
Definition: DVecDVecSubExpr.h:232
DVecDVecSubExpr< VT1, VT2, TF > This
Type of this DVecDVecSubExpr instance.
Definition: DVecDVecSubExpr.h:155
ElementType * PointerType
Pointer return type.
Definition: DVecDVecSubExpr.h:183
Constraint on the data type.
friend const ConstIterator operator-(const ConstIterator &it, size_t dec)
Subtraction between a ConstIterator and an integral value.
Definition: DVecDVecSubExpr.h:411
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:90
ConstIterator & operator+=(size_t inc)
Addition assignment operator.
Definition: DVecDVecSubExpr.h:219
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:259
IntrinsicTrait< ElementType >::Type IntrinsicType
Resulting intrinsic element type.
Definition: DVecDVecSubExpr.h:159
ReturnType operator*() const
Direct access to the element at the current iterator position.
Definition: DVecDVecSubExpr.h:288
bool operator<=(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:353
Compile time type selection.The SelectType class template selects one of the two given types T1 and T...
Definition: SelectType.h:59
bool operator>=(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:364
Header file for the IsTemporary type trait class.
IteratorCategory iterator_category
The iterator category.
Definition: DVecDVecSubExpr.h:188
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two iterators.
Definition: DVecDVecSubExpr.h:375
VT1::CompositeType CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:100
ConstIterator & operator--()
Pre-decrement operator.
Definition: DVecDVecSubExpr.h:266
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
std::random_access_iterator_tag IteratorCategory
The iterator category.
Definition: DVecDVecSubExpr.h:181
Iterator over the elements of the dense vector.
Definition: DVecDVecSubExpr.h:177
ResultType::ElementType ElementType
Resulting element type.
Definition: DVecDVecSubExpr.h:158
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
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:320
DifferenceType difference_type
Difference between two iterators.
Definition: DVecDVecSubExpr.h:192
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
const size_t SMP_DVECDVECSUB_THRESHOLD
SMP dense vector/dense vector subtraction threshold.This threshold specifies when a dense vector/dens...
Definition: Thresholds.h:253
SelectType< useAssign, const ResultType, const DVecDVecSubExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: DVecDVecSubExpr.h:165
Constraint on the data type.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DVecDVecSubExpr.h:157
#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
Constraint on the data type.
VT2::ElementType ET2
Element type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:103
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:104
LeftIteratorType left_
Iterator to the current left-hand side element.
Definition: DVecDVecSubExpr.h:418
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
friend const ConstIterator operator+(size_t inc, const ConstIterator &it)
Addition between an integral value and a ConstIterator.
Definition: DVecDVecSubExpr.h:399
Header file for the VecVecSubExpr base class.
Header file for the SelectType class template.
ConstIterator end() const
Returns an iterator just past the last non-zero element of the dense vector.
Definition: DVecDVecSubExpr.h:491
Header file for all forward declarations for expression class templates.
ConstIterator begin() const
Returns an iterator to the first non-zero element of the dense vector.
Definition: DVecDVecSubExpr.h:481
ElementType ValueType
Type of the underlying elements.
Definition: DVecDVecSubExpr.h:182
Header file for the EnableIf class template.
Header file for the serial shim.
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DVecDVecSubExpr.h:566
ConstIterator(LeftIteratorType left, RightIteratorType right)
Constructor for the ConstIterator class.
Definition: DVecDVecSubExpr.h:207
bool operator<(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:331
VT1::ConstIterator LeftIteratorType
ConstIterator type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:195
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
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.
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
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DVecDVecSubExpr.h:162
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
ElementType & ReferenceType
Reference return type.
Definition: DVecDVecSubExpr.h:184
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DVecDVecSubExpr.h:546
ConstIterator & operator++()
Pre-increment operator.
Definition: DVecDVecSubExpr.h:244
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DVecDVecSubExpr.h:156
VT2::ConstIterator RightIteratorType
ConstIterator type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:198
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:309
IntrinsicType load() const
Access to the intrinsic elements of the vector.
Definition: DVecDVecSubExpr.h:298
RightOperand rhs_
Right-hand side dense vector of the subtraction expression.
Definition: DVecDVecSubExpr.h:575
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DVecDVecSubExpr.h:116
Header file for all intrinsic functionality.
Constraint on the data type.
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a dense, N-dimensional vector type...
Definition: DenseVector.h:79
VT2::CompositeType CT2
Composite type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:101
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DVecDVecSubExpr.h:533
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_VECVECSUBEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/vector ...
Definition: VecVecSubExpr.h:165
Header file for the IsComputation type trait class.
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 vector operand.
Definition: DVecDVecSubExpr.h:511
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:108
#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
VT2::ResultType RT2
Result type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:97
Header file for basic type definitions.
IntrinsicType load(size_t index) const
Access to the intrinsic elements of the vector.
Definition: DVecDVecSubExpr.h:466
Header file for the SubvectorExprTrait class template.
Base template for the SubTrait class.
Definition: SubTrait.h:142
VT1::ResultType RT1
Result type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:96
RightOperand rightOperand() const
Returns the right-hand side dense vector operand.
Definition: DVecDVecSubExpr.h:521
Base class for all vector/vector subtraction expression templates.The VecVecSubExpr class serves as a...
Definition: VecVecSubExpr.h:65
LeftOperand lhs_
Left-hand side dense vector of the subtraction expression.
Definition: DVecDVecSubExpr.h:574
const ConstIterator operator--(int)
Post-decrement operator.
Definition: DVecDVecSubExpr.h:278
Header file for the SubExprTrait class template.
size_t size() const
Returns the current size/dimension of the vector.
Definition: DVecDVecSubExpr.h:501
Expression object for dense vector-dense vector subtractions.The DVecDVecSubExpr class represents the...
Definition: DVecDVecSubExpr.h:90
EnableIf< IsDenseVector< VT1 > >::Type smpMultAssign(Vector< VT1, TF1 > &lhs, const Vector< VT2, TF2 > &rhs)
Default implementation of the SMP multiplication assignment of a vector to a dense vector...
Definition: DenseVector.h:189
Header file for the Size type trait.
#define BLAZE_CONSTRAINT_MUST_BE_VECTOR_WITH_TRANSPOSE_FLAG(T, TF)
Constraint on the data type.In case the given data type T is not a dense or sparse vector type and in...
Definition: TransposeFlag.h:238
Header file for the thresholds for matrix/vector and matrix/matrix multiplications.
bool operator>(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:342
#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