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>
71 #include <blaze/util/Assert.h>
74 #include <blaze/util/InvalidType.h>
75 #include <blaze/util/mpl/And.h>
76 #include <blaze/util/mpl/If.h>
77 #include <blaze/util/Types.h>
78 
79 
80 namespace blaze {
81 
82 //=================================================================================================
83 //
84 // CLASS SMATEVALEXPR
85 //
86 //=================================================================================================
87 
88 //*************************************************************************************************
95 template< typename MT // Type of the sparse matrix
96  , bool SO > // Storage order
97 class SMatEvalExpr : public SparseMatrix< SMatEvalExpr<MT,SO>, SO >
98  , private MatEvalExpr
99  , private Computation
100 {
101  public:
102  //**Type definitions****************************************************************************
109 
111  typedef const ResultType CompositeType;
112 
114  typedef If_< IsExpression<MT>, const MT, const MT& > Operand;
115  //**********************************************************************************************
116 
117  //**Compilation flags***************************************************************************
119  enum : bool { smpAssignable = MT::smpAssignable };
120  //**********************************************************************************************
121 
122  //**Constructor*********************************************************************************
127  explicit inline SMatEvalExpr( const MT& sm ) noexcept
128  : sm_( sm ) // Sparse matrix of the evaluation expression
129  {}
130  //**********************************************************************************************
131 
132  //**Access operator*****************************************************************************
139  inline ReturnType operator()( size_t i, size_t j ) const {
140  BLAZE_INTERNAL_ASSERT( i < sm_.rows() , "Invalid row access index" );
141  BLAZE_INTERNAL_ASSERT( j < sm_.columns(), "Invalid column access index" );
142  return sm_(i,j);
143  }
144  //**********************************************************************************************
145 
146  //**At function*********************************************************************************
154  inline ReturnType at( size_t i, size_t j ) const {
155  if( i >= sm_.rows() ) {
156  BLAZE_THROW_OUT_OF_RANGE( "Invalid row access index" );
157  }
158  if( j >= sm_.columns() ) {
159  BLAZE_THROW_OUT_OF_RANGE( "Invalid column access index" );
160  }
161  return (*this)(i,j);
162  }
163  //**********************************************************************************************
164 
165  //**Rows function*******************************************************************************
170  inline size_t rows() const noexcept {
171  return sm_.rows();
172  }
173  //**********************************************************************************************
174 
175  //**Columns function****************************************************************************
180  inline size_t columns() const noexcept {
181  return sm_.columns();
182  }
183  //**********************************************************************************************
184 
185  //**NonZeros function***************************************************************************
190  inline size_t nonZeros() const {
191  return sm_.nonZeros();
192  }
193  //**********************************************************************************************
194 
195  //**NonZeros function***************************************************************************
201  inline size_t nonZeros( size_t i ) const {
202  return sm_.nonZeros(i);
203  }
204  //**********************************************************************************************
205 
206  //**Operand access******************************************************************************
211  inline Operand 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****************************************************************************
252  Operand sm_;
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 SMatEvalExpr& 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 SMatEvalExpr& 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 SMatEvalExpr& 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 SMatEvalExpr& 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****************************************************
371  template< typename MT2 // Type of the target dense matrix
372  , bool SO2 > // Storage order of the target dense matrix
373  friend inline void subAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
374  {
376 
377  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
378  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
379 
380  subAssign( ~lhs, rhs.sm_ );
381  }
383  //**********************************************************************************************
384 
385  //**Subtraction assignment to sparse matrices***************************************************
397  template< typename MT2 // Type of the target sparse matrix
398  , bool SO2 > // Storage order of the target sparse matrix
399  friend inline void subAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
400  {
402 
403  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
404  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
405 
406  subAssign( ~lhs, rhs.sm_ );
407  }
409  //**********************************************************************************************
410 
411  //**Multiplication assignment to dense matrices*************************************************
423  template< typename MT2 // Type of the target dense matrix
424  , bool SO2 > // Storage order of the target dense matrix
425  friend inline void multAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
426  {
428 
429  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
430  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
431 
432  multAssign( ~lhs, rhs.sm_ );
433  }
435  //**********************************************************************************************
436 
437  //**Multiplication assignment to sparse matrices************************************************
449  template< typename MT2 // Type of the target sparse matrix
450  , bool SO2 > // Storage order of the target sparse matrix
451  friend inline void multAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
452  {
454 
455  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
456  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
457 
458  multAssign( ~lhs, rhs.sm_ );
459  }
461  //**********************************************************************************************
462 
463  //**SMP assignment to dense matrices************************************************************
475  template< typename MT2 // Type of the target dense matrix
476  , bool SO2 > // Storage order of the target dense matrix
477  friend inline void smpAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
478  {
480 
481  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
482  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
483 
484  smpAssign( ~lhs, rhs.sm_ );
485  }
487  //**********************************************************************************************
488 
489  //**SMP assignment to sparse matrices***********************************************************
501  template< typename MT2 // Type of the target sparse matrix
502  , bool SO2 > // Storage order of the target sparse matrix
503  friend inline void smpAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
504  {
506 
507  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
508  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
509 
510  smpAssign( ~lhs, rhs.sm_ );
511  }
513  //**********************************************************************************************
514 
515  //**SMP addition assignment to dense matrices***************************************************
527  template< typename MT2 // Type of the target dense matrix
528  , bool SO2 > // Storage order of the target dense matrix
529  friend inline void smpAddAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
530  {
532 
533  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
534  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
535 
536  smpAddAssign( ~lhs, rhs.sm_ );
537  }
539  //**********************************************************************************************
540 
541  //**SMP addition assignment to sparse matrices**************************************************
553  template< typename MT2 // Type of the target sparse matrix
554  , bool SO2 > // Storage order of the target sparse matrix
555  friend inline void smpAddAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
556  {
558 
559  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
560  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
561 
562  smpAddAssign( ~lhs, rhs.sm_ );
563  }
565  //**********************************************************************************************
566 
567  //**SMP subtraction assignment to dense matrices************************************************
579  template< typename MT2 // Type of the target dense matrix
580  , bool SO2 > // Storage order of the target dense matrix
581  friend inline void smpSubAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
582  {
584 
585  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
586  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
587 
588  smpSubAssign( ~lhs, rhs.sm_ );
589  }
591  //**********************************************************************************************
592 
593  //**SMP subtraction assignment to sparse matrices***********************************************
605  template< typename MT2 // Type of the target sparse matrix
606  , bool SO2 > // Storage order of the target sparse matrix
607  friend inline void smpSubAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
608  {
610 
611  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
612  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
613 
614  smpSubAssign( ~lhs, rhs.sm_ );
615  }
617  //**********************************************************************************************
618 
619  //**SMP multiplication assignment to dense matrices*********************************************
632  template< typename MT2 // Type of the target dense matrix
633  , bool SO2 > // Storage order of the target dense matrix
634  friend inline void smpMultAssign( DenseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
635  {
637 
638  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
639  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
640 
641  smpMultAssign( ~lhs, rhs.sm_ );
642  }
644  //**********************************************************************************************
645 
646  //**SMP multiplication assignment to sparse matrices********************************************
659  template< typename MT2 // Type of the target sparse matrix
660  , bool SO2 > // Storage order of the target sparse matrix
661  friend inline void smpMultAssign( SparseMatrix<MT2,SO2>& lhs, const SMatEvalExpr& rhs )
662  {
664 
665  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
666  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
667 
668  smpMultAssign( ~lhs, rhs.sm_ );
669  }
671  //**********************************************************************************************
672 
673  //**Compile time checks*************************************************************************
678  //**********************************************************************************************
679 };
680 //*************************************************************************************************
681 
682 
683 
684 
685 //=================================================================================================
686 //
687 // GLOBAL FUNCTIONS
688 //
689 //=================================================================================================
690 
691 //*************************************************************************************************
708 template< typename MT // Type of the sparse matrix
709  , bool SO > // Storage order
711 {
713 
714  return SMatEvalExpr<MT,SO>( ~sm );
715 }
716 //*************************************************************************************************
717 
718 
719 
720 
721 //=================================================================================================
722 //
723 // GLOBAL RESTRUCTURING FUNCTIONS
724 //
725 //=================================================================================================
726 
727 //*************************************************************************************************
738 template< typename MT // Type of the sparse matrix
739  , bool SO > // Storage order
740 inline const SMatEvalExpr<MT,SO> eval( const SMatEvalExpr<MT,SO>& sm )
741 {
742  return sm;
743 }
745 //*************************************************************************************************
746 
747 
748 
749 
750 //=================================================================================================
751 //
752 // ROWS SPECIALIZATIONS
753 //
754 //=================================================================================================
755 
756 //*************************************************************************************************
758 template< typename MT, bool SO >
759 struct Rows< SMatEvalExpr<MT,SO> > : public Rows<MT>
760 {};
762 //*************************************************************************************************
763 
764 
765 
766 
767 //=================================================================================================
768 //
769 // COLUMNS SPECIALIZATIONS
770 //
771 //=================================================================================================
772 
773 //*************************************************************************************************
775 template< typename MT, bool SO >
776 struct Columns< SMatEvalExpr<MT,SO> > : public Columns<MT>
777 {};
779 //*************************************************************************************************
780 
781 
782 
783 
784 //=================================================================================================
785 //
786 // ISSYMMETRIC SPECIALIZATIONS
787 //
788 //=================================================================================================
789 
790 //*************************************************************************************************
792 template< typename MT, bool SO >
793 struct IsSymmetric< SMatEvalExpr<MT,SO> >
794  : public BoolConstant< IsSymmetric<MT>::value >
795 {};
797 //*************************************************************************************************
798 
799 
800 
801 
802 //=================================================================================================
803 //
804 // ISHERMITIAN SPECIALIZATIONS
805 //
806 //=================================================================================================
807 
808 //*************************************************************************************************
810 template< typename MT, bool SO >
811 struct IsHermitian< SMatEvalExpr<MT,SO> >
812  : public BoolConstant< IsHermitian<MT>::value >
813 {};
815 //*************************************************************************************************
816 
817 
818 
819 
820 //=================================================================================================
821 //
822 // ISLOWER SPECIALIZATIONS
823 //
824 //=================================================================================================
825 
826 //*************************************************************************************************
828 template< typename MT, bool SO >
829 struct IsLower< SMatEvalExpr<MT,SO> >
830  : public BoolConstant< IsLower<MT>::value >
831 {};
833 //*************************************************************************************************
834 
835 
836 
837 
838 //=================================================================================================
839 //
840 // ISUNILOWER SPECIALIZATIONS
841 //
842 //=================================================================================================
843 
844 //*************************************************************************************************
846 template< typename MT, bool SO >
847 struct IsUniLower< SMatEvalExpr<MT,SO> >
848  : public BoolConstant< IsUniLower<MT>::value >
849 {};
851 //*************************************************************************************************
852 
853 
854 
855 
856 //=================================================================================================
857 //
858 // ISSTRICTLYLOWER SPECIALIZATIONS
859 //
860 //=================================================================================================
861 
862 //*************************************************************************************************
864 template< typename MT, bool SO >
865 struct IsStrictlyLower< SMatEvalExpr<MT,SO> >
866  : public BoolConstant< IsStrictlyLower<MT>::value >
867 {};
869 //*************************************************************************************************
870 
871 
872 
873 
874 //=================================================================================================
875 //
876 // ISUPPER SPECIALIZATIONS
877 //
878 //=================================================================================================
879 
880 //*************************************************************************************************
882 template< typename MT, bool SO >
883 struct IsUpper< SMatEvalExpr<MT,SO> >
884  : public BoolConstant< IsUpper<MT>::value >
885 {};
887 //*************************************************************************************************
888 
889 
890 
891 
892 //=================================================================================================
893 //
894 // ISUNIUPPER SPECIALIZATIONS
895 //
896 //=================================================================================================
897 
898 //*************************************************************************************************
900 template< typename MT, bool SO >
901 struct IsUniUpper< SMatEvalExpr<MT,SO> >
902  : public BoolConstant< IsUniUpper<MT>::value >
903 {};
905 //*************************************************************************************************
906 
907 
908 
909 
910 //=================================================================================================
911 //
912 // ISSTRICTLYUPPER SPECIALIZATIONS
913 //
914 //=================================================================================================
915 
916 //*************************************************************************************************
918 template< typename MT, bool SO >
919 struct IsStrictlyUpper< SMatEvalExpr<MT,SO> >
920  : public BoolConstant< IsStrictlyUpper<MT>::value >
921 {};
923 //*************************************************************************************************
924 
925 
926 
927 
928 //=================================================================================================
929 //
930 // EXPRESSION TRAIT SPECIALIZATIONS
931 //
932 //=================================================================================================
933 
934 //*************************************************************************************************
936 template< typename MT >
937 struct SMatEvalExprTrait< SMatEvalExpr<MT,false> >
938 {
939  public:
940  //**********************************************************************************************
943  , INVALID_TYPE >;
944  //**********************************************************************************************
945 };
947 //*************************************************************************************************
948 
949 
950 //*************************************************************************************************
952 template< typename MT >
953 struct TSMatEvalExprTrait< SMatEvalExpr<MT,true> >
954 {
955  public:
956  //**********************************************************************************************
959  , INVALID_TYPE >;
960  //**********************************************************************************************
961 };
963 //*************************************************************************************************
964 
965 
966 //*************************************************************************************************
968 template< typename MT, bool SO, bool AF >
969 struct SubmatrixExprTrait< SMatEvalExpr<MT,SO>, AF >
970 {
971  public:
972  //**********************************************************************************************
974  //**********************************************************************************************
975 };
977 //*************************************************************************************************
978 
979 
980 //*************************************************************************************************
982 template< typename MT, bool SO >
983 struct RowExprTrait< SMatEvalExpr<MT,SO> >
984 {
985  public:
986  //**********************************************************************************************
988  //**********************************************************************************************
989 };
991 //*************************************************************************************************
992 
993 
994 //*************************************************************************************************
996 template< typename MT, bool SO >
997 struct ColumnExprTrait< SMatEvalExpr<MT,SO> >
998 {
999  public:
1000  //**********************************************************************************************
1002  //**********************************************************************************************
1003 };
1005 //*************************************************************************************************
1006 
1007 } // namespace blaze
1008 
1009 #endif
Header file for auxiliary alias declarations.
Header file for the Rows type trait.
Header file for the IsUniUpper type trait.
Header file for basic type definitions.
Evaluation of the expression type of a sparse matrix evaluation operation.Via this type trait it is p...
Definition: SMatEvalExprTrait.h:74
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: SMatEvalExpr.h:180
Evaluation of the expression type of a sparse matrix evaluation operation.Via this type trait it is p...
Definition: TSMatEvalExprTrait.h:74
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:160
Header file for the IsSparseMatrix type trait.
#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
Header file for the ColumnExprTrait class template.
Header file for the TSMatEvalExprTrait class template.
Header file for the IsColumnMajorMatrix type trait.
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 And class template.
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.
typename EvalExprTrait< T >::Type EvalExprTrait_
Auxiliary alias declaration for the EvalExprTrait class template.The EvalExprTrait_ alias declaration...
Definition: EvalExprTrait.h:142
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.
typename T::ResultType ResultType_
Alias declaration for nested ResultType type definitions.The ResultType_ alias declaration provides a...
Definition: Aliases.h:323
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:129
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:71
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:119
typename T::ReturnType ReturnType_
Alias declaration for nested ReturnType type definitions.The ReturnType_ alias declaration provides a...
Definition: Aliases.h:343
SMatEvalExpr< MT, SO > This
Type of this SMatEvalExpr instance.
Definition: SMatEvalExpr.h:103
Header file for the SparseMatrix base class.
SMatEvalExpr(const MT &sm) noexcept
Constructor for the SMatEvalExpr class.
Definition: SMatEvalExpr.h:127
Constraint on the data type.
Header file for the SMatEvalExprTrait class template.
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.
Compile time check for row-major matrix types.This type trait tests whether or not the given template...
Definition: IsRowMajorMatrix.h:83
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SMatEvalExpr.h:223
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatEvalExpr.h:235
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: SMatEvalExpr.h:154
EnableIf_< IsDenseMatrix< MT1 > > smpAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:98
#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
Evaluation of the expression type type of a submatrix operation.Via this type trait it is possible to...
Definition: SubmatrixExprTrait.h:80
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: SMatEvalExpr.h:245
Header file for the IsLower type trait.
Operand operand() const noexcept
Returns the sparse matrix operand.
Definition: SMatEvalExpr.h:211
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:201
Evaluation of the expression type type of a row operation.Via this type trait it is possible to evalu...
Definition: RowExprTrait.h:79
Header file for the RowExprTrait class template.
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:252
Header file for the SubmatrixExprTrait class template.
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SMatEvalExpr.h:190
Header file for run time assertion macros.
Compile time check for column-major matrix types.This type trait tests whether or not the given templ...
Definition: IsColumnMajorMatrix.h:83
TransposeType_< MT > TransposeType
Transpose type for expression template evaluations.
Definition: SMatEvalExpr.h:106
Utility type for generic codes.
ElementType_< MT > ElementType
Resulting element type.
Definition: SMatEvalExpr.h:107
ReturnType_< MT > ReturnType
Return type for expression template evaluations.
Definition: SMatEvalExpr.h:108
Header file for the EvalExprTrait class template.
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:160
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:93
Compile time check for Hermitian matrices.This type trait tests whether or not the given template par...
Definition: IsHermitian.h:85
ResultType_< MT > ResultType
Result type for expression template evaluations.
Definition: SMatEvalExpr.h:104
If_< IsExpression< MT >, const MT, const MT &> Operand
Composite data type of the sparse matrix expression.
Definition: SMatEvalExpr.h:114
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SMatEvalExpr.h:139
const DMatEvalExpr< MT, SO > eval(const DenseMatrix< MT, SO > &dm)
Forces the evaluation of the given dense matrix expression dm.
Definition: DMatEvalExpr.h:703
typename T::OppositeType OppositeType_
Alias declaration for nested OppositeType type definitions.The OppositeType_ alias declaration provid...
Definition: Aliases.h:243
Compile time check for strictly lower triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyLower.h:86
const ResultType CompositeType
Data type for composite expression templates.
Definition: SMatEvalExpr.h:111
Header file for the IsRowMajorMatrix type trait.
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:76
Expression object for the forced evaluation of sparse matrices.The SMatEvalExpr class represents the ...
Definition: Forward.h:103
Compile time evaluation of the number of rows of a matrix.The Rows type trait evaluates the number of...
Definition: Rows.h:76
OppositeType_< MT > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatEvalExpr.h:105
typename T::TransposeType TransposeType_
Alias declaration for nested TransposeType type definitions.The TransposeType_ alias declaration prov...
Definition: Aliases.h:403
Header file for the IsUpper type trait.
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: SMatEvalExpr.h:170
Header file for the IsHermitian type trait.
Evaluation of the expression type type of a column operation.Via this type trait it is possible to ev...
Definition: ColumnExprTrait.h:78
#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.