All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
DMatTDMatSubExpr.h
Go to the documentation of this file.
1 //=================================================================================================
33 //=================================================================================================
34 
35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATTDMATSUBEXPR_H_
37 
38 
39 //*************************************************************************************************
40 // Includes
41 //*************************************************************************************************
42 
43 #include <stdexcept>
66 #include <blaze/util/Assert.h>
68 #include <blaze/util/DisableIf.h>
69 #include <blaze/util/EnableIf.h>
71 #include <blaze/util/mpl/Max.h>
72 #include <blaze/util/SelectType.h>
73 #include <blaze/util/Types.h>
75 
76 
77 namespace blaze {
78 
79 //=================================================================================================
80 //
81 // CLASS DMATTDMATSUBEXPR
82 //
83 //=================================================================================================
84 
85 //*************************************************************************************************
92 template< typename MT1 // Type of the left-hand side dense matrix
93  , typename MT2 > // Type of the right-hand side dense matrix
94 class DMatTDMatSubExpr : public DenseMatrix< DMatTDMatSubExpr<MT1,MT2>, false >
95  , private MatMatSubExpr
96  , private Computation
97 {
98  private:
99  //**Type definitions****************************************************************************
100  typedef typename MT1::ResultType RT1;
101  typedef typename MT2::ResultType RT2;
102  typedef typename MT1::ReturnType RN1;
103  typedef typename MT2::ReturnType RN2;
104  typedef typename MT1::CompositeType CT1;
105  typedef typename MT2::CompositeType CT2;
106  //**********************************************************************************************
107 
108  //**Return type evaluation**********************************************************************
110 
115  enum { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
116 
119  //**********************************************************************************************
120 
121  //**Serial evaluation strategy******************************************************************
123 
129  enum { useAssign = RequiresEvaluation<MT1>::value || RequiresEvaluation<MT2>::value || !returnExpr };
130 
132  template< typename MT >
134  struct UseAssign {
135  enum { value = useAssign };
136  };
138  //**********************************************************************************************
139 
140  //**Parallel evaluation strategy****************************************************************
142 
148  template< typename MT >
149  struct UseSMPAssign {
150  enum { value = ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign };
151  };
153  //**********************************************************************************************
154 
155  public:
156  //**Type definitions****************************************************************************
162 
165 
167  typedef const ResultType CompositeType;
168 
170  typedef typename SelectType< IsExpression<MT1>::value, const MT1, const MT1& >::Type LeftOperand;
171 
173  typedef typename SelectType< IsExpression<MT2>::value, const MT2, const MT2& >::Type RightOperand;
174  //**********************************************************************************************
175 
176  //**Compilation flags***************************************************************************
178  enum { vectorizable = 0 };
179 
181  enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
182  //**********************************************************************************************
183 
184  //**Constructor*********************************************************************************
190  explicit inline DMatTDMatSubExpr( const MT1& lhs, const MT2& rhs )
191  : lhs_( lhs ) // Left-hand side dense matrix of the subtraction expression
192  , rhs_( rhs ) // Right-hand side dense matrix of the subtraction expression
193  {
194  BLAZE_INTERNAL_ASSERT( lhs.rows() == rhs.rows() , "Invalid number of rows" );
195  BLAZE_INTERNAL_ASSERT( lhs.columns() == rhs.columns(), "Invalid number of columns" );
196  }
197  //**********************************************************************************************
198 
199  //**Access operator*****************************************************************************
206  inline ReturnType operator()( size_t i, size_t j ) const {
207  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
208  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
209  return lhs_(i,j) - rhs_(i,j);
210  }
211  //**********************************************************************************************
212 
213  //**Rows function*******************************************************************************
218  inline size_t rows() const {
219  return lhs_.rows();
220  }
221  //**********************************************************************************************
222 
223  //**Columns function****************************************************************************
228  inline size_t columns() const {
229  return lhs_.columns();
230  }
231  //**********************************************************************************************
232 
233  //**Left operand access*************************************************************************
238  inline LeftOperand leftOperand() const {
239  return lhs_;
240  }
241  //**********************************************************************************************
242 
243  //**Right operand access************************************************************************
248  inline RightOperand rightOperand() const {
249  return rhs_;
250  }
251  //**********************************************************************************************
252 
253  //**********************************************************************************************
259  template< typename T >
260  inline bool canAlias( const T* alias ) const {
261  return ( IsExpression<MT1>::value && ( RequiresEvaluation<MT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
262  ( IsExpression<MT2>::value && ( RequiresEvaluation<MT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
263  }
264  //**********************************************************************************************
265 
266  //**********************************************************************************************
272  template< typename T >
273  inline bool isAliased( const T* alias ) const {
274  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
275  }
276  //**********************************************************************************************
277 
278  //**********************************************************************************************
283  inline bool isAligned() const {
284  return lhs_.isAligned() && rhs_.isAligned();
285  }
286  //**********************************************************************************************
287 
288  //**********************************************************************************************
293  inline bool canSMPAssign() const {
294  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
296  }
297  //**********************************************************************************************
298 
299  private:
300  //**Member variables****************************************************************************
303  //**********************************************************************************************
304 
305  //**Assignment to dense matrices****************************************************************
319  template< typename MT // Type of the target dense matrix
320  , bool SO2 > // Storage order of the target dense matrix
321  friend inline typename DisableIf< UseAssign<MT> >::Type
322  assign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
323  {
325 
326  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
327  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
328 
329  const size_t m( rhs.rows() );
330  const size_t n( rhs.columns() );
331  const size_t block( 16UL );
332 
333  for( size_t ii=0UL; ii<m; ii+=block ) {
334  const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
335  for( size_t jj=0UL; jj<n; jj+=block ) {
336  const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
337  for( size_t i=ii; i<iend; ++i ) {
338  for( size_t j=jj; j<jend; ++j ) {
339  (~lhs)(i,j) = rhs.lhs_(i,j) - rhs.rhs_(i,j);
340  }
341  }
342  }
343  }
344  }
346  //**********************************************************************************************
347 
348  //**Assignment to dense matrices****************************************************************
362  template< typename MT // Type of the target dense matrix
363  , bool SO2 > // Storage order of the target dense matrix
364  friend inline typename EnableIf< UseAssign<MT> >::Type
365  assign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
366  {
368 
369  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
370  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
371 
372  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
373  subAssign( ~lhs, rhs.rhs_ );
374  }
375  else {
376  assign ( ~lhs, rhs.lhs_ );
377  subAssign( ~lhs, rhs.rhs_ );
378  }
379  }
381  //**********************************************************************************************
382 
383  //**Assignment to sparse matrices***************************************************************
395  template< typename MT // Type of the target sparse matrix
396  , bool SO2 > // Storage order of the target sparse matrix
397  friend inline void assign( SparseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
398  {
400 
401  typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
402 
409 
410  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
411  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
412 
413  const TmpType tmp( serial( rhs ) );
414  assign( ~lhs, tmp );
415  }
417  //**********************************************************************************************
418 
419  //**Addition assignment to dense matrices*******************************************************
434  template< typename MT // Type of the target dense matrix
435  , bool SO2 > // Storage order of the target dense matrix
436  friend inline typename DisableIf< UseAssign<MT> >::Type
437  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
438  {
440 
441  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
442  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
443 
444  const size_t m( rhs.rows() );
445  const size_t n( rhs.columns() );
446  const size_t block( 16UL );
447 
448  for( size_t ii=0UL; ii<m; ii+=block ) {
449  const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
450  for( size_t jj=0UL; jj<n; jj+=block ) {
451  const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
452  for( size_t i=ii; i<iend; ++i ) {
453  for( size_t j=jj; j<jend; ++j ) {
454  (~lhs)(i,j) += rhs.lhs_(i,j) - rhs.rhs_(i,j);
455  }
456  }
457  }
458  }
459  }
461  //**********************************************************************************************
462 
463  //**Addition assignment to dense matrices*******************************************************
478  template< typename MT // Type of the target dense matrix
479  , bool SO2 > // Storage order of the target dense matrix
480  friend inline typename EnableIf< UseAssign<MT> >::Type
481  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
482  {
484 
485  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
486  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
487 
488  addAssign( ~lhs, rhs.lhs_ );
489  subAssign( ~lhs, rhs.rhs_ );
490  }
492  //**********************************************************************************************
493 
494  //**Addition assignment to sparse matrices******************************************************
495  // No special implementation for the addition assignment to sparse matrices.
496  //**********************************************************************************************
497 
498  //**Subtraction assignment to dense matrices****************************************************
513  template< typename MT // Type of the target dense matrix
514  , bool SO2 > // Storage order of the target dense matrix
515  friend inline typename DisableIf< UseAssign<MT> >::Type
516  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
517  {
519 
520  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
521  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
522 
523  const size_t m( rhs.rows() );
524  const size_t n( rhs.columns() );
525  const size_t block( 16UL );
526 
527  for( size_t ii=0UL; ii<m; ii+=block ) {
528  const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
529  for( size_t jj=0UL; jj<n; jj+=block ) {
530  const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
531  for( size_t i=ii; i<iend; ++i ) {
532  for( size_t j=jj; j<jend; ++j ) {
533  (~lhs)(i,j) -= rhs.lhs_(i,j) - rhs.rhs_(i,j);
534  }
535  }
536  }
537  }
538  }
540  //**********************************************************************************************
541 
542  //**Subtraction assignment to dense matrices****************************************************
557  template< typename MT // Type of the target dense matrix
558  , bool SO2 > // Storage order of the target dense matrix
559  friend inline typename EnableIf< UseAssign<MT> >::Type
560  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
561  {
563 
564  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
565  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
566 
567  subAssign( ~lhs, rhs.lhs_ );
568  addAssign( ~lhs, rhs.rhs_ );
569  }
571  //**********************************************************************************************
572 
573  //**Subtraction assignment to sparse matrices***************************************************
574  // No special implementation for the subtraction assignment to sparse matrices.
575  //**********************************************************************************************
576 
577  //**Multiplication assignment to dense matrices*************************************************
578  // No special implementation for the multiplication assignment to dense matrices.
579  //**********************************************************************************************
580 
581  //**Multiplication assignment to sparse matrices************************************************
582  // No special implementation for the multiplication assignment to sparse matrices.
583  //**********************************************************************************************
584 
585  //**SMP assignment to dense matrices************************************************************
599  template< typename MT // Type of the target dense matrix
600  , bool SO2 > // Storage order of the target dense matrix
601  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
602  smpAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
603  {
605 
606  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
607  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
608 
609  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
610  smpSubAssign( ~lhs, rhs.rhs_ );
611  }
612  else {
613  smpAssign ( ~lhs, rhs.lhs_ );
614  smpSubAssign( ~lhs, rhs.rhs_ );
615  }
616  }
618  //**********************************************************************************************
619 
620  //**SMP assignment to sparse matrices***********************************************************
634  template< typename MT // Type of the target sparse matrix
635  , bool SO2 > // Storage order of the target sparse matrix
636  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
637  smpAssign( SparseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
638  {
640 
641  typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
642 
649 
650  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
651  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
652 
653  const TmpType tmp( rhs );
654  smpAssign( ~lhs, tmp );
655  }
657  //**********************************************************************************************
658 
659  //**SMP addition assignment to dense matrices***************************************************
674  template< typename MT // Type of the target dense matrix
675  , bool SO2 > // Storage order of the target dense matrix
676  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
677  smpAddAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
678  {
680 
681  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
682  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
683 
684  smpAddAssign( ~lhs, rhs.lhs_ );
685  smpSubAssign( ~lhs, rhs.rhs_ );
686  }
688  //**********************************************************************************************
689 
690  //**SMP addition assignment to sparse matrices**************************************************
691  // No special implementation for the SMP addition assignment to sparse matrices.
692  //**********************************************************************************************
693 
694  //**SMP subtraction assignment to dense matrices************************************************
709  template< typename MT // Type of the target dense matrix
710  , bool SO2 > // Storage order of the target dense matrix
711  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
712  smpSubAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
713  {
715 
716  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
717  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
718 
719  smpSubAssign( ~lhs, rhs.lhs_ );
720  smpAddAssign( ~lhs, rhs.rhs_ );
721  }
723  //**********************************************************************************************
724 
725  //**SMP subtraction assignment to sparse matrices***********************************************
726  // No special implementation for the SMP subtraction assignment to sparse matrices.
727  //**********************************************************************************************
728 
729  //**SMP multiplication assignment to dense matrices*********************************************
730  // No special implementation for the SMP multiplication assignment to dense matrices.
731  //**********************************************************************************************
732 
733  //**SMP multiplication assignment to sparse matrices********************************************
734  // No special implementation for the SMP multiplication assignment to sparse matrices.
735  //**********************************************************************************************
736 
737  //**Compile time checks*************************************************************************
744  //**********************************************************************************************
745 };
746 //*************************************************************************************************
747 
748 
749 
750 
751 //=================================================================================================
752 //
753 // GLOBAL BINARY ARITHMETIC OPERATORS
754 //
755 //=================================================================================================
756 
757 //*************************************************************************************************
786 template< typename T1 // Type of the left-hand side dense matrix
787  , typename T2 > // Type of the right-hand side dense matrix
788 inline typename EnableIf< And< Not< IsSymmetric<T1> >, Not< IsSymmetric<T2> > >
789  , const DMatTDMatSubExpr<T1,T2> >::Type
791 {
793 
794  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
795  throw std::invalid_argument( "Matrix sizes do not match" );
796 
797  return DMatTDMatSubExpr<T1,T2>( ~lhs, ~rhs );
798 }
799 //*************************************************************************************************
800 
801 
802 //*************************************************************************************************
831 template< typename T1 // Type of the left-hand side dense matrix
832  , typename T2 > // Type of the right-hand side dense matrix
833 inline typename EnableIf< And< Not< IsSymmetric<T1> >, Not< IsSymmetric<T2> > >
834  , const DMatTDMatSubExpr<T1,T2> >::Type
836 {
838 
839  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
840  throw std::invalid_argument( "Matrix sizes do not match" );
841 
842  return DMatTDMatSubExpr<T1,T2>( ~lhs, ~rhs );
843 }
844 //*************************************************************************************************
845 
846 
847 
848 
849 //=================================================================================================
850 //
851 // GLOBAL RESTRUCTURING BINARY ARITHMETIC OPERATORS
852 //
853 //=================================================================================================
854 
855 //*************************************************************************************************
869 template< typename T1 // Type of the left-hand side dense matrix
870  , typename T2 > // Type of the right-hand side dense matrix
871 inline typename EnableIf< And< IsSymmetric<T1>, Not< IsSymmetric<T2> > >
872  , const typename SubExprTrait<T1,T2>::Type >::Type
873  operator-( const DenseMatrix<T1,false>& lhs, const DenseMatrix<T2,true>& rhs )
874 {
876 
877  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
878  throw std::invalid_argument( "Matrix sizes do not match" );
879 
880  return trans( ~lhs ) - ~rhs;
881 }
883 //*************************************************************************************************
884 
885 
886 //*************************************************************************************************
900 template< typename T1 // Type of the left-hand side dense matrix
901  , typename T2 > // Type of the right-hand side dense matrix
902 inline typename EnableIf< IsSymmetric<T2>
903  , const typename SubExprTrait<T1,T2>::Type >::Type
904  operator-( const DenseMatrix<T1,false>& lhs, const DenseMatrix<T2,true>& rhs )
905 {
907 
908  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
909  throw std::invalid_argument( "Matrix sizes do not match" );
910 
911  return (~lhs) - trans( ~rhs );
912 }
914 //*************************************************************************************************
915 
916 
917 //*************************************************************************************************
931 template< typename T1 // Type of the left-hand side dense matrix
932  , typename T2 > // Type of the right-hand side dense matrix
933 inline typename EnableIf< And< Not< IsSymmetric<T1> >, IsSymmetric<T2> >
934  , const typename SubExprTrait<T1,T2>::Type >::Type
935  operator-( const DenseMatrix<T1,true>& lhs, const DenseMatrix<T2,false>& rhs )
936 {
938 
939  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
940  throw std::invalid_argument( "Matrix sizes do not match" );
941 
942  return (~lhs) - trans( ~rhs );
943 }
945 //*************************************************************************************************
946 
947 
948 //*************************************************************************************************
962 template< typename T1 // Type of the left-hand side dense matrix
963  , typename T2 > // Type of the right-hand side dense matrix
964 inline typename EnableIf< IsSymmetric<T1>
965  , const typename SubExprTrait<T1,T2>::Type >::Type
966  operator-( const DenseMatrix<T1,true>& lhs, const DenseMatrix<T2,false>& rhs )
967 {
969 
970  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
971  throw std::invalid_argument( "Matrix sizes do not match" );
972 
973  return trans( ~lhs ) - (~rhs);
974 }
976 //*************************************************************************************************
977 
978 
979 
980 
981 //=================================================================================================
982 //
983 // ROWS SPECIALIZATIONS
984 //
985 //=================================================================================================
986 
987 //*************************************************************************************************
989 template< typename MT1, typename MT2 >
990 struct Rows< DMatTDMatSubExpr<MT1,MT2> >
991  : public Max< Rows<MT1>, Rows<MT2> >::Type
992 {};
994 //*************************************************************************************************
995 
996 
997 
998 
999 //=================================================================================================
1000 //
1001 // COLUMNS SPECIALIZATIONS
1002 //
1003 //=================================================================================================
1004 
1005 //*************************************************************************************************
1007 template< typename MT1, typename MT2 >
1008 struct Columns< DMatTDMatSubExpr<MT1,MT2> >
1009  : public Max< Columns<MT1>, Columns<MT2> >::Type
1010 {};
1012 //*************************************************************************************************
1013 
1014 
1015 
1016 
1017 //=================================================================================================
1018 //
1019 // ISSYMMETRIC SPECIALIZATIONS
1020 //
1021 //=================================================================================================
1022 
1023 //*************************************************************************************************
1025 template< typename MT1, typename MT2 >
1026 struct IsSymmetric< DMatTDMatSubExpr<MT1,MT2> >
1027  : public IsTrue< IsSymmetric<MT1>::value && IsSymmetric<MT2>::value >
1028 {};
1030 //*************************************************************************************************
1031 
1032 
1033 
1034 
1035 //=================================================================================================
1036 //
1037 // ISLOWER SPECIALIZATIONS
1038 //
1039 //=================================================================================================
1040 
1041 //*************************************************************************************************
1043 template< typename MT1, typename MT2 >
1044 struct IsLower< DMatTDMatSubExpr<MT1,MT2> >
1045  : public IsTrue< IsLower<MT1>::value && IsLower<MT2>::value >
1046 {};
1048 //*************************************************************************************************
1049 
1050 
1051 
1052 
1053 //=================================================================================================
1054 //
1055 // ISUPPER SPECIALIZATIONS
1056 //
1057 //=================================================================================================
1058 
1059 //*************************************************************************************************
1061 template< typename MT1, typename MT2 >
1062 struct IsUpper< DMatTDMatSubExpr<MT1,MT2> >
1063  : public IsTrue< IsUpper<MT1>::value && IsUpper<MT2>::value >
1064 {};
1066 //*************************************************************************************************
1067 
1068 
1069 
1070 
1071 //=================================================================================================
1072 //
1073 // EXPRESSION TRAIT SPECIALIZATIONS
1074 //
1075 //=================================================================================================
1076 
1077 //*************************************************************************************************
1079 template< typename MT1, typename MT2, bool AF >
1080 struct SubmatrixExprTrait< DMatTDMatSubExpr<MT1,MT2>, AF >
1081 {
1082  public:
1083  //**********************************************************************************************
1084  typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
1085  , typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
1086  //**********************************************************************************************
1087 };
1089 //*************************************************************************************************
1090 
1091 
1092 //*************************************************************************************************
1094 template< typename MT1, typename MT2 >
1095 struct RowExprTrait< DMatTDMatSubExpr<MT1,MT2> >
1096 {
1097  public:
1098  //**********************************************************************************************
1099  typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
1100  , typename RowExprTrait<const MT2>::Type >::Type Type;
1101  //**********************************************************************************************
1102 };
1104 //*************************************************************************************************
1105 
1106 
1107 //*************************************************************************************************
1109 template< typename MT1, typename MT2 >
1110 struct ColumnExprTrait< DMatTDMatSubExpr<MT1,MT2> >
1111 {
1112  public:
1113  //**********************************************************************************************
1114  typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
1115  , typename ColumnExprTrait<const MT2>::Type >::Type Type;
1116  //**********************************************************************************************
1117 };
1119 //*************************************************************************************************
1120 
1121 } // namespace blaze
1122 
1123 #endif
Constraint on the data type.
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:103
Header file for the Max class template.
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTDMatSubExpr.h:260
Compile time check whether the given type is a temporary vector or matrix type.This type trait class ...
Definition: IsTemporary.h:87
Header file for the Rows type trait.
Header file for the subtraction trait.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:160
const size_t SMP_DMATTDMATSUB_THRESHOLD
SMP row-major dense matrix/column-major dense matrix subtraction threshold.This threshold specifies w...
Definition: Thresholds.h:788
Efficient implementation of a compressed matrix.The CompressedMatrix class template is the represent...
Definition: CompressedMatrix.h:205
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is not a dense, N-dimensional matrix type...
Definition: DenseMatrix.h:79
BLAZE_ALWAYS_INLINE bool isSame(const Matrix< MT1, SO1 > &a, const Matrix< MT2, SO2 > &b)
Returns whether the two given matrices represent the same observable state.
Definition: Matrix.h:946
Header file for the ColumnExprTrait class template.
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatTDMatSubExpr.h:118
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2478
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:257
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:301
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:695
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:104
Header file for the Computation base class.
Header file for the RequiresEvaluation type trait.
MT1::ReturnType RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:102
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:70
Constraint on the data type.
const SelectType< returnExpr, ExprReturnType, ElementType >::Type ReturnType
Return type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:164
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:101
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:90
Expression object for dense matrix-transpose dense matrix subtractions.The DMatTDMatSubExpr class rep...
Definition: DMatTDMatSubExpr.h:94
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTDMatSubExpr.h:167
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 DisableIf class template.
Header file for the IsTemporary type trait class.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:100
Header file for the IsSymmetric type trait.
Base class for all matrix/matrix subtraction expression templates.The MatMatSubExpr class serves as a...
Definition: MatMatSubExpr.h:65
#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: StorageOrder.h:161
bool isAliased(const T *alias) const
Returns whether the expression is aliased with the given address alias.
Definition: DMatTDMatSubExpr.h:273
Header file for the DenseMatrix base class.
BLAZE_ALWAYS_INLINE void assign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the assignment of a matrix to a matrix.
Definition: Matrix.h:635
Header file for the Columns type trait.
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTDMatSubExpr.h:161
Header file for the MatMatSubExpr base class.
DMatTDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTDMatSubExpr class.
Definition: DMatTDMatSubExpr.h:190
Header file for the IsLower type trait.
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:302
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:158
#define BLAZE_CONSTRAINT_MUST_BE_REFERENCE_TYPE(T)
Constraint on the data type.In case the given data type T is not a reference type, a compilation error is created.
Definition: Reference.h:78
Constraints on the storage order of matrix types.
Evaluation of the return type of a subtraction expression.Via this type trait it is possible to evalu...
Definition: SubExprTrait.h:104
Header file for the SelectType class template.
Header file for the RowExprTrait class template.
Header file for all forward declarations for expression class templates.
Header file for the EnableIf class template.
Header file for the serial shim.
#define BLAZE_CONSTRAINT_MATRICES_MUST_HAVE_DIFFERENT_STORAGE_ORDER(T1, T2)
Constraint on the data type.In case either of the two given data types T1 or T2 is not a matrix type ...
Definition: StorageOrder.h:325
bool isAligned() const
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatTDMatSubExpr.h:283
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:105
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.
#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: StorageOrder.h:81
const Type & ReturnType
Return type for expression template evaluations.
Definition: CompressedMatrix.h:2477
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
BLAZE_ALWAYS_INLINE void addAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the addition assignment of a matrix to a matrix.
Definition: Matrix.h:742
DMatTDMatSubExpr< MT1, MT2 > This
Type of this DMatTDMatSubExpr instance.
Definition: DMatTDMatSubExpr.h:157
const DenseIterator< Type > operator-(const DenseIterator< Type > &it, ptrdiff_t inc)
Subtraction between a DenseIterator and an integral value.
Definition: DenseIterator.h:585
Substitution Failure Is Not An Error (SFINAE) class.The EnableIf class template is an auxiliary tool ...
Definition: EnableIf.h:184
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTDMatSubExpr.h:206
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTDMatSubExpr.h:218
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:170
Substitution Failure Is Not An Error (SFINAE) class.The DisableIf class template is an auxiliary tool...
Definition: DisableIf.h:184
#define BLAZE_CONSTRAINT_MATRICES_MUST_HAVE_SAME_STORAGE_ORDER(T1, T2)
Constraint on the data type.In case either of the two given data types T1 or T2 is not a matrix type ...
Definition: StorageOrder.h:283
SelectType< IsExpression< MT2 >::value, const MT2, const MT2 & >::Type RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:173
const DMatTransExpr< MT,!SO > trans(const DenseMatrix< MT, SO > &dm)
Calculation of the transpose of the given dense matrix.
Definition: DMatTransExpr.h:932
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
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
LeftOperand leftOperand() const
Returns the left-hand side dense matrix operand.
Definition: DMatTDMatSubExpr.h:238
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2473
Header file for the IsTrue value trait.
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: DMatTDMatSubExpr.h:248
Header file for basic type definitions.
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_MATMATSUBEXPR(T1, T2)
Constraint on the data type.In case the given data types T1 and T2 do not form a valid matrix/matrix ...
Definition: MatMatSubExpr.h:165
Base template for the SubTrait class.
Definition: SubTrait.h:142
Header file for the IsUpper type trait.
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTDMatSubExpr.h:159
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTDMatSubExpr.h:293
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTDMatSubExpr.h:228
Header file for the SubExprTrait class template.
Header file for the thresholds for matrix/vector and matrix/matrix multiplications.
#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
Compile time check whether the given type is an expression template.This type trait class tests wheth...
Definition: IsExpression.h:87
Header file for the IsExpression type trait class.
Header file for the FunctionTrace class.
BLAZE_ALWAYS_INLINE void subAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the subtraction assignment of a matrix to matrix.
Definition: Matrix.h:849