SMatSerialExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_SMATSERIALEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SMATSERIALEXPR_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 SMATSERIALEXPR
75 //
76 //=================================================================================================
77 
78 //*************************************************************************************************
85 template< typename MT // Type of the sparse matrix
86  , bool SO > // Storage order
87 class SMatSerialExpr
88  : public MatSerialExpr< SparseMatrix< SMatSerialExpr<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 SMatSerialExpr( const MT& sm ) noexcept
118  : sm_( sm ) // Sparse matrix of the serial 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  //**Conversion operator*************************************************************************
211  inline operator Operand() const noexcept {
212  return sm_;
213  }
214  //**********************************************************************************************
215 
216  //**********************************************************************************************
222  template< typename T >
223  inline bool canAlias( const T* alias ) const noexcept {
224  return sm_.canAlias( alias );
225  }
226  //**********************************************************************************************
227 
228  //**********************************************************************************************
234  template< typename T >
235  inline bool isAliased( const T* alias ) const noexcept {
236  return sm_.isAliased( alias );
237  }
238  //**********************************************************************************************
239 
240  //**********************************************************************************************
245  inline bool canSMPAssign() const noexcept {
246  return sm_.canSMPAssign();
247  }
248  //**********************************************************************************************
249 
250  private:
251  //**Member variables****************************************************************************
253  //**********************************************************************************************
254 
255  //**Assignment to dense matrices****************************************************************
267  template< typename MT2 // Type of the target dense matrix
268  , bool SO2 > // Storage order of the target dense matrix
269  friend inline void assign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
270  {
272 
273  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
274  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
275 
276  assign( ~lhs, rhs.sm_ );
277  }
279  //**********************************************************************************************
280 
281  //**Assignment to sparse matrices***************************************************************
293  template< typename MT2 // Type of the target sparse matrix
294  , bool SO2 > // Storage order of the target sparse matrix
295  friend inline void assign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
296  {
298 
299  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
300  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
301 
302  assign( ~lhs, rhs.sm_ );
303  }
305  //**********************************************************************************************
306 
307  //**Addition assignment to dense matrices*******************************************************
319  template< typename MT2 // Type of the target dense matrix
320  , bool SO2 > // Storage order of the target dense matrix
321  friend inline void addAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
322  {
324 
325  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
326  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
327 
328  addAssign( ~lhs, rhs.sm_ );
329  }
331  //**********************************************************************************************
332 
333  //**Addition assignment to sparse matrices******************************************************
345  template< typename MT2 // Type of the target sparse matrix
346  , bool SO2 > // Storage order of the target sparse matrix
347  friend inline void addAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
348  {
350 
351  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
352  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
353 
354  addAssign( ~lhs, rhs.sm_ );
355  }
357  //**********************************************************************************************
358 
359  //**Subtraction assignment to dense matrices****************************************************
372  template< typename MT2 // Type of the target dense matrix
373  , bool SO2 > // Storage order of the target dense matrix
374  friend inline void subAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
375  {
377 
378  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
379  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
380 
381  subAssign( ~lhs, rhs.sm_ );
382  }
384  //**********************************************************************************************
385 
386  //**Subtraction assignment to sparse matrices***************************************************
399  template< typename MT2 // Type of the target sparse matrix
400  , bool SO2 > // Storage order of the target sparse matrix
401  friend inline void subAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
402  {
404 
405  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
406  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
407 
408  subAssign( ~lhs, rhs.sm_ );
409  }
411  //**********************************************************************************************
412 
413  //**Schur product assignment to dense matrices**************************************************
426  template< typename MT2 // Type of the target dense matrix
427  , bool SO2 > // Storage order of the target dense matrix
428  friend inline void schurAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
429  {
431 
432  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
433  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
434 
435  schurAssign( ~lhs, rhs.sm_ );
436  }
438  //**********************************************************************************************
439 
440  //**Schur product assignment to sparse matrices*************************************************
453  template< typename MT2 // Type of the target sparse matrix
454  , bool SO2 > // Storage order of the target sparse matrix
455  friend inline void schurAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
456  {
458 
459  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
460  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
461 
462  schurAssign( ~lhs, rhs.sm_ );
463  }
465  //**********************************************************************************************
466 
467  //**Multiplication assignment to dense matrices*************************************************
480  template< typename MT2 // Type of the target dense matrix
481  , bool SO2 > // Storage order of the target dense matrix
482  friend inline void multAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
483  {
485 
486  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
487  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
488 
489  multAssign( ~lhs, rhs.sm_ );
490  }
492  //**********************************************************************************************
493 
494  //**Multiplication assignment to sparse matrices************************************************
507  template< typename MT2 // Type of the target sparse matrix
508  , bool SO2 > // Storage order of the target sparse matrix
509  friend inline void multAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
510  {
512 
513  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
514  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
515 
516  multAssign( ~lhs, rhs.sm_ );
517  }
519  //**********************************************************************************************
520 
521  //**SMP assignment to dense matrices************************************************************
533  template< typename MT2 // Type of the target dense matrix
534  , bool SO2 > // Storage order of the target dense matrix
535  friend inline void smpAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
536  {
538 
539  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
540  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
541 
542  assign( ~lhs, rhs.sm_ );
543  }
545  //**********************************************************************************************
546 
547  //**SMP assignment to sparse matrices***********************************************************
559  template< typename MT2 // Type of the target sparse matrix
560  , bool SO2 > // Storage order of the target sparse matrix
561  friend inline void smpAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
562  {
564 
565  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
566  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
567 
568  assign( ~lhs, rhs.sm_ );
569  }
571  //**********************************************************************************************
572 
573  //**SMP addition assignment to dense matrices***************************************************
586  template< typename MT2 // Type of the target dense matrix
587  , bool SO2 > // Storage order of the target dense matrix
588  friend inline void smpAddAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
589  {
591 
592  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
593  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
594 
595  addAssign( ~lhs, rhs.sm_ );
596  }
598  //**********************************************************************************************
599 
600  //**SMP addition assignment to sparse matrices**************************************************
613  template< typename MT2 // Type of the target sparse matrix
614  , bool SO2 > // Storage order of the target sparse matrix
615  friend inline void smpAddAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
616  {
618 
619  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
620  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
621 
622  addAssign( ~lhs, rhs.sm_ );
623  }
625  //**********************************************************************************************
626 
627  //**SMP subtraction assignment to dense matrices************************************************
640  template< typename MT2 // Type of the target dense matrix
641  , bool SO2 > // Storage order of the target dense matrix
642  friend inline void smpSubAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
643  {
645 
646  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
647  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
648 
649  subAssign( ~lhs, rhs.sm_ );
650  }
652  //**********************************************************************************************
653 
654  //**SMP subtraction assignment to sparse matrices***********************************************
667  template< typename MT2 // Type of the target sparse matrix
668  , bool SO2 > // Storage order of the target sparse matrix
669  friend inline void smpSubAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
670  {
672 
673  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
674  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
675 
676  subAssign( ~lhs, rhs.sm_ );
677  }
679  //**********************************************************************************************
680 
681  //**SMP Schur product assignment to dense matrices**********************************************
694  template< typename MT2 // Type of the target dense matrix
695  , bool SO2 > // Storage order of the target dense matrix
696  friend inline void smpSchurAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
697  {
699 
700  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
701  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
702 
703  schurAssign( ~lhs, rhs.sm_ );
704  }
706  //**********************************************************************************************
707 
708  //**SMP Schur product assignment to sparse matrices*********************************************
721  template< typename MT2 // Type of the target sparse matrix
722  , bool SO2 > // Storage order of the target sparse matrix
723  friend inline void smpSchurAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
724  {
726 
727  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
728  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
729 
730  schurAssign( ~lhs, rhs.sm_ );
731  }
733  //**********************************************************************************************
734 
735  //**SMP multiplication assignment to dense matrices*********************************************
748  template< typename MT2 // Type of the target dense matrix
749  , bool SO2 > // Storage order of the target dense matrix
750  friend inline void smpMultAssign( DenseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
751  {
753 
754  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
755  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
756 
757  multAssign( ~lhs, rhs.sm_ );
758  }
760  //**********************************************************************************************
761 
762  //**SMP multiplication assignment to sparse matrices********************************************
775  template< typename MT2 // Type of the target sparse matrix
776  , bool SO2 > // Storage order of the target sparse matrix
777  friend inline void smpMultAssign( SparseMatrix<MT2,SO2>& lhs, const SMatSerialExpr& rhs )
778  {
780 
781  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
782  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
783 
784  multAssign( ~lhs, rhs.sm_ );
785  }
787  //**********************************************************************************************
788 
789  //**Compile time checks*************************************************************************
794  //**********************************************************************************************
795 };
796 //*************************************************************************************************
797 
798 
799 
800 
801 //=================================================================================================
802 //
803 // GLOBAL FUNCTIONS
804 //
805 //=================================================================================================
806 
807 //*************************************************************************************************
824 template< typename MT // Type of the sparse matrix
825  , bool SO > // Storage order
826 inline decltype(auto) serial( const SparseMatrix<MT,SO>& sm )
827 {
829 
830  using ReturnType = const SMatSerialExpr<MT,SO>;
831  return ReturnType( ~sm );
832 }
833 //*************************************************************************************************
834 
835 
836 
837 
838 //=================================================================================================
839 //
840 // GLOBAL RESTRUCTURING FUNCTIONS
841 //
842 //=================================================================================================
843 
844 //*************************************************************************************************
855 template< typename MT // Type of the sparse matrix
856  , bool SO > // Storage order
857 inline decltype(auto) serial( const SMatSerialExpr<MT,SO>& sm )
858 {
859  return sm;
860 }
862 //*************************************************************************************************
863 
864 
865 
866 
867 //=================================================================================================
868 //
869 // ROWS SPECIALIZATIONS
870 //
871 //=================================================================================================
872 
873 //*************************************************************************************************
875 template< typename MT, bool SO >
876 struct Rows< SMatSerialExpr<MT,SO> >
877  : public Rows<MT>
878 {};
880 //*************************************************************************************************
881 
882 
883 
884 
885 //=================================================================================================
886 //
887 // COLUMNS SPECIALIZATIONS
888 //
889 //=================================================================================================
890 
891 //*************************************************************************************************
893 template< typename MT, bool SO >
894 struct Columns< SMatSerialExpr<MT,SO> >
895  : public Columns<MT>
896 {};
898 //*************************************************************************************************
899 
900 
901 
902 
903 //=================================================================================================
904 //
905 // ISSYMMETRIC SPECIALIZATIONS
906 //
907 //=================================================================================================
908 
909 //*************************************************************************************************
911 template< typename MT, bool SO >
912 struct IsSymmetric< SMatSerialExpr<MT,SO> >
913  : public BoolConstant< IsSymmetric<MT>::value >
914 {};
916 //*************************************************************************************************
917 
918 
919 
920 
921 //=================================================================================================
922 //
923 // ISHERMITIAN SPECIALIZATIONS
924 //
925 //=================================================================================================
926 
927 //*************************************************************************************************
929 template< typename MT, bool SO >
930 struct IsHermitian< SMatSerialExpr<MT,SO> >
931  : public BoolConstant< IsHermitian<MT>::value >
932 {};
934 //*************************************************************************************************
935 
936 
937 
938 
939 //=================================================================================================
940 //
941 // ISLOWER SPECIALIZATIONS
942 //
943 //=================================================================================================
944 
945 //*************************************************************************************************
947 template< typename MT, bool SO >
948 struct IsLower< SMatSerialExpr<MT,SO> >
949  : public BoolConstant< IsLower<MT>::value >
950 {};
952 //*************************************************************************************************
953 
954 
955 
956 
957 //=================================================================================================
958 //
959 // ISUNILOWER SPECIALIZATIONS
960 //
961 //=================================================================================================
962 
963 //*************************************************************************************************
965 template< typename MT, bool SO >
966 struct IsUniLower< SMatSerialExpr<MT,SO> >
967  : public BoolConstant< IsUniLower<MT>::value >
968 {};
970 //*************************************************************************************************
971 
972 
973 
974 
975 //=================================================================================================
976 //
977 // ISSTRICTLYLOWER SPECIALIZATIONS
978 //
979 //=================================================================================================
980 
981 //*************************************************************************************************
983 template< typename MT, bool SO >
984 struct IsStrictlyLower< SMatSerialExpr<MT,SO> >
985  : public BoolConstant< IsStrictlyLower<MT>::value >
986 {};
988 //*************************************************************************************************
989 
990 
991 
992 
993 //=================================================================================================
994 //
995 // ISUPPER SPECIALIZATIONS
996 //
997 //=================================================================================================
998 
999 //*************************************************************************************************
1001 template< typename MT, bool SO >
1002 struct IsUpper< SMatSerialExpr<MT,SO> >
1003  : public BoolConstant< IsUpper<MT>::value >
1004 {};
1006 //*************************************************************************************************
1007 
1008 
1009 
1010 
1011 //=================================================================================================
1012 //
1013 // ISUNIUPPER SPECIALIZATIONS
1014 //
1015 //=================================================================================================
1016 
1017 //*************************************************************************************************
1019 template< typename MT, bool SO >
1020 struct IsUniUpper< SMatSerialExpr<MT,SO> >
1021  : public BoolConstant< IsUniUpper<MT>::value >
1022 {};
1024 //*************************************************************************************************
1025 
1026 
1027 
1028 
1029 //=================================================================================================
1030 //
1031 // ISSTRICTLYUPPER SPECIALIZATIONS
1032 //
1033 //=================================================================================================
1034 
1035 //*************************************************************************************************
1037 template< typename MT, bool SO >
1038 struct IsStrictlyUpper< SMatSerialExpr<MT,SO> >
1039  : public BoolConstant< IsStrictlyUpper<MT>::value >
1040 {};
1042 //*************************************************************************************************
1043 
1044 } // namespace blaze
1045 
1046 #endif
Header file for auxiliary alias declarations.
SMatSerialExpr(const MT &sm) noexcept
Constructor for the SMatSerialExpr class.
Definition: SMatSerialExpr.h:117
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatSerialExpr.h:235
ResultType_< MT > ResultType
Result type for expression template evaluations.
Definition: SMatSerialExpr.h:94
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.
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
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SMatSerialExpr.h:180
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 Computation base class.
Compile time check for upper triangular matrices.This type trait tests whether or not the given templ...
Definition: IsUpper.h:88
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: SMatSerialExpr.h:144
If_< IsExpression< MT >, const MT, const MT &> Operand
Composite data type of the sparse matrix expression.
Definition: SMatSerialExpr.h:104
Header file for the IsUniLower type trait.
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
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.
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
Operand operand() const noexcept
Returns the sparse matrix operand.
Definition: SMatSerialExpr.h:201
Header file for the If class template.
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
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: SMatSerialExpr.h:160
Header file for the Columns type trait.
ElementType_< MT > ElementType
Resulting element type.
Definition: SMatSerialExpr.h:97
typename T::ElementType ElementType_
Alias declaration for nested ElementType type definitions.The ElementType_ alias declaration provides...
Definition: Aliases.h:163
Header file for the IsLower type trait.
TransposeType_< MT > TransposeType
Transpose type for expression template evaluations.
Definition: SMatSerialExpr.h:96
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
Header file for all forward declarations for expression class templates.
Header file for the MatSerialExpr base class.
Header file for the IsStrictlyLower type trait.
Operand sm_
Sparse matrix of the serial evaluation expression.
Definition: SMatSerialExpr.h:252
Compile time check for lower unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniLower.h:86
Header file for run time assertion macros.
Utility type for generic codes.
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: SMatSerialExpr.h:245
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SMatSerialExpr.h:223
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
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row.
Definition: SMatSerialExpr.h:191
decltype(auto) serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:819
typename T::OppositeType OppositeType_
Alias declaration for nested OppositeType type definitions.The OppositeType_ alias declaration provid...
Definition: Aliases.h:263
ReturnType_< MT > ReturnType
Return type for expression template evaluations.
Definition: SMatSerialExpr.h:98
OppositeType_< MT > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatSerialExpr.h:95
const ResultType CompositeType
Data type for composite expression templates.
Definition: SMatSerialExpr.h:101
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
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
Compile time evaluation of the number of rows of a matrix.The Rows type trait evaluates the number of...
Definition: Rows.h:75
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SMatSerialExpr.h:129
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 columns() const noexcept
Returns the current number of columns of the matrix.
Definition: SMatSerialExpr.h:170
Expression object for the forced serial evaluation of sparse matrices.The SMatSerialExpr class repres...
Definition: Forward.h:115
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.