SMatTransExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_SMATTRANSEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SMATTRANSEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <iterator>
44 #include <blaze/math/Aliases.h>
48 #include <blaze/math/Exception.h>
68 #include <blaze/util/Assert.h>
69 #include <blaze/util/EnableIf.h>
71 #include <blaze/util/InvalidType.h>
72 #include <blaze/util/mpl/If.h>
73 #include <blaze/util/Types.h>
76 
77 
78 namespace blaze {
79 
80 //=================================================================================================
81 //
82 // CLASS SMATTRANSEXPR
83 //
84 //=================================================================================================
85 
86 //*************************************************************************************************
93 template< typename MT // Type of the sparse matrix
94  , bool SO > // Storage order
95 class SMatTransExpr
96  : public MatTransExpr< SparseMatrix< SMatTransExpr<MT,SO>, SO > >
97  , private If< IsComputation<MT>, Computation, Transformation >::Type
98 {
99  private:
100  //**Type definitions****************************************************************************
101  using RT = ResultType_<MT>;
103  //**********************************************************************************************
104 
105  //**Serial evaluation strategy******************************************************************
107 
113  enum : bool { useAssign = RequiresEvaluation<MT>::value };
114 
116  template< typename MT2 >
118  struct UseAssign {
119  enum : bool { value = useAssign };
120  };
122  //**********************************************************************************************
123 
124  //**Parallel evaluation strategy****************************************************************
126 
131  template< typename MT2 >
132  struct UseSMPAssign {
133  enum : bool { value = MT2::smpAssignable && useAssign };
134  };
136  //**********************************************************************************************
137 
138  //**********************************************************************************************
140 
144  template< typename MT2 >
145  struct GetConstIterator {
147  struct Success { using Type = typename MT2::ConstIterator; };
148  struct Failure { using Type = INVALID_TYPE; };
149  using Type = typename If_< HasConstIterator<MT2>, Success, Failure >::Type;
150  };
152  //**********************************************************************************************
153 
154  public:
155  //**Type definitions****************************************************************************
162 
165 
167  using ConstIterator = typename GetConstIterator<MT>::Type;
168 
170  using Operand = If_< IsExpression<MT>, const MT, const MT& >;
171  //**********************************************************************************************
172 
173  //**Compilation flags***************************************************************************
175  enum : bool { smpAssignable = MT::smpAssignable };
176  //**********************************************************************************************
177 
178  //**Constructor*********************************************************************************
183  explicit inline SMatTransExpr( const MT& sm ) noexcept
184  : sm_( sm )
185  {}
186  //**********************************************************************************************
187 
188  //**Access operator*****************************************************************************
195  inline ReturnType operator()( size_t i, size_t j ) const {
196  BLAZE_INTERNAL_ASSERT( i < sm_.columns(), "Invalid row access index" );
197  BLAZE_INTERNAL_ASSERT( j < sm_.rows() , "Invalid column access index" );
198  return sm_(j,i);
199  }
200  //**********************************************************************************************
201 
202  //**At function*********************************************************************************
210  inline ReturnType at( size_t i, size_t j ) const {
211  if( i >= sm_.columns() ) {
212  BLAZE_THROW_OUT_OF_RANGE( "Invalid row access index" );
213  }
214  if( j >= sm_.rows() ) {
215  BLAZE_THROW_OUT_OF_RANGE( "Invalid column access index" );
216  }
217  return (*this)(i,j);
218  }
219  //**********************************************************************************************
220 
221  //**Begin function******************************************************************************
227  inline ConstIterator begin( size_t i ) const {
228  return sm_.begin(i);
229  }
230  //**********************************************************************************************
231 
232  //**End function********************************************************************************
238  inline ConstIterator end( size_t i ) const {
239  return sm_.end(i);
240  }
241  //**********************************************************************************************
242 
243  //**Rows function*******************************************************************************
248  inline size_t rows() const noexcept {
249  return sm_.columns();
250  }
251  //**********************************************************************************************
252 
253  //**Columns function****************************************************************************
258  inline size_t columns() const noexcept {
259  return sm_.rows();
260  }
261  //**********************************************************************************************
262 
263  //**NonZeros function***************************************************************************
268  inline size_t nonZeros() const {
269  return sm_.nonZeros();
270  }
271  //**********************************************************************************************
272 
273  //**NonZeros function***************************************************************************
279  inline size_t nonZeros( size_t i ) const {
280  return sm_.nonZeros( i );
281  }
282  //**********************************************************************************************
283 
284  //**Find function*******************************************************************************
291  inline ConstIterator find( size_t i, size_t j ) const {
293  return sm_.find( j, i );
294  }
295  //**********************************************************************************************
296 
297  //**LowerBound function*************************************************************************
304  inline ConstIterator lowerBound( size_t i, size_t j ) const {
306  return sm_.lowerBound( j, i );
307  }
308  //**********************************************************************************************
309 
310  //**UpperBound function*************************************************************************
317  inline ConstIterator upperBound( size_t i, size_t j ) const {
319  return sm_.upperBound( j, i );
320  }
321  //**********************************************************************************************
322 
323  //**Operand access******************************************************************************
328  inline Operand operand() const noexcept {
329  return sm_;
330  }
331  //**********************************************************************************************
332 
333  //**********************************************************************************************
339  template< typename T >
340  inline bool canAlias( const T* alias ) const noexcept {
341  return sm_.isAliased( alias );
342  }
343  //**********************************************************************************************
344 
345  //**********************************************************************************************
351  template< typename T >
352  inline bool isAliased( const T* alias ) const noexcept {
353  return sm_.isAliased( alias );
354  }
355  //**********************************************************************************************
356 
357  //**********************************************************************************************
362  inline bool canSMPAssign() const noexcept {
363  return sm_.canSMPAssign();
364  }
365  //**********************************************************************************************
366 
367  private:
368  //**Member variables****************************************************************************
370  //**********************************************************************************************
371 
372  //**Assignment to dense matrices****************************************************************
386  template< typename MT2 // Type of the target dense matrix
387  , bool SO2 > // Storage order of the target dense matrix
388  friend inline EnableIf_< UseAssign<MT2> >
389  assign( DenseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
390  {
392 
393  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
394  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
395 
396  DMatTransposer<MT2,!SO2> tmp( ~lhs );
397  assign( tmp, rhs.sm_ );
398  }
400  //**********************************************************************************************
401 
402  //**Assignment to sparse matrices***************************************************************
416  template< typename MT2 // Type of the target sparse matrix
417  , bool SO2 > // Storage order of the target sparse matrix
418  friend inline EnableIf_< UseAssign<MT2> >
419  assign( SparseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
420  {
422 
423  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
424  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
425 
426  SMatTransposer<MT2,!SO2> tmp( ~lhs );
427  assign( tmp, rhs.sm_ );
428  }
430  //**********************************************************************************************
431 
432  //**Addition assignment to dense matrices*******************************************************
446  template< typename MT2 // Type of the target dense matrix
447  , bool SO2 > // Storage order of the target dense matrix
448  friend inline EnableIf_< UseAssign<MT2> >
449  addAssign( DenseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
450  {
452 
453  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
454  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
455 
456  DMatTransposer<MT2,!SO2> tmp( ~lhs );
457  addAssign( tmp, rhs.sm_ );
458  }
460  //**********************************************************************************************
461 
462  //**Addition assignment to sparse matrices******************************************************
463  // No special implementation for the addition assignment to sparse matrices.
464  //**********************************************************************************************
465 
466  //**Subtraction 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 EnableIf_< UseAssign<MT2> >
483  subAssign( DenseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
484  {
486 
487  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
488  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
489 
490  DMatTransposer<MT2,!SO2> tmp( ~lhs );
491  subAssign( tmp, rhs.sm_ );
492  }
494  //**********************************************************************************************
495 
496  //**Subtraction assignment to sparse matrices***************************************************
497  // No special implementation for the subtraction assignment to sparse matrices.
498  //**********************************************************************************************
499 
500  //**Schur product assignment to dense matrices**************************************************
514  template< typename MT2 // Type of the target dense matrix
515  , bool SO2 > // Storage order of the target dense matrix
516  friend inline EnableIf_< UseAssign<MT2> >
517  schurAssign( DenseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
518  {
520 
521  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
522  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
523 
524  DMatTransposer<MT2,!SO2> tmp( ~lhs );
525  schurAssign( tmp, rhs.sm_ );
526  }
528  //**********************************************************************************************
529 
530  //**Schur product assignment to sparse matrices*************************************************
531  // No special implementation for the Schur product assignment to sparse matrices.
532  //**********************************************************************************************
533 
534  //**Multiplication assignment to dense matrices*************************************************
535  // No special implementation for the multiplication assignment to dense matrices.
536  //**********************************************************************************************
537 
538  //**Multiplication assignment to sparse matrices************************************************
539  // No special implementation for the multiplication assignment to sparse matrices.
540  //**********************************************************************************************
541 
542  //**SMP assignment to dense matrices************************************************************
556  template< typename MT2 // Type of the target dense matrix
557  , bool SO2 > // Storage order of the target dense matrix
558  friend inline EnableIf_< UseSMPAssign<MT2> >
559  smpAssign( DenseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
560  {
562 
563  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
564  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
565 
566  DMatTransposer<MT2,!SO2> tmp( ~lhs );
567  smpAssign( tmp, rhs.sm_ );
568  }
570  //**********************************************************************************************
571 
572  //**SMP assignment to sparse matrices***********************************************************
586  template< typename MT2 // Type of the target sparse matrix
587  , bool SO2 > // Storage order of the target sparse matrix
588  friend inline EnableIf_< UseSMPAssign<MT2> >
589  smpAssign( SparseMatrix<MT2,SO2>& lhs, const SMatTransExpr& rhs )
590  {
592 
593  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
594  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
595 
596  SMatTransposer<MT2,!SO2> tmp( ~lhs );
597  smpAssign( tmp, rhs.sm_ );
598  }
600  //**********************************************************************************************
601 
602  //**SMP addition assignment to dense matrices***************************************************
616  template< typename MT2 // Type of the target dense matrix
617  , bool SO2 > // Storage order of the target dense matrix
618  friend inline EnableIf_< UseSMPAssign<MT2> >
620  {
622 
623  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
624  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
625 
626  DMatTransposer<MT2,!SO2> tmp( ~lhs );
627  smpAddAssign( tmp, rhs.sm_ );
628  }
630  //**********************************************************************************************
631 
632  //**SMP addition assignment to sparse matrices**************************************************
633  // No special implementation for the SMP addition assignment to sparse matrices.
634  //**********************************************************************************************
635 
636  //**SMP subtraction assignment to dense matrices************************************************
651  template< typename MT2 // Type of the target dense matrix
652  , bool SO2 > // Storage order of the target dense matrix
653  friend inline EnableIf_< UseSMPAssign<MT2> >
655  {
657 
658  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
659  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
660 
661  DMatTransposer<MT2,!SO2> tmp( ~lhs );
662  smpSubAssign( tmp, rhs.sm_ );
663  }
665  //**********************************************************************************************
666 
667  //**SMP subtraction assignment to sparse matrices***********************************************
668  // No special implementation for the SMP subtraction assignment to sparse matrices.
669  //**********************************************************************************************
670 
671  //**SMP Schur product assignment to dense matrices**********************************************
686  template< typename MT2 // Type of the target dense matrix
687  , bool SO2 > // Storage order of the target dense matrix
688  friend inline EnableIf_< UseSMPAssign<MT2> >
690  {
692 
693  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
694  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
695 
696  DMatTransposer<MT2,!SO2> tmp( ~lhs );
697  smpSchurAssign( tmp, rhs.sm_ );
698  }
700  //**********************************************************************************************
701 
702  //**SMP Schur product assignment to sparse matrices*********************************************
703  // No special implementation for the SMP Schur product assignment to sparse matrices.
704  //**********************************************************************************************
705 
706  //**SMP multiplication assignment to dense matrices*********************************************
707  // No special implementation for the SMP multiplication assignment to dense matrices.
708  //**********************************************************************************************
709 
710  //**SMP multiplication assignment to sparse matrices********************************************
711  // No special implementation for the SMP multiplication assignment to sparse matrices.
712  //**********************************************************************************************
713 
714  //**Compile time checks*************************************************************************
719  //**********************************************************************************************
720 };
721 //*************************************************************************************************
722 
723 
724 
725 
726 //=================================================================================================
727 //
728 // GLOBAL OPERATORS
729 //
730 //=================================================================================================
731 
732 //*************************************************************************************************
751 template< typename MT // Type of the sparse matrix
752  , bool SO > // Storage order
753 inline decltype(auto) trans( const SparseMatrix<MT,SO>& sm )
754 {
756 
757  using ReturnType = const SMatTransExpr<MT,!SO>;
758  return ReturnType( ~sm );
759 }
760 //*************************************************************************************************
761 
762 
763 
764 
765 //=================================================================================================
766 //
767 // GLOBAL RESTRUCTURING FUNCTIONS
768 //
769 //=================================================================================================
770 
771 //*************************************************************************************************
791 template< typename MT // Type of the dense matrix
792  , bool SO > // Storage order
793 inline decltype(auto) trans( const SMatTransExpr<MT,SO>& sm )
794 {
796 
797  return sm.operand();
798 }
800 //*************************************************************************************************
801 
802 
803 //*************************************************************************************************
815 template< typename MT // Type of the left-hand side sparse matrix
816  , typename ST // Type of the right-hand side scalar value
817  , bool SO > // Storage order
818 inline decltype(auto) trans( const SMatScalarMultExpr<MT,ST,SO>& sm )
819 {
821 
822  return trans( sm.leftOperand() ) * sm.rightOperand();
823 }
825 //*************************************************************************************************
826 
827 
828 
829 
830 //=================================================================================================
831 //
832 // SIZE SPECIALIZATIONS
833 //
834 //=================================================================================================
835 
836 //*************************************************************************************************
838 template< typename MT, bool SO >
839 struct Size< SMatTransExpr<MT,SO>, 0UL >
840  : public Size<MT,1UL>
841 {};
842 
843 template< typename MT, bool SO >
844 struct Size< SMatTransExpr<MT,SO>, 1UL >
845  : public Size<MT,0UL>
846 {};
848 //*************************************************************************************************
849 
850 
851 
852 
853 //=================================================================================================
854 //
855 // ISSYMMETRIC SPECIALIZATIONS
856 //
857 //=================================================================================================
858 
859 //*************************************************************************************************
861 template< typename MT, bool SO >
862 struct IsSymmetric< SMatTransExpr<MT,SO> >
863  : public IsSymmetric<MT>
864 {};
866 //*************************************************************************************************
867 
868 
869 
870 
871 //=================================================================================================
872 //
873 // ISHERMITIAN SPECIALIZATIONS
874 //
875 //=================================================================================================
876 
877 //*************************************************************************************************
879 template< typename MT, bool SO >
880 struct IsHermitian< SMatTransExpr<MT,SO> >
881  : public IsHermitian<MT>
882 {};
884 //*************************************************************************************************
885 
886 
887 
888 
889 //=================================================================================================
890 //
891 // ISLOWER SPECIALIZATIONS
892 //
893 //=================================================================================================
894 
895 //*************************************************************************************************
897 template< typename MT, bool SO >
898 struct IsLower< SMatTransExpr<MT,SO> >
899  : public IsUpper<MT>
900 {};
902 //*************************************************************************************************
903 
904 
905 
906 
907 //=================================================================================================
908 //
909 // ISUNILOWER SPECIALIZATIONS
910 //
911 //=================================================================================================
912 
913 //*************************************************************************************************
915 template< typename MT, bool SO >
916 struct IsUniLower< SMatTransExpr<MT,SO> >
917  : public IsUniUpper<MT>
918 {};
920 //*************************************************************************************************
921 
922 
923 
924 
925 //=================================================================================================
926 //
927 // ISSTRICTLYLOWER SPECIALIZATIONS
928 //
929 //=================================================================================================
930 
931 //*************************************************************************************************
933 template< typename MT, bool SO >
934 struct IsStrictlyLower< SMatTransExpr<MT,SO> >
935  : public IsStrictlyUpper<MT>
936 {};
938 //*************************************************************************************************
939 
940 
941 
942 
943 //=================================================================================================
944 //
945 // ISUPPER SPECIALIZATIONS
946 //
947 //=================================================================================================
948 
949 //*************************************************************************************************
951 template< typename MT, bool SO >
952 struct IsUpper< SMatTransExpr<MT,SO> >
953  : public IsLower<MT>
954 {};
956 //*************************************************************************************************
957 
958 
959 
960 
961 //=================================================================================================
962 //
963 // ISUNIUPPER SPECIALIZATIONS
964 //
965 //=================================================================================================
966 
967 //*************************************************************************************************
969 template< typename MT, bool SO >
970 struct IsUniUpper< SMatTransExpr<MT,SO> >
971  : public IsUniLower<MT>
972 {};
974 //*************************************************************************************************
975 
976 
977 
978 
979 //=================================================================================================
980 //
981 // ISSTRICTLYUPPER SPECIALIZATIONS
982 //
983 //=================================================================================================
984 
985 //*************************************************************************************************
987 template< typename MT, bool SO >
988 struct IsStrictlyUpper< SMatTransExpr<MT,SO> >
989  : public IsStrictlyLower<MT>
990 {};
992 //*************************************************************************************************
993 
994 } // namespace blaze
995 
996 #endif
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row/column.
Definition: SMatTransExpr.h:279
TransposeType_< MT > ResultType
Result type for expression template evaluations.
Definition: SMatTransExpr.h:157
Header file for the sparse matrix transposer.
Header file for auxiliary alias declarations.
IfTrue_< useAssign, const ResultType, const SMatTransExpr &> CompositeType
Data type for composite expression templates.
Definition: SMatTransExpr.h:164
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
RightOperand rightOperand() const noexcept
Returns the right-hand side scalar operand.
Definition: SMatScalarMultExpr.h:475
Header file for basic type definitions.
ResultType_< MT > TransposeType
Transpose type for expression template evaluations.
Definition: SMatTransExpr.h:159
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SMatTransExpr.h:268
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
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: SMatTransExpr.h:210
#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
SMatTransExpr(const MT &sm) noexcept
Constructor for the SMatTransExpr class.
Definition: SMatTransExpr.h:183
Header file for the MatTransExpr base class.
If_< IsExpression< MT >, const MT, const MT &> Operand
Composite data type of the sparse matrix expression.
Definition: SMatTransExpr.h:170
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SMatTransExpr.h:195
ElementType_< MT > ElementType
Resulting element type.
Definition: SMatTransExpr.h:160
LeftOperand leftOperand() const noexcept
Returns the left-hand side sparse matrix operand.
Definition: SMatScalarMultExpr.h:465
Compile time check for lower triangular matrices.This type trait tests whether or not the given templ...
Definition: IsLower.h:87
#define BLAZE_CREATE_HAS_TYPE_MEMBER_TYPE_TRAIT(TYPE_TRAIT_NAME, MEMBER_NAME)
Macro for the creation of a type trait for compile time checks for member types.This macro creates th...
Definition: HasMember.h:182
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:87
Header file for the RequiresEvaluation type trait.
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:80
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:129
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatTransExpr.h:352
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: SMatTransExpr.h:362
typename IfTrue< Condition, T1, T2 >::Type IfTrue_
Auxiliary alias declaration for the IfTrue class template.The IfTrue_ alias declaration provides a co...
Definition: If.h:109
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 to query the requirement to evaluate an expression.Via this type trait it is possi...
Definition: RequiresEvaluation.h:71
typename T::CompositeType CompositeType_
Alias declaration for nested CompositeType type definitions.The CompositeType_ alias declaration prov...
Definition: Aliases.h:83
Compile time check for upper unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniUpper.h:86
ResultType_< MT > RT
Result type of the sparse matrix expression.
Definition: SMatTransExpr.h:101
Operand sm_
Sparse matrix of the transposition expression.
Definition: SMatTransExpr.h:369
Header file for the IsStrictlyUpper type trait.
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: SMatTransExpr.h:258
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row/column i.
Definition: SMatTransExpr.h:227
Header file for the IsSymmetric type trait.
OppositeType_< ResultType > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatTransExpr.h:158
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
Header file for the If class template.
Header file for the Transformation base class.
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
CompositeType_< MT > CT
Composite type of the sparse matrix expression.
Definition: SMatTransExpr.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
ReturnType_< MT > ReturnType
Return type for expression template evaluations.
Definition: SMatTransExpr.h:161
const Element * ConstIterator
Iterator over constant elements.
Definition: CompressedMatrix.h:3085
typename T::ElementType ElementType_
Alias declaration for nested ElementType type definitions.The ElementType_ alias declaration provides...
Definition: Aliases.h:163
Expression object for the transposition of a dense matrix.The DMatTransposer class is a wrapper objec...
Definition: DMatTransposer.h:77
ConstIterator upperBound(size_t i, size_t j) const
Returns an iterator to the first index greater then the given index.
Definition: SMatTransExpr.h:317
Header file for the IsLower type trait.
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
Constraint on the data type.
Expression object for sparse matrix-scalar multiplications.The SMatScalarMult class represents the co...
Definition: Forward.h:114
Header file for all forward declarations for expression class templates.
Header file for the EnableIf class template.
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
Header file for the dense matrix transposer.
Expression object for the transposition of a sparse matrix.The SMatTransposer class is a wrapper obje...
Definition: Forward.h:124
typename GetConstIterator< MT >::Type ConstIterator
Iterator over the elements of the dense matrix.
Definition: SMatTransExpr.h:167
ConstIterator find(size_t i, size_t j) const
Searches for a specific matrix element.
Definition: SMatTransExpr.h:291
Header file for run time assertion macros.
Utility type for generic codes.
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
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row/column i.
Definition: SMatTransExpr.h:238
#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
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SMatTransExpr.h:340
Header file for the HasMember type traits.
#define BLAZE_CONSTRAINT_MUST_NOT_REQUIRE_EVALUATION(T)
Constraint on the data type.In case the given data type T requires an intermediate evaluation within ...
Definition: RequiresEvaluation.h:81
Header file for the RemoveReference type trait.
typename EnableIf< Condition, T >::Type EnableIf_
Auxiliary alias declaration for the EnableIf class template.The EnableIf_ alias declaration provides ...
Definition: EnableIf.h:224
typename T::OppositeType OppositeType_
Alias declaration for nested OppositeType type definitions.The OppositeType_ alias declaration provid...
Definition: Aliases.h:263
Operand operand() const noexcept
Returns the sparse matrix operand.
Definition: SMatTransExpr.h:328
Compile time check for strictly lower triangular matrices.This type trait tests whether or not the gi...
Definition: IsStrictlyLower.h:86
ConstIterator lowerBound(size_t i, size_t j) const
Returns an iterator to the first index not less then the given index.
Definition: SMatTransExpr.h:304
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:3080
decltype(auto) trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:789
Header file for the IsComputation type trait class.
Compile time evaluation of the size of vectors and matrices.The Size type trait evaluates the size of...
Definition: Size.h:80
Expression object for sparse matrix transpositions.The SMatTransExpr class represents the compile tim...
Definition: Forward.h:123
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: SMatTransExpr.h:248
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.
Header file for the IsHermitian type trait.
Header file for the Size 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.