All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
SVecTDVecMultExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_SVECTDVECMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SVECTDVECMULTEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <iterator>
55 #include <blaze/math/Intrinsics.h>
71 #include <blaze/util/Assert.h>
73 #include <blaze/util/EnableIf.h>
75 #include <blaze/util/SelectType.h>
76 #include <blaze/util/Types.h>
79 #include <blaze/util/Unused.h>
80 
81 
82 namespace blaze {
83 
84 //=================================================================================================
85 //
86 // CLASS SVECTDVECMULTEXPR
87 //
88 //=================================================================================================
89 
90 //*************************************************************************************************
97 template< typename VT1 // Type of the left-hand side sparse vector
98  , typename VT2 > // Type of the right-hand side dense vector
99 class SVecTDVecMultExpr : public SparseMatrix< SVecTDVecMultExpr<VT1,VT2>, true >
100  , private VecTVecMultExpr
101  , private Computation
102 {
103  private:
104  //**Type definitions****************************************************************************
105  typedef typename VT1::ResultType RT1;
106  typedef typename VT2::ResultType RT2;
107  typedef typename VT1::ReturnType RN1;
108  typedef typename VT2::ReturnType RN2;
109  typedef typename VT1::CompositeType CT1;
110  typedef typename VT2::CompositeType CT2;
111  typedef typename VT1::ElementType ET1;
112  typedef typename VT2::ElementType ET2;
113  //**********************************************************************************************
114 
115  //**Return type evaluation**********************************************************************
117 
122  enum { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
123 
126  //**********************************************************************************************
127 
128  //**Evaluation strategy*************************************************************************
130 
136  enum { useAssign = ( IsComputation<VT1>::value || !IsNumeric<ET1>::value ||
138 
140  template< typename MT >
142  struct UseAssign {
143  enum { value = useAssign };
144  };
146  //**********************************************************************************************
147 
148  //**********************************************************************************************
150 
153  template< typename T1, typename T2, typename T3 >
154  struct UseVectorizedKernel {
155  enum { value = T1::vectorizable && T3::vectorizable &&
159  };
161  //**********************************************************************************************
162 
163  //**********************************************************************************************
165 
168  template< typename T1, typename T2, typename T3 >
169  struct UseDefaultKernel {
170  enum { value = !UseVectorizedKernel<T1,T2,T3>::value };
171  };
173  //**********************************************************************************************
174 
175  public:
176  //**Type definitions****************************************************************************
182 
185 
188 
190  typedef typename SelectType< IsExpression<VT1>::value, const VT1, const VT1& >::Type LeftOperand;
191 
193  typedef typename SelectType< IsExpression<VT2>::value, const VT2, const VT2& >::Type RightOperand;
194 
196  typedef typename SelectType< IsComputation<VT1>::value, const RT1, CT1 >::Type LT;
197 
199  typedef typename SelectType< IsComputation<VT2>::value, const RT2, CT2 >::Type RT;
200  //**********************************************************************************************
201 
202  //**Compilation flags***************************************************************************
204  enum { smpAssignable = 0 };
205  //**********************************************************************************************
206 
207  //**ConstIterator class definition**************************************************************
211  {
212  public:
213  //**Type definitions*************************************************************************
216 
219 
221  typedef ET2 RightElement;
222 
223  typedef std::forward_iterator_tag IteratorCategory;
224  typedef Element ValueType;
228 
229  // STL iterator requirements
235  //*******************************************************************************************
236 
237  //**Constructor******************************************************************************
241  : it_( it ) // Iterator over the elements of the left-hand side sparse vector expression
242  , v_ ( v ) // Element of the right-hand side dense vector expression.
243  {}
244  //*******************************************************************************************
245 
246  //**Prefix increment operator****************************************************************
252  ++it_;
253  return *this;
254  }
255  //*******************************************************************************************
256 
257  //**Element access operator******************************************************************
262  inline const Element operator*() const {
263  return Element( it_->value() * v_, it_->index() );
264  }
265  //*******************************************************************************************
266 
267  //**Element access operator******************************************************************
272  inline const ConstIterator* operator->() const {
273  return this;
274  }
275  //*******************************************************************************************
276 
277  //**Value function***************************************************************************
282  inline ReturnType value() const {
283  return it_->value() * v_;
284  }
285  //*******************************************************************************************
286 
287  //**Index function***************************************************************************
292  inline size_t index() const {
293  return it_->index();
294  }
295  //*******************************************************************************************
296 
297  //**Equality operator************************************************************************
303  inline bool operator==( const ConstIterator& rhs ) const {
304  return it_ == rhs.it_;
305  }
306  //*******************************************************************************************
307 
308  //**Inequality operator**********************************************************************
314  inline bool operator!=( const ConstIterator& rhs ) const {
315  return it_ != rhs.it_;
316  }
317  //*******************************************************************************************
318 
319  //**Subtraction operator*********************************************************************
325  inline DifferenceType operator-( const ConstIterator& rhs ) const {
326  return it_ - rhs.it_;
327  }
328  //*******************************************************************************************
329 
330  private:
331  //**Member variables*************************************************************************
334  //*******************************************************************************************
335  };
336  //**********************************************************************************************
337 
338  //**Constructor*********************************************************************************
344  explicit inline SVecTDVecMultExpr( const VT1& lhs, const VT2& rhs )
345  : lhs_( lhs ) // Left-hand side sparse vector of the multiplication expression
346  , rhs_( rhs ) // Right-hand side dense vector of the multiplication expression
347  {}
348  //**********************************************************************************************
349 
350  //**Access operator*****************************************************************************
357  inline ReturnType operator()( size_t i, size_t j ) const {
358  BLAZE_INTERNAL_ASSERT( i < lhs_.size(), "Invalid row access index" );
359  BLAZE_INTERNAL_ASSERT( j < rhs_.size(), "Invalid column access index" );
360 
361  return lhs_[i] * rhs_[j];
362  }
363  //**********************************************************************************************
364 
365  //**Begin function******************************************************************************
371  inline ConstIterator begin( size_t i ) const {
372  return ConstIterator( lhs_.begin(), rhs_[i] );
373  }
374  //**********************************************************************************************
375 
376  //**End function********************************************************************************
382  inline ConstIterator end( size_t i ) const {
383  return ConstIterator( lhs_.end(), rhs_[i] );
384  }
385  //**********************************************************************************************
386 
387  //**Rows function*******************************************************************************
392  inline size_t rows() const {
393  return lhs_.size();
394  }
395  //**********************************************************************************************
396 
397  //**Columns function****************************************************************************
402  inline size_t columns() const {
403  return rhs_.size();
404  }
405  //**********************************************************************************************
406 
407  //**NonZeros function***************************************************************************
412  inline size_t nonZeros() const {
413  return lhs_.nonZeros() * rhs_.size();
414  }
415  //**********************************************************************************************
416 
417  //**NonZeros function***************************************************************************
423  inline size_t nonZeros( size_t i ) const {
424  UNUSED_PARAMETER( i );
425  return lhs_.nonZeros();
426  }
427  //**********************************************************************************************
428 
429  //**Find function*******************************************************************************
436  inline ConstIterator find( size_t i, size_t j ) const {
438  return ConstIterator( lhs_.find( i ), rhs_[j] );
439  }
440  //**********************************************************************************************
441 
442  //**LowerBound function*************************************************************************
449  inline ConstIterator lowerBound( size_t i, size_t j ) const {
451  return ConstIterator( lhs_.lowerBound( i ), rhs_[j] );
452  }
453  //**********************************************************************************************
454 
455  //**UpperBound function*************************************************************************
462  inline ConstIterator upperBound( size_t i, size_t j ) const {
464  return ConstIterator( lhs_.upperBound( i ), rhs_[j] );
465  }
466  //**********************************************************************************************
467 
468  //**Left operand access*************************************************************************
473  inline LeftOperand leftOperand() const {
474  return lhs_;
475  }
476  //**********************************************************************************************
477 
478  //**Right operand access************************************************************************
483  inline RightOperand rightOperand() const {
484  return rhs_;
485  }
486  //**********************************************************************************************
487 
488  //**********************************************************************************************
494  template< typename T >
495  inline bool canAlias( const T* alias ) const {
496  return ( lhs_.canAlias( alias ) || rhs_.canAlias( alias ) );
497  }
498  //**********************************************************************************************
499 
500  //**********************************************************************************************
506  template< typename T >
507  inline bool isAliased( const T* alias ) const {
508  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
509  }
510  //**********************************************************************************************
511 
512  private:
513  //**Member variables****************************************************************************
514  LeftOperand lhs_;
515  RightOperand rhs_;
516  //**********************************************************************************************
517 
518  //**Assignment to row-major dense matrices******************************************************
530  template< typename MT > // Type of the target dense matrix
531  friend inline void assign( DenseMatrix<MT,false>& lhs, const SVecTDVecMultExpr& rhs )
532  {
534 
536 
537  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
538  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
539 
540  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
541  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
542 
543  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
544  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
545  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
546  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
547 
548  SVecTDVecMultExpr::selectAssignKernel( ~lhs, x, y );
549  }
551  //**********************************************************************************************
552 
553  //**Default assignment to row-major dense matrices**********************************************
567  template< typename MT // Type of the left-hand side target matrix
568  , typename VT3 // Type of the left-hand side vector operand
569  , typename VT4 > // Type of the right-hand side vector operand
570  static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
571  selectAssignKernel( DenseMatrix<MT,false>& A, const VT3& x, const VT4& y )
572  {
574 
575  const ConstIterator end( x.end() );
576 
577  for( ConstIterator element=x.begin(); element!=end; ++element ) {
578  if( !isDefault( element->value() ) ) {
579  for( size_t j=0UL; j<y.size(); ++j ) {
580  (~A)(element->index(),j) = element->value() * y[j];
581  }
582  }
583  }
584  }
586  //**********************************************************************************************
587 
588  //**Vectorized assignment to row-major dense matrices*******************************************
602  template< typename MT // Type of the left-hand side target matrix
603  , typename VT3 // Type of the left-hand side vector operand
604  , typename VT4 > // Type of the right-hand side vector operand
605  static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
606  selectAssignKernel( DenseMatrix<MT,false>& A, const VT3& x, const VT4& y )
607  {
609 
610  typedef IntrinsicTrait<ElementType> IT;
611  typedef typename IT::Type IntrinsicType;
612 
613  const ConstIterator begin( x.begin() );
614  const ConstIterator end ( x.end() );
615 
616  for( ConstIterator element=begin; element!=end; ++element )
617  {
618  const IntrinsicType x1( set( element->value() ) );
619 
620  for( size_t j=0UL; j<(~A).columns(); j+=IT::size ) {
621  (~A).store( element->index(), j, x1 * y.load(j) );
622  }
623  }
624  }
626  //**********************************************************************************************
627 
628  //**Assignment to column-major dense matrices***************************************************
644  template< typename MT > // Type of the target dense matrix
645  friend inline typename EnableIf< UseAssign<MT> >::Type
646  assign( DenseMatrix<MT,true>& lhs, const SVecTDVecMultExpr& rhs )
647  {
649 
650  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
651  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
652 
654 
655  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
656  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
657 
658  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
659  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
660  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
661  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
662 
663  const ConstIterator end( x.end() );
664 
665  for( size_t i=0UL; i<y.size(); ++i ) {
666  if( !isDefault( y[i] ) ) {
667  for( ConstIterator element=x.begin(); element!=end; ++element ) {
668  (~lhs)(element->index(),i) = element->value() * y[i];
669  }
670  }
671  }
672  }
674  //**********************************************************************************************
675 
676  //**Assignment to row-major sparse matrices*****************************************************
688  template< typename MT > // Type of the target sparse matrix
689  friend inline void assign( SparseMatrix<MT,false>& lhs, const SVecTDVecMultExpr& rhs )
690  {
692 
694 
695  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
696  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
697 
699 
700  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
701  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
702 
703  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
704  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
705  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
706  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
707 
708  const ConstIterator end( x.end() );
709 
710  for( ConstIterator element=x.begin(); element!=end; ++element ) {
711  if( !isDefault( element->value() ) ) {
712  (~lhs).reserve( element->index(), y.size() );
713  for( size_t i=0UL; i<y.size(); ++i ) {
714  (~lhs).append( element->index(), i, element->value() * y[i] );
715  }
716  }
717  }
718  }
720  //**********************************************************************************************
721 
722  //**Assignment to column-major sparse matrices**************************************************
738  template< typename MT > // Type of the target sparse matrix
739  friend inline typename EnableIf< UseAssign<MT> >::Type
740  assign( SparseMatrix<MT,true>& lhs, const SVecTDVecMultExpr& rhs )
741  {
743 
744  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
745  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
746 
748 
749  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
750  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
751 
752  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
753  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
754  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
755  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
756 
757  const ConstIterator end( x.end() );
758 
759  for( size_t i=0UL; i<y.size(); ++i ) {
760  if( !isDefault( y[i] ) ) {
761  (~lhs).reserve( i, x.nonZeros() );
762  for( ConstIterator element=x.begin(); element!=end; ++element ) {
763  (~lhs).append( element->index(), i, element->value() * y[i] );
764  }
765  }
766  }
767  }
769  //**********************************************************************************************
770 
771  //**Addition assignment to row-major dense matrices*********************************************
784  template< typename MT > // Type of the target dense matrix
785  friend inline void addAssign( DenseMatrix<MT,false>& lhs, const SVecTDVecMultExpr& rhs )
786  {
788 
790 
791  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
792  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
793 
794  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
795  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
796 
797  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
798  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
799  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
800  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
801 
802  SVecTDVecMultExpr::selectAddAssignKernel( ~lhs, x, y );
803  }
805  //**********************************************************************************************
806 
807  //**Default addition assignment to row-major dense matrices*************************************
821  template< typename MT // Type of the left-hand side target matrix
822  , typename VT3 // Type of the left-hand side vector operand
823  , typename VT4 > // Type of the right-hand side vector operand
824  static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
825  selectAddAssignKernel( DenseMatrix<MT,false>& A, const VT3& x, const VT4& y )
826  {
828 
829  const ConstIterator end( x.end() );
830 
831  for( ConstIterator element=x.begin(); element!=end; ++element ) {
832  if( !isDefault( element->value() ) ) {
833  for( size_t i=0UL; i<y.size(); ++i ) {
834  (~A)(element->index(),i) += element->value() * y[i];
835  }
836  }
837  }
838  }
840  //**********************************************************************************************
841 
842  //**Vectorized addition assignment to row-major dense matrices**********************************
856  template< typename MT // Type of the left-hand side target matrix
857  , typename VT3 // Type of the left-hand side vector operand
858  , typename VT4 > // Type of the right-hand side vector operand
859  static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
860  selectAddAssignKernel( DenseMatrix<MT,false>& A, const VT3& x, const VT4& y )
861  {
863 
864  typedef IntrinsicTrait<ElementType> IT;
865  typedef typename IT::Type IntrinsicType;
866 
867  const ConstIterator begin( x.begin() );
868  const ConstIterator end ( x.end() );
869 
870  for( ConstIterator element=begin; element!=end; ++element )
871  {
872  const IntrinsicType x1( set( element->value() ) );
873 
874  for( size_t j=0UL; j<(~A).columns(); j+=IT::size ) {
875  (~A).store( element->index(), j, (~A).load(element->index(),j) + x1 * y.load(j) );
876  }
877  }
878  }
880  //**********************************************************************************************
881 
882  //**Addition assignment to column-major dense matrices******************************************
898  template< typename MT > // Type of the target dense matrix
899  friend inline typename EnableIf< UseAssign<MT> >::Type
900  addAssign( DenseMatrix<MT,true>& lhs, const SVecTDVecMultExpr& rhs )
901  {
903 
904  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
905  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
906 
908 
909  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
910  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
911 
912  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
913  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
914  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
915  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
916 
917  const ConstIterator end( x.end() );
918 
919  for( size_t i=0UL; i<y.size(); ++i ) {
920  if( !isDefault( y[i] ) ) {
921  for( ConstIterator element=x.begin(); element!=end; ++element ) {
922  (~lhs)(element->index(),i) += element->value() * y[i];
923  }
924  }
925  }
926  }
928  //**********************************************************************************************
929 
930  //**Addition assignment to sparse matrices******************************************************
931  // No special implementation for the addition assignment to sparse matrices.
932  //**********************************************************************************************
933 
934  //**Subtraction assignment to row-major dense matrices******************************************
947  template< typename MT > // Type of the target dense matrix
948  friend inline void subAssign( DenseMatrix<MT,false>& lhs, const SVecTDVecMultExpr& rhs )
949  {
951 
953 
954  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
955  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
956 
957  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
958  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
959 
960  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
961  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
962  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
963  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
964 
965  SVecTDVecMultExpr::selectSubAssignKernel( ~lhs, x, y );
966  }
968  //**********************************************************************************************
969 
970  //**Default subtraction assignment to row-major dense matrices**********************************
984  template< typename MT // Type of the left-hand side target matrix
985  , typename VT3 // Type of the left-hand side vector operand
986  , typename VT4 > // Type of the right-hand side vector operand
987  static inline typename EnableIf< UseDefaultKernel<MT,VT3,VT4> >::Type
988  selectSubAssignKernel( DenseMatrix<MT,false>& A, const VT3& x, const VT4& y )
989  {
991 
992  const ConstIterator end( x.end() );
993 
994  for( ConstIterator element=x.begin(); element!=end; ++element ) {
995  if( !isDefault( element->value() ) ) {
996  for( size_t i=0UL; i<y.size(); ++i ) {
997  (~A)(element->index(),i) -= element->value() * y[i];
998  }
999  }
1000  }
1001  }
1003  //**********************************************************************************************
1004 
1005  //**Vectorized subtraction assignment to row-major dense matrices*******************************
1019  template< typename MT // Type of the left-hand side target matrix
1020  , typename VT3 // Type of the left-hand side vector operand
1021  , typename VT4 > // Type of the right-hand side vector operand
1022  static inline typename EnableIf< UseVectorizedKernel<MT,VT3,VT4> >::Type
1023  selectSubAssignKernel( DenseMatrix<MT,false>& A, const VT3& x, const VT4& y )
1024  {
1026 
1027  typedef IntrinsicTrait<ElementType> IT;
1028  typedef typename IT::Type IntrinsicType;
1029 
1030  const ConstIterator begin( x.begin() );
1031  const ConstIterator end ( x.end() );
1032 
1033  for( ConstIterator element=begin; element!=end; ++element )
1034  {
1035  const IntrinsicType x1( set( element->value() ) );
1036 
1037  for( size_t j=0UL; j<(~A).columns(); j+=IT::size ) {
1038  (~A).store( element->index(), j, (~A).load(element->index(),j) - x1 * y.load(j) );
1039  }
1040  }
1041  }
1043  //**********************************************************************************************
1044 
1045  //**Subtraction assignment to column-major dense matrices***************************************
1061  template< typename MT > // Type of the target dense matrix
1062  friend inline typename EnableIf< UseAssign<MT> >::Type
1063  subAssign( DenseMatrix<MT,true>& lhs, const SVecTDVecMultExpr& rhs )
1064  {
1066 
1067  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
1068  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
1069 
1071 
1072  LT x( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse vector operand
1073  RT y( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side dense vector operand
1074 
1075  BLAZE_INTERNAL_ASSERT( x.size() == rhs.lhs_.size() , "Invalid vector size" );
1076  BLAZE_INTERNAL_ASSERT( y.size() == rhs.rhs_.size() , "Invalid vector size" );
1077  BLAZE_INTERNAL_ASSERT( x.size() == (~lhs).rows() , "Invalid vector size" );
1078  BLAZE_INTERNAL_ASSERT( y.size() == (~lhs).columns(), "Invalid vector size" );
1079 
1080  const ConstIterator end( x.end() );
1081 
1082  for( size_t i=0UL; i<y.size(); ++i ) {
1083  if( !isDefault( y[i] ) ) {
1084  for( ConstIterator element=x.begin(); element!=end; ++element ) {
1085  (~lhs)(element->index(),i) -= element->value() * y[i];
1086  }
1087  }
1088  }
1089  }
1091  //**********************************************************************************************
1092 
1093  //**Subtraction assignment to sparse matrices***************************************************
1094  // No special implementation for the subtraction assignment to sparse matrices.
1095  //**********************************************************************************************
1096 
1097  //**Multiplication assignment to dense matrices*************************************************
1098  // No special implementation for the multiplication assignment to dense matrices.
1099  //**********************************************************************************************
1100 
1101  //**Multiplication assignment to sparse matrices************************************************
1102  // No special implementation for the multiplication assignment to sparse matrices.
1103  //**********************************************************************************************
1104 
1105  //**Compile time checks*************************************************************************
1113  //**********************************************************************************************
1114 };
1115 //*************************************************************************************************
1116 
1117 
1118 
1119 
1120 //=================================================================================================
1121 //
1122 // GLOBAL BINARY ARITHMETIC OPERATORS
1123 //
1124 //=================================================================================================
1125 
1126 //*************************************************************************************************
1155 template< typename T1 // Type of the left-hand side sparse vector
1156  , typename T2 > // Type of the right-hand side dense vector
1157 inline const SVecTDVecMultExpr<T1,T2>
1159 {
1161 
1162  return SVecTDVecMultExpr<T1,T2>( ~lhs, ~rhs );
1163 }
1164 //*************************************************************************************************
1165 
1166 
1167 
1168 
1169 //=================================================================================================
1170 //
1171 // ROWS SPECIALIZATIONS
1172 //
1173 //=================================================================================================
1174 
1175 //*************************************************************************************************
1177 template< typename VT1, typename VT2 >
1178 struct Rows< SVecTDVecMultExpr<VT1,VT2> > : public Size<VT1>
1179 {};
1181 //*************************************************************************************************
1182 
1183 
1184 
1185 
1186 //=================================================================================================
1187 //
1188 // COLUMNS SPECIALIZATIONS
1189 //
1190 //=================================================================================================
1191 
1192 //*************************************************************************************************
1194 template< typename VT1, typename VT2 >
1195 struct Columns< SVecTDVecMultExpr<VT1,VT2> > : public Size<VT2>
1196 {};
1198 //*************************************************************************************************
1199 
1200 
1201 
1202 
1203 //=================================================================================================
1204 //
1205 // EXPRESSION TRAIT SPECIALIZATIONS
1206 //
1207 //=================================================================================================
1208 
1209 //*************************************************************************************************
1211 template< typename VT1, typename VT2, bool AF >
1212 struct SubmatrixExprTrait< SVecTDVecMultExpr<VT1,VT2>, AF >
1213 {
1214  public:
1215  //**********************************************************************************************
1216  typedef typename MultExprTrait< typename SubvectorExprTrait<const VT1,AF>::Type
1217  , typename SubvectorExprTrait<const VT2,AF>::Type >::Type Type;
1218  //**********************************************************************************************
1219 };
1221 //*************************************************************************************************
1222 
1223 
1224 //*************************************************************************************************
1226 template< typename VT1, typename VT2 >
1227 struct RowExprTrait< SVecTDVecMultExpr<VT1,VT2> >
1228 {
1229  public:
1230  //**********************************************************************************************
1231  typedef typename MultExprTrait< typename VT1::ReturnType, VT2 >::Type Type;
1232  //**********************************************************************************************
1233 };
1235 //*************************************************************************************************
1236 
1237 
1238 //*************************************************************************************************
1240 template< typename VT1, typename VT2 >
1241 struct ColumnExprTrait< SVecTDVecMultExpr<VT1,VT2> >
1242 {
1243  public:
1244  //**********************************************************************************************
1245  typedef typename MultExprTrait< VT1, typename VT2::ReturnType >::Type Type;
1246  //**********************************************************************************************
1247 };
1249 //*************************************************************************************************
1250 
1251 } // namespace blaze
1252 
1253 #endif
std::forward_iterator_tag IteratorCategory
The iterator category.
Definition: SVecTDVecMultExpr.h:223
Pointer difference type of the Blaze library.
ValueType value_type
Type of the underlying pointers.
Definition: SVecTDVecMultExpr.h:231
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
Compile time check for numeric types.This type trait tests whether or not the given template paramete...
Definition: IsNumeric.h:98
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.
Header file for the UNUSED_PARAMETER function template.
const DMatDMatMultExpr< T1, T2 > operator*(const DenseMatrix< T1, false > &lhs, const DenseMatrix< T2, false > &rhs)
Multiplication operator for the multiplication of two row-major dense matrices ( ).
Definition: DMatDMatMultExpr.h:4838
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SVecTDVecMultExpr.h:179
Iterator over the elements of the sparse vector-dense vector outer product expression.
Definition: SVecTDVecMultExpr.h:210
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: SVecTDVecMultExpr.h:495
RightOperand rightOperand() const
Returns the right-hand side dense vector operand.
Definition: SVecTDVecMultExpr.h:483
BLAZE_ALWAYS_INLINE size_t size(const Vector< VT, TF > &vector)
Returns the current size/dimension of the vector.
Definition: Vector.h:258
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a column dense or sparse vector type...
Definition: TransposeFlag.h:159
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:205
SelectType< IsExpression< VT2 >::value, const VT2, const VT2 & >::Type RightOperand
Composite type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:193
Header file for the ColumnExprTrait class template.
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of column i.
Definition: SVecTDVecMultExpr.h:371
ConstIterator lowerBound(size_t i, size_t j) const
Returns an iterator to the first index not less then the given index.
Definition: SVecTDVecMultExpr.h:449
Expression object for sparse vector-dense vector outer products.The SVecTDVecMultExpr class represent...
Definition: Forward.h:121
VT1::ResultType RT1
Result type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:105
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
size_t index() const
Access to the current index of the sparse element.
Definition: SVecTDVecMultExpr.h:292
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
BLAZE_ALWAYS_INLINE EnableIf< IsIntegral< T >, Load< T, sizeof(T)> >::Type::Type load(const T *address)
Loads a vector of integral values.
Definition: Load.h:224
void UNUSED_PARAMETER(const T1 &)
Suppression of unused parameter warnings.
Definition: Unused.h:81
SelectType< IsComputation< VT1 >::value, const RT1, CT1 >::Type LT
Type for the assignment of the left-hand side dense vector operand.
Definition: SVecTDVecMultExpr.h:196
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of column i.
Definition: SVecTDVecMultExpr.h:382
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
ConstIterator(IteratorType it, RightElement v)
Constructor for the ConstIterator class.
Definition: SVecTDVecMultExpr.h:240
const Element operator*() const
Direct access to the sparse matrix element at the current iterator position.
Definition: SVecTDVecMultExpr.h:262
RightOperand rhs_
Right-hand side dense vector of the multiplication expression.
Definition: SVecTDVecMultExpr.h:515
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Constraint on the data type.
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: SVecTDVecMultExpr.h:227
DifferenceType difference_type
Difference between two iterators.
Definition: SVecTDVecMultExpr.h:234
Header file for the SparseMatrix base class.
Constraint on the data type.
Header file for the MultExprTrait class template.
IteratorType it_
Iterator over the elements of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:332
SVecTDVecMultExpr(const VT1 &lhs, const VT2 &rhs)
Constructor for the SVecTDVecMultExpr class.
Definition: SVecTDVecMultExpr.h:344
Header file for the ValueIndexPair class.
SelectType< IsComputation< VT2 >::value, const RT2, CT2 >::Type RT
Type for the assignment of the right-hand side dense vector operand.
Definition: SVecTDVecMultExpr.h:199
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.
Header file for the multiplication trait.
RemoveReference< LeftOperand >::Type::ConstIterator IteratorType
Iterator type of the sparse vector expression.
Definition: SVecTDVecMultExpr.h:218
size_t columns() const
Returns the current number of columns of the matrix.
Definition: SVecTDVecMultExpr.h:402
VT2::CompositeType CT2
Composite type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:110
bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: SVecTDVecMultExpr.h:314
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2482
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_VECTVECMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/vector ...
Definition: VecTVecMultExpr.h:161
MultExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: SVecTDVecMultExpr.h:125
Element ValueType
Type of the underlying pointers.
Definition: SVecTDVecMultExpr.h:224
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.
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
#define BLAZE_CONSTRAINT_MUST_BE_SPARSE_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a sparse, N-dimensional vector type...
Definition: SparseVector.h:79
ET2 RightElement
Element type of the dense vector expression.
Definition: SVecTDVecMultExpr.h:221
Constraint on the data type.
VT2::ResultType RT2
Result type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:106
size_t rows() const
Returns the current number of rows of the matrix.
Definition: SVecTDVecMultExpr.h:392
RightElement v_
Element of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:333
Constraint on the data type.
SelectType< useAssign, const ResultType, const SVecTDVecMultExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: SVecTDVecMultExpr.h:187
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2476
Constraint on the data type.
SVecTDVecMultExpr< VT1, VT2 > This
Type of this SVecTDVecMultExpr instance.
Definition: SVecTDVecMultExpr.h:177
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Header file for the VecTVecMultExpr base class.
VT1::CompositeType CT1
Composite type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:109
Header file for the EnableIf class template.
Header file for the serial shim.
ReferenceType reference
Reference return type.
Definition: SVecTDVecMultExpr.h:233
ConstIterator find(size_t i, size_t j) const
Searches for a specific matrix element.
Definition: SVecTDVecMultExpr.h:436
bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: SVecTDVecMultExpr.h:303
Header file for the IsNumeric type trait.
VT2::ReturnType RN2
Return type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:108
BLAZE_ALWAYS_INLINE EnableIf< IsIntegral< T >, Set< T, sizeof(T)> >::Type::Type set(T value)
Sets all values in the vector to the given integral value.
Definition: Set.h:211
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: SVecTDVecMultExpr.h:184
#define BLAZE_CONSTRAINT_MUST_NOT_BE_SYMMETRIC_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is a symmetric matrix type, a compilation error is created.
Definition: Symmetric.h:116
Header file for the SubmatrixExprTrait class template.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2477
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
IteratorCategory iterator_category
The iterator category.
Definition: SVecTDVecMultExpr.h:230
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 MultTrait class.
Definition: MultTrait.h:142
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
ReturnType value() const
Access to the current value of the sparse element.
Definition: SVecTDVecMultExpr.h:282
LeftOperand leftOperand() const
Returns the left-hand side sparse vector operand.
Definition: SVecTDVecMultExpr.h:473
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
VT2::ElementType ET2
Element type of the right-hand side dense vector expression.
Definition: SVecTDVecMultExpr.h:112
ResultType::ElementType ElementType
Resulting element type.
Definition: SVecTDVecMultExpr.h:181
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: SVecTDVecMultExpr.h:507
Header file for the isDefault shim.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: SVecTDVecMultExpr.h:180
BLAZE_ALWAYS_INLINE bool isDefault(const NonNumericProxy< MT > &proxy)
Returns whether the represented element is in default state.
Definition: NonNumericProxy.h:874
Constraint on the data type.
ValueType * PointerType
Pointer return type.
Definition: SVecTDVecMultExpr.h:225
ConstIterator upperBound(size_t i, size_t j) const
Returns an iterator to the first index greater then the given index.
Definition: SVecTDVecMultExpr.h:462
VT1::ElementType ET1
Element type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:111
ValueIndexPair< ElementType > Element
Element type of the sparse matrix expression.
Definition: SVecTDVecMultExpr.h:215
#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:118
Constraint on the data type.
Header file for the RemoveReference type trait.
LeftOperand lhs_
Left-hand side sparse vector of the multiplication expression.
Definition: SVecTDVecMultExpr.h:514
Header file for all intrinsic functionality.
PointerType pointer
Pointer return type.
Definition: SVecTDVecMultExpr.h:232
DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two expression iterators.
Definition: SVecTDVecMultExpr.h:325
#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
Index-value-pair for sparse vectors and matrices.The ValueIndexPair class represents a single index-v...
Definition: ValueIndexPair.h:70
Header file for the IsComputation type trait class.
const ConstIterator * operator->() const
Direct access to the sparse matrix element at the current iterator position.
Definition: SVecTDVecMultExpr.h:272
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
Header file for basic type definitions.
ConstIterator & operator++()
Pre-increment operator.
Definition: SVecTDVecMultExpr.h:251
#define BLAZE_CONSTRAINT_MUST_BE_ROW_VECTOR_TYPE(T)
Constraint on the data type.In case the given data type T is not a row dense or sparse vector type (i...
Definition: TransposeFlag.h:81
Header file for the SubvectorExprTrait class template.
ValueType & ReferenceType
Reference return type.
Definition: SVecTDVecMultExpr.h:226
MultTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: SVecTDVecMultExpr.h:178
VT1::ReturnType RN1
Return type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:107
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified column.
Definition: SVecTDVecMultExpr.h:423
BLAZE_ALWAYS_INLINE EnableIf< IsIntegral< T > >::Type store(T *address, const typename Store< T, sizeof(T)>::Type &value)
Aligned store of a vector of integral values.
Definition: Store.h:225
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SVecTDVecMultExpr.h:357
Header file for the Size type trait.
Evaluation of the resulting expression type of a multiplication.Via this type trait it is possible to...
Definition: MultExprTrait.h:137
SelectType< IsExpression< VT1 >::value, const VT1, const VT1 & >::Type LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: SVecTDVecMultExpr.h:190
#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
Header file for the IsExpression type trait class.
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SVecTDVecMultExpr.h:412
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