SMatSMatAddExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_SMATSMATADDEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_SMATSMATADDEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <blaze/math/Aliases.h>
49 #include <blaze/math/Exception.h>
75 #include <blaze/util/Assert.h>
76 #include <blaze/util/DisableIf.h>
77 #include <blaze/util/EnableIf.h>
80 #include <blaze/util/mpl/And.h>
81 #include <blaze/util/mpl/If.h>
82 #include <blaze/util/mpl/Max.h>
83 #include <blaze/util/mpl/Or.h>
84 #include <blaze/util/Types.h>
86 
87 
88 namespace blaze {
89 
90 //=================================================================================================
91 //
92 // CLASS SMATSMATADDEXPR
93 //
94 //=================================================================================================
95 
96 //*************************************************************************************************
103 template< typename MT1 // Type of the left-hand side sparse matrix
104  , typename MT2 > // Type of the right-hand side sparse matrix
105 class SMatSMatAddExpr : public SparseMatrix< SMatSMatAddExpr<MT1,MT2>, false >
106  , private MatMatAddExpr
107  , private Computation
108 {
109  private:
110  //**Type definitions****************************************************************************
117  //**********************************************************************************************
118 
119  //**Return type evaluation**********************************************************************
121 
126  enum : bool { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
127 
130  //**********************************************************************************************
131 
132  //**Serial evaluation strategy******************************************************************
134 
139  template< typename T1, typename T2, typename T3 >
140  struct UseSymmetricKernel {
142  enum : bool { value = IsSymmetric<T2>::value && IsSymmetric<T3>::value };
143  };
145  //**********************************************************************************************
146 
147  //**Parallel evaluation strategy****************************************************************
149 
154  template< typename MT >
155  struct UseSMPAssign {
156  enum : bool { value = MT::smpAssignable };
157  };
159  //**********************************************************************************************
160 
161  public:
162  //**Type definitions****************************************************************************
168 
171 
173  typedef const ResultType CompositeType;
174 
176  typedef If_< IsExpression<MT1>, const MT1, const MT1& > LeftOperand;
177 
179  typedef If_< IsExpression<MT2>, const MT2, const MT2& > RightOperand;
180  //**********************************************************************************************
181 
182  //**Compilation flags***************************************************************************
184  enum : bool { smpAssignable = false };
185  //**********************************************************************************************
186 
187  //**Constructor*********************************************************************************
193  explicit inline SMatSMatAddExpr( const MT1& lhs, const MT2& rhs ) noexcept
194  : lhs_( lhs ) // Left-hand side sparse matrix of the addition expression
195  , rhs_( rhs ) // Right-hand side sparse matrix of the addition expression
196  {
197  BLAZE_INTERNAL_ASSERT( lhs.rows() == rhs.rows() , "Invalid number of rows" );
198  BLAZE_INTERNAL_ASSERT( lhs.columns() == rhs.columns(), "Invalid number of columns" );
199  }
200  //**********************************************************************************************
201 
202  //**Access operator*****************************************************************************
209  inline ReturnType operator()( size_t i, size_t j ) const {
210  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
211  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
212  return lhs_(i,j) + rhs_(i,j);
213  }
214  //**********************************************************************************************
215 
216  //**At function*********************************************************************************
224  inline ReturnType at( size_t i, size_t j ) const {
225  if( i >= lhs_.rows() ) {
226  BLAZE_THROW_OUT_OF_RANGE( "Invalid row access index" );
227  }
228  if( j >= lhs_.columns() ) {
229  BLAZE_THROW_OUT_OF_RANGE( "Invalid column access index" );
230  }
231  return (*this)(i,j);
232  }
233  //**********************************************************************************************
234 
235  //**Rows function*******************************************************************************
240  inline size_t rows() const noexcept {
241  return lhs_.rows();
242  }
243  //**********************************************************************************************
244 
245  //**Columns function****************************************************************************
250  inline size_t columns() const noexcept {
251  return lhs_.columns();
252  }
253  //**********************************************************************************************
254 
255  //**NonZeros function***************************************************************************
260  inline size_t nonZeros() const {
261  return lhs_.nonZeros() + rhs_.nonZeros();
262  }
263  //**********************************************************************************************
264 
265  //**NonZeros function***************************************************************************
271  inline size_t nonZeros( size_t i ) const {
272  return lhs_.nonZeros(i) + rhs_.nonZeros(i);
273  }
274  //**********************************************************************************************
275 
276  //**Left operand access*************************************************************************
281  inline LeftOperand leftOperand() const noexcept {
282  return lhs_;
283  }
284  //**********************************************************************************************
285 
286  //**Right operand access************************************************************************
291  inline RightOperand rightOperand() const noexcept {
292  return rhs_;
293  }
294  //**********************************************************************************************
295 
296  //**********************************************************************************************
302  template< typename T >
303  inline bool canAlias( const T* alias ) const noexcept {
304  return ( lhs_.canAlias( alias ) || rhs_.canAlias( alias ) );
305  }
306  //**********************************************************************************************
307 
308  //**********************************************************************************************
314  template< typename T >
315  inline bool isAliased( const T* alias ) const noexcept {
316  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
317  }
318  //**********************************************************************************************
319 
320  private:
321  //**Member variables****************************************************************************
322  LeftOperand lhs_;
323  RightOperand rhs_;
324  //**********************************************************************************************
325 
326  //**Assignment to dense matrices****************************************************************
338  template< typename MT // Type of the target dense matrix
339  , bool SO > // Storage order of the target dense matrix
340  friend inline void assign( DenseMatrix<MT,SO>& lhs, const SMatSMatAddExpr& rhs )
341  {
343 
344  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
345  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
346 
347  typedef ConstIterator_< RemoveReference_<CT2> > RightIterator;
348 
349  assign( ~lhs, rhs.lhs_ );
350 
351  if( !IsResizable< ElementType_<MT> >::value ) {
352  addAssign( ~lhs, rhs.rhs_ );
353  }
354  else
355  {
356  CT2 B( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side sparse matrix operand
357 
358  BLAZE_INTERNAL_ASSERT( B.rows() == rhs.rhs_.rows() , "Invalid number of rows" );
359  BLAZE_INTERNAL_ASSERT( B.columns() == rhs.rhs_.columns(), "Invalid number of columns" );
360  BLAZE_INTERNAL_ASSERT( B.rows() == (~lhs).rows() , "Invalid number of rows" );
361  BLAZE_INTERNAL_ASSERT( B.columns() == (~lhs).columns() , "Invalid number of columns" );
362 
363  for( size_t i=0UL; i<(~lhs).rows(); ++i ) {
364  const RightIterator end( B.end(i) );
365  for( RightIterator element=B.begin(i); element!=end; ++element ) {
366  if( isDefault( (~lhs)(i,element->index()) ) )
367  (~lhs)(i,element->index()) = element->value();
368  else
369  (~lhs)(i,element->index()) += element->value();
370  }
371  }
372  }
373  }
375  //**********************************************************************************************
376 
377  //**Assignment to row-major sparse matrices*****************************************************
389  template< typename MT > // Type of the target sparse matrix
390  friend inline void assign( SparseMatrix<MT,false>& lhs, const SMatSMatAddExpr& rhs )
391  {
393 
394  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
395  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
396 
397  typedef ConstIterator_< RemoveReference_<CT1> > LeftIterator;
398  typedef ConstIterator_< RemoveReference_<CT2> > RightIterator;
399 
400  CT1 A( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse matrix operand
401  CT2 B( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side sparse matrix operand
402 
403  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.lhs_.rows() , "Invalid number of rows" );
404  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.lhs_.columns(), "Invalid number of columns" );
405  BLAZE_INTERNAL_ASSERT( B.rows() == rhs.rhs_.rows() , "Invalid number of rows" );
406  BLAZE_INTERNAL_ASSERT( B.columns() == rhs.rhs_.columns(), "Invalid number of columns" );
407  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).rows() , "Invalid number of rows" );
408  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).columns() , "Invalid number of columns" );
409 
410  // Final memory allocation (based on the evaluated operands)
411  (~lhs).reserve( A.nonZeros() + B.nonZeros() );
412 
413  // Performing the matrix addition
414  for( size_t i=0UL; i<(~lhs).rows(); ++i )
415  {
416  const LeftIterator lend( A.end(i) );
417  const RightIterator rend( B.end(i) );
418 
419  LeftIterator l( A.begin(i) );
420  RightIterator r( B.begin(i) );
421 
422  while( l != lend && r != rend )
423  {
424  if( l->index() < r->index() ) {
425  (~lhs).append( i, l->index(), l->value() );
426  ++l;
427  }
428  else if( l->index() > r->index() ) {
429  (~lhs).append( i, r->index(), r->value() );
430  ++r;
431  }
432  else {
433  (~lhs).append( i, l->index(), l->value()+r->value() );
434  ++l;
435  ++r;
436  }
437  }
438 
439  while( l != lend ) {
440  (~lhs).append( i, l->index(), l->value() );
441  ++l;
442  }
443 
444  while( r != rend ) {
445  (~lhs).append( i, r->index(), r->value() );
446  ++r;
447  }
448 
449  (~lhs).finalize( i );
450  }
451  }
453  //**********************************************************************************************
454 
455  //**Assignment to column-major sparse matrices**************************************************
467  template< typename MT > // Type of the target sparse matrix
469  assign( SparseMatrix<MT,true>& lhs, const SMatSMatAddExpr& rhs )
470  {
472 
474 
475  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
476  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
477 
478  typedef ConstIterator_< RemoveReference_<CT1> > LeftIterator;
479  typedef ConstIterator_< RemoveReference_<CT2> > RightIterator;
480 
481  CT1 A( serial( rhs.lhs_ ) ); // Evaluation of the left-hand side sparse matrix operand
482  CT2 B( serial( rhs.rhs_ ) ); // Evaluation of the right-hand side sparse matrix operand
483 
484  BLAZE_INTERNAL_ASSERT( A.rows() == rhs.lhs_.rows() , "Invalid number of rows" );
485  BLAZE_INTERNAL_ASSERT( A.columns() == rhs.lhs_.columns(), "Invalid number of columns" );
486  BLAZE_INTERNAL_ASSERT( B.rows() == rhs.rhs_.rows() , "Invalid number of rows" );
487  BLAZE_INTERNAL_ASSERT( B.columns() == rhs.rhs_.columns(), "Invalid number of columns" );
488  BLAZE_INTERNAL_ASSERT( A.rows() == (~lhs).rows() , "Invalid number of rows" );
489  BLAZE_INTERNAL_ASSERT( A.columns() == (~lhs).columns() , "Invalid number of columns" );
490 
491  const size_t m( A.rows() );
492  const size_t n( A.columns() );
493 
494  // Counting the number of elements per column
495  std::vector<size_t> nonzeros( n, 0UL );
496  for( size_t i=0UL; i<m; ++i )
497  {
498  const LeftIterator lend( A.end(i) );
499  const RightIterator rend( B.end(i) );
500 
501  LeftIterator l( A.begin(i) );
502  RightIterator r( B.begin(i) );
503 
504  while( l != lend && r != rend )
505  {
506  if( l->index() < r->index() ) {
507  ++nonzeros[l->index()];
508  ++l;
509  }
510  else if( l->index() > r->index() ) {
511  ++nonzeros[r->index()];
512  ++r;
513  }
514  else {
515  ++nonzeros[l->index()];
516  ++l;
517  ++r;
518  }
519  }
520 
521  while( l != lend ) {
522  ++nonzeros[l->index()];
523  ++l;
524  }
525 
526  while( r != rend ) {
527  ++nonzeros[r->index()];
528  ++r;
529  }
530  }
531 
532  // Resizing the left-hand side sparse matrix
533  for( size_t j=0UL; j<n; ++j ) {
534  (~lhs).reserve( j, nonzeros[j] );
535  }
536 
537  // Performing the matrix addition
538  for( size_t i=0UL; i<m; ++i )
539  {
540  const LeftIterator lend( A.end(i) );
541  const RightIterator rend( B.end(i) );
542 
543  LeftIterator l( A.begin(i) );
544  RightIterator r( B.begin(i) );
545 
546  while( l != lend && r != rend )
547  {
548  if( l->index() < r->index() ) {
549  (~lhs).append( i, l->index(), l->value() );
550  ++l;
551  }
552  else if( l->index() > r->index() ) {
553  (~lhs).append( i, r->index(), r->value() );
554  ++r;
555  }
556  else {
557  (~lhs).append( i, l->index(), l->value()+r->value() );
558  ++l;
559  ++r;
560  }
561  }
562 
563  while( l != lend ) {
564  (~lhs).append( i, l->index(), l->value() );
565  ++l;
566  }
567 
568  while( r != rend ) {
569  (~lhs).append( i, r->index(), r->value() );
570  ++r;
571  }
572  }
573  }
575  //**********************************************************************************************
576 
577  //**Assignment to column-major sparse matrices**************************************************
589  template< typename MT > // Type of the target sparse matrix
591  assign( SparseMatrix<MT,true>& lhs, const SMatSMatAddExpr& rhs )
592  {
594 
596 
597  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
598  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
599 
600  assign( ~lhs, trans( rhs.lhs_ ) + trans( rhs.rhs_ ) );
601  }
603  //**********************************************************************************************
604 
605  //**Addition assignment to dense matrices*******************************************************
617  template< typename MT // Type of the target dense matrix
618  , bool SO > // Storage order of the target dense matrix
619  friend inline void addAssign( DenseMatrix<MT,SO>& lhs, const SMatSMatAddExpr& rhs )
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  addAssign( ~lhs, rhs.lhs_ );
627  addAssign( ~lhs, rhs.rhs_ );
628  }
630  //**********************************************************************************************
631 
632  //**Addition assignment to sparse matrices******************************************************
633  // No special implementation for the addition assignment to sparse matrices.
634  //**********************************************************************************************
635 
636  //**Subtraction assignment to dense matrices****************************************************
648  template< typename MT // Type of the target dense matrix
649  , bool SO > // Storage order of the target dense matrix
650  friend inline void subAssign( DenseMatrix<MT,SO>& lhs, const SMatSMatAddExpr& rhs )
651  {
653 
654  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
655  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
656 
657  subAssign( ~lhs, rhs.lhs_ );
658  subAssign( ~lhs, rhs.rhs_ );
659  }
661  //**********************************************************************************************
662 
663  //**Subtraction assignment to sparse matrices***************************************************
664  // No special implementation for the subtraction assignment to sparse matrices.
665  //**********************************************************************************************
666 
667  //**Multiplication assignment to dense matrices*************************************************
668  // No special implementation for the multiplication assignment to dense matrices.
669  //**********************************************************************************************
670 
671  //**Multiplication assignment to sparse matrices************************************************
672  // No special implementation for the multiplication assignment to sparse matrices.
673  //**********************************************************************************************
674 
675  //**SMP assignment to dense matrices************************************************************
676  // No special implementation for the SMP assignment to dense matrices.
677  //**********************************************************************************************
678 
679  //**SMP assignment to sparse matrices***********************************************************
680  // No special implementation for the SMP assignment to sparse matrices.
681  //**********************************************************************************************
682 
683  //**SMP addition assignment to dense matrices***************************************************
697  template< typename MT // Type of the target dense matrix
698  , bool SO > // Storage order of the target dense matrix
699  friend inline EnableIf_< UseSMPAssign<MT> >
701  {
703 
704  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
705  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
706 
707  smpAddAssign( ~lhs, rhs.lhs_ );
708  smpAddAssign( ~lhs, rhs.rhs_ );
709  }
711  //**********************************************************************************************
712 
713  //**SMP addition assignment to sparse matrices**************************************************
714  // No special implementation for the SMP addition assignment to sparse matrices.
715  //**********************************************************************************************
716 
717  //**SMP subtraction assignment to dense matrices************************************************
731  template< typename MT // Type of the target dense matrix
732  , bool SO > // Storage order of the target dense matrix
733  friend inline EnableIf_< UseSMPAssign<MT> >
735  {
737 
738  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
739  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
740 
741  smpSubAssign( ~lhs, rhs.lhs_ );
742  smpSubAssign( ~lhs, rhs.rhs_ );
743  }
745  //**********************************************************************************************
746 
747  //**SMP subtraction assignment to sparse matrices***********************************************
748  // No special implementation for the SMP subtraction assignment to sparse matrices.
749  //**********************************************************************************************
750 
751  //**SMP multiplication assignment to dense matrices*********************************************
752  // No special implementation for the SMP multiplication assignment to dense matrices.
753  //**********************************************************************************************
754 
755  //**SMP multiplication assignment to sparse matrices********************************************
756  // No special implementation for the SMP multiplication assignment to sparse matrices.
757  //**********************************************************************************************
758 
759  //**Compile time checks*************************************************************************
767  //**********************************************************************************************
768 };
769 //*************************************************************************************************
770 
771 
772 
773 
774 //=================================================================================================
775 //
776 // GLOBAL BINARY ARITHMETIC OPERATORS
777 //
778 //=================================================================================================
779 
780 //*************************************************************************************************
806 template< typename T1 // Type of the left-hand side sparse matrix
807  , typename T2 > // Type of the right-hand side sparse matrix
808 inline const SMatSMatAddExpr<T1,T2>
810 {
812 
813  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() ) {
814  BLAZE_THROW_INVALID_ARGUMENT( "Matrix sizes do not match" );
815  }
816 
817  return SMatSMatAddExpr<T1,T2>( ~lhs, ~rhs );
818 }
819 //*************************************************************************************************
820 
821 
822 
823 
824 //=================================================================================================
825 //
826 // ROWS SPECIALIZATIONS
827 //
828 //=================================================================================================
829 
830 //*************************************************************************************************
832 template< typename MT1, typename MT2 >
833 struct Rows< SMatSMatAddExpr<MT1,MT2> >
834  : public Max< Rows<MT1>, Rows<MT2> >
835 {};
837 //*************************************************************************************************
838 
839 
840 
841 
842 //=================================================================================================
843 //
844 // COLUMNS SPECIALIZATIONS
845 //
846 //=================================================================================================
847 
848 //*************************************************************************************************
850 template< typename MT1, typename MT2 >
851 struct Columns< SMatSMatAddExpr<MT1,MT2> >
852  : public Max< Columns<MT1>, Columns<MT2> >
853 {};
855 //*************************************************************************************************
856 
857 
858 
859 
860 //=================================================================================================
861 //
862 // ISSYMMETRIC SPECIALIZATIONS
863 //
864 //=================================================================================================
865 
866 //*************************************************************************************************
868 template< typename MT1, typename MT2 >
869 struct IsSymmetric< SMatSMatAddExpr<MT1,MT2> >
870  : public BoolConstant< IsSymmetric<MT1>::value && IsSymmetric<MT2>::value >
871 {};
873 //*************************************************************************************************
874 
875 
876 
877 
878 //=================================================================================================
879 //
880 // ISHERMITIAN SPECIALIZATIONS
881 //
882 //=================================================================================================
883 
884 //*************************************************************************************************
886 template< typename MT1, typename MT2 >
887 struct IsHermitian< SMatSMatAddExpr<MT1,MT2> >
888  : public BoolConstant< IsHermitian<MT1>::value && IsHermitian<MT2>::value >
889 {};
891 //*************************************************************************************************
892 
893 
894 
895 
896 //=================================================================================================
897 //
898 // ISLOWER SPECIALIZATIONS
899 //
900 //=================================================================================================
901 
902 //*************************************************************************************************
904 template< typename MT1, typename MT2 >
905 struct IsLower< SMatSMatAddExpr<MT1,MT2> >
906  : public BoolConstant< And< IsLower<MT1>, IsLower<MT2> >::value >
907 {};
909 //*************************************************************************************************
910 
911 
912 
913 
914 //=================================================================================================
915 //
916 // ISUNILOWER SPECIALIZATIONS
917 //
918 //=================================================================================================
919 
920 //*************************************************************************************************
922 template< typename MT1, typename MT2 >
923 struct IsUniLower< SMatSMatAddExpr<MT1,MT2> >
924  : public BoolConstant< Or< And< IsUniLower<MT1>, IsStrictlyLower<MT2> >
925  , And< IsUniLower<MT2>, IsStrictlyLower<MT1> > >::value >
926 {};
928 //*************************************************************************************************
929 
930 
931 
932 
933 //=================================================================================================
934 //
935 // ISSTRICTLYLOWER SPECIALIZATIONS
936 //
937 //=================================================================================================
938 
939 //*************************************************************************************************
941 template< typename MT1, typename MT2 >
942 struct IsStrictlyLower< SMatSMatAddExpr<MT1,MT2> >
943  : public BoolConstant< And< IsStrictlyLower<MT1>, IsStrictlyLower<MT2> >::value >
944 {};
946 //*************************************************************************************************
947 
948 
949 
950 
951 //=================================================================================================
952 //
953 // ISUPPER SPECIALIZATIONS
954 //
955 //=================================================================================================
956 
957 //*************************************************************************************************
959 template< typename MT1, typename MT2 >
960 struct IsUpper< SMatSMatAddExpr<MT1,MT2> >
961  : public BoolConstant< And< IsUpper<MT1>, IsUpper<MT2> >::value >
962 {};
964 //*************************************************************************************************
965 
966 
967 
968 
969 //=================================================================================================
970 //
971 // ISUNIUPPER SPECIALIZATIONS
972 //
973 //=================================================================================================
974 
975 //*************************************************************************************************
977 template< typename MT1, typename MT2 >
978 struct IsUniUpper< SMatSMatAddExpr<MT1,MT2> >
979  : public BoolConstant< Or< And< IsUniUpper<MT1>, IsStrictlyUpper<MT2> >
980  , And< IsUniUpper<MT2>, IsStrictlyUpper<MT1> > >::value >
981 {};
983 //*************************************************************************************************
984 
985 
986 
987 
988 //=================================================================================================
989 //
990 // ISSTRICTLYUPPER SPECIALIZATIONS
991 //
992 //=================================================================================================
993 
994 //*************************************************************************************************
996 template< typename MT1, typename MT2 >
997 struct IsStrictlyUpper< SMatSMatAddExpr<MT1,MT2> >
998  : public BoolConstant< And< IsStrictlyUpper<MT1>, IsStrictlyUpper<MT2> >::value >
999 {};
1001 //*************************************************************************************************
1002 
1003 
1004 
1005 
1006 //=================================================================================================
1007 //
1008 // EXPRESSION TRAIT SPECIALIZATIONS
1009 //
1010 //=================================================================================================
1011 
1012 //*************************************************************************************************
1014 template< typename MT1, typename MT2, bool AF >
1015 struct SubmatrixExprTrait< SMatSMatAddExpr<MT1,MT2>, AF >
1016 {
1017  public:
1018  //**********************************************************************************************
1021  //**********************************************************************************************
1022 };
1024 //*************************************************************************************************
1025 
1026 
1027 //*************************************************************************************************
1029 template< typename MT1, typename MT2 >
1030 struct RowExprTrait< SMatSMatAddExpr<MT1,MT2> >
1031 {
1032  public:
1033  //**********************************************************************************************
1036  //**********************************************************************************************
1037 };
1039 //*************************************************************************************************
1040 
1041 
1042 //*************************************************************************************************
1044 template< typename MT1, typename MT2 >
1045 struct ColumnExprTrait< SMatSMatAddExpr<MT1,MT2> >
1046 {
1047  public:
1048  //**********************************************************************************************
1051  //**********************************************************************************************
1052 };
1054 //*************************************************************************************************
1055 
1056 } // namespace blaze
1057 
1058 #endif
typename SubmatrixExprTrait< MT, AF >::Type SubmatrixExprTrait_
Auxiliary alias declaration for the SubmatrixExprTrait type trait.The SubmatrixExprTrait_ alias decla...
Definition: SubmatrixExprTrait.h:134
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.This macro encapsulates the default way o...
Definition: Exception.h:235
Header file for auxiliary alias declarations.
Header file for the Max class template.
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:70
Header file for the Rows type trait.
Header file for the IsUniUpper type trait.
Header file for basic type definitions.
EnableIf_< IsDenseMatrix< MT1 > > smpSubAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:160
Header file for the serial shim.
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.
const DenseIterator< Type, AF > operator+(const DenseIterator< Type, AF > &it, ptrdiff_t inc) noexcept
Addition between a DenseIterator and an integral value.
Definition: DenseIterator.h:699
typename DisableIf< Condition, T >::Type DisableIf_
Auxiliary type for the DisableIf class template.The DisableIf_ alias declaration provides a convenien...
Definition: DisableIf.h:223
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 AddExprTrait class template.
size_t nonZeros(size_t i) const
Returns the number of non-zero elements in the specified row.
Definition: SMatSMatAddExpr.h:271
SMatSMatAddExpr< MT1, MT2 > This
Type of this SMatSMatAddExpr instance.
Definition: SMatSMatAddExpr.h:163
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:721
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: SMatSMatAddExpr.h:224
Header file for the Computation base class.
Compile time check for upper triangular matrices.This type trait tests whether or not the given templ...
Definition: IsUpper.h:88
Constraints on the storage order of matrix types.
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
SMatSMatAddExpr(const MT1 &lhs, const MT2 &rhs) noexcept
Constructor for the SMatSMatAddExpr class.
Definition: SMatSMatAddExpr.h:193
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
ElementType_< ResultType > ElementType
Resulting element type.
Definition: SMatSMatAddExpr.h:167
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:343
RightOperand rightOperand() const noexcept
Returns the right-hand side sparse matrix operand.
Definition: SMatSMatAddExpr.h:291
Header file for the SparseMatrix base class.
Constraint on the data type.
If_< IsExpression< MT1 >, const MT1, const MT1 &> LeftOperand
Composite type of the left-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:176
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
Header file for the DisableIf class template.
typename RowExprTrait< MT >::Type RowExprTrait_
Auxiliary alias declaration for the RowExprTrait type trait.The RowExprTrait_ alias declaration provi...
Definition: RowExprTrait.h:134
Header file for the IsTemporary type trait class.
Header file for the IsStrictlyUpper type trait.
Header file for the IsSymmetric type trait.
typename AddExprTrait< T1, T2 >::Type AddExprTrait_
Auxiliary alias declaration for the AddExprTrait class template.The AddExprTrait_ alias declaration p...
Definition: AddExprTrait.h:219
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
Header file for the If class template.
#define BLAZE_CONSTRAINT_MUST_BE_COLUMN_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a column-major dense or sparse matri...
Definition: ColumnMajorMatrix.h:61
Header file for the Or class template.
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exception.This macro encapsulates the default way of Bl...
Definition: Exception.h:331
size_t nonZeros() const
Returns the number of non-zero elements in the sparse matrix.
Definition: SMatSMatAddExpr.h:260
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
CompositeType_< MT2 > CT2
Composite type of the right-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:114
Evaluation of the expression type type of a submatrix operation.Via this type trait it is possible to...
Definition: SubmatrixExprTrait.h:80
Header file for the IsLower type trait.
typename ColumnExprTrait< MT >::Type ColumnExprTrait_
Auxiliary alias declaration for the ColumnExprTrait type trait.The ColumnExprTrait_ alias declaration...
Definition: ColumnExprTrait.h:133
LeftOperand leftOperand() const noexcept
Returns the left-hand side sparse matrix operand.
Definition: SMatSMatAddExpr.h:281
RightOperand rhs_
Right-hand side sparse matrix of the addition expression.
Definition: SMatSMatAddExpr.h:323
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
OppositeType_< ResultType > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: SMatSMatAddExpr.h:165
Evaluation of the expression type type of a row operation.Via this type trait it is possible to evalu...
Definition: RowExprTrait.h:79
BLAZE_ALWAYS_INLINE MT::Iterator end(Matrix< MT, SO > &matrix, size_t i)
Returns an iterator just past the last element of row/column i.
Definition: Matrix.h:260
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_MATMATADDEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid matrix/matrix ...
Definition: MatMatAddExpr.h:109
Header file for the RowExprTrait class template.
CompositeType_< MT1 > CT1
Composite type of the left-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:113
Header file for all forward declarations for expression class templates.
Header file for the EnableIf class template.
Header file for the IsStrictlyLower type trait.
ResultType_< MT1 > RT1
Result type of the left-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:111
Compile time check for lower unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniLower.h:86
AddTrait_< RT1, RT2 > ResultType
Result type for expression template evaluations.
Definition: SMatSMatAddExpr.h:164
Compile time check for resizable data types.This type trait tests whether the given data type is a re...
Definition: IsResizable.h:75
#define BLAZE_CONSTRAINT_MUST_NOT_BE_SYMMETRIC_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is a symmetric matrix type, a compilation error is created.
Definition: Symmetric.h:79
Header file for the SubmatrixExprTrait class template.
ResultType_< MT2 > RT2
Result type of the right-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:112
#define BLAZE_CONSTRAINT_MUST_BE_ROW_MAJOR_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a row-major dense or sparse matrix t...
Definition: RowMajorMatrix.h:61
Header file for run time assertion macros.
Compile time value evaluation.The Max class template selects the larger of the two given template arg...
Definition: Max.h:72
Header file for the addition trait.
LeftOperand lhs_
Left-hand side sparse matrix of the addition expression.
Definition: SMatSMatAddExpr.h:322
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
TransposeType_< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: SMatSMatAddExpr.h:166
Constraint on the data type.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: SMatSMatAddExpr.h:303
#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
Header file for the isDefault shim.
Compile time check for Hermitian matrices.This type trait tests whether or not the given template par...
Definition: IsHermitian.h:85
Constraint on the data type.
Constraints on the storage order of matrix types.
If_< IsExpression< MT2 >, const MT2, const MT2 &> RightOperand
Composite type of the right-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:179
Header file for the MatMatAddExpr base class.
AddExprTrait_< RN1, RN2 > ExprReturnType
Expression return type for the subscript operator.
Definition: SMatSMatAddExpr.h:129
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:223
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
ReturnType_< MT1 > RN1
Return type of the left-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:115
typename T::ConstIterator ConstIterator_
Alias declaration for nested ConstIterator type definitions.The ConstIterator_ alias declaration prov...
Definition: Aliases.h:103
ReturnType_< MT2 > RN2
Return type of the right-hand side sparse matrix expression.
Definition: SMatSMatAddExpr.h:116
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: SMatSMatAddExpr.h:240
const ResultType CompositeType
Data type for composite expression templates.
Definition: SMatSMatAddExpr.h:173
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:733
Header file for the IsComputation type trait class.
Expression object for sparse matrix-sparse matrix additions.The SMatSMatAddExpr class represents the ...
Definition: Forward.h:108
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: SMatSMatAddExpr.h:315
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
const IfTrue_< returnExpr, ExprReturnType, ElementType > ReturnType
Return type for expression template evaluations.
Definition: SMatSMatAddExpr.h:170
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: SMatSMatAddExpr.h:209
Compile time evaluation of the number of rows of a matrix.The Rows type trait evaluates the number of...
Definition: Rows.h:76
bool isDefault(const DiagonalProxy< MT > &proxy)
Returns whether the represented element is in default state.
Definition: DiagonalProxy.h:573
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.
Header file for the IsHermitian type trait.
Header file for the IsResizable type trait.
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: SMatSMatAddExpr.h:250
typename AddTrait< T1, T2 >::Type AddTrait_
Auxiliary alias declaration for the AddTrait class template.The AddTrait_ alias declaration provides ...
Definition: AddTrait.h:245
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.