SMatEvalExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_SMATEVALEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SMATEVALEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <blaze/math/Aliases.h>
46 #include <blaze/math/Exception.h>
62 #include <blaze/util/Assert.h>
65 #include <blaze/util/InvalidType.h>
66 #include <blaze/util/mpl/If.h>
67 #include <blaze/util/Types.h>
68 
69 
70 namespace blaze {
71 
72 //=================================================================================================
73 //
74 // CLASS SMATEVALEXPR
75 //
76 //=================================================================================================
77 
78 //*************************************************************************************************
85 template< typename MT // Type of the sparse matrix
86  , bool SO > // Storage order
87 class SMatEvalExpr
88  : public MatEvalExpr< SparseMatrix< SMatEvalExpr<MT,SO>, SO > >
89  , private Computation
90 {
91  public:
92  //**Type definitions****************************************************************************
99 
101  using CompositeType = const ResultType;
102 
104  using Operand = If_< IsExpression<MT>, const MT, const MT& >;
105  //**********************************************************************************************
106 
107  //**Compilation flags***************************************************************************
109  enum : bool { smpAssignable = MT::smpAssignable };
110  //**********************************************************************************************
111 
112  //**Constructor*********************************************************************************
117  explicit inline SMatEvalExpr( const MT& sm ) noexcept
118  : sm_( sm ) // Sparse matrix of the evaluation expression
119  {}
120  //**********************************************************************************************
121 
122  //**Access operator*****************************************************************************
129  inline ReturnType operator()( size_t i, size_t j ) const {
130  BLAZE_INTERNAL_ASSERT( i < sm_.rows() , "Invalid row access index" );
131  BLAZE_INTERNAL_ASSERT( j < sm_.columns(), "Invalid column access index" );
132  return sm_(i,j);
133  }
134  //**********************************************************************************************
135 
136  //**At function*********************************************************************************
144  inline ReturnType at( size_t i, size_t j ) const {
145  if( i >= sm_.rows() ) {
146  BLAZE_THROW_OUT_OF_RANGE( "Invalid row access index" );
147  }
148  if( j >= sm_.columns() ) {
149  BLAZE_THROW_OUT_OF_RANGE( "Invalid column access index" );
150  }
151  return (*this)(i,j);
152  }
153  //**********************************************************************************************
154 
155  //**Rows function*******************************************************************************
160  inline size_t rows() const noexcept {
161  return sm_.rows();
162  }
163  //**********************************************************************************************
164 
165  //**Columns function****************************************************************************
170  inline size_t columns() const noexcept {
171  return sm_.columns();
172  }
173  //**********************************************************************************************
174 
175  //**NonZeros function***************************************************************************
180  inline size_t nonZeros() const {
181  return sm_.nonZeros();
182  }
183  //**********************************************************************************************
184 
185  //**NonZeros function***************************************************************************
191  inline size_t nonZeros( size_t i ) const {
192  return sm_.nonZeros(i);
193  }
194  //**********************************************************************************************
195 
196  //**Operand access******************************************************************************
201  inline Operand operand() const noexcept {
202  return sm_;
203  }
204  //**********************************************************************************************
205 
206  //**********************************************************************************************
212  template< typename T >
213  inline bool canAlias( const T* alias ) const noexcept {
214  return sm_.canAlias( alias );
215  }
216  //**********************************************************************************************
217 
218  //**********************************************************************************************
224  template< typename T >
225  inline bool isAliased( const T* alias ) const noexcept {
226  return sm_.isAliased( alias );
227  }
228  //**********************************************************************************************
229 
230  //**********************************************************************************************
235  inline bool canSMPAssign() const noexcept {
236  return sm_.canSMPAssign();
237  }
238  //**********************************************************************************************
239 
240  private:
241  //**Member variables****************************************************************************
243  //**********************************************************************************************
244 
245  //**Assignment to dense matrices****************************************************************
257  template< typename MT2 // Type of the target dense matrix
258  , bool SO2 > // Storage order of the target dense matrix
259  friend inline void assign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
260  {
262 
263  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
264  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
265 
266  assign( ~lhs, rhs.sm_ );
267  }
269  //**********************************************************************************************
270 
271  //**Assignment to sparse matrices***************************************************************
283  template< typename MT2 // Type of the target sparse matrix
284  , bool SO2 > // Storage order of the target sparse matrix
285  friend inline void assign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
286  {
288 
289  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
290  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
291 
292  assign( ~lhs, rhs.sm_ );
293  }
295  //**********************************************************************************************
296 
297  //**Addition assignment to dense matrices*******************************************************
309  template< typename MT2 // Type of the target dense matrix
310  , bool SO2 > // Storage order of the target dense matrix
311  friend inline void addAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
312  {
314 
315  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
316  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
317 
318  addAssign( ~lhs, rhs.sm_ );
319  }
321  //**********************************************************************************************
322 
323  //**Addition assignment to sparse matrices******************************************************
335  template< typename MT2 // Type of the target sparse matrix
336  , bool SO2 > // Storage order of the target sparse matrix
337  friend inline void addAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
338  {
340 
341  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
342  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
343 
344  addAssign( ~lhs, rhs.sm_ );
345  }
347  //**********************************************************************************************
348 
349  //**Subtraction assignment to dense matrices****************************************************
361  template< typename MT2 // Type of the target dense matrix
362  , bool SO2 > // Storage order of the target dense matrix
363  friend inline void subAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
364  {
366 
367  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
368  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
369 
370  subAssign( ~lhs, rhs.sm_ );
371  }
373  //**********************************************************************************************
374 
375  //**Subtraction assignment to sparse matrices***************************************************
387  template< typename MT2 // Type of the target sparse matrix
388  , bool SO2 > // Storage order of the target sparse matrix
389  friend inline void subAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
390  {
392 
393  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
394  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
395 
396  subAssign( ~lhs, rhs.sm_ );
397  }
399  //**********************************************************************************************
400 
401  //**Schur product assignment to dense matrices**************************************************
413  template< typename MT2 // Type of the target dense matrix
414  , bool SO2 > // Storage order of the target dense matrix
415  friend inline void schurAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
416  {
418 
419  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
420  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
421 
422  schurAssign( ~lhs, rhs.sm_ );
423  }
425  //**********************************************************************************************
426 
427  //**Schur product assignment to sparse matrices*************************************************
439  template< typename MT2 // Type of the target sparse matrix
440  , bool SO2 > // Storage order of the target sparse matrix
441  friend inline void schurAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
442  {
444 
445  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
446  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
447 
448  schurAssign( ~lhs, rhs.sm_ );
449  }
451  //**********************************************************************************************
452 
453  //**Multiplication assignment to dense matrices*************************************************
465  template< typename MT2 // Type of the target dense matrix
466  , bool SO2 > // Storage order of the target dense matrix
467  friend inline void multAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
468  {
470 
471  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
472  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
473 
474  multAssign( ~lhs, rhs.sm_ );
475  }
477  //**********************************************************************************************
478 
479  //**Multiplication assignment to sparse matrices************************************************
491  template< typename MT2 // Type of the target sparse matrix
492  , bool SO2 > // Storage order of the target sparse matrix
493  friend inline void multAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
494  {
496 
497  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
498  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
499 
500  multAssign( ~lhs, rhs.sm_ );
501  }
503  //**********************************************************************************************
504 
505  //**SMP assignment to dense matrices************************************************************
517  template< typename MT2 // Type of the target dense matrix
518  , bool SO2 > // Storage order of the target dense matrix
519  friend inline void smpAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
520  {
522 
523  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
524  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
525 
526  smpAssign( ~lhs, rhs.sm_ );
527  }
529  //**********************************************************************************************
530 
531  //**SMP assignment to sparse matrices***********************************************************
543  template< typename MT2 // Type of the target sparse matrix
544  , bool SO2 > // Storage order of the target sparse matrix
545  friend inline void smpAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
546  {
548 
549  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
550  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
551 
552  smpAssign( ~lhs, rhs.sm_ );
553  }
555  //**********************************************************************************************
556 
557  //**SMP addition assignment to dense matrices***************************************************
569  template< typename MT2 // Type of the target dense matrix
570  , bool SO2 > // Storage order of the target dense matrix
571  friend inline void smpAddAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
572  {
574 
575  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
576  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
577 
578  smpAddAssign( ~lhs, rhs.sm_ );
579  }
581  //**********************************************************************************************
582 
583  //**SMP addition assignment to sparse matrices**************************************************
595  template< typename MT2 // Type of the target sparse matrix
596  , bool SO2 > // Storage order of the target sparse matrix
597  friend inline void smpAddAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
598  {
600 
601  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
602  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
603 
604  smpAddAssign( ~lhs, rhs.sm_ );
605  }
607  //**********************************************************************************************
608 
609  //**SMP subtraction assignment to dense matrices************************************************
621  template< typename MT2 // Type of the target dense matrix
622  , bool SO2 > // Storage order of the target dense matrix
623  friend inline void smpSubAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
624  {
626 
627  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
628  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
629 
630  smpSubAssign( ~lhs, rhs.sm_ );
631  }
633  //**********************************************************************************************
634 
635  //**SMP subtraction assignment to sparse matrices***********************************************
647  template< typename MT2 // Type of the target sparse matrix
648  , bool SO2 > // Storage order of the target sparse matrix
649  friend inline void smpSubAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
650  {
652 
653  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
654  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
655 
656  smpSubAssign( ~lhs, rhs.sm_ );
657  }
659  //**********************************************************************************************
660 
661  //**SMP Schur product assignment to dense matrices**********************************************
674  template< typename MT2 // Type of the target dense matrix
675  , bool SO2 > // Storage order of the target dense matrix
676  friend inline void smpSchurAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
677  {
679 
680  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
681  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
682 
683  smpSchurAssign( ~lhs, rhs.sm_ );
684  }
686  //**********************************************************************************************
687 
688  //**SMP Schur product assignment to sparse matrices*********************************************
701  template< typename MT2 // Type of the target sparse matrix
702  , bool SO2 > // Storage order of the target sparse matrix
703  friend inline void smpSchurAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
704  {
706 
707  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
708  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
709 
710  smpSchurAssign( ~lhs, rhs.sm_ );
711  }
713  //**********************************************************************************************
714 
715  //**SMP multiplication assignment to dense matrices*********************************************
728  template< typename MT2 // Type of the target dense matrix
729  , bool SO2 > // Storage order of the target dense matrix
730  friend inline void smpMultAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
731  {
733 
734  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
735  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
736 
737  smpMultAssign( ~lhs, rhs.sm_ );
738  }
740  //**********************************************************************************************
741 
742  //**SMP multiplication assignment to sparse matrices********************************************
755  template< typename MT2 // Type of the target sparse matrix
756  , bool SO2 > // Storage order of the target sparse matrix
757  friend inline void smpMultAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
758  {
760 
761  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
762  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
763 
764  smpMultAssign( ~lhs, rhs.sm_ );
765  }
767  //**********************************************************************************************
768 
769  //**Compile time checks*************************************************************************
774  //**********************************************************************************************
775 };
776 //*************************************************************************************************
777 
778 
779 
780 
781 //=================================================================================================
782 //
783 // GLOBAL FUNCTIONS
784 //
785 //=================================================================================================
786 
787 //*************************************************************************************************
804 template< typename MT // Type of the sparse matrix
805  , bool SO > // Storage order
806 inline decltype(auto) eval( const SparseMatrix<MT,SO>& sm )
807 {
809 
810  using ReturnType = const SMatEvalExpr<MT,SO>;
811  return ReturnType( ~sm );
812 }
813 //*************************************************************************************************
814 
815 
816 
817 
818 //=================================================================================================
819 //
820 // GLOBAL RESTRUCTURING FUNCTIONS
821 //
822 //=================================================================================================
823 
824 //*************************************************************************************************
835 template< typename MT // Type of the sparse matrix
836  , bool SO > // Storage order
837 inline decltype(auto) eval( const SMatEvalExpr<MT,SO>& sm )
838 {
839  return sm;
840 }
842 //*************************************************************************************************
843 
844 
845 
846 
847 //=================================================================================================
848 //
849 // ROWS SPECIALIZATIONS
850 //
851 //=================================================================================================
852 
853 //*************************************************************************************************
855 template< typename MT, bool SO >
856 struct Rows< SMatEvalExpr<MT,SO> >
857  : public Rows<MT>
858 {};
860 //*************************************************************************************************
861 
862 
863 
864 
865 //=================================================================================================
866 //
867 // COLUMNS SPECIALIZATIONS
868 //
869 //=================================================================================================
870 
871 //*************************************************************************************************
873 template< typename MT, bool SO >
874 struct Columns< SMatEvalExpr<MT,SO> >
875  : public Columns<MT>
876 {};
878 //*************************************************************************************************
879 
880 
881 
882 
883 //=================================================================================================
884 //
885 // ISSYMMETRIC SPECIALIZATIONS
886 //
887 //=================================================================================================
888 
889 //*************************************************************************************************
891 template< typename MT, bool SO >
892 struct IsSymmetric< SMatEvalExpr<MT,SO> >
893  : public BoolConstant< IsSymmetric<MT>::value >
894 {};
896 //*************************************************************************************************
897 
898 
899 
900 
901 //=================================================================================================
902 //
903 // ISHERMITIAN SPECIALIZATIONS
904 //
905 //=================================================================================================
906 
907 //*************************************************************************************************
909 template< typename MT, bool SO >
910 struct IsHermitian< SMatEvalExpr<MT,SO> >
911  : public BoolConstant< IsHermitian<MT>::value >
912 {};
914 //*************************************************************************************************
915 
916 
917 
918 
919 //=================================================================================================
920 //
921 // ISLOWER SPECIALIZATIONS
922 //
923 //=================================================================================================
924 
925 //*************************************************************************************************
927 template< typename MT, bool SO >
928 struct IsLower< SMatEvalExpr<MT,SO> >
929  : public BoolConstant< IsLower<MT>::value >
930 {};
932 //*************************************************************************************************
933 
934 
935 
936 
937 //=================================================================================================
938 //
939 // ISUNILOWER SPECIALIZATIONS
940 //
941 //=================================================================================================
942 
943 //*************************************************************************************************
945 template< typename MT, bool SO >
946 struct IsUniLower< SMatEvalExpr<MT,SO> >
947  : public BoolConstant< IsUniLower<MT>::value >
948 {};
950 //*************************************************************************************************
951 
952 
953 
954 
955 //=================================================================================================
956 //
957 // ISSTRICTLYLOWER SPECIALIZATIONS
958 //
959 //=================================================================================================
960 
961 //*************************************************************************************************
963 template< typename MT, bool SO >
964 struct IsStrictlyLower< SMatEvalExpr<MT,SO> >
965  : public BoolConstant< IsStrictlyLower<MT>::value >
966 {};
968 //*************************************************************************************************
969 
970 
971 
972 
973 //=================================================================================================
974 //
975 // ISUPPER SPECIALIZATIONS
976 //
977 //=================================================================================================
978 
979 //*************************************************************************************************
981 template< typename MT, bool SO >
982 struct IsUpper< SMatEvalExpr<MT,SO> >
983  : public BoolConstant< IsUpper<MT>::value >
984 {};
986 //*************************************************************************************************
987 
988 
989 
990 
991 //=================================================================================================
992 //
993 // ISUNIUPPER SPECIALIZATIONS
994 //
995 //=================================================================================================
996 
997 //*************************************************************************************************
999 template< typename MT, bool SO >
1000 struct IsUniUpper< SMatEvalExpr<MT,SO> >
1001  : public BoolConstant< IsUniUpper<MT>::value >
1002 {};
1004 //*************************************************************************************************
1005 
1006 
1007 
1008 
1009 //=================================================================================================
1010 //
1011 // ISSTRICTLYUPPER SPECIALIZATIONS
1012 //
1013 //=================================================================================================
1014 
1015 //*************************************************************************************************
1017 template< typename MT, bool SO >
1018 struct IsStrictlyUpper< SMatEvalExpr<MT,SO> >
1019  : public BoolConstant< IsStrictlyUpper<MT>::value >
1020 {};
1022 //*************************************************************************************************
1023 
1024 } // namespace blaze
1025 
1026 #endif
Header file for auxiliary alias declarations.
Header file for the Rows type trait.
Header file for the IsUniUpper type trait.
EnableIf_< IsDenseMatrix< MT1 > > smpSchurAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP Schur product assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:196
Header file for basic type definitions.
const ResultType CompositeType
Data type for composite expression templates.
Definition: SMatEvalExpr.h:101
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: SMatEvalExpr.h:170
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
#define BLAZE_CONSTRAINT_MUST_BE_MATRIX_WITH_STORAGE_ORDER(T, SO)
Constraint on the data type.In case the given data type T is not a dense or sparse matrix type and in...
Definition: StorageOrder.h:63
Generic wrapper for a compile time constant integral value.The IntegralConstant class template repres...
Definition: IntegralConstant.h:71
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
Compile time check for lower triangular matrices.This type trait tests whether or not the given templ...
Definition: IsLower.h:88
Header file for the MatEvalExpr base class.
Header file for the Computation base class.
Compile time check for upper triangular matrices.This type trait tests whether or not the given templ...
Definition: IsUpper.h:88
Header file for the IsUniLower type trait.
ElementType_< MT > ElementType
Resulting element type.
Definition: SMatEvalExpr.h:97
typename T::ResultType ResultType_
Alias declaration for nested ResultType type definitions.The ResultType_ alias declaration provides a...
Definition: Aliases.h:343
EnableIf_< IsDenseMatrix< MT1 > > smpAddAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP addition assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:133
ReturnType_< MT > ReturnType
Return type for expression template evaluations.
Definition: SMatEvalExpr.h:98
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:78
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:129
typename T::ReturnType ReturnType_
Alias declaration for nested ReturnType type definitions.The ReturnType_ alias declaration provides a...
Definition: Aliases.h:363
Header file for the SparseMatrix base class.
SMatEvalExpr(const MT &sm) noexcept
Constructor for the SMatEvalExpr class.
Definition: SMatEvalExpr.h:117
Constraint on the data type.
Compile time check for upper unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniUpper.h:86
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
Header file for the If class template.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SMatEvalExpr.h:213
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatEvalExpr.h:225
If_< IsExpression< MT >, const MT, const MT &> Operand
Composite data type of the sparse matrix expression.
Definition: SMatEvalExpr.h:104
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: SMatEvalExpr.h:144
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
#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
Header file for the Columns type trait.
typename T::ElementType ElementType_
Alias declaration for nested ElementType type definitions.The ElementType_ alias declaration provides...
Definition: Aliases.h:163
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: SMatEvalExpr.h:235
Header file for the IsLower type trait.
Operand operand() const noexcept
Returns the sparse matrix operand.
Definition: SMatEvalExpr.h:201
decltype(auto) eval(const DenseMatrix< MT, SO > &dm)
Forces the evaluation of the given dense matrix expression dm.
Definition: DMatEvalExpr.h:797
Constraints on the storage order of matrix types.
Compile time check for symmetric matrices.This type trait tests whether or not the given template par...
Definition: IsSymmetric.h:85
Header file for the exception macros of the math module.
Compile time check for strictly upper triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyUpper.h:86
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row.
Definition: SMatEvalExpr.h:191
Header file for all forward declarations for expression class templates.
Header file for the IsStrictlyLower type trait.
Compile time check for lower unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniLower.h:86
Operand sm_
Sparse matrix of the evaluation expression.
Definition: SMatEvalExpr.h:242
OppositeType_< MT > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatEvalExpr.h:95
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SMatEvalExpr.h:180
Header file for run time assertion macros.
Utility type for generic codes.
ResultType_< MT > ResultType
Result type for expression template evaluations.
Definition: SMatEvalExpr.h:94
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
#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
Compile time check for Hermitian matrices.This type trait tests whether or not the given template par...
Definition: IsHermitian.h:85
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SMatEvalExpr.h:129
typename T::OppositeType OppositeType_
Alias declaration for nested OppositeType type definitions.The OppositeType_ alias declaration provid...
Definition: Aliases.h:263
Compile time check for strictly lower triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyLower.h:86
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3082
TransposeType_< MT > TransposeType
Transpose type for expression template evaluations.
Definition: SMatEvalExpr.h:96
Header file for the IntegralConstant class template.
Compile time evaluation of the number of columns of a matrix.The Columns type trait evaluates the num...
Definition: Columns.h:75
Expression object for the forced evaluation of sparse matrices.The SMatEvalExpr class represents the ...
Definition: Forward.h:111
Compile time evaluation of the number of rows of a matrix.The Rows type trait evaluates the number of...
Definition: Rows.h:75
typename T::TransposeType TransposeType_
Alias declaration for nested TransposeType type definitions.The TransposeType_ alias declaration prov...
Definition: Aliases.h:423
Header file for the IsUpper type trait.
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: SMatEvalExpr.h:160
Header file for the IsHermitian type trait.
#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.