TSVecTDMatMultExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_TSVECTDMATMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_TSVECTDMATMULTEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
54 #include <blaze/math/shims/Reset.h>
73 #include <blaze/util/Assert.h>
75 #include <blaze/util/DisableIf.h>
76 #include <blaze/util/EnableIf.h>
77 #include <blaze/util/Exception.h>
79 #include <blaze/util/mpl/Or.h>
80 #include <blaze/util/SelectType.h>
81 #include <blaze/util/Types.h>
83 
84 
85 namespace blaze {
86 
87 //=================================================================================================
88 //
89 // CLASS TSVECDMATMULTEXPR
90 //
91 //=================================================================================================
92 
93 //*************************************************************************************************
100 template< typename VT // Type of the left-hand side sparse vector
101  , typename MT > // Type of the right-hand side dense matrix
102 class TSVecTDMatMultExpr : public DenseVector< TSVecTDMatMultExpr<VT,MT>, true >
103  , private TVecMatMultExpr
104  , private Computation
105 {
106  private:
107  //**Type definitions****************************************************************************
108  typedef typename VT::ResultType VRT;
109  typedef typename MT::ResultType MRT;
110  typedef typename VT::CompositeType VCT;
111  typedef typename MT::CompositeType MCT;
112  //**********************************************************************************************
113 
114  //**********************************************************************************************
116  enum { evaluateVector = IsComputation<VT>::value || RequiresEvaluation<VT>::value };
117  //**********************************************************************************************
118 
119  //**********************************************************************************************
121  enum { evaluateMatrix = RequiresEvaluation<MT>::value };
122  //**********************************************************************************************
123 
124  //**********************************************************************************************
126 
132  enum { useAssign = ( evaluateVector || evaluateMatrix ) };
133  //**********************************************************************************************
134 
135  //**********************************************************************************************
137  template< typename VT2 >
139  struct UseAssign {
140  enum { value = useAssign };
141  };
143  //**********************************************************************************************
144 
145  //**********************************************************************************************
147 
151  template< typename T1 >
152  struct UseSMPAssign {
153  enum { value = useAssign };
154  };
156  //**********************************************************************************************
157 
158  public:
159  //**Type definitions****************************************************************************
164  typedef const ElementType ReturnType;
165 
168 
170  typedef typename SelectType< IsExpression<VT>::value, const VT, const VT& >::Type LeftOperand;
171 
173  typedef typename SelectType< IsExpression<MT>::value, const MT, const MT& >::Type RightOperand;
174 
177 
180  //**********************************************************************************************
181 
182  //**Compilation flags***************************************************************************
184  enum { vectorizable = 0 };
185 
187  enum { smpAssignable = !evaluateVector && VT::smpAssignable &&
188  !evaluateMatrix && MT::smpAssignable };
189  //**********************************************************************************************
190 
191  //**Constructor*********************************************************************************
197  explicit inline TSVecTDMatMultExpr( const VT& vec, const MT& mat )
198  : vec_( vec ) // Left-hand side sparse vector of the multiplication expression
199  , mat_( mat ) // Right-hand side dense matrix of the multiplication expression
200  {
201  BLAZE_INTERNAL_ASSERT( vec_.size() == mat_.rows(), "Invalid vector and matrix sizes" );
202  }
203  //**********************************************************************************************
204 
205  //**Subscript operator**************************************************************************
211  inline ReturnType operator[]( size_t index ) const {
212  BLAZE_INTERNAL_ASSERT( index < mat_.columns(), "Invalid vector access index" );
213 
215 
216  VCT x( vec_ ); // Evaluation of the left-hand side sparse vector operand
217 
218  BLAZE_INTERNAL_ASSERT( x.size() == vec_.size(), "Invalid vector size" );
219 
220  const ConstIterator end( ( IsUpper<MT>::value )
221  ?( IsStrictlyUpper<MT>::value ? x.lowerBound( index ) : x.upperBound( index ) )
222  :( x.end() ) );
223  ConstIterator element( ( IsLower<MT>::value )
224  ?( IsStrictlyLower<MT>::value ? x.upperBound( index ) : x.lowerBound( index ) )
225  :( x.begin() ) );
226 
227  ElementType res = ElementType();
228 
229  if( element != end ) {
230  res = element->value() * mat_( element->index(), index );
231  ++element;
232  for( ; element!=end; ++element ) {
233  res += element->value() * mat_( element->index(), index );
234  }
235  }
236 
237  return res;
238  }
239  //**********************************************************************************************
240 
241  //**At function*********************************************************************************
248  inline ReturnType at( size_t index ) const {
249  if( index >= mat_.columns() ) {
250  BLAZE_THROW_OUT_OF_RANGE( "Invalid vector access index" );
251  }
252  return (*this)[index];
253  }
254  //**********************************************************************************************
255 
256  //**Size function*******************************************************************************
261  inline size_t size() const {
262  return mat_.columns();
263  }
264  //**********************************************************************************************
265 
266  //**Left operand access*************************************************************************
271  inline LeftOperand leftOperand() const {
272  return vec_;
273  }
274  //**********************************************************************************************
275 
276  //**Right operand access************************************************************************
281  inline RightOperand rightOperand() const {
282  return mat_;
283  }
284  //**********************************************************************************************
285 
286  //**********************************************************************************************
292  template< typename T >
293  inline bool canAlias( const T* alias ) const {
294  return vec_.isAliased( alias ) || mat_.isAliased( alias );
295  }
296  //**********************************************************************************************
297 
298  //**********************************************************************************************
304  template< typename T >
305  inline bool isAliased( const T* alias ) const {
306  return vec_.isAliased( alias ) || mat_.isAliased( alias );
307  }
308  //**********************************************************************************************
309 
310  //**********************************************************************************************
315  inline bool isAligned() const {
316  return mat_.isAligned();
317  }
318  //**********************************************************************************************
319 
320  //**********************************************************************************************
325  inline bool canSMPAssign() const {
326  return ( size() > SMP_TSVECTDMATMULT_THRESHOLD );
327  }
328  //**********************************************************************************************
329 
330  private:
331  //**Member variables****************************************************************************
332  LeftOperand vec_;
333  RightOperand mat_;
334  //**********************************************************************************************
335 
336  //**Assignment to dense vectors*****************************************************************
352  template< typename VT2 > // Type of the target dense vector
353  friend inline typename EnableIf< UseAssign<VT2> >::Type
354  assign( DenseVector<VT2,true>& lhs, const TSVecTDMatMultExpr& rhs )
355  {
357 
358  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
359 
361 
362  // Evaluation of the left-hand side sparse vector operand
363  LT x( serial( rhs.vec_ ) );
364  if( x.nonZeros() == 0UL ) {
365  reset( ~lhs );
366  return;
367  }
368 
369  // Evaluation of the right-hand side dense matrix operand
370  RT A( serial( rhs.mat_ ) );
371 
372  // Checking the evaluated operands
373  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
374  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
375  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
376  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).size() , "Invalid vector size" );
377 
378  // Performing the sparse vector-dense matrix multiplication
379  assign( ~lhs, x * A );
380  }
382  //**********************************************************************************************
383 
384  //**Assignment to sparse vectors****************************************************************
400  template< typename VT2 > // Type of the target sparse vector
401  friend inline typename EnableIf< UseAssign<VT2> >::Type
402  assign( SparseVector<VT2,true>& lhs, const TSVecTDMatMultExpr& rhs )
403  {
405 
409 
410  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
411 
412  const ResultType tmp( serial( rhs ) );
413  assign( ~lhs, tmp );
414  }
416  //**********************************************************************************************
417 
418  //**Addition assignment to dense vectors********************************************************
433  template< typename VT2 > // Type of the target dense vector
434  friend inline typename EnableIf< UseAssign<VT2> >::Type
436  {
438 
439  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
440 
442 
443  // Evaluation of the left-hand side sparse vector operand
444  LT x( serial( rhs.vec_ ) );
445  if( x.nonZeros() == 0UL ) return;
446 
447  // Evaluation of the right-hand side dense matrix operand
448  RT A( serial( rhs.mat_ ) );
449 
450  // Checking the evaluated operands
451  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
452  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
453  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
454  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).size() , "Invalid vector size" );
455 
456  // Performing the sparse vector-dense matrix multiplication
457  addAssign( ~lhs, x * A );
458  }
459  //**********************************************************************************************
460 
461  //**Addition assignment to sparse vectors*******************************************************
462  // No special implementation for the addition assignment to sparse vectors.
463  //**********************************************************************************************
464 
465  //**Subtraction assignment to dense vectors*****************************************************
480  template< typename VT2 > // Type of the target dense vector
481  friend inline typename EnableIf< UseAssign<VT2> >::Type
483  {
485 
486  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
487 
489 
490  // Evaluation of the left-hand side sparse vector operand
491  LT x( serial( rhs.vec_ ) );
492  if( x.nonZeros() == 0UL ) return;
493 
494  // Evaluation of the right-hand side dense matrix operand
495  RT A( serial( rhs.mat_ ) );
496 
497  // Checking the evaluated operands
498  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
499  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
500  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
501  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).size() , "Invalid vector size" );
502 
503  // Performing the sparse vector-dense matrix multiplication
504  subAssign( ~lhs, x * A );
505  }
506  //**********************************************************************************************
507 
508  //**Subtraction assignment to sparse vectors****************************************************
509  // No special implementation for the subtraction assignment to sparse vectors.
510  //**********************************************************************************************
511 
512  //**Multiplication assignment to dense vectors**************************************************
527  template< typename VT2 > // Type of the target dense vector
528  friend inline typename EnableIf< UseAssign<VT2> >::Type
530  {
532 
536 
537  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
538 
539  const ResultType tmp( serial( rhs ) );
540  multAssign( ~lhs, tmp );
541  }
542  //**********************************************************************************************
543 
544  //**Multiplication assignment to sparse vectors*************************************************
545  // No special implementation for the multiplication assignment to sparse vectors.
546  //**********************************************************************************************
547 
548  //**SMP assignment to dense vectors*************************************************************
563  template< typename VT2 > // Type of the target dense vector
564  friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
566  {
568 
569  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
570 
572 
573  // Evaluation of the left-hand side sparse vector operand
574  LT x( rhs.vec_ );
575  if( x.nonZeros() == 0UL ) {
576  reset( ~lhs );
577  return;
578  }
579 
580  // Evaluation of the right-hand side dense matrix operand
581  RT A( rhs.mat_ );
582 
583  // Checking the evaluated operands
584  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
585  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
586  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
587  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).size() , "Invalid vector size" );
588 
589  // Performing the sparse vector-dense matrix multiplication
590  smpAssign( ~lhs, x * A );
591  }
593  //**********************************************************************************************
594 
595  //**SMP assignment to sparse vectors************************************************************
610  template< typename VT2 > // Type of the target sparse vector
611  friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
612  smpAssign( SparseVector<VT2,true>& lhs, const TSVecTDMatMultExpr& rhs )
613  {
615 
619 
620  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
621 
622  const ResultType tmp( rhs );
623  smpAssign( ~lhs, tmp );
624  }
626  //**********************************************************************************************
627 
628  //**SMP addition assignment to dense vectors****************************************************
642  template< typename VT2 > // Type of the target dense vector
643  friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
645  {
647 
648  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
649 
651 
652  // Evaluation of the left-hand side sparse vector operand
653  LT x( rhs.vec_ );
654  if( x.nonZeros() == 0UL ) return;
655 
656  // Evaluation of the right-hand side dense matrix operand
657  RT A( rhs.mat_ );
658 
659  // Checking the evaluated operands
660  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
661  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
662  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
663  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).size() , "Invalid vector size" );
664 
665  // Performing the sparse vector-dense matrix multiplication
666  smpAddAssign( ~lhs, x * A );
667  }
668  //**********************************************************************************************
669 
670  //**SMP addition assignment to sparse vectors***************************************************
671  // No special implementation for the SMP addition assignment to sparse vectors.
672  //**********************************************************************************************
673 
674  //**SMP subtraction assignment to dense vectors*************************************************
688  template< typename VT2 > // Type of the target dense vector
689  friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
691  {
693 
694  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
695 
697 
698  // Evaluation of the left-hand side sparse vector operand
699  LT x( rhs.vec_ );
700  if( x.nonZeros() == 0UL ) return;
701 
702  // Evaluation of the right-hand side dense matrix operand
703  RT A( rhs.mat_ );
704 
705  // Checking the evaluated operands
706  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
707  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
708  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
709  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).size() , "Invalid vector size" );
710 
711  // Performing the sparse vector-dense matrix multiplication
712  smpSubAssign( ~lhs, x * A );
713  }
714  //**********************************************************************************************
715 
716  //**SMP subtraction assignment to sparse vectors************************************************
717  // No special implementation for the SMP subtraction assignment to sparse vectors.
718  //**********************************************************************************************
719 
720  //**SMP multiplication assignment to dense vectors**********************************************
735  template< typename VT2 > // Type of the target dense vector
736  friend inline typename EnableIf< UseSMPAssign<VT2> >::Type
738  {
740 
744 
745  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
746 
747  const ResultType tmp( rhs );
748  smpMultAssign( ~lhs, tmp );
749  }
750  //**********************************************************************************************
751 
752  //**SMP multiplication assignment to sparse vectors*********************************************
753  // No special implementation for the SMP multiplication assignment to sparse vectors.
754  //**********************************************************************************************
755 
756  //**Compile time checks*************************************************************************
765  //**********************************************************************************************
766 };
767 //*************************************************************************************************
768 
769 
770 
771 
772 //=================================================================================================
773 //
774 // GLOBAL BINARY ARITHMETIC OPERATORS
775 //
776 //=================================================================================================
777 
778 //*************************************************************************************************
809 template< typename T1, typename T2 >
810 inline const typename DisableIf< Or< IsSymmetric<T2>, IsMatMatMultExpr<T2> >
811  , TSVecTDMatMultExpr<T1,T2> >::Type
813 {
815 
816  if( (~vec).size() != (~mat).rows() ) {
817  BLAZE_THROW_INVALID_ARGUMENT( "Vector and matrix sizes do not match" );
818  }
819 
820  return TSVecTDMatMultExpr<T1,T2>( ~vec, ~mat );
821 }
822 //*************************************************************************************************
823 
824 
825 
826 
827 //=================================================================================================
828 //
829 // GLOBAL RESTRUCTURING BINARY ARITHMETIC OPERATORS
830 //
831 //=================================================================================================
832 
833 //*************************************************************************************************
848 template< typename T1 // Type of the left-hand side sparse vector
849  , typename T2 > // Type of the right-hand side dense matrix
850 inline const typename EnableIf< IsSymmetric<T2>, typename MultExprTrait<T1,T2>::Type >::Type
851  operator*( const SparseVector<T1,true>& vec, const DenseMatrix<T2,true>& mat )
852 {
854 
856 
857  if( (~vec).size() != (~mat).rows() ) {
858  BLAZE_THROW_INVALID_ARGUMENT( "Vector and matrix sizes do not match" );
859  }
860 
861  return (~vec) * trans( ~mat );
862 }
864 //*************************************************************************************************
865 
866 
867 
868 
869 //=================================================================================================
870 //
871 // SIZE SPECIALIZATIONS
872 //
873 //=================================================================================================
874 
875 //*************************************************************************************************
877 template< typename VT, typename MT >
878 struct Size< TSVecTDMatMultExpr<VT,MT> > : public Columns<MT>
879 {};
881 //*************************************************************************************************
882 
883 
884 
885 
886 //=================================================================================================
887 //
888 // ISALIGNED SPECIALIZATIONS
889 //
890 //=================================================================================================
891 
892 //*************************************************************************************************
894 template< typename VT, typename MT >
895 struct IsAligned< TSVecTDMatMultExpr<VT,MT> > : public IsTrue< IsAligned<MT>::value >
896 {};
898 //*************************************************************************************************
899 
900 
901 
902 
903 //=================================================================================================
904 //
905 // EXPRESSION TRAIT SPECIALIZATIONS
906 //
907 //=================================================================================================
908 
909 //*************************************************************************************************
911 template< typename VT, typename MT, bool AF >
912 struct SubvectorExprTrait< TSVecTDMatMultExpr<VT,MT>, AF >
913 {
914  public:
915  //**********************************************************************************************
916  typedef typename MultExprTrait< typename SubvectorExprTrait<const VT,AF>::Type
917  , typename SubmatrixExprTrait<const MT,AF>::Type >::Type Type;
918  //**********************************************************************************************
919 };
921 //*************************************************************************************************
922 
923 } // namespace blaze
924 
925 #endif
LeftOperand vec_
Left-hand side sparse vector of the multiplication expression.
Definition: TSVecTDMatMultExpr.h:332
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exceptionThis macro encapsulates the default way of...
Definition: Exception.h:187
const ElementType ReturnType
Return type for expression template evaluations.
Definition: TSVecTDMatMultExpr.h:164
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:89
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:7820
Header file for basic type definitions.
BLAZE_ALWAYS_INLINE MT::Iterator end(Matrix< MT, SO > &matrix, size_t i)
Returns an iterator just past the last element of row/column i.
Definition: Matrix.h:250
SelectType< evaluateVector, const VRT, VCT >::Type LT
Type for the assignment of the left-hand side sparse vector operand.
Definition: TSVecTDMatMultExpr.h:176
size_t size() const
Returns the current size/dimension of the vector.
Definition: TSVecTDMatMultExpr.h:261
SelectType< evaluateMatrix, const MRT, MCT >::Type RT
Type for the assignment of the left-hand side dense matrix operand.
Definition: TSVecTDMatMultExpr.h:179
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:207
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a dense, N-dimensional matrix type...
Definition: DenseMatrix.h:79
void reset(const DiagonalProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: DiagonalProxy.h:507
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2588
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:259
BLAZE_ALWAYS_INLINE size_t rows(const Matrix< MT, SO > &matrix)
Returns the current number of rows of the matrix.
Definition: Matrix.h:308
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: TSVecTDMatMultExpr.h:281
Header file for the DenseVector base class.
Compile time check for lower triangular matrices.This type trait tests whether or not the given templ...
Definition: IsLower.h:90
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:721
Header file for the Computation base class.
MT::ResultType MRT
Result type of the right-hand side dense matrix expression.
Definition: TSVecTDMatMultExpr.h:109
Compile time check for upper triangular matrices.This type trait tests whether or not the given templ...
Definition: IsUpper.h:90
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: TSVecTDMatMultExpr.h:248
Header file for the RequiresEvaluation type trait.
ResultType::ElementType ElementType
Resulting element type.
Definition: TSVecTDMatMultExpr.h:163
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.
VT::CompositeType VCT
Composite type of the left-hand side sparse vector expression.
Definition: TSVecTDMatMultExpr.h:110
Constraint on the data type.
friend EnableIf< UseSMPAssign< VT2 > >::Type smpMultAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
SMP multiplication assignment of a transpose sparse vector-transpose dense matrix multiplication to a...
Definition: TSVecTDMatMultExpr.h:737
Header file for the MultExprTrait class template.
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:90
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:261
friend EnableIf< UseSMPAssign< VT2 > >::Type smpSubAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
SMP subtraction assignment of a transpose sparse vector-transpose dense matrix multiplication to a de...
Definition: TSVecTDMatMultExpr.h:690
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 DisableIf class template.
Header file for the multiplication trait.
Header file for the IsStrictlyUpper type trait.
Header file for the IsSymmetric type trait.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a column-major dense or sparse matri...
Definition: ColumnMajorMatrix.h:79
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: TSVecTDMatMultExpr.h:315
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:2592
Header file for the Or class template.
friend EnableIf< UseAssign< VT2 > >::Type multAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
Multiplication assignment of a transpose sparse vector-transpose dense matrix multiplication to a den...
Definition: TSVecTDMatMultExpr.h:529
Header file for the IsMatMatMultExpr type trait class.
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exceptionThis macro encapsulates the default way of Bla...
Definition: Exception.h:331
Header file for the Columns type trait.
SelectType< useAssign, const ResultType, const TSVecTDMatMultExpr & >::Type CompositeType
Data type for composite expression templates.
Definition: TSVecTDMatMultExpr.h:167
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
Header file for the IsLower type trait.
#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
Header file for the IsAligned type trait.
#define BLAZE_CONSTRAINT_MUST_NOT_BE_MATMATMULTEXPR_TYPE(T)
Constraint on the data type.In case the given data type T is a matrix/matrix multiplication expressio...
Definition: MatMatMultExpr.h:126
#define BLAZE_CONSTRAINT_MUST_BE_REFERENCE_TYPE(T)
Constraint on the data type.In case the given data type T is not a reference type, a compilation error is created.
Definition: Reference.h:78
Constraint on the data type.
Compile time check for strictly upper triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyUpper.h:86
Header file for the SelectType class template.
Header file for all forward declarations for expression class templates.
Constraint on the data type.
Header file for the EnableIf class template.
Header file for the IsStrictlyLower type trait.
Header file for the serial shim.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: TSVecTDMatMultExpr.h:293
TSVecTDMatMultExpr< VT, MT > This
Type of this TSVecTDMatMultExpr instance.
Definition: TSVecTDMatMultExpr.h:160
#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.
Removal of reference modifiers.The RemoveCV type trait removes any reference modifiers from the given...
Definition: RemoveReference.h:69
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: TSVecTDMatMultExpr.h:162
Header file for run time assertion macros.
EnableIf< IsDenseMatrix< MT1 > >::Type smpAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:98
friend EnableIf< UseAssign< VT2 > >::Type subAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
Subtraction assignment of a transpose sparse vector-transpose dense matrix multiplication to a dense ...
Definition: TSVecTDMatMultExpr.h:482
MultTrait< VRT, MRT >::Type ResultType
Result type for expression template evaluations.
Definition: TSVecTDMatMultExpr.h:161
friend EnableIf< UseAssign< VT2 > >::Type addAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
Addition assignment of a transpose sparse vector-transpose dense matrix multiplication to a dense vec...
Definition: TSVecTDMatMultExpr.h:435
Base template for the MultTrait class.
Definition: MultTrait.h:138
MT::CompositeType MCT
Composite type of the right-hand side dense matrix expression.
Definition: TSVecTDMatMultExpr.h:111
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
Header file for the reset shim.
Constraint on the data type.
TSVecTDMatMultExpr(const VT &vec, const MT &mat)
Constructor for the TSVecTDMatMultExpr class.
Definition: TSVecTDMatMultExpr.h:197
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: TSVecTDMatMultExpr.h:305
Header file for the TVecMatMultExpr base class.
Constraints on the storage order of matrix types.
RightOperand mat_
Right-hand side dense matrix of the multiplication expression.
Definition: TSVecTDMatMultExpr.h:333
Header file for the RemoveReference type trait.
Expression object for transpose sparse vector-transpose dense matrix multiplications.The TSVecTDMatMultExpr class represents the compile time expression for multiplications between transpose sparse vectors and column-major dense matrices.
Definition: Forward.h:162
LeftOperand leftOperand() const
Returns the left-hand side sparse vector operand.
Definition: TSVecTDMatMultExpr.h:271
Compile time check for strictly lower triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyLower.h:86
#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
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_TVECMATMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid vector/matrix ...
Definition: TVecMatMultExpr.h:166
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: TSVecTDMatMultExpr.h:211
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:944
Header file for the IsComputation type trait class.
friend EnableIf< UseSMPAssign< VT2 > >::Type smpAddAssign(DenseVector< VT2, true > &lhs, const TSVecTDMatMultExpr &rhs)
SMP addition assignment of a transpose sparse vector-transpose dense matrix multiplication to a dense...
Definition: TSVecTDMatMultExpr.h:644
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:118
#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:2583
#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: RowVector.h:79
Header file for the SubvectorExprTrait class template.
Header file for the IsUpper type trait.
Header file for exception macros.
Constraint on the data type.
SelectType< IsExpression< MT >::value, const MT, const MT & >::Type RightOperand
Composite type of the right-hand side sparse matrix expression.
Definition: TSVecTDMatMultExpr.h:173
SelectType< IsExpression< VT >::value, const VT, const VT & >::Type LeftOperand
Composite type of the left-hand side sparse vector expression.
Definition: TSVecTDMatMultExpr.h:170
Header file for the Size type trait.
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
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: TSVecTDMatMultExpr.h:325
Constraint on the transpose flag of vector types.
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
VT::ResultType VRT
Result type of the left-hand side sparse vector expression.
Definition: TSVecTDMatMultExpr.h:108