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>
60 #include <blaze/util/Assert.h>
63 #include <blaze/util/SelectType.h>
64 #include <blaze/util/Types.h>
65 
66 
67 namespace blaze {
68 
69 //=================================================================================================
70 //
71 // CLASS DMATTDMATSUBEXPR
72 //
73 //=================================================================================================
74 
75 //*************************************************************************************************
82 template< typename MT1 // Type of the left-hand side dense matrix
83  , typename MT2 > // Type of the right-hand side dense matrix
84 class DMatTDMatSubExpr : public DenseMatrix< DMatTDMatSubExpr<MT1,MT2>, false >
85  , private MatMatSubExpr
86  , private Computation
87 {
88  private:
89  //**Type definitions****************************************************************************
90  typedef typename MT1::ResultType RT1;
91  typedef typename MT2::ResultType RT2;
92  typedef typename MT1::ReturnType RN1;
93  typedef typename MT2::ReturnType RN2;
94  typedef typename MT1::CompositeType CT1;
95  typedef typename MT2::CompositeType CT2;
96  //**********************************************************************************************
97 
98  //**Return type evaluation**********************************************************************
100 
105  enum { returnExpr = !IsTemporary<RN1>::value && !IsTemporary<RN2>::value };
106 
109  //**********************************************************************************************
110 
111  //**Serial evaluation strategy******************************************************************
113 
119  enum { useAssign = RequiresEvaluation<MT1>::value || RequiresEvaluation<MT2>::value || !returnExpr };
120 
122  template< typename MT >
124  struct UseAssign {
125  enum { value = useAssign };
126  };
128  //**********************************************************************************************
129 
130  //**Parallel evaluation strategy****************************************************************
132 
138  template< typename MT >
139  struct UseSMPAssign {
140  enum { value = ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign };
141  };
143  //**********************************************************************************************
144 
145  public:
146  //**Type definitions****************************************************************************
152 
155 
157  typedef const ResultType CompositeType;
158 
160  typedef typename SelectType< IsExpression<MT1>::value, const MT1, const MT1& >::Type LeftOperand;
161 
163  typedef typename SelectType< IsExpression<MT2>::value, const MT2, const MT2& >::Type RightOperand;
164  //**********************************************************************************************
165 
166  //**Compilation flags***************************************************************************
168  enum { vectorizable = 0 };
169 
171  enum { smpAssignable = MT1::smpAssignable && MT2::smpAssignable };
172  //**********************************************************************************************
173 
174  //**Constructor*********************************************************************************
180  explicit inline DMatTDMatSubExpr( const MT1& lhs, const MT2& rhs )
181  : lhs_( lhs ) // Left-hand side dense matrix of the subtraction expression
182  , rhs_( rhs ) // Right-hand side dense matrix of the subtraction expression
183  {
184  BLAZE_INTERNAL_ASSERT( lhs.rows() == rhs.rows() , "Invalid number of rows" );
185  BLAZE_INTERNAL_ASSERT( lhs.columns() == rhs.columns(), "Invalid number of columns" );
186  }
187  //**********************************************************************************************
188 
189  //**Access operator*****************************************************************************
196  inline ReturnType operator()( size_t i, size_t j ) const {
197  BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
198  BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
199  return lhs_(i,j) - rhs_(i,j);
200  }
201  //**********************************************************************************************
202 
203  //**Rows function*******************************************************************************
208  inline size_t rows() const {
209  return lhs_.rows();
210  }
211  //**********************************************************************************************
212 
213  //**Columns function****************************************************************************
218  inline size_t columns() const {
219  return lhs_.columns();
220  }
221  //**********************************************************************************************
222 
223  //**Left operand access*************************************************************************
228  inline LeftOperand leftOperand() const {
229  return lhs_;
230  }
231  //**********************************************************************************************
232 
233  //**Right operand access************************************************************************
238  inline RightOperand rightOperand() const {
239  return rhs_;
240  }
241  //**********************************************************************************************
242 
243  //**********************************************************************************************
249  template< typename T >
250  inline bool canAlias( const T* alias ) const {
251  return ( IsExpression<MT1>::value && ( RequiresEvaluation<MT1>::value ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
252  ( IsExpression<MT2>::value && ( RequiresEvaluation<MT2>::value ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
253  }
254  //**********************************************************************************************
255 
256  //**********************************************************************************************
262  template< typename T >
263  inline bool isAliased( const T* alias ) const {
264  return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
265  }
266  //**********************************************************************************************
267 
268  //**********************************************************************************************
273  inline bool isAligned() const {
274  return lhs_.isAligned() && rhs_.isAligned();
275  }
276  //**********************************************************************************************
277 
278  //**********************************************************************************************
283  inline bool canSMPAssign() const {
284  return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
286  }
287  //**********************************************************************************************
288 
289  private:
290  //**Member variables****************************************************************************
293  //**********************************************************************************************
294 
295  //**Assignment to dense matrices****************************************************************
309  template< typename MT // Type of the target dense matrix
310  , bool SO2 > // Storage order of the target dense matrix
311  friend inline typename DisableIf< UseAssign<MT> >::Type
312  assign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
313  {
315 
316  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
317  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
318 
319  const size_t m( rhs.rows() );
320  const size_t n( rhs.columns() );
321  const size_t block( 16UL );
322 
323  for( size_t ii=0UL; ii<m; ii+=block ) {
324  const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
325  for( size_t jj=0UL; jj<n; jj+=block ) {
326  const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
327  for( size_t i=ii; i<iend; ++i ) {
328  for( size_t j=jj; j<jend; ++j ) {
329  (~lhs)(i,j) = rhs.lhs_(i,j) - rhs.rhs_(i,j);
330  }
331  }
332  }
333  }
334  }
336  //**********************************************************************************************
337 
338  //**Assignment to dense matrices****************************************************************
352  template< typename MT // Type of the target dense matrix
353  , bool SO2 > // Storage order of the target dense matrix
354  friend inline typename EnableIf< UseAssign<MT> >::Type
355  assign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
356  {
358 
359  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
360  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
361 
362  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
363  subAssign( ~lhs, rhs.rhs_ );
364  }
365  else {
366  assign ( ~lhs, rhs.lhs_ );
367  subAssign( ~lhs, rhs.rhs_ );
368  }
369  }
371  //**********************************************************************************************
372 
373  //**Assignment to sparse matrices***************************************************************
385  template< typename MT // Type of the target sparse matrix
386  , bool SO2 > // Storage order of the target sparse matrix
387  friend inline void assign( SparseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
388  {
390 
391  typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
392 
399 
400  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
401  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
402 
403  const TmpType tmp( serial( rhs ) );
404  assign( ~lhs, tmp );
405  }
407  //**********************************************************************************************
408 
409  //**Addition assignment to dense matrices*******************************************************
424  template< typename MT // Type of the target dense matrix
425  , bool SO2 > // Storage order of the target dense matrix
426  friend inline typename DisableIf< UseAssign<MT> >::Type
427  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
428  {
430 
431  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
432  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
433 
434  const size_t m( rhs.rows() );
435  const size_t n( rhs.columns() );
436  const size_t block( 16UL );
437 
438  for( size_t ii=0UL; ii<m; ii+=block ) {
439  const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
440  for( size_t jj=0UL; jj<n; jj+=block ) {
441  const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
442  for( size_t i=ii; i<iend; ++i ) {
443  for( size_t j=jj; j<jend; ++j ) {
444  (~lhs)(i,j) += rhs.lhs_(i,j) - rhs.rhs_(i,j);
445  }
446  }
447  }
448  }
449  }
451  //**********************************************************************************************
452 
453  //**Addition assignment to dense matrices*******************************************************
468  template< typename MT // Type of the target dense matrix
469  , bool SO2 > // Storage order of the target dense matrix
470  friend inline typename EnableIf< UseAssign<MT> >::Type
471  addAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
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  addAssign( ~lhs, rhs.lhs_ );
479  subAssign( ~lhs, rhs.rhs_ );
480  }
482  //**********************************************************************************************
483 
484  //**Addition assignment to sparse matrices******************************************************
485  // No special implementation for the addition assignment to sparse matrices.
486  //**********************************************************************************************
487 
488  //**Subtraction assignment to dense matrices****************************************************
503  template< typename MT // Type of the target dense matrix
504  , bool SO2 > // Storage order of the target dense matrix
505  friend inline typename DisableIf< UseAssign<MT> >::Type
506  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
507  {
509 
510  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
511  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
512 
513  const size_t m( rhs.rows() );
514  const size_t n( rhs.columns() );
515  const size_t block( 16UL );
516 
517  for( size_t ii=0UL; ii<m; ii+=block ) {
518  const size_t iend( ( m < ii+block )?( m ):( ii+block ) );
519  for( size_t jj=0UL; jj<n; jj+=block ) {
520  const size_t jend( ( n < jj+block )?( n ):( jj+block ) );
521  for( size_t i=ii; i<iend; ++i ) {
522  for( size_t j=jj; j<jend; ++j ) {
523  (~lhs)(i,j) -= rhs.lhs_(i,j) - rhs.rhs_(i,j);
524  }
525  }
526  }
527  }
528  }
530  //**********************************************************************************************
531 
532  //**Subtraction assignment to dense matrices****************************************************
547  template< typename MT // Type of the target dense matrix
548  , bool SO2 > // Storage order of the target dense matrix
549  friend inline typename EnableIf< UseAssign<MT> >::Type
550  subAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
551  {
553 
554  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
555  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
556 
557  subAssign( ~lhs, rhs.lhs_ );
558  addAssign( ~lhs, rhs.rhs_ );
559  }
561  //**********************************************************************************************
562 
563  //**Subtraction assignment to sparse matrices***************************************************
564  // No special implementation for the subtraction assignment to sparse matrices.
565  //**********************************************************************************************
566 
567  //**Multiplication assignment to dense matrices*************************************************
568  // No special implementation for the multiplication assignment to dense matrices.
569  //**********************************************************************************************
570 
571  //**Multiplication assignment to sparse matrices************************************************
572  // No special implementation for the multiplication assignment to sparse matrices.
573  //**********************************************************************************************
574 
575  //**SMP assignment to dense matrices************************************************************
589  template< typename MT // Type of the target dense matrix
590  , bool SO2 > // Storage order of the target dense matrix
591  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
592  smpAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
593  {
595 
596  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
597  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
598 
599  if( !IsExpression<MT1>::value && isSame( ~lhs, rhs.lhs_ ) ) {
600  smpSubAssign( ~lhs, rhs.rhs_ );
601  }
602  else {
603  smpAssign ( ~lhs, rhs.lhs_ );
604  smpSubAssign( ~lhs, rhs.rhs_ );
605  }
606  }
608  //**********************************************************************************************
609 
610  //**SMP assignment to sparse matrices***********************************************************
624  template< typename MT // Type of the target sparse matrix
625  , bool SO2 > // Storage order of the target sparse matrix
626  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
627  smpAssign( SparseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
628  {
630 
631  typedef typename SelectType< SO2, OppositeType, ResultType >::Type TmpType;
632 
639 
640  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
641  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
642 
643  const TmpType tmp( rhs );
644  smpAssign( ~lhs, tmp );
645  }
647  //**********************************************************************************************
648 
649  //**SMP addition assignment to dense matrices***************************************************
664  template< typename MT // Type of the target dense matrix
665  , bool SO2 > // Storage order of the target dense matrix
666  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
667  smpAddAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
668  {
670 
671  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
672  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
673 
674  smpAddAssign( ~lhs, rhs.lhs_ );
675  smpSubAssign( ~lhs, rhs.rhs_ );
676  }
678  //**********************************************************************************************
679 
680  //**SMP addition assignment to sparse matrices**************************************************
681  // No special implementation for the SMP addition assignment to sparse matrices.
682  //**********************************************************************************************
683 
684  //**SMP subtraction assignment to dense matrices************************************************
699  template< typename MT // Type of the target dense matrix
700  , bool SO2 > // Storage order of the target dense matrix
701  friend inline typename EnableIf< UseSMPAssign<MT> >::Type
702  smpSubAssign( DenseMatrix<MT,SO2>& lhs, const DMatTDMatSubExpr& rhs )
703  {
705 
706  BLAZE_INTERNAL_ASSERT( (~lhs).rows() == rhs.rows() , "Invalid number of rows" );
707  BLAZE_INTERNAL_ASSERT( (~lhs).columns() == rhs.columns(), "Invalid number of columns" );
708 
709  smpSubAssign( ~lhs, rhs.lhs_ );
710  smpAddAssign( ~lhs, rhs.rhs_ );
711  }
713  //**********************************************************************************************
714 
715  //**SMP subtraction assignment to sparse matrices***********************************************
716  // No special implementation for the SMP subtraction assignment to sparse matrices.
717  //**********************************************************************************************
718 
719  //**SMP multiplication assignment to dense matrices*********************************************
720  // No special implementation for the SMP multiplication assignment to dense matrices.
721  //**********************************************************************************************
722 
723  //**SMP multiplication assignment to sparse matrices********************************************
724  // No special implementation for the SMP multiplication assignment to sparse matrices.
725  //**********************************************************************************************
726 
727  //**Compile time checks*************************************************************************
733  //**********************************************************************************************
734 };
735 //*************************************************************************************************
736 
737 
738 
739 
740 //=================================================================================================
741 //
742 // GLOBAL BINARY ARITHMETIC OPERATORS
743 //
744 //=================================================================================================
745 
746 //*************************************************************************************************
775 template< typename T1 // Type of the left-hand side dense matrix
776  , typename T2 > // Type of the right-hand side dense matrix
777 inline const DMatTDMatSubExpr<T1,T2>
779 {
781 
782  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
783  throw std::invalid_argument( "Matrix sizes do not match" );
784 
785  return DMatTDMatSubExpr<T1,T2>( ~lhs, ~rhs );
786 }
787 //*************************************************************************************************
788 
789 
790 //*************************************************************************************************
819 template< typename T1 // Type of the left-hand side dense matrix
820  , typename T2 > // Type of the right-hand side dense matrix
821 inline const DMatTDMatSubExpr<T1,T2>
823 {
825 
826  if( (~lhs).rows() != (~rhs).rows() || (~lhs).columns() != (~rhs).columns() )
827  throw std::invalid_argument( "Matrix sizes do not match" );
828 
829  return DMatTDMatSubExpr<T1,T2>( ~lhs, ~rhs );
830 }
831 //*************************************************************************************************
832 
833 
834 
835 
836 //=================================================================================================
837 //
838 // EXPRESSION TRAIT SPECIALIZATIONS
839 //
840 //=================================================================================================
841 
842 //*************************************************************************************************
844 template< typename MT1, typename MT2, bool AF >
845 struct SubmatrixExprTrait< DMatTDMatSubExpr<MT1,MT2>, AF >
846 {
847  public:
848  //**********************************************************************************************
849  typedef typename SubExprTrait< typename SubmatrixExprTrait<const MT1,AF>::Type
850  , typename SubmatrixExprTrait<const MT2,AF>::Type >::Type Type;
851  //**********************************************************************************************
852 };
854 //*************************************************************************************************
855 
856 
857 //*************************************************************************************************
859 template< typename MT1, typename MT2 >
860 struct RowExprTrait< DMatTDMatSubExpr<MT1,MT2> >
861 {
862  public:
863  //**********************************************************************************************
864  typedef typename SubExprTrait< typename RowExprTrait<const MT1>::Type
865  , typename RowExprTrait<const MT2>::Type >::Type Type;
866  //**********************************************************************************************
867 };
869 //*************************************************************************************************
870 
871 
872 //*************************************************************************************************
874 template< typename MT1, typename MT2 >
875 struct ColumnExprTrait< DMatTDMatSubExpr<MT1,MT2> >
876 {
877  public:
878  //**********************************************************************************************
879  typedef typename SubExprTrait< typename ColumnExprTrait<const MT1>::Type
880  , typename ColumnExprTrait<const MT2>::Type >::Type Type;
881  //**********************************************************************************************
882 };
884 //*************************************************************************************************
885 
886 } // namespace blaze
887 
888 #endif
MT2::ReturnType RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:93
bool canAlias(const T *alias) const
Returns whether the expression can alias with the given address alias.
Definition: DMatTDMatSubExpr.h:250
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 subtraction trait.
ResultType::TransposeType TransposeType
Transpose type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:150
void smpSubAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:152
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:199
#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
Header file for the ColumnExprTrait class template.
SubExprTrait< RN1, RN2 >::Type ExprReturnType
Expression return type for the subscript operator.
Definition: DMatTDMatSubExpr.h:108
const This & CompositeType
Data type for composite expression templates.
Definition: CompressedMatrix.h:2408
Type ElementType
Type of the sparse matrix elements.
Definition: CompressedMatrix.h:251
LeftOperand lhs_
Left-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:291
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:690
MT1::CompositeType CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:94
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:92
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:409
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:154
MT2::ResultType RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:91
Constraint on the data type.
void smpAddAssign(DenseMatrix< 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:122
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:84
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatTDMatSubExpr.h:157
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 IsTemporary type trait class.
MT1::ResultType RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:90
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:263
Header file for the DenseMatrix base class.
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:271
ResultType::ElementType ElementType
Resulting element type.
Definition: DMatTDMatSubExpr.h:151
Header file for the MatMatSubExpr base class.
DMatTDMatSubExpr(const MT1 &lhs, const MT2 &rhs)
Constructor for the DMatTDMatSubExpr class.
Definition: DMatTDMatSubExpr.h:180
RightOperand rhs_
Right-hand side dense matrix of the subtraction expression.
Definition: DMatTDMatSubExpr.h:292
SubTrait< RT1, RT2 >::Type ResultType
Result type for expression template evaluations.
Definition: DMatTDMatSubExpr.h:148
#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:103
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 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:273
MT2::CompositeType CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:95
void smpAssign(DenseMatrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs)
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:92
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:2407
Header file for run time assertion macros.
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:301
DMatTDMatSubExpr< MT1, MT2 > This
Type of this DMatTDMatSubExpr instance.
Definition: DMatTDMatSubExpr.h:147
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
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:331
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatTDMatSubExpr.h:196
size_t rows() const
Returns the current number of rows of the matrix.
Definition: DMatTDMatSubExpr.h:208
SelectType< IsExpression< MT1 >::value, const MT1, const MT1 & >::Type LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatTDMatSubExpr.h:160
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:163
Base class for all compute expression templates.The Computation class serves as a tag for all computa...
Definition: Computation.h:59
#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:228
This ResultType
Result type for expression template evaluations.
Definition: CompressedMatrix.h:2403
RightOperand rightOperand() const
Returns the right-hand side transpose dense matrix operand.
Definition: DMatTDMatSubExpr.h:238
Header file for basic type definitions.
Base template for the SubTrait class.
Definition: SubTrait.h:141
ResultType::OppositeType OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatTDMatSubExpr.h:149
bool canSMPAssign() const
Returns whether the expression can be used in SMP assignments.
Definition: DMatTDMatSubExpr.h:283
size_t columns() const
Returns the current number of columns of the matrix.
Definition: DMatTDMatSubExpr.h:218
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.