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>
64 #include <blaze/math/views/Check.h>
66 #include <blaze/util/Assert.h>
67 #include <blaze/util/DisableIf.h>
68 #include <blaze/util/EnableIf.h>
70 #include <blaze/util/mpl/If.h>
71 #include <blaze/util/Types.h>
73 
74 
75 namespace blaze {
76 
77 //=================================================================================================
78 //
79 // CLASS SMATDVECMULTEXPR
80 //
81 //=================================================================================================
82 
83 //*************************************************************************************************
90 template< typename MT // Type of the left-hand side sparse matrix
91  , typename VT > // Type of the right-hand side dense vector
92 class SMatDVecMultExpr
93  : public MatVecMultExpr< DenseVector< SMatDVecMultExpr<MT,VT>, false > >
94  , private Computation
95 {
96  private:
97  //**Type definitions****************************************************************************
98  using MRT = ResultType_<MT>;
99  using VRT = ResultType_<VT>;
102  //**********************************************************************************************
103 
104  //**********************************************************************************************
106  enum : bool { evaluateMatrix = RequiresEvaluation<MT>::value };
107  //**********************************************************************************************
108 
109  //**********************************************************************************************
111  enum : bool { evaluateVector = IsComputation<VT>::value || RequiresEvaluation<VT>::value };
112  //**********************************************************************************************
113 
114  //**********************************************************************************************
116 
122  enum : bool { useAssign = evaluateMatrix || evaluateVector };
123  //**********************************************************************************************
124 
125  //**********************************************************************************************
127  template< typename VT2 >
129  struct UseAssign {
130  enum : bool { value = useAssign };
131  };
133  //**********************************************************************************************
134 
135  //**********************************************************************************************
137 
141  template< typename T1 >
142  struct UseSMPAssign {
143  enum : bool { value = useAssign };
144  };
146  //**********************************************************************************************
147 
148  public:
149  //**Type definitions****************************************************************************
154  using ReturnType = const ElementType;
155 
158 
160  using LeftOperand = If_< IsExpression<MT>, const MT, const MT& >;
161 
163  using RightOperand = If_< IsExpression<VT>, const VT, const VT& >;
164 
167 
170  //**********************************************************************************************
171 
172  //**Compilation flags***************************************************************************
174  enum : bool { simdEnabled = false };
175 
177  enum : bool { smpAssignable = !evaluateMatrix && MT::smpAssignable &&
178  !evaluateVector && VT::smpAssignable };
179  //**********************************************************************************************
180 
181  //**Constructor*********************************************************************************
187  explicit inline SMatDVecMultExpr( const MT& mat, const VT& vec ) noexcept
188  : mat_( mat ) // Left-hand side sparse matrix of the multiplication expression
189  , vec_( vec ) // Right-hand side dense vector of the multiplication expression
190  {
191  BLAZE_INTERNAL_ASSERT( mat_.columns() == vec_.size(), "Invalid matrix and vector sizes" );
192  }
193  //**********************************************************************************************
194 
195  //**Subscript operator**************************************************************************
201  inline ReturnType operator[]( size_t index ) const {
202  BLAZE_INTERNAL_ASSERT( index < mat_.rows(), "Invalid vector access index" );
203  return row( mat_, index, unchecked ) * vec_;
204  }
205  //**********************************************************************************************
206 
207  //**At function*********************************************************************************
214  inline ReturnType at( size_t index ) const {
215  if( index >= mat_.rows() ) {
216  BLAZE_THROW_OUT_OF_RANGE( "Invalid vector access index" );
217  }
218  return (*this)[index];
219  }
220  //**********************************************************************************************
221 
222  //**Size function*******************************************************************************
227  inline size_t size() const noexcept {
228  return mat_.rows();
229  }
230  //**********************************************************************************************
231 
232  //**Left operand access*************************************************************************
237  inline LeftOperand leftOperand() const noexcept {
238  return mat_;
239  }
240  //**********************************************************************************************
241 
242  //**Right operand access************************************************************************
247  inline RightOperand rightOperand() const noexcept {
248  return vec_;
249  }
250  //**********************************************************************************************
251 
252  //**********************************************************************************************
258  template< typename T >
259  inline bool canAlias( const T* alias ) const noexcept {
260  return ( mat_.isAliased( alias ) || vec_.isAliased( alias ) );
261  }
262  //**********************************************************************************************
263 
264  //**********************************************************************************************
270  template< typename T >
271  inline bool isAliased( const T* alias ) const noexcept {
272  return ( mat_.isAliased( alias ) || vec_.isAliased( alias ) );
273  }
274  //**********************************************************************************************
275 
276  //**********************************************************************************************
281  inline bool isAligned() const noexcept {
282  return vec_.isAligned();
283  }
284  //**********************************************************************************************
285 
286  //**********************************************************************************************
291  inline bool canSMPAssign() const noexcept {
292  return ( size() > SMP_SMATDVECMULT_THRESHOLD );
293  }
294  //**********************************************************************************************
295 
296  private:
297  //**Member variables****************************************************************************
300  //**********************************************************************************************
301 
302  //**Assignment to dense vectors*****************************************************************
318  template< typename VT1 > // Type of the target dense vector
319  friend inline EnableIf_< UseAssign<VT1> >
320  assign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
321  {
323 
324  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
325 
326  if( rhs.mat_.columns() == 0UL ) {
327  reset( ~lhs );
328  return;
329  }
330 
331  LT A( serial( rhs.mat_ ) ); // Evaluation of the left-hand side sparse matrix operand
332  RT x( serial( rhs.vec_ ) ); // Evaluation of the right-hand side dense vector operand
333 
334  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
335  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
336  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
337  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
338 
339  assign( ~lhs, A * x );
340  }
342  //**********************************************************************************************
343 
344  //**Assignment to sparse vectors****************************************************************
360  template< typename VT1 > // Type of the target sparse vector
361  friend inline EnableIf_< UseAssign<VT1> >
362  assign( SparseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
363  {
365 
369 
370  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
371 
372  const ResultType tmp( serial( rhs ) );
373  assign( ~lhs, tmp );
374  }
376  //**********************************************************************************************
377 
378  //**Addition assignment to dense vectors********************************************************
394  template< typename VT1 > // Type of the target dense vector
395  friend inline EnableIf_< UseAssign<VT1> >
396  addAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
397  {
399 
400  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
401 
402  if( rhs.mat_.columns() == 0UL ) {
403  return;
404  }
405 
406  LT A( serial( rhs.mat_ ) ); // Evaluation of the left-hand side sparse matrix operand
407  RT x( serial( rhs.vec_ ) ); // Evaluation of the right-hand side dense vector operand
408 
409  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
410  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
411  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
412  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
413 
414  addAssign( ~lhs, A * x );
415  }
417  //**********************************************************************************************
418 
419  //**Addition assignment to sparse vectors*******************************************************
420  // No special implementation for the addition assignment to sparse vectors.
421  //**********************************************************************************************
422 
423  //**Subtraction assignment to dense vectors*****************************************************
439  template< typename VT1 > // Type of the target dense vector
440  friend inline EnableIf_< UseAssign<VT1> >
441  subAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
442  {
444 
445  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
446 
447  if( rhs.mat_.columns() == 0UL ) {
448  return;
449  }
450 
451  LT A( serial( rhs.mat_ ) ); // Evaluation of the left-hand side sparse matrix operand
452  RT x( serial( rhs.vec_ ) ); // Evaluation of the right-hand side dense vector operand
453 
454  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
455  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
456  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
457  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
458 
459  subAssign( ~lhs, A * x );
460  }
462  //**********************************************************************************************
463 
464  //**Subtraction assignment to sparse vectors****************************************************
465  // No special implementation for the subtraction assignment to sparse vectors.
466  //**********************************************************************************************
467 
468  //**Multiplication assignment to dense vectors**************************************************
484  template< typename VT1 > // Type of the target dense vector
485  friend inline EnableIf_< UseAssign<VT1> >
486  multAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
487  {
489 
493 
494  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
495 
496  const ResultType tmp( serial( rhs ) );
497  multAssign( ~lhs, tmp );
498  }
500  //**********************************************************************************************
501 
502  //**Multiplication assignment to sparse vectors*************************************************
503  // No special implementation for the multiplication assignment to sparse vectors.
504  //**********************************************************************************************
505 
506  //**Division assignment to dense vectors********************************************************
522  template< typename VT1 > // Type of the target dense vector
523  friend inline EnableIf_< UseAssign<VT1> >
524  divAssign( DenseVector<VT1,false>& lhs, const SMatDVecMultExpr& rhs )
525  {
527 
531 
532  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
533 
534  const ResultType tmp( serial( rhs ) );
535  divAssign( ~lhs, tmp );
536  }
538  //**********************************************************************************************
539 
540  //**Division assignment to sparse vectors*******************************************************
541  // No special implementation for the division assignment to sparse vectors.
542  //**********************************************************************************************
543 
544  //**SMP assignment to dense vectors*************************************************************
559  template< typename VT1 > // Type of the target dense vector
560  friend inline EnableIf_< UseSMPAssign<VT1> >
562  {
564 
565  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
566 
567  if( rhs.mat_.columns() == 0UL ) {
568  reset( ~lhs );
569  return;
570  }
571 
572  LT A( rhs.mat_ ); // Evaluation of the left-hand side sparse matrix operand
573  RT x( rhs.vec_ ); // Evaluation of the right-hand side dense vector operand
574 
575  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
576  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
577  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
578  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
579 
580  smpAssign( ~lhs, A * x );
581  }
583  //**********************************************************************************************
584 
585  //**SMP assignment to sparse vectors************************************************************
600  template< typename VT1 > // Type of the target sparse vector
601  friend inline EnableIf_< UseSMPAssign<VT1> >
603  {
605 
609 
610  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
611 
612  const ResultType tmp( rhs );
613  smpAssign( ~lhs, tmp );
614  }
616  //**********************************************************************************************
617 
618  //**SMP addition assignment to dense vectors****************************************************
633  template< typename VT1 > // Type of the target dense vector
634  friend inline EnableIf_< UseSMPAssign<VT1> >
636  {
638 
639  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
640 
641  if( rhs.mat_.columns() == 0UL ) {
642  return;
643  }
644 
645  LT A( rhs.mat_ ); // Evaluation of the left-hand side sparse matrix operand
646  RT x( rhs.vec_ ); // Evaluation of the right-hand side dense vector operand
647 
648  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
649  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
650  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
651  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
652 
653  smpAddAssign( ~lhs, A * x );
654  }
656  //**********************************************************************************************
657 
658  //**SMP addition assignment to sparse vectors***************************************************
659  // No special implementation for the SMP addition assignment to sparse vectors.
660  //**********************************************************************************************
661 
662  //**SMP subtraction assignment to dense vectors*************************************************
677  template< typename VT1 > // Type of the target dense vector
678  friend inline EnableIf_< UseSMPAssign<VT1> >
680  {
682 
683  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
684 
685  if( rhs.mat_.columns() == 0UL ) {
686  return;
687  }
688 
689  LT A( rhs.mat_ ); // Evaluation of the left-hand side sparse matrix operand
690  RT x( rhs.vec_ ); // Evaluation of the right-hand side dense vector operand
691 
692  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.mat_.rows() , "Invalid number of rows" );
693  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.mat_.columns(), "Invalid number of columns" );
694  BLAZE_INTERNAL_ASSERT( x.size() == rhs.vec_.size() , "Invalid vector size" );
695  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).size() , "Invalid vector size" );
696 
697  smpSubAssign( ~lhs, A * x );
698  }
700  //**********************************************************************************************
701 
702  //**SMP subtraction assignment to sparse vectors************************************************
703  // No special implementation for the SMP subtraction assignment to sparse vectors.
704  //**********************************************************************************************
705 
706  //**SMP multiplication assignment to dense vectors**********************************************
721  template< typename VT1 > // Type of the target dense vector
722  friend inline EnableIf_< UseSMPAssign<VT1> >
724  {
726 
730 
731  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
732 
733  const ResultType tmp( rhs );
734  smpMultAssign( ~lhs, tmp );
735  }
737  //**********************************************************************************************
738 
739  //**SMP multiplication assignment to sparse vectors*********************************************
740  // No special implementation for the SMP multiplication assignment to sparse vectors.
741  //**********************************************************************************************
742 
743  //**SMP division assignment to dense vectors****************************************************
758  template< typename VT1 > // Type of the target dense vector
759  friend inline EnableIf_< UseSMPAssign<VT1> >
761  {
763 
767 
768  BLAZE_INTERNAL_ASSERT( (~lhs).size() == rhs.size(), "Invalid vector sizes" );
769 
770  const ResultType tmp( rhs );
771  smpDivAssign( ~lhs, tmp );
772  }
774  //**********************************************************************************************
775 
776  //**SMP division assignment to sparse vectors***************************************************
777  // No special implementation for the SMP division assignment to sparse vectors.
778  //**********************************************************************************************
779 
780  //**Compile time checks*************************************************************************
788  //**********************************************************************************************
789 };
790 //*************************************************************************************************
791 
792 
793 
794 
795 //=================================================================================================
796 //
797 // GLOBAL BINARY ARITHMETIC OPERATORS
798 //
799 //=================================================================================================
800 
801 //*************************************************************************************************
832 template< typename MT // Type of the left-hand side sparse matrix
833  , typename VT > // Type of the right-hand side dense vector
834 inline decltype(auto)
835  operator*( const SparseMatrix<MT,false>& mat, const DenseVector<VT,false>& vec )
836 {
838 
840 
841  if( (~mat).columns() != (~vec).size() ) {
842  BLAZE_THROW_INVALID_ARGUMENT( "Matrix and vector sizes do not match" );
843  }
844 
845  using ReturnType = const SMatDVecMultExpr<MT,VT>;
846  return ReturnType( ~mat, ~vec );
847 }
848 //*************************************************************************************************
849 
850 
851 
852 
853 //=================================================================================================
854 //
855 // SIZE SPECIALIZATIONS
856 //
857 //=================================================================================================
858 
859 //*************************************************************************************************
861 template< typename MT, typename VT >
862 struct Size< SMatDVecMultExpr<MT,VT>, 0UL >
863  : public Size<MT,0UL>
864 {};
866 //*************************************************************************************************
867 
868 
869 
870 
871 //=================================================================================================
872 //
873 // ISALIGNED SPECIALIZATIONS
874 //
875 //=================================================================================================
876 
877 //*************************************************************************************************
879 template< typename MT, typename VT >
880 struct IsAligned< SMatDVecMultExpr<MT,VT> >
881  : public IsAligned<VT>
882 {};
884 //*************************************************************************************************
885 
886 } // namespace blaze
887 
888 #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.
Header file for the blaze::checked and blaze::unchecked instances.
Compile time check whether the given type is a computational expression template.This type trait clas...
Definition: IsComputation.h:71
Header file for basic type definitions.
ElementType_< ResultType > ElementType
Resulting element type.
Definition: SMatDVecMultExpr.h:153
#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.
CompositeType_< MT > MCT
Composite type of the left-hand side sparse matrix expression.
Definition: SMatDVecMultExpr.h:100
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:588
ReturnType at(size_t index) const
Checked access to the vector elements.
Definition: SMatDVecMultExpr.h:214
constexpr Unchecked unchecked
Global Unchecked instance.The blaze::unchecked instance is an optional token for the creation of view...
Definition: Check.h:138
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:259
typename MultTrait< T1, T2 >::Type MultTrait_
Auxiliary alias declaration for the MultTrait class template.The MultTrait_ alias declaration provide...
Definition: MultTrait.h:291
size_t size() const noexcept
Returns the current size/dimension of the vector.
Definition: SMatDVecMultExpr.h:227
Header file for the Computation base class.
SMatDVecMultExpr(const MT &mat, const VT &vec) noexcept
Constructor for the SMatDVecMultExpr class.
Definition: SMatDVecMultExpr.h:187
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:166
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:154
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:201
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.
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:71
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:98
If_< IsExpression< MT >, const MT, const MT &> LeftOperand
Composite type of the left-hand side sparse matrix expression.
Definition: SMatDVecMultExpr.h:160
Header file for the DisableIf class template.
Header file for the multiplication trait.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatDVecMultExpr.h:271
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:281
IfTrue_< evaluateVector, const VRT, VCT > RT
Type for the assignment of the right-hand side dense matrix operand.
Definition: SMatDVecMultExpr.h:169
TransposeType_< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: SMatDVecMultExpr.h:152
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:151
RightOperand vec_
Right-hand side dense vector of the multiplication expression.
Definition: SMatDVecMultExpr.h:299
#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:298
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:76
CompositeType_< VT > VCT
Composite type of the right-hand side dense vector expression.
Definition: SMatDVecMultExpr.h:101
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:506
#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:87
#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:108
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:291
#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:237
RightOperand rightOperand() const noexcept
Returns the right-hand side dense vector operand.
Definition: SMatDVecMultExpr.h:247
If_< IsExpression< VT >, const VT, const VT &> RightOperand
Composite type of the right-hand side dense vector expression.
Definition: SMatDVecMultExpr.h:163
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
decltype(auto) row(Matrix< MT, SO > &, RRAs...)
Creating a view on a specific row of the given matrix.
Definition: Row.h:131
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:816
#define BLAZE_CONSTRAINT_MUST_NOT_REQUIRE_EVALUATION(T)
Constraint on the data type.In case the given data type T requires an intermediate evaluation within ...
Definition: RequiresEvaluation.h:81
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:99
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3080
#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 vectors and matrices.The Size type trait evaluates the size of...
Definition: Size.h:80
Base class for sparse vectors.The SparseVector class is a base class for all arbitrarily sized (N-dim...
Definition: Forward.h:130
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:157
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.