SMatDVecMultExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_SMATDVECMULTEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SMATDVECMULTEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <blaze/math/Aliases.h>
51 #include <blaze/math/Exception.h>
56 #include <blaze/math/shims/Reset.h>
67 #include <blaze/util/Assert.h>
68 #include <blaze/util/DisableIf.h>
69 #include <blaze/util/EnableIf.h>
72 #include <blaze/util/mpl/If.h>
73 #include <blaze/util/Types.h>
75 
76 
77 namespace blaze {
78 
79 //=================================================================================================
80 //
81 // CLASS SMATDVECMULTEXPR
82 //
83 //=================================================================================================
84 
85 //*************************************************************************************************
92 template< typename MT // Type of the left-hand side sparse matrix
93  , typename VT > // Type of the right-hand side dense vector
94 class SMatDVecMultExpr
95  : public MatVecMultExpr< DenseVector< SMatDVecMultExpr<MT,VT>, false > >
96  , private Computation
97 {
98  private:
99  //**Type definitions****************************************************************************
104  //**********************************************************************************************
105 
106  //**********************************************************************************************
108  enum : bool { evaluateMatrix = RequiresEvaluation<MT>::value };
109  //**********************************************************************************************
110 
111  //**********************************************************************************************
113  enum : bool { evaluateVector = IsComputation<VT>::value || RequiresEvaluation<VT>::value };
114  //**********************************************************************************************
115 
116  //**********************************************************************************************
118 
124  enum : bool { useAssign = evaluateMatrix || evaluateVector };
125  //**********************************************************************************************
126 
127  //**********************************************************************************************
129  template< typename VT2 >
131  struct UseAssign {
132  enum : bool { value = useAssign };
133  };
135  //**********************************************************************************************
136 
137  //**********************************************************************************************
139 
143  template< typename T1 >
144  struct UseSMPAssign {
145  enum : bool { value = useAssign };
146  };
148  //**********************************************************************************************
149 
150  public:
151  //**Type definitions****************************************************************************
156  using ReturnType = const ElementType;
157 
160 
162  using LeftOperand = If_< IsExpression<MT>, const MT, const MT& >;
163 
165  using RightOperand = If_< IsExpression<VT>, const VT, const VT& >;
166 
169 
172  //**********************************************************************************************
173 
174  //**Compilation flags***************************************************************************
176  enum : bool { simdEnabled = false };
177 
179  enum : bool { smpAssignable = !evaluateMatrix && MT::smpAssignable &&
180  !evaluateVector && VT::smpAssignable };
181  //**********************************************************************************************
182 
183  //**Constructor*********************************************************************************
189  explicit inline SMatDVecMultExpr( const MT& mat, const VT& vec ) noexcept
190  : mat_( mat ) // Left-hand side sparse matrix of the multiplication expression
191  , vec_( vec ) // Right-hand side dense vector of the multiplication expression
192  {
193  BLAZE_INTERNAL_ASSERT( mat_.columns() == vec_.size(), "Invalid matrix and vector sizes" );
194  }
195  //**********************************************************************************************
196 
197  //**Subscript operator**************************************************************************
203  inline ReturnType operator[]( size_t index ) const {
204  BLAZE_INTERNAL_ASSERT( index < mat_.rows(), "Invalid vector access index" );
205  return row( mat_, index ) * vec_;
206  }
207  //**********************************************************************************************
208 
209  //**At function*********************************************************************************
216  inline ReturnType at( size_t index ) const {
217  if( index >= mat_.rows() ) {
218  BLAZE_THROW_OUT_OF_RANGE( "Invalid vector access index" );
219  }
220  return (*this)[index];
221  }
222  //**********************************************************************************************
223 
224  //**Size function*******************************************************************************
229  inline size_t size() const noexcept {
230  return mat_.rows();
231  }
232  //**********************************************************************************************
233 
234  //**Left operand access*************************************************************************
239  inline LeftOperand leftOperand() const noexcept {
240  return mat_;
241  }
242  //**********************************************************************************************
243 
244  //**Right operand access************************************************************************
249  inline RightOperand rightOperand() const noexcept {
250  return vec_;
251  }
252  //**********************************************************************************************
253 
254  //**********************************************************************************************
260  template< typename T >
261  inline bool canAlias( const T* alias ) const noexcept {
262  return ( mat_.isAliased( alias ) || vec_.isAliased( alias ) );
263  }
264  //**********************************************************************************************
265 
266  //**********************************************************************************************
272  template< typename T >
273  inline bool isAliased( const T* alias ) const noexcept {
274  return ( mat_.isAliased( alias ) || vec_.isAliased( alias ) );
275  }
276  //**********************************************************************************************
277 
278  //**********************************************************************************************
283  inline bool isAligned() const noexcept {
284  return vec_.isAligned();
285  }
286  //**********************************************************************************************
287 
288  //**********************************************************************************************
293  inline bool canSMPAssign() const noexcept {
294  return ( size() > SMP_SMATDVECMULT_THRESHOLD );
295  }
296  //**********************************************************************************************
297 
298  private:
299  //**Member variables****************************************************************************
302  //**********************************************************************************************
303 
304  //**Assignment to dense vectors*****************************************************************
320  template< typename VT1 > // Type of the target dense vector
321  friend inline EnableIf_< UseAssign<VT1> >
322  assign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
323  {
325 
326  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
327 
328  if( rhs.mat_.columns() == 0UL ) {
329  reset( ~lhs );
330  return;
331  }
332 
333  LT A( serial( rhs.mat_ ) ); // Evaluation of the left-hand side sparse matrix operand
334  RT x( serial( rhs.vec_ ) ); // Evaluation of the right-hand side dense vector operand
335 
336  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
337  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
338  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
339  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
340 
341  assign( ~lhs, A * x );
342  }
344  //**********************************************************************************************
345 
346  //**Assignment to sparse vectors****************************************************************
362  template< typename VT1 > // Type of the target sparse vector
363  friend inline EnableIf_< UseAssign<VT1> >
364  assign( SparseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
365  {
367 
371 
372  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
373 
374  const ResultType tmp( serial( rhs ) );
375  assign( ~lhs, tmp );
376  }
378  //**********************************************************************************************
379 
380  //**Addition assignment to dense vectors********************************************************
396  template< typename VT1 > // Type of the target dense vector
397  friend inline EnableIf_< UseAssign<VT1> >
398  addAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
399  {
401 
402  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
403 
404  if( rhs.mat_.columns() == 0UL ) {
405  return;
406  }
407 
408  LT A( serial( rhs.mat_ ) ); // Evaluation of the left-hand side sparse matrix operand
409  RT x( serial( rhs.vec_ ) ); // Evaluation of the right-hand side dense vector operand
410 
411  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
412  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
413  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
414  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
415 
416  addAssign( ~lhs, A * x );
417  }
419  //**********************************************************************************************
420 
421  //**Addition assignment to sparse vectors*******************************************************
422  // No special implementation for the addition assignment to sparse vectors.
423  //**********************************************************************************************
424 
425  //**Subtraction assignment to dense vectors*****************************************************
441  template< typename VT1 > // Type of the target dense vector
442  friend inline EnableIf_< UseAssign<VT1> >
443  subAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
444  {
446 
447  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
448 
449  if( rhs.mat_.columns() == 0UL ) {
450  return;
451  }
452 
453  LT A( serial( rhs.mat_ ) ); // Evaluation of the left-hand side sparse matrix operand
454  RT x( serial( rhs.vec_ ) ); // Evaluation of the right-hand side dense vector operand
455 
456  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
457  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
458  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
459  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
460 
461  subAssign( ~lhs, A * x );
462  }
464  //**********************************************************************************************
465 
466  //**Subtraction assignment to sparse vectors****************************************************
467  // No special implementation for the subtraction assignment to sparse vectors.
468  //**********************************************************************************************
469 
470  //**Multiplication assignment to dense vectors**************************************************
486  template< typename VT1 > // Type of the target dense vector
487  friend inline EnableIf_< UseAssign<VT1> >
488  multAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
489  {
491 
495 
496  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
497 
498  const ResultType tmp( serial( rhs ) );
499  multAssign( ~lhs, tmp );
500  }
502  //**********************************************************************************************
503 
504  //**Multiplication assignment to sparse vectors*************************************************
505  // No special implementation for the multiplication assignment to sparse vectors.
506  //**********************************************************************************************
507 
508  //**Division assignment to dense vectors********************************************************
524  template< typename VT1 > // Type of the target dense vector
525  friend inline EnableIf_< UseAssign<VT1> >
526  divAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
527  {
529 
533 
534  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
535 
536  const ResultType tmp( serial( rhs ) );
537  divAssign( ~lhs, tmp );
538  }
540  //**********************************************************************************************
541 
542  //**Division assignment to sparse vectors*******************************************************
543  // No special implementation for the division assignment to sparse vectors.
544  //**********************************************************************************************
545 
546  //**SMP assignment to dense vectors*************************************************************
561  template< typename VT1 > // Type of the target dense vector
562  friend inline EnableIf_< UseSMPAssign<VT1> >
564  {
566 
567  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
568 
569  if( rhs.mat_.columns() == 0UL ) {
570  reset( ~lhs );
571  return;
572  }
573 
574  LT A( rhs.mat_ ); // Evaluation of the left-hand side sparse matrix operand
575  RT x( rhs.vec_ ); // Evaluation of the right-hand side dense vector operand
576 
577  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
578  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
579  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
580  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
581 
582  smpAssign( ~lhs, A * x );
583  }
585  //**********************************************************************************************
586 
587  //**SMP assignment to sparse vectors************************************************************
602  template< typename VT1 > // Type of the target sparse vector
603  friend inline EnableIf_< UseSMPAssign<VT1> >
605  {
607 
611 
612  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
613 
614  const ResultType tmp( rhs );
615  smpAssign( ~lhs, tmp );
616  }
618  //**********************************************************************************************
619 
620  //**SMP addition assignment to dense vectors****************************************************
635  template< typename VT1 > // Type of the target dense vector
636  friend inline EnableIf_< UseSMPAssign<VT1> >
638  {
640 
641  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
642 
643  if( rhs.mat_.columns() == 0UL ) {
644  return;
645  }
646 
647  LT A( rhs.mat_ ); // Evaluation of the left-hand side sparse matrix operand
648  RT x( rhs.vec_ ); // Evaluation of the right-hand side dense vector operand
649 
650  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
651  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
652  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
653  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
654 
655  smpAddAssign( ~lhs, A * x );
656  }
658  //**********************************************************************************************
659 
660  //**SMP addition assignment to sparse vectors***************************************************
661  // No special implementation for the SMP addition assignment to sparse vectors.
662  //**********************************************************************************************
663 
664  //**SMP subtraction assignment to dense vectors*************************************************
679  template< typename VT1 > // Type of the target dense vector
680  friend inline EnableIf_< UseSMPAssign<VT1> >
682  {
684 
685  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
686 
687  if( rhs.mat_.columns() == 0UL ) {
688  return;
689  }
690 
691  LT A( rhs.mat_ ); // Evaluation of the left-hand side sparse matrix operand
692  RT x( rhs.vec_ ); // Evaluation of the right-hand side dense vector operand
693 
694  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
695  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
696  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
697  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
698 
699  smpSubAssign( ~lhs, A * x );
700  }
702  //**********************************************************************************************
703 
704  //**SMP subtraction assignment to sparse vectors************************************************
705  // No special implementation for the SMP subtraction assignment to sparse vectors.
706  //**********************************************************************************************
707 
708  //**SMP multiplication assignment to dense vectors**********************************************
723  template< typename VT1 > // Type of the target dense vector
724  friend inline EnableIf_< UseSMPAssign<VT1> >
726  {
728 
732 
733  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
734 
735  const ResultType tmp( rhs );
736  smpMultAssign( ~lhs, tmp );
737  }
739  //**********************************************************************************************
740 
741  //**SMP multiplication assignment to sparse vectors*********************************************
742  // No special implementation for the SMP multiplication assignment to sparse vectors.
743  //**********************************************************************************************
744 
745  //**SMP division assignment to dense vectors****************************************************
760  template< typename VT1 > // Type of the target dense vector
761  friend inline EnableIf_< UseSMPAssign<VT1> >
763  {
765 
769 
770  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
771 
772  const ResultType tmp( rhs );
773  smpDivAssign( ~lhs, tmp );
774  }
776  //**********************************************************************************************
777 
778  //**SMP division assignment to sparse vectors***************************************************
779  // No special implementation for the SMP division assignment to sparse vectors.
780  //**********************************************************************************************
781 
782  //**Compile time checks*************************************************************************
790  //**********************************************************************************************
791 };
792 //*************************************************************************************************
793 
794 
795 
796 
797 //=================================================================================================
798 //
799 // GLOBAL BINARY ARITHMETIC OPERATORS
800 //
801 //=================================================================================================
802 
803 //*************************************************************************************************
834 template< typename MT // Type of the left-hand side sparse matrix
835  , typename VT > // Type of the right-hand side dense vector
836 inline decltype(auto)
837  operator*( const SparseMatrix<MT,false>& mat, const DenseVector<VT,false>& vec )
838 {
840 
842 
843  if( (~mat).columns() != (~vec).size() ) {
844  BLAZE_THROW_INVALID_ARGUMENT( "Matrix and vector sizes do not match" );
845  }
846 
847  using ReturnType = const SMatDVecMultExpr<MT,VT>;
848  return ReturnType( ~mat, ~vec );
849 }
850 //*************************************************************************************************
851 
852 
853 
854 
855 //=================================================================================================
856 //
857 // SIZE SPECIALIZATIONS
858 //
859 //=================================================================================================
860 
861 //*************************************************************************************************
863 template< typename MT, typename VT >
864 struct Size< SMatDVecMultExpr<MT,VT> >
865  : public Rows<MT>
866 {};
868 //*************************************************************************************************
869 
870 
871 
872 
873 //=================================================================================================
874 //
875 // ISALIGNED SPECIALIZATIONS
876 //
877 //=================================================================================================
878 
879 //*************************************************************************************************
881 template< typename MT, typename VT >
882 struct IsAligned< SMatDVecMultExpr<MT,VT> >
883  : public BoolConstant< IsAligned<VT>::value >
884 {};
886 //*************************************************************************************************
887 
888 } // namespace blaze
889 
890 #endif
#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
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:72
Header file for the Rows type trait.
Header file for basic type definitions.
ElementType_< ResultType > ElementType
Resulting element type.
Definition: SMatDVecMultExpr.h:155
#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: ColumnVector.h:61
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
Header file for the serial shim.
Header file for the IsDiagonal type trait.
Generic wrapper for a compile time constant integral value.The IntegralConstant class template repres...
Definition: IntegralConstant.h:71
CompositeType_< MT > MCT
Composite type of the left-hand side sparse matrix expression.
Definition: SMatDVecMultExpr.h:102
Expression object for sparse matrix-dense vector multiplications.The SMatDVecMultExpr class represent...
Definition: Forward.h:110
void reset(const DiagonalProxy< MT > &proxy)
Resetting the represented element to the default initial values.
Definition: DiagonalProxy.h:560
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: SMatDVecMultExpr.h:216
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 DenseVector base class.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SMatDVecMultExpr.h:261
typename MultTrait< T1, T2 >::Type MultTrait_
Auxiliary alias declaration for the MultTrait class template.The MultTrait_ alias declaration provide...
Definition: MultTrait.h:250
size_t size() const noexcept
Returns the current size/dimension of the vector.
Definition: SMatDVecMultExpr.h:229
Header file for the Computation base class.
SMatDVecMultExpr(const MT &mat, const VT &vec) noexcept
Constructor for the SMatDVecMultExpr class.
Definition: SMatDVecMultExpr.h:189
Constraints on the storage order of matrix types.
Header file for the RequiresEvaluation type trait.
IfTrue_< evaluateMatrix, const MRT, MCT > LT
Type for the assignment of the left-hand side sparse matrix operand.
Definition: SMatDVecMultExpr.h:168
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
const ElementType ReturnType
Return type for expression template evaluations.
Definition: SMatDVecMultExpr.h:156
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:129
ReturnType operator[](size_t index) const
Subscript operator for the direct access to the vector elements.
Definition: SMatDVecMultExpr.h:203
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
Constraint on the transpose flag of vector types.
Row< MT > row(Matrix< MT, SO > &matrix, size_t index)
Creating a view on a specific row of the given matrix.
Definition: Row.h:124
Compile time check for the alignment of data types.This type trait tests whether the given data type ...
Definition: IsAligned.h:87
Constraint on the data type.
Compile time check to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:72
typename T::CompositeType CompositeType_
Alias declaration for nested CompositeType type definitions.The CompositeType_ alias declaration prov...
Definition: Aliases.h:83
ResultType_< MT > MRT
Result type of the left-hand side sparse matrix expression.
Definition: SMatDVecMultExpr.h:100
If_< IsExpression< MT >, const MT, const MT &> LeftOperand
Composite type of the left-hand side sparse matrix expression.
Definition: SMatDVecMultExpr.h:162
Header file for the DisableIf class template.
Header file for the multiplication trait.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatDVecMultExpr.h:273
Header file for the If class template.
bool isAligned() const noexcept
Returns whether the operands of the expression are properly aligned in memory.
Definition: SMatDVecMultExpr.h:283
IfTrue_< evaluateVector, const VRT, VCT > RT
Type for the assignment of the right-hand side dense matrix operand.
Definition: SMatDVecMultExpr.h:171
TransposeType_< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: SMatDVecMultExpr.h:154
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
MultTrait_< MRT, VRT > ResultType
Result type for expression template evaluations.
Definition: SMatDVecMultExpr.h:153
RightOperand vec_
Right-hand side dense vector of the multiplication expression.
Definition: SMatDVecMultExpr.h:301
#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
LeftOperand mat_
Left-hand side sparse matrix of the multiplication expression.
Definition: SMatDVecMultExpr.h:300
typename T::ElementType ElementType_
Alias declaration for nested ElementType type definitions.The ElementType_ alias declaration provides...
Definition: Aliases.h:163
Base class for N-dimensional dense vectors.The DenseVector class is a base class for all arbitrarily ...
Definition: DenseVector.h:70
CompositeType_< VT > VCT
Composite type of the right-hand side dense vector expression.
Definition: SMatDVecMultExpr.h:103
Header file for the IsAligned type trait.
BLAZE_ALWAYS_INLINE size_t columns(const Matrix< MT, SO > &matrix) noexcept
Returns the current number of columns of the matrix.
Definition: Matrix.h:340
#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:88
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_MATVECMULTEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid matrix/vector ...
Definition: MatVecMultExpr.h:109
Constraint on the data type.
Header file for the exception macros of the math module.
Constraint on the data type.
Header file for all forward declarations for expression class templates.
Header file for the EnableIf class template.
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: SMatDVecMultExpr.h:293
#define BLAZE_CONSTRAINT_MUST_BE_ROW_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a row-major dense or sparse matrix t...
Definition: RowMajorMatrix.h:61
Header file for run time assertion macros.
LeftOperand leftOperand() const noexcept
Returns the left-hand side sparse matrix operand.
Definition: SMatDVecMultExpr.h:239
RightOperand rightOperand() const noexcept
Returns the right-hand side dense vector operand.
Definition: SMatDVecMultExpr.h:249
If_< IsExpression< VT >, const VT, const VT &> RightOperand
Composite type of the right-hand side dense vector expression.
Definition: SMatDVecMultExpr.h:165
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
Header file for the reset shim.
#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
decltype(auto) serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:819
#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
Header file for the RemoveReference type trait.
typename EnableIf< Condition, T >::Type EnableIf_
Auxiliary alias declaration for the EnableIf class template.The EnableIf_ alias declaration provides ...
Definition: EnableIf.h:224
ResultType_< VT > VRT
Result type of the right-hand side dense vector expression.
Definition: SMatDVecMultExpr.h:101
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3082
#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
Header file for the IsComputation type trait class.
Compile time evaluation of the size of a vector.The Size type trait evaluates the size of the given v...
Definition: Size.h:74
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:130
Header file for the IntegralConstant class template.
Compile time evaluation of the number of rows of a matrix.The Rows type trait evaluates the number of...
Definition: Rows.h:75
Constraint on the data type.
typename T::TransposeType TransposeType_
Alias declaration for nested TransposeType type definitions.The TransposeType_ alias declaration prov...
Definition: Aliases.h:423
Header file for the MatVecMultExpr base class.
Constraint on the data type.
IfTrue_< useAssign, const ResultType, const SMatDVecMultExpr &> CompositeType
Data type for composite expression templates.
Definition: SMatDVecMultExpr.h:159
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
#define BLAZE_CONSTRAINT_MUST_BE_SPARSE_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a sparse, N-dimensional matrix type...
Definition: SparseMatrix.h:61
Header file for the IsExpression type trait class.
Header file for the function trace functionality.