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 <cmath>
70 #include <blaze/util/Assert.h>
72 #include <blaze/util/Exception.h>
73 #include <blaze/util/InvalidType.h>
75 #include <blaze/util/SelectType.h>
76 #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****************************************************************************
104  typedef typename MT::ResultType ResultType;
105  typedef typename MT::OppositeType OppositeType;
106  typedef typename MT::TransposeType TransposeType;
107  typedef typename MT::ElementType ElementType;
108  typedef typename MT::ReturnType ReturnType;
109 
111  typedef const ResultType CompositeType;
112 
114  typedef typename SelectType< IsExpression<MT>::value, const MT, const MT& >::Type Operand;
115  //**********************************************************************************************
116 
117  //**Compilation flags***************************************************************************
119  enum { smpAssignable = MT::smpAssignable };
120  //**********************************************************************************************
121 
122  //**Constructor*********************************************************************************
127  explicit inline SMatEvalExpr( const MT& sm )
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 {
171  return sm_.rows();
172  }
173  //**********************************************************************************************
174 
175  //**Columns function****************************************************************************
180  inline size_t columns() const {
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 {
212  return sm_;
213  }
214  //**********************************************************************************************
215 
216  //**********************************************************************************************
222  template< typename T >
223  inline bool canAlias( const T* alias ) const {
224  return sm_.canAlias( alias );
225  }
226  //**********************************************************************************************
227 
228  //**********************************************************************************************
234  template< typename T >
235  inline bool isAliased( const T* alias ) const {
236  return sm_.isAliased( alias );
237  }
238  //**********************************************************************************************
239 
240  //**********************************************************************************************
245  inline bool canSMPAssign() const {
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> > : public IsTrue< IsSymmetric<MT>::value >
794 {};
796 //*************************************************************************************************
797 
798 
799 
800 
801 //=================================================================================================
802 //
803 // ISHERMITIAN SPECIALIZATIONS
804 //
805 //=================================================================================================
806 
807 //*************************************************************************************************
809 template< typename MT, bool SO >
810 struct IsHermitian< SMatEvalExpr<MT,SO> > : public IsTrue< IsHermitian<MT>::value >
811 {};
813 //*************************************************************************************************
814 
815 
816 
817 
818 //=================================================================================================
819 //
820 // ISLOWER SPECIALIZATIONS
821 //
822 //=================================================================================================
823 
824 //*************************************************************************************************
826 template< typename MT, bool SO >
827 struct IsLower< SMatEvalExpr<MT,SO> > : public IsTrue< IsLower<MT>::value >
828 {};
830 //*************************************************************************************************
831 
832 
833 
834 
835 //=================================================================================================
836 //
837 // ISUNILOWER SPECIALIZATIONS
838 //
839 //=================================================================================================
840 
841 //*************************************************************************************************
843 template< typename MT, bool SO >
844 struct IsUniLower< SMatEvalExpr<MT,SO> > : public IsTrue< IsUniLower<MT>::value >
845 {};
847 //*************************************************************************************************
848 
849 
850 
851 
852 //=================================================================================================
853 //
854 // ISSTRICTLYLOWER SPECIALIZATIONS
855 //
856 //=================================================================================================
857 
858 //*************************************************************************************************
860 template< typename MT, bool SO >
861 struct IsStrictlyLower< SMatEvalExpr<MT,SO> > : public IsTrue< IsStrictlyLower<MT>::value >
862 {};
864 //*************************************************************************************************
865 
866 
867 
868 
869 //=================================================================================================
870 //
871 // ISUPPER SPECIALIZATIONS
872 //
873 //=================================================================================================
874 
875 //*************************************************************************************************
877 template< typename MT, bool SO >
878 struct IsUpper< SMatEvalExpr<MT,SO> > : public IsTrue< IsUpper<MT>::value >
879 {};
881 //*************************************************************************************************
882 
883 
884 
885 
886 //=================================================================================================
887 //
888 // ISUNIUPPER SPECIALIZATIONS
889 //
890 //=================================================================================================
891 
892 //*************************************************************************************************
894 template< typename MT, bool SO >
895 struct IsUniUpper< SMatEvalExpr<MT,SO> > : public IsTrue< IsUniUpper<MT>::value >
896 {};
898 //*************************************************************************************************
899 
900 
901 
902 
903 //=================================================================================================
904 //
905 // ISSTRICTLYUPPER SPECIALIZATIONS
906 //
907 //=================================================================================================
908 
909 //*************************************************************************************************
911 template< typename MT, bool SO >
912 struct IsStrictlyUpper< SMatEvalExpr<MT,SO> > : public IsTrue< IsStrictlyUpper<MT>::value >
913 {};
915 //*************************************************************************************************
916 
917 
918 
919 
920 //=================================================================================================
921 //
922 // EXPRESSION TRAIT SPECIALIZATIONS
923 //
924 //=================================================================================================
925 
926 //*************************************************************************************************
928 template< typename MT >
929 struct SMatEvalExprTrait< SMatEvalExpr<MT,false> >
930 {
931  public:
932  //**********************************************************************************************
933  typedef typename SelectType< IsSparseMatrix<MT>::value && IsRowMajorMatrix<MT>::value
934  , SMatEvalExpr<MT,false>
935  , INVALID_TYPE >::Type Type;
936  //**********************************************************************************************
937 };
939 //*************************************************************************************************
940 
941 
942 //*************************************************************************************************
944 template< typename MT >
945 struct TSMatEvalExprTrait< SMatEvalExpr<MT,true> >
946 {
947  public:
948  //**********************************************************************************************
949  typedef typename SelectType< IsSparseMatrix<MT>::value && IsColumnMajorMatrix<MT>::value
950  , SMatEvalExpr<MT,true>
951  , INVALID_TYPE >::Type Type;
952  //**********************************************************************************************
953 };
955 //*************************************************************************************************
956 
957 
958 //*************************************************************************************************
960 template< typename MT, bool SO, bool AF >
961 struct SubmatrixExprTrait< SMatEvalExpr<MT,SO>, AF >
962 {
963  public:
964  //**********************************************************************************************
965  typedef typename EvalExprTrait< typename SubmatrixExprTrait<const MT,AF>::Type >::Type Type;
966  //**********************************************************************************************
967 };
969 //*************************************************************************************************
970 
971 
972 //*************************************************************************************************
974 template< typename MT, bool SO >
975 struct RowExprTrait< SMatEvalExpr<MT,SO> >
976 {
977  public:
978  //**********************************************************************************************
979  typedef typename EvalExprTrait< typename RowExprTrait<const MT>::Type >::Type Type;
980  //**********************************************************************************************
981 };
983 //*************************************************************************************************
984 
985 
986 //*************************************************************************************************
988 template< typename MT, bool SO >
989 struct ColumnExprTrait< SMatEvalExpr<MT,SO> >
990 {
991  public:
992  //**********************************************************************************************
993  typedef typename EvalExprTrait< typename ColumnExprTrait<const MT>::Type >::Type Type;
994  //**********************************************************************************************
995 };
997 //*************************************************************************************************
998 
999 } // namespace blaze
1000 
1001 #endif
MT::ReturnType ReturnType
Return type for expression template evaluations.
Definition: SMatEvalExpr.h:108
Header file for the Rows type trait.
Header file for the IsUniUpper type trait.
Header file for basic type definitions.
size_t columns() const
Returns the current number of columns of the matrix.
Definition: SMatEvalExpr.h:180
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:81
Header file for the ColumnExprTrait class template.
Header file for the TSMatEvalExprTrait class template.
Header file for the IsColumnMajorMatrix type trait.
Header file for the MatEvalExpr base class.
Header file for the Computation base class.
Operand operand() const
Returns the sparse matrix operand.
Definition: SMatEvalExpr.h:211
Header file for the IsUniLower type trait.
CompressedMatrix< Type, false > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: CompressedMatrix.h:2584
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:117
Constraint on the data type.
SMatEvalExpr< MT, SO > This
Type of this SMatEvalExpr instance.
Definition: SMatEvalExpr.h:103
Header file for the SparseMatrix base class.
Constraint on the data type.
Header file for the SMatEvalExprTrait class template.
Compile time type selection.The SelectType class template selects one of the two given types T1 and T...
Definition: SelectType.h:59
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
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exceptionThis macro encapsulates the default way of Bla...
Definition: Exception.h:331
Header file for the Columns type trait.
Header file for the IsLower type trait.
SelectType< IsExpression< MT >::value, const MT, const MT & >::Type 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
CompressedMatrix< Type, false > TransposeType
Transpose type for expression template evaluations.
Definition: CompressedMatrix.h:2585
Constraints on the storage order of matrix types.
size_t rows() const
Returns the current number of rows of the matrix.
Definition: SMatEvalExpr.h:170
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:2586
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row.
Definition: SMatEvalExpr.h:201
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: SMatEvalExpr.h:154
Header file for the IsStrictlyLower type trait.
Operand sm_
Sparse matrix of the evaluation expression.
Definition: SMatEvalExpr.h:252
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SMatEvalExpr.h:190
EnableIf< IsDenseMatrix< MT1 > >::Type 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 SubmatrixExprTrait class template.
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2587
Header file for run time assertion macros.
EnableIf< IsDenseMatrix< MT1 > >::Type 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
MT::ElementType ElementType
Resulting element type.
Definition: SMatEvalExpr.h:107
Utility type for generic codes.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: SMatEvalExpr.h:223
Header file for the EvalExprTrait class template.
const DMatEvalExpr< MT, SO > eval(const DenseMatrix< MT, SO > &dm)
Forces the evaluation of the given dense matrix expression dm.
Definition: DMatEvalExpr.h:703
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: SMatEvalExpr.h:245
const ResultType CompositeType
Data type for composite expression templates.
Definition: SMatEvalExpr.h:111
Header file for the IsRowMajorMatrix type trait.
EnableIf< IsDenseMatrix< MT1 > >::Type 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
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:157
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2583
Header file for the IsTrue value trait.
Expression object for the forced evaluation of sparse matrices.The SMatEvalExpr class represents the ...
Definition: Forward.h:100
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: SMatEvalExpr.h:235
SMatEvalExpr(const MT &sm)
Constructor for the SMatEvalExpr class.
Definition: SMatEvalExpr.h:127
Header file for the IsUpper type trait.
Header file for exception macros.
Header file for the IsHermitian type trait.
MT::ResultType ResultType
Result type for expression template evaluations.
Definition: SMatEvalExpr.h:104
EnableIf< IsDenseVector< VT1 > >::Type 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:189
#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
MT::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: SMatEvalExpr.h:106
#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:79
MT::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatEvalExpr.h:105
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.