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 <iterator>
44 #include <blaze/math/Aliases.h>
49 #include <blaze/math/Exception.h>
55 #include <blaze/math/SIMD.h>
66 #include <blaze/system/Inline.h>
68 #include <blaze/util/Assert.h>
69 #include <blaze/util/EnableIf.h>
71 #include <blaze/util/mpl/And.h>
72 #include <blaze/util/mpl/If.h>
73 #include <blaze/util/mpl/Maximum.h>
74 #include <blaze/util/Types.h>
75 
76 
77 namespace blaze {
78 
79 //=================================================================================================
80 //
81 // CLASS DVECDVECSUBEXPR
82 //
83 //=================================================================================================
84 
85 //*************************************************************************************************
92 template< typename VT1 // Type of the left-hand side dense vector
93  , typename VT2 // Type of the right-hand side dense vector
94  , bool TF > // Transpose flag
96  : public VecVecSubExpr< DenseVector< DVecDVecSubExpr<VT1,VT2,TF>, TF > >
97  , private Computation
98 {
99  private:
100  //**Type definitions****************************************************************************
109  //**********************************************************************************************
110 
111  //**Return type evaluation**********************************************************************
113 
118  enum : bool { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
119 
122  //**********************************************************************************************
123 
124  //**Serial evaluation strategy******************************************************************
126 
132  enum : bool { useAssign = ( RequiresEvaluation<VT1>::value || RequiresEvaluation<VT2>::value || !returnExpr ) };
133 
135  template< typename VT >
137  struct UseAssign {
138  enum : bool { value = useAssign };
139  };
141  //**********************************************************************************************
142 
143  //**Parallel evaluation strategy****************************************************************
145 
151  template< typename VT >
152  struct UseSMPAssign {
153  enum : bool { value = ( !VT1::smpAssignable || !VT2::smpAssignable ) && useAssign };
154  };
156  //**********************************************************************************************
157 
158  public:
159  //**Type definitions****************************************************************************
164 
167 
170 
172  using LeftOperand = If_< IsExpression<VT1>, const VT1, const VT1& >;
173 
175  using RightOperand = If_< IsExpression<VT2>, const VT2, const VT2& >;
176  //**********************************************************************************************
177 
178  //**ConstIterator class definition**************************************************************
182  {
183  public:
184  //**Type definitions*************************************************************************
185  using IteratorCategory = std::random_access_iterator_tag;
190 
191  // STL iterator requirements
197 
200 
203  //*******************************************************************************************
204 
205  //**Constructor******************************************************************************
211  explicit inline ConstIterator( LeftIteratorType left, RightIteratorType right )
212  : left_ ( left ) // Iterator to the current left-hand side element
213  , right_( right ) // Iterator to the current right-hand side element
214  {}
215  //*******************************************************************************************
216 
217  //**Addition assignment operator*************************************************************
223  inline ConstIterator& operator+=( size_t inc ) {
224  left_ += inc;
225  right_ += inc;
226  return *this;
227  }
228  //*******************************************************************************************
229 
230  //**Subtraction assignment operator**********************************************************
236  inline ConstIterator& operator-=( size_t dec ) {
237  left_ -= dec;
238  right_ -= dec;
239  return *this;
240  }
241  //*******************************************************************************************
242 
243  //**Prefix increment operator****************************************************************
249  ++left_;
250  ++right_;
251  return *this;
252  }
253  //*******************************************************************************************
254 
255  //**Postfix increment operator***************************************************************
260  inline const ConstIterator operator++( int ) {
261  return ConstIterator( left_++, right_++ );
262  }
263  //*******************************************************************************************
264 
265  //**Prefix decrement operator****************************************************************
271  --left_;
272  --right_;
273  return *this;
274  }
275  //*******************************************************************************************
276 
277  //**Postfix decrement operator***************************************************************
282  inline const ConstIterator operator--( int ) {
283  return ConstIterator( left_--, right_-- );
284  }
285  //*******************************************************************************************
286 
287  //**Element access operator******************************************************************
292  inline ReturnType operator*() const {
293  return (*left_) - (*right_);
294  }
295  //*******************************************************************************************
296 
297  //**Load function****************************************************************************
302  inline auto load() const noexcept {
303  return left_.load() - right_.load();
304  }
305  //*******************************************************************************************
306 
307  //**Equality operator************************************************************************
313  inline bool operator==( const ConstIterator& rhs ) const {
314  return left_ == rhs.left_;
315  }
316  //*******************************************************************************************
317 
318  //**Inequality operator**********************************************************************
324  inline bool operator!=( const ConstIterator& rhs ) const {
325  return left_ != rhs.left_;
326  }
327  //*******************************************************************************************
328 
329  //**Less-than operator***********************************************************************
335  inline bool operator<( const ConstIterator& rhs ) const {
336  return left_ < rhs.left_;
337  }
338  //*******************************************************************************************
339 
340  //**Greater-than operator********************************************************************
346  inline bool operator>( const ConstIterator& rhs ) const {
347  return left_ > rhs.left_;
348  }
349  //*******************************************************************************************
350 
351  //**Less-or-equal-than operator**************************************************************
357  inline bool operator<=( const ConstIterator& rhs ) const {
358  return left_ <= rhs.left_;
359  }
360  //*******************************************************************************************
361 
362  //**Greater-or-equal-than operator***********************************************************
368  inline bool operator>=( const ConstIterator& rhs ) const {
369  return left_ >= rhs.left_;
370  }
371  //*******************************************************************************************
372 
373  //**Subtraction operator*********************************************************************
379  inline DifferenceType operator-( const ConstIterator& rhs ) const {
380  return left_ - rhs.left_;
381  }
382  //*******************************************************************************************
383 
384  //**Addition operator************************************************************************
391  friend inline const ConstIterator operator+( const ConstIterator& it, size_t inc ) {
392  return ConstIterator( it.left_ + inc, it.right_ + inc );
393  }
394  //*******************************************************************************************
395 
396  //**Addition operator************************************************************************
403  friend inline const ConstIterator operator+( size_t inc, const ConstIterator& it ) {
404  return ConstIterator( it.left_ + inc, it.right_ + inc );
405  }
406  //*******************************************************************************************
407 
408  //**Subtraction operator*********************************************************************
415  friend inline const ConstIterator operator-( const ConstIterator& it, size_t dec ) {
416  return ConstIterator( it.left_ - dec, it.right_ - dec );
417  }
418  //*******************************************************************************************
419 
420  private:
421  //**Member variables*************************************************************************
424  //*******************************************************************************************
425  };
426  //**********************************************************************************************
427 
428  //**Compilation flags***************************************************************************
430  enum : bool { simdEnabled = VT1::simdEnabled && VT2::simdEnabled &&
432 
434  enum : bool { smpAssignable = VT1::smpAssignable && VT2::smpAssignable };
435  //**********************************************************************************************
436 
437  //**SIMD properties*****************************************************************************
439  enum : size_t { SIMDSIZE = SIMDTrait<ElementType>::size };
440  //**********************************************************************************************
441 
442  //**Constructor*********************************************************************************
448  explicit inline DVecDVecSubExpr( const VT1& lhs, const VT2& rhs ) noexcept
449  : lhs_( lhs ) // Left-hand side dense vector of the subtraction expression
450  , rhs_( rhs ) // Right-hand side dense vector of the subtraction expression
451  {
452  BLAZE_INTERNAL_ASSERT( lhs.size() == rhs.size(), "Invalid vector sizes" );
453  }
454  //**********************************************************************************************
455 
456  //**Subscript operator**************************************************************************
462  inline ReturnType operator[]( size_t index ) const {
463  BLAZE_INTERNAL_ASSERT( index < lhs_.size(), "Invalid vector access index" );
464  return lhs_[index] - rhs_[index];
465  }
466  //**********************************************************************************************
467 
468  //**At function*********************************************************************************
475  inline ReturnType at( size_t index ) const {
476  if( index >= lhs_.size() ) {
477  BLAZE_THROW_OUT_OF_RANGE( "Invalid vector access index" );
478  }
479  return (*this)[index];
480  }
481  //**********************************************************************************************
482 
483  //**Load function*******************************************************************************
489  BLAZE_ALWAYS_INLINE auto load( size_t index ) const noexcept {
490  BLAZE_INTERNAL_ASSERT( index < lhs_.size() , "Invalid vector access index" );
491  BLAZE_INTERNAL_ASSERT( index % SIMDSIZE == 0UL, "Invalid vector access index" );
492  return lhs_.load( index ) - rhs_.load( index );
493  }
494  //**********************************************************************************************
495 
496  //**Begin function******************************************************************************
501  inline ConstIterator begin() const {
502  return ConstIterator( lhs_.begin(), rhs_.begin() );
503  }
504  //**********************************************************************************************
505 
506  //**End function********************************************************************************
511  inline ConstIterator end() const {
512  return ConstIterator( lhs_.end(), rhs_.end() );
513  }
514  //**********************************************************************************************
515 
516  //**Size function*******************************************************************************
521  inline size_t size() const noexcept {
522  return lhs_.size();
523  }
524  //**********************************************************************************************
525 
526  //**Left operand access*************************************************************************
531  inline LeftOperand leftOperand() const noexcept {
532  return lhs_;
533  }
534  //**********************************************************************************************
535 
536  //**Right operand access************************************************************************
541  inline RightOperand rightOperand() const noexcept {
542  return rhs_;
543  }
544  //**********************************************************************************************
545 
546  //**********************************************************************************************
552  template< typename T >
553  inline bool canAlias( const T* alias ) const noexcept {
554  return ( IsExpression<VT1>::value && ( RequiresEvaluation<VT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
555  ( IsExpression<VT2>::value && ( RequiresEvaluation<VT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
556  }
557  //**********************************************************************************************
558 
559  //**********************************************************************************************
565  template< typename T >
566  inline bool isAliased( const T* alias ) const noexcept {
567  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
568  }
569  //**********************************************************************************************
570 
571  //**********************************************************************************************
576  inline bool isAligned() const noexcept {
577  return lhs_.isAligned() && rhs_.isAligned();
578  }
579  //**********************************************************************************************
580 
581  //**********************************************************************************************
586  inline bool canSMPAssign() const noexcept {
587  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
588  ( size() > SMP_DVECDVECSUB_THRESHOLD );
589  }
590  //**********************************************************************************************
591 
592  private:
593  //**Member variables****************************************************************************
596  //**********************************************************************************************
597 
598  //**Assignment to dense vectors*****************************************************************
612  template< typename VT > // Type of the target dense vector
613  friend inline EnableIf_< UseAssign<VT> >
614  assign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
615  {
617 
618  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
619 
620  if( !IsComputation<VT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
621  subAssign( ~lhs, rhs.rhs_ );
622  }
623  else {
624  assign ( ~lhs, rhs.lhs_ );
625  subAssign( ~lhs, rhs.rhs_ );
626  }
627  }
629  //**********************************************************************************************
630 
631  //**Assignment to sparse vectors****************************************************************
645  template< typename VT > // Type of the target sparse vector
646  friend inline EnableIf_< UseAssign<VT> >
647  assign( SparseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
648  {
650 
654 
655  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
656 
657  const ResultType tmp( serial( rhs ) );
658  assign( ~lhs, tmp );
659  }
661  //**********************************************************************************************
662 
663  //**Addition assignment to dense vectors********************************************************
677  template< typename VT > // Type of the target dense vector
678  friend inline EnableIf_< UseAssign<VT> >
679  addAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
680  {
682 
683  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
684 
685  addAssign( ~lhs, rhs.lhs_ );
686  subAssign( ~lhs, rhs.rhs_ );
687  }
689  //**********************************************************************************************
690 
691  //**Addition assignment to sparse vectors*******************************************************
692  // No special implementation for the addition assignment to sparse vectors.
693  //**********************************************************************************************
694 
695  //**Subtraction assignment to dense vectors*****************************************************
709  template< typename VT > // Type of the target dense vector
710  friend inline EnableIf_< UseAssign<VT> >
711  subAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
712  {
714 
715  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
716 
717  subAssign( ~lhs, rhs.lhs_ );
718  addAssign( ~lhs, rhs.rhs_ );
719  }
721  //**********************************************************************************************
722 
723  //**Subtraction assignment to sparse vectors****************************************************
724  // No special implementation for the subtraction assignment to sparse vectors.
725  //**********************************************************************************************
726 
727  //**Multiplication assignment to dense vectors**************************************************
741  template< typename VT > // Type of the target dense vector
742  friend inline EnableIf_< UseAssign<VT> >
743  multAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
744  {
746 
750 
751  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
752 
753  const ResultType tmp( serial( rhs ) );
754  multAssign( ~lhs, tmp );
755  }
757  //**********************************************************************************************
758 
759  //**Multiplication assignment to sparse vectors*************************************************
760  // No special implementation for the multiplication assignment to sparse vectors.
761  //**********************************************************************************************
762 
763  //**Division assignment to dense vectors********************************************************
777  template< typename VT > // Type of the target dense vector
778  friend inline EnableIf_< UseAssign<VT> >
779  divAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
780  {
782 
786 
787  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
788 
789  const ResultType tmp( serial( rhs ) );
790  divAssign( ~lhs, tmp );
791  }
793  //**********************************************************************************************
794 
795  //**Division assignment to sparse vectors*******************************************************
796  // No special implementation for the division assignment to sparse vectors.
797  //**********************************************************************************************
798 
799  //**SMP assignment to dense vectors*************************************************************
813  template< typename VT > // Type of the target dense vector
814  friend inline EnableIf_< UseSMPAssign<VT> >
815  smpAssign( DenseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
816  {
818 
819  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
820 
821  if( !IsComputation<VT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
822  smpSubAssign( ~lhs, rhs.rhs_ );
823  }
824  else {
825  smpAssign ( ~lhs, rhs.lhs_ );
826  smpSubAssign( ~lhs, rhs.rhs_ );
827  }
828  }
830  //**********************************************************************************************
831 
832  //**SMP assignment to sparse vectors************************************************************
846  template< typename VT > // Type of the target sparse vector
847  friend inline EnableIf_< UseSMPAssign<VT> >
848  smpAssign( SparseVector<VT,TF>& lhs, const DVecDVecSubExpr& rhs )
849  {
851 
855 
856  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
857 
858  const ResultType tmp( rhs );
859  smpAssign( ~lhs, tmp );
860  }
862  //**********************************************************************************************
863 
864  //**SMP addition assignment to dense vectors****************************************************
878  template< typename VT > // Type of the target dense vector
879  friend inline EnableIf_< UseSMPAssign<VT> >
881  {
883 
884  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
885 
886  smpAddAssign( ~lhs, rhs.lhs_ );
887  smpSubAssign( ~lhs, rhs.rhs_ );
888  }
890  //**********************************************************************************************
891 
892  //**SMP addition assignment to sparse vectors***************************************************
893  // No special implementation for the SMP addition assignment to sparse vectors.
894  //**********************************************************************************************
895 
896  //**SMP subtraction assignment to dense vectors*************************************************
910  template< typename VT > // Type of the target dense vector
911  friend inline EnableIf_< UseSMPAssign<VT> >
913  {
915 
916  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
917 
918  smpSubAssign( ~lhs, rhs.lhs_ );
919  smpAddAssign( ~lhs, rhs.rhs_ );
920  }
922  //**********************************************************************************************
923 
924  //**SMP subtraction assignment to sparse vectors************************************************
925  // No special implementation for the SMP subtraction assignment to sparse vectors.
926  //**********************************************************************************************
927 
928  //**SMP multiplication assignment to dense vectors**********************************************
942  template< typename VT > // Type of the target dense vector
943  friend inline EnableIf_< UseSMPAssign<VT> >
945  {
947 
951 
952  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
953 
954  const ResultType tmp( rhs );
955  smpMultAssign( ~lhs, tmp );
956  }
958  //**********************************************************************************************
959 
960  //**SMP multiplication assignment to sparse vectors*********************************************
961  // No special implementation for the SMP multiplication assignment to sparse vectors.
962  //**********************************************************************************************
963 
964  //**SMP division assignment to dense vectors****************************************************
978  template< typename VT > // Type of the target dense vector
979  friend inline EnableIf_< UseSMPAssign<VT> >
981  {
983 
987 
988  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
989 
990  const ResultType tmp( rhs );
991  smpDivAssign( ~lhs, tmp );
992  }
994  //**********************************************************************************************
995 
996  //**SMP division assignment to sparse vectors***************************************************
997  // No special implementation for the SMP division assignment to sparse vectors.
998  //**********************************************************************************************
999 
1000  //**Compile time checks*************************************************************************
1008  //**********************************************************************************************
1009 };
1010 //*************************************************************************************************
1011 
1012 
1013 
1014 
1015 //=================================================================================================
1016 //
1017 // GLOBAL BINARY ARITHMETIC OPERATORS
1018 //
1019 //=================================================================================================
1020 
1021 //*************************************************************************************************
1045 template< typename VT1 // Type of the left-hand side dense vector
1046  , typename VT2 // Type of the right-hand side dense vector
1047  , bool TF > // Transpose flag
1048 inline decltype(auto)
1049  operator-( const DenseVector<VT1,TF>& lhs, const DenseVector<VT2,TF>& rhs )
1050 {
1052 
1053  if( (~lhs).size() != (~rhs).size() ) {
1054  BLAZE_THROW_INVALID_ARGUMENT( "Vector sizes do not match" );
1055  }
1056 
1058  return ReturnType( ~lhs, ~rhs );
1059 }
1060 //*************************************************************************************************
1061 
1062 
1063 
1064 
1065 //=================================================================================================
1066 //
1067 // SIZE SPECIALIZATIONS
1068 //
1069 //=================================================================================================
1070 
1071 //*************************************************************************************************
1073 template< typename VT1, typename VT2, bool TF >
1074 struct Size< DVecDVecSubExpr<VT1,VT2,TF>, 0UL >
1075  : public Maximum< Size<VT1,0UL>, Size<VT2,0UL> >
1076 {};
1078 //*************************************************************************************************
1079 
1080 
1081 
1082 
1083 //=================================================================================================
1084 //
1085 // ISALIGNED SPECIALIZATIONS
1086 //
1087 //=================================================================================================
1088 
1089 //*************************************************************************************************
1091 template< typename VT1, typename VT2, bool TF >
1092 struct IsAligned< DVecDVecSubExpr<VT1,VT2,TF> >
1093  : public And< IsAligned<VT1>, IsAligned<VT2> >
1094 {};
1096 //*************************************************************************************************
1097 
1098 
1099 
1100 
1101 //=================================================================================================
1102 //
1103 // ISPADDED SPECIALIZATIONS
1104 //
1105 //=================================================================================================
1106 
1107 //*************************************************************************************************
1109 template< typename VT1, typename VT2, bool TF >
1110 struct IsPadded< DVecDVecSubExpr<VT1,VT2,TF> >
1111  : public And< IsPadded<VT1>, IsPadded<VT2> >
1112 {};
1114 //*************************************************************************************************
1115 
1116 } // namespace blaze
1117 
1118 #endif
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:324
RightIteratorType right_
Iterator to the current right-hand side element.
Definition: DVecDVecSubExpr.h:423
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.This macro encapsulates the default way o...
Definition: Exception.h:235
Pointer difference type of the Blaze library.
Header file for auxiliary alias declarations.
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:71
Availability of a SIMD subtraction for the given data types.Depending on the available instruction se...
Definition: HasSIMDSub.h:171
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:69
bool isAligned() const noexcept
Returns whether the operands of the expression are properly aligned in memory.
Definition: DVecDVecSubExpr.h:576
If_< IsExpression< VT1 >, const VT1, const VT1 &> LeftOperand
Composite type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:172
Header file for the subtraction trait.
LeftOperand leftOperand() const noexcept
Returns the left-hand side dense vector operand.
Definition: DVecDVecSubExpr.h:531
BLAZE_ALWAYS_INLINE bool isSame(const Matrix< MT1, SO1 > &a, const Matrix< MT2, SO2 > &b) noexcept
Returns whether the two given matrices represent the same observable state.
Definition: Matrix.h:949
Header file for basic type definitions.
ElementType_< VT1 > ET1
Element type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:107
bool operator>(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:346
const ConstIterator operator++(int)
Post-increment operator.
Definition: DVecDVecSubExpr.h:260
std::random_access_iterator_tag IteratorCategory
The iterator category.
Definition: DVecDVecSubExpr.h:185
EnableIf_< IsDenseMatrix< MT1 > > 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:164
ElementType_< VT2 > ET2
Element type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:108
Header file for the serial shim.
ConstIterator begin() const
Returns an iterator to the first non-zero element of the dense vector.
Definition: DVecDVecSubExpr.h:501
ConstIterator end() const
Returns an iterator just past the last non-zero element of the dense vector.
Definition: DVecDVecSubExpr.h:511
friend const ConstIterator operator+(const ConstIterator &it, size_t inc)
Addition between a ConstIterator and an integral value.
Definition: DVecDVecSubExpr.h:391
IteratorCategory iterator_category
The iterator category.
Definition: DVecDVecSubExpr.h:192
SubTrait_< RT1, RT2 > ResultType
Result type for expression template evaluations.
Definition: DVecDVecSubExpr.h:161
EnableIf_< IsDenseVector< VT1 > > 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:193
Header file for the And class template.
Compile time value evaluation.The Maximum alias declaration selects the larger of the two given templ...
Definition: Maximum.h:73
Header file for the DenseVector base class.
ValueType value_type
Type of the underlying elements.
Definition: DVecDVecSubExpr.h:193
ResultType_< VT1 > RT1
Result type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:101
Header file for the Computation base class.
Header file for the RequiresEvaluation type trait.
ConstIterator & operator-=(size_t dec)
Subtraction assignment operator.
Definition: DVecDVecSubExpr.h:236
auto load() const noexcept
Access to the SIMD elements of the vector.
Definition: DVecDVecSubExpr.h:302
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:313
typename T::ResultType ResultType_
Alias declaration for nested ResultType type definitions.The ResultType_ alias declaration provides a...
Definition: Aliases.h:343
EnableIf_< IsDenseMatrix< MT1 > > 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:133
ReturnType_< VT1 > RN1
Return type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:103
TransposeType_< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: DVecDVecSubExpr.h:162
typename IfTrue< Condition, T1, T2 >::Type IfTrue_
Auxiliary alias declaration for the IfTrue class template.The IfTrue_ alias declaration provides a co...
Definition: If.h:109
typename T::ReturnType ReturnType_
Alias declaration for nested ReturnType type definitions.The ReturnType_ alias declaration provides a...
Definition: Aliases.h:363
ElementType * PointerType
Pointer return type.
Definition: DVecDVecSubExpr.h:187
CompositeType_< VT1 > CT1
Composite type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:105
friend const ConstIterator operator-(const ConstIterator &it, size_t dec)
Subtraction between a ConstIterator and an integral value.
Definition: DVecDVecSubExpr.h:415
Compile time check for the alignment of data types.This type trait tests whether the given data type ...
Definition: IsAligned.h:87
typename SubExprTrait< T1, T2 >::Type SubExprTrait_
Auxiliary alias declaration for the SubExprTrait class template.The SubExprTrait_ alias declaration p...
Definition: SubExprTrait.h:112
Header file for the Maximum class template.
bool operator>=(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:368
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:71
ConstIterator & operator+=(size_t inc)
Addition assignment operator.
Definition: DVecDVecSubExpr.h:223
typename T::CompositeType CompositeType_
Alias declaration for nested CompositeType type definitions.The CompositeType_ alias declaration prov...
Definition: Aliases.h:83
Header file for the IsTemporary type trait class.
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two iterators.
Definition: DVecDVecSubExpr.h:379
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
#define BLAZE_ALWAYS_INLINE
Platform dependent setup of an enforced inline keyword.
Definition: Inline.h:85
Header file for the If class template.
ConstIterator & operator--()
Pre-decrement operator.
Definition: DVecDVecSubExpr.h:270
Compile time check for data types with padding.This type trait tests whether the given data type empl...
Definition: IsPadded.h:76
EnableIf_< IsDenseMatrix< MT1 > > 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:102
Iterator over the elements of the dense vector.
Definition: DVecDVecSubExpr.h:181
RightOperand rightOperand() const noexcept
Returns the right-hand side dense vector operand.
Definition: DVecDVecSubExpr.h:541
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exception.This macro encapsulates the default way of Bl...
Definition: Exception.h:331
ConstIterator_< VT1 > LeftIteratorType
ConstIterator type of the left-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:199
ReturnType operator*() const
Direct access to the element at the current iterator position.
Definition: DVecDVecSubExpr.h:292
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:3085
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: DVecDVecSubExpr.h:553
typename T::ElementType ElementType_
Alias declaration for nested ElementType type definitions.The ElementType_ alias declaration provides...
Definition: Aliases.h:163
Header file for all SIMD functionality.
ElementType_< ResultType > ElementType
Resulting element type.
Definition: DVecDVecSubExpr.h:163
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:76
SubExprTrait_< RN1, RN2 > ExprReturnType
Expression return type for the subscript operator.
Definition: DVecDVecSubExpr.h:121
ElementType & ReferenceType
Reference return type.
Definition: DVecDVecSubExpr.h:188
PointerType pointer
Pointer return type.
Definition: DVecDVecSubExpr.h:194
const IfTrue_< returnExpr, ExprReturnType, ElementType > ReturnType
Return type for expression template evaluations.
Definition: DVecDVecSubExpr.h:166
Header file for the IsAligned type trait.
bool operator<=(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:357
Constraint on the data type.
Constraint on the data type.
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: DVecDVecSubExpr.h:475
Header file for the exception macros of the math module.
LeftIteratorType left_
Iterator to the current left-hand side element.
Definition: DVecDVecSubExpr.h:422
friend const ConstIterator operator+(size_t inc, const ConstIterator &it)
Addition between an integral value and a ConstIterator.
Definition: DVecDVecSubExpr.h:403
Header file for the VecVecSubExpr base class.
Constraint on the data type.
Header file for all forward declarations for expression class templates.
Header file for the EnableIf class template.
Header file for the IsPadded type trait.
ConstIterator(LeftIteratorType left, RightIteratorType right)
Constructor for the ConstIterator class.
Definition: DVecDVecSubExpr.h:211
size_t size() const noexcept
Returns the current size/dimension of the vector.
Definition: DVecDVecSubExpr.h:521
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: DVecDVecSubExpr.h:566
ElementType ValueType
Type of the underlying elements.
Definition: DVecDVecSubExpr.h:186
Header file for run time assertion macros.
ConstIterator_< VT2 > RightIteratorType
ConstIterator type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:202
typename If< T1, T2, T3 >::Type If_
Auxiliary alias declaration for the If class template.The If_ alias declaration provides a convenient...
Definition: If.h:154
EnableIf_< IsDenseVector< VT1 > > smpDivAssign(Vector< VT1, TF1 > &lhs, const Vector< VT2, TF2 > &rhs)
Default implementation of the SMP division assignment of a vector to a dense vector.
Definition: DenseVector.h:222
DVecDVecSubExpr(const VT1 &lhs, const VT2 &rhs) noexcept
Constructor for the DVecDVecSubExpr class.
Definition: DVecDVecSubExpr.h:448
ConstIterator & operator++()
Pre-increment operator.
Definition: DVecDVecSubExpr.h:248
SIMD characteristics of data types.The SIMDTrait class template provides the SIMD characteristics of ...
Definition: SIMDTrait.h:296
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:94
Header file for the HasSIMDSub type trait.
decltype(auto) serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:816
#define BLAZE_CONSTRAINT_MUST_NOT_REQUIRE_EVALUATION(T)
Constraint on the data type.In case the given data type T requires an intermediate evaluation within ...
Definition: RequiresEvaluation.h:81
BLAZE_ALWAYS_INLINE auto load(size_t index) const noexcept
Access to the SIMD elements of the vector.
Definition: DVecDVecSubExpr.h:489
typename EnableIf< Condition, T >::Type EnableIf_
Auxiliary alias declaration for the EnableIf class template.The EnableIf_ alias declaration provides ...
Definition: EnableIf.h:224
RightOperand rhs_
Right-hand side dense vector of the subtraction expression.
Definition: DVecDVecSubExpr.h:595
ReferenceType reference
Reference return type.
Definition: DVecDVecSubExpr.h:195
Constraint on the data type.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3080
typename T::ConstIterator ConstIterator_
Alias declaration for nested ConstIterator type definitions.The ConstIterator_ alias declaration prov...
Definition: Aliases.h:103
#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:61
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: DVecDVecSubExpr.h:586
#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:108
Header file for the IsComputation type trait class.
IfTrue_< useAssign, const ResultType, const DVecDVecSubExpr &> CompositeType
Data type for composite expression templates.
Definition: DVecDVecSubExpr.h:169
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:66
typename SubTrait< T1, T2 >::Type SubTrait_
Auxiliary alias declaration for the SubTrait class template.The SubTrait_ alias declaration provides ...
Definition: SubTrait.h:291
Compile time evaluation of the size of vectors and matrices.The Size type trait evaluates the size of...
Definition: Size.h:80
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:130
bool operator<(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DVecDVecSubExpr.h:335
If_< IsExpression< VT2 >, const VT2, const VT2 &> RightOperand
Composite type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:175
ReturnType_< VT2 > RN2
Return type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:104
Compile time logical &#39;and&#39; evaluation.The And alias declaration performs at compile time a logical &#39;a...
Definition: And.h:76
typename T::TransposeType TransposeType_
Alias declaration for nested TransposeType type definitions.The TransposeType_ alias declaration prov...
Definition: Aliases.h:423
CompositeType_< VT2 > CT2
Composite type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:106
Base class for all vector/vector subtraction expression templates.The VecVecSubExpr class serves as a...
Definition: VecVecSubExpr.h:66
LeftOperand lhs_
Left-hand side dense vector of the subtraction expression.
Definition: DVecDVecSubExpr.h:594
const ConstIterator operator--(int)
Post-decrement operator.
Definition: DVecDVecSubExpr.h:282
Header file for the SubExprTrait class template.
Expression object for dense vector-dense vector subtractions.The DVecDVecSubExpr class represents the...
Definition: DVecDVecSubExpr.h:95
System settings for the inline keywords.
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:63
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:95
Header file for the IsExpression type trait class.
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: DVecDVecSubExpr.h:462
Header file for the function trace functionality.
ResultType_< VT2 > RT2
Result type of the right-hand side dense vector expression.
Definition: DVecDVecSubExpr.h:102