Blaze 3.9
DMatDMatSchurExpr.h
Go to the documentation of this file.
1//=================================================================================================
33//=================================================================================================
34
35#ifndef _BLAZE_MATH_EXPRESSIONS_DMATDMATSCHUREXPR_H_
36#define _BLAZE_MATH_EXPRESSIONS_DMATDMATSCHUREXPR_H_
37
38
39//*************************************************************************************************
40// Includes
41//*************************************************************************************************
42
43#include <iterator>
44#include <utility>
45#include <blaze/math/Aliases.h>
58#include <blaze/math/SIMD.h>
76#include <blaze/system/Inline.h>
78#include <blaze/util/Assert.h>
79#include <blaze/util/EnableIf.h>
83#include <blaze/util/mpl/If.h>
84#include <blaze/util/Types.h>
85
86
87namespace blaze {
88
89//=================================================================================================
90//
91// CLASS DMATDMATSCHUREXPR
92//
93//=================================================================================================
94
95//*************************************************************************************************
102template< typename MT1 // Type of the left-hand side dense matrix
103 , typename MT2 // Type of the right-hand side dense matrix
104 , bool SO > // Storage order
106 : public SchurExpr< DenseMatrix< DMatDMatSchurExpr<MT1,MT2,SO>, SO > >
107 , private Computation
108{
109 private:
110 //**Type definitions****************************************************************************
119 //**********************************************************************************************
120
121 //**Return type evaluation**********************************************************************
123
128 static constexpr bool returnExpr = ( !IsTemporary_v<RN1> && !IsTemporary_v<RN2> );
129
131 using ExprReturnType = decltype( std::declval<RN1>() * std::declval<RN2>() );
132 //**********************************************************************************************
133
134 //**Serial evaluation strategy******************************************************************
136
142 static constexpr bool useAssign =
143 ( RequiresEvaluation_v<MT1> || RequiresEvaluation_v<MT2> || !returnExpr );
144
147 template< typename MT >
148 static constexpr bool UseAssign_v = useAssign;
150 //**********************************************************************************************
151
152 //**Parallel evaluation strategy****************************************************************
154
160 template< typename MT >
161 static constexpr bool UseSMPAssign_v =
162 ( ( !MT1::smpAssignable || !MT2::smpAssignable ) && useAssign );
164 //**********************************************************************************************
165
166 public:
167 //**Type definitions****************************************************************************
170
173
178
181
184
186 using LeftOperand = If_t< IsExpression_v<MT1>, const MT1, const MT1& >;
187
189 using RightOperand = If_t< IsExpression_v<MT2>, const MT2, const MT2& >;
190 //**********************************************************************************************
191
192 //**ConstIterator class definition**************************************************************
196 {
197 public:
198 //**Type definitions*************************************************************************
199 using IteratorCategory = std::random_access_iterator_tag;
203 using DifferenceType = ptrdiff_t;
204
205 // STL iterator requirements
211
214
217 //*******************************************************************************************
218
219 //**Constructor******************************************************************************
226 : left_ ( left ) // Iterator to the current left-hand side element
227 , right_( right ) // Iterator to the current right-hand side element
228 {}
229 //*******************************************************************************************
230
231 //**Addition assignment operator*************************************************************
238 left_ += inc;
239 right_ += inc;
240 return *this;
241 }
242 //*******************************************************************************************
243
244 //**Subtraction assignment operator**********************************************************
251 left_ -= dec;
252 right_ -= dec;
253 return *this;
254 }
255 //*******************************************************************************************
256
257 //**Prefix increment operator****************************************************************
263 ++left_;
264 ++right_;
265 return *this;
266 }
267 //*******************************************************************************************
268
269 //**Postfix increment operator***************************************************************
275 return ConstIterator( left_++, right_++ );
276 }
277 //*******************************************************************************************
278
279 //**Prefix decrement operator****************************************************************
285 --left_;
286 --right_;
287 return *this;
288 }
289 //*******************************************************************************************
290
291 //**Postfix decrement operator***************************************************************
297 return ConstIterator( left_--, right_-- );
298 }
299 //*******************************************************************************************
300
301 //**Element access operator******************************************************************
306 inline ReturnType operator*() const {
307 return (*left_) * (*right_);
308 }
309 //*******************************************************************************************
310
311 //**Load function****************************************************************************
316 inline auto load() const noexcept {
317 return left_.load() * right_.load();
318 }
319 //*******************************************************************************************
320
321 //**Equality operator************************************************************************
327 inline BLAZE_DEVICE_CALLABLE bool operator==( const ConstIterator& rhs ) const {
328 return left_ == rhs.left_;
329 }
330 //*******************************************************************************************
331
332 //**Inequality operator**********************************************************************
338 inline BLAZE_DEVICE_CALLABLE bool operator!=( const ConstIterator& rhs ) const {
339 return left_ != rhs.left_;
340 }
341 //*******************************************************************************************
342
343 //**Less-than operator***********************************************************************
349 inline BLAZE_DEVICE_CALLABLE bool operator<( const ConstIterator& rhs ) const {
350 return left_ < rhs.left_;
351 }
352 //*******************************************************************************************
353
354 //**Greater-than operator********************************************************************
360 inline BLAZE_DEVICE_CALLABLE bool operator>( const ConstIterator& rhs ) const {
361 return left_ > rhs.left_;
362 }
363 //*******************************************************************************************
364
365 //**Less-or-equal-than operator**************************************************************
371 inline BLAZE_DEVICE_CALLABLE bool operator<=( const ConstIterator& rhs ) const {
372 return left_ <= rhs.left_;
373 }
374 //*******************************************************************************************
375
376 //**Greater-or-equal-than operator***********************************************************
382 inline BLAZE_DEVICE_CALLABLE bool operator>=( const ConstIterator& rhs ) const {
383 return left_ >= rhs.left_;
384 }
385 //*******************************************************************************************
386
387 //**Subtraction operator*********************************************************************
394 return left_ - rhs.left_;
395 }
396 //*******************************************************************************************
397
398 //**Addition operator************************************************************************
405 friend inline const ConstIterator operator+( const ConstIterator& it, size_t inc ) {
406 return ConstIterator( it.left_ + inc, it.right_ + inc );
407 }
408 //*******************************************************************************************
409
410 //**Addition operator************************************************************************
417 friend inline const ConstIterator operator+( size_t inc, const ConstIterator& it ) {
418 return ConstIterator( it.left_ + inc, it.right_ + inc );
419 }
420 //*******************************************************************************************
421
422 //**Subtraction operator*********************************************************************
429 friend inline BLAZE_DEVICE_CALLABLE const ConstIterator operator-( const ConstIterator& it, size_t dec ) {
430 return ConstIterator( it.left_ - dec, it.right_ - dec );
431 }
432 //*******************************************************************************************
433
434 private:
435 //**Member variables*************************************************************************
438 //*******************************************************************************************
439 };
440 //**********************************************************************************************
441
442 //**Compilation flags***************************************************************************
444 static constexpr bool simdEnabled =
445 ( MT1::simdEnabled && MT2::simdEnabled && HasSIMDMult_v<ET1,ET2> );
446
448 static constexpr bool smpAssignable = ( MT1::smpAssignable && MT2::smpAssignable );
449 //**********************************************************************************************
450
451 //**SIMD properties*****************************************************************************
453 static constexpr size_t SIMDSIZE = SIMDTrait<ElementType>::size;
454 //**********************************************************************************************
455
456 //**Constructor*********************************************************************************
462 inline DMatDMatSchurExpr( const MT1& lhs, const MT2& rhs ) noexcept
463 : lhs_( lhs ) // Left-hand side dense matrix of the Schur product expression
464 , rhs_( rhs ) // Right-hand side dense matrix of the Schur product expression
465 {
466 BLAZE_INTERNAL_ASSERT( lhs.rows() == rhs.rows() , "Invalid number of rows" );
467 BLAZE_INTERNAL_ASSERT( lhs.columns() == rhs.columns(), "Invalid number of columns" );
468 }
469 //**********************************************************************************************
470
471 //**Access operator*****************************************************************************
478 inline ReturnType operator()( size_t i, size_t j ) const {
479 BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
480 BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
481 return lhs_(i,j) * rhs_(i,j);
482 }
483 //**********************************************************************************************
484
485 //**At function*********************************************************************************
493 inline ReturnType at( size_t i, size_t j ) const {
494 if( i >= lhs_.rows() ) {
495 BLAZE_THROW_OUT_OF_RANGE( "Invalid row access index" );
496 }
497 if( j >= lhs_.columns() ) {
498 BLAZE_THROW_OUT_OF_RANGE( "Invalid column access index" );
499 }
500 return (*this)(i,j);
501 }
502 //**********************************************************************************************
503
504 //**Load function*******************************************************************************
511 BLAZE_ALWAYS_INLINE auto load( size_t i, size_t j ) const noexcept {
512 BLAZE_INTERNAL_ASSERT( i < lhs_.rows() , "Invalid row access index" );
513 BLAZE_INTERNAL_ASSERT( j < lhs_.columns(), "Invalid column access index" );
514 BLAZE_INTERNAL_ASSERT( !SO || ( i % SIMDSIZE == 0UL ), "Invalid row access index" );
515 BLAZE_INTERNAL_ASSERT( SO || ( j % SIMDSIZE == 0UL ), "Invalid column access index" );
516 return lhs_.load(i,j) * rhs_.load(i,j);
517 }
518 //**********************************************************************************************
519
520 //**Begin function******************************************************************************
526 inline ConstIterator begin( size_t i ) const {
527 return ConstIterator( lhs_.begin(i), rhs_.begin(i) );
528 }
529 //**********************************************************************************************
530
531 //**End function********************************************************************************
537 inline ConstIterator end( size_t i ) const {
538 return ConstIterator( lhs_.end(i), rhs_.end(i) );
539 }
540 //**********************************************************************************************
541
542 //**Rows function*******************************************************************************
547 inline size_t rows() const noexcept {
548 return lhs_.rows();
549 }
550 //**********************************************************************************************
551
552 //**Columns function****************************************************************************
557 inline size_t columns() const noexcept {
558 return lhs_.columns();
559 }
560 //**********************************************************************************************
561
562 //**Left operand access*************************************************************************
567 inline LeftOperand leftOperand() const noexcept {
568 return lhs_;
569 }
570 //**********************************************************************************************
571
572 //**Right operand access************************************************************************
577 inline RightOperand rightOperand() const noexcept {
578 return rhs_;
579 }
580 //**********************************************************************************************
581
582 //**********************************************************************************************
588 template< typename T >
589 inline bool canAlias( const T* alias ) const noexcept {
590 return ( IsExpression_v<MT1> && ( RequiresEvaluation_v<MT1> ? lhs_.isAliased( alias ) : lhs_.canAlias( alias ) ) ) ||
591 ( IsExpression_v<MT2> && ( RequiresEvaluation_v<MT2> ? rhs_.isAliased( alias ) : rhs_.canAlias( alias ) ) );
592 }
593 //**********************************************************************************************
594
595 //**********************************************************************************************
601 template< typename T >
602 inline bool isAliased( const T* alias ) const noexcept {
603 return ( lhs_.isAliased( alias ) || rhs_.isAliased( alias ) );
604 }
605 //**********************************************************************************************
606
607 //**********************************************************************************************
612 inline bool isAligned() const noexcept {
613 return lhs_.isAligned() && rhs_.isAligned();
614 }
615 //**********************************************************************************************
616
617 //**********************************************************************************************
622 inline bool canSMPAssign() const noexcept {
623 return lhs_.canSMPAssign() || rhs_.canSMPAssign() ||
624 ( rows() * columns() >= SMP_DMATDMATSCHUR_THRESHOLD );
625 }
626 //**********************************************************************************************
627
628 private:
629 //**Member variables****************************************************************************
632 //**********************************************************************************************
633
634 //**Assignment to dense matrices****************************************************************
649 template< typename MT // Type of the target dense matrix
650 , bool SO2 > // Storage order of the target dense matrix
651 friend inline auto assign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
652 -> EnableIf_t< UseAssign_v<MT> && !IsCommutative_v<MT1,MT2> >
653 {
655
656 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
657 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
658
659 if( !IsOperation_v<MT1> && isSame( *lhs, rhs.lhs_ ) ) {
660 schurAssign( *lhs, rhs.rhs_ );
661 }
662 else {
663 CT1 A( serial( rhs.lhs_ ) );
664 CT2 B( serial( rhs.rhs_ ) );
665 assign( *lhs, A % B );
666 }
667 }
669 //**********************************************************************************************
670
671 //**Assignment to dense matrices****************************************************************
685 template< typename MT // Type of the target dense matrix
686 , bool SO2 > // Storage order of the target dense matrix
687 friend inline auto assign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
688 -> EnableIf_t< UseAssign_v<MT> && IsCommutative_v<MT1,MT2> >
689 {
691
692 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
693 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
694
695 if( !IsOperation_v<MT1> && isSame( *lhs, rhs.lhs_ ) ) {
696 schurAssign( *lhs, rhs.rhs_ );
697 }
698 else if( !IsOperation_v<MT2> && isSame( *lhs, rhs.rhs_ ) ) {
699 schurAssign( *lhs, rhs.lhs_ );
700 }
701 else if( !RequiresEvaluation_v<MT2> ) {
702 assign ( *lhs, rhs.rhs_ );
703 schurAssign( *lhs, rhs.lhs_ );
704 }
705 else {
706 assign ( *lhs, rhs.lhs_ );
707 schurAssign( *lhs, rhs.rhs_ );
708 }
709 }
711 //**********************************************************************************************
712
713 //**Assignment to sparse matrices***************************************************************
727 template< typename MT // Type of the target sparse matrix
728 , bool SO2 > // Storage order of the target sparse matrix
729 friend inline auto assign( SparseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
730 -> EnableIf_t< UseAssign_v<MT> >
731 {
733
734 using TmpType = If_t< SO == SO2, ResultType, OppositeType >;
735
742
743 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
744 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
745
746 const TmpType tmp( serial( rhs ) );
747 assign( *lhs, tmp );
748 }
750 //**********************************************************************************************
751
752 //**Addition assignment to dense matrices*******************************************************
766 template< typename MT // Type of the target dense matrix
767 , bool SO2 > // Storage order of the target dense matrix
768 friend inline auto addAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
769 -> EnableIf_t< UseAssign_v<MT> >
770 {
772
776
777 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
778 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
779
780 const ResultType tmp( serial( rhs ) );
781 addAssign( *lhs, tmp );
782 }
784 //**********************************************************************************************
785
786 //**Addition assignment to sparse matrices******************************************************
787 // No special implementation for the addition assignment to sparse matrices.
788 //**********************************************************************************************
789
790 //**Subtraction assignment to dense matrices****************************************************
804 template< typename MT // Type of the target dense matrix
805 , bool SO2 > // Storage order of the target dense matrix
806 friend inline auto subAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
807 -> EnableIf_t< UseAssign_v<MT> >
808 {
810
814
815 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
816 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
817
818 const ResultType tmp( serial( rhs ) );
819 subAssign( *lhs, tmp );
820 }
822 //**********************************************************************************************
823
824 //**Subtraction assignment to sparse matrices***************************************************
825 // No special implementation for the subtraction assignment to sparse matrices.
826 //**********************************************************************************************
827
828 //**Schur product assignment to dense matrices**************************************************
843 template< typename MT // Type of the target dense matrix
844 , bool SO2 > // Storage order of the target dense matrix
845 friend inline auto schurAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
846 -> EnableIf_t< UseAssign_v<MT> && !IsCommutative_v<MT1,MT2> >
847 {
849
853
854 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
855 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
856
857 const ResultType tmp( serial( rhs ) );
858 schurAssign( *lhs, tmp );
859 }
861 //**********************************************************************************************
862
863 //**Schur product assignment to dense matrices**************************************************
878 template< typename MT // Type of the target dense matrix
879 , bool SO2 > // Storage order of the target dense matrix
880 friend inline auto schurAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
881 -> EnableIf_t< UseAssign_v<MT> && IsCommutative_v<MT1,MT2> >
882 {
884
885 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
886 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
887
888 if( !RequiresEvaluation_v<MT2> ) {
889 schurAssign( *lhs, rhs.rhs_ );
890 schurAssign( *lhs, rhs.lhs_ );
891 }
892 else {
893 schurAssign( *lhs, rhs.lhs_ );
894 schurAssign( *lhs, rhs.rhs_ );
895 }
896 }
898 //**********************************************************************************************
899
900 //**Schur product assignment to sparse matrices*************************************************
901 // No special implementation for the Schur product assignment to sparse matrices.
902 //**********************************************************************************************
903
904 //**Multiplication assignment to dense matrices*************************************************
905 // No special implementation for the multiplication assignment to dense matrices.
906 //**********************************************************************************************
907
908 //**Multiplication assignment to sparse matrices************************************************
909 // No special implementation for the multiplication assignment to sparse matrices.
910 //**********************************************************************************************
911
912 //**SMP assignment to dense matrices************************************************************
927 template< typename MT // Type of the target dense matrix
928 , bool SO2 > // Storage order of the target dense matrix
929 friend inline auto smpAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
930 -> EnableIf_t< UseSMPAssign_v<MT> && !IsCommutative_v<MT1,MT2> >
931 {
933
934 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
935 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
936
937 if( !IsOperation_v<MT1> && isSame( *lhs, rhs.lhs_ ) ) {
938 smpSchurAssign( *lhs, rhs.rhs_ );
939 }
940 else {
941 CT1 A( rhs.lhs_ );
942 CT2 B( rhs.rhs_ );
943 smpAssign( *lhs, A % B );
944 }
945 }
947 //**********************************************************************************************
948
949 //**SMP assignment to dense matrices************************************************************
964 template< typename MT // Type of the target dense matrix
965 , bool SO2 > // Storage order of the target dense matrix
966 friend inline auto smpAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
967 -> EnableIf_t< UseSMPAssign_v<MT> && IsCommutative_v<MT1,MT2> >
968 {
970
971 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
972 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
973
974 if( !IsOperation_v<MT1> && isSame( *lhs, rhs.lhs_ ) ) {
975 smpSchurAssign( *lhs, rhs.rhs_ );
976 }
977 else if( !IsOperation_v<MT2> && isSame( *lhs, rhs.rhs_ ) ) {
978 smpSchurAssign( *lhs, rhs.lhs_ );
979 }
980 else if( !RequiresEvaluation_v<MT2> ) {
981 smpAssign ( *lhs, rhs.rhs_ );
982 smpSchurAssign( *lhs, rhs.lhs_ );
983 }
984 else {
985 smpAssign ( *lhs, rhs.lhs_ );
986 smpSchurAssign( *lhs, rhs.rhs_ );
987 }
988 }
990 //**********************************************************************************************
991
992 //**SMP assignment to sparse matrices***********************************************************
1006 template< typename MT // Type of the target sparse matrix
1007 , bool SO2 > // Storage order of the target sparse matrix
1008 friend inline auto smpAssign( SparseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
1009 -> EnableIf_t< UseSMPAssign_v<MT> >
1010 {
1012
1013 using TmpType = If_t< SO == SO2, ResultType, OppositeType >;
1014
1021
1022 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
1023 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
1024
1025 const TmpType tmp( rhs );
1026 smpAssign( *lhs, tmp );
1027 }
1029 //**********************************************************************************************
1030
1031 //**SMP addition assignment to dense matrices***************************************************
1045 template< typename MT // Type of the target dense matrix
1046 , bool SO2 > // Storage order of the target dense matrix
1047 friend inline auto smpAddAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
1048 -> EnableIf_t< UseAssign_v<MT> >
1049 {
1051
1055
1056 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
1057 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
1058
1059 const ResultType tmp( rhs );
1060 smpAddAssign( *lhs, tmp );
1061 }
1063 //**********************************************************************************************
1064
1065 //**SMP addition assignment to sparse matrices**************************************************
1066 // No special implementation for the SMP addition assignment to sparse matrices.
1067 //**********************************************************************************************
1068
1069 //**SMP subtraction assignment to dense matrices************************************************
1084 template< typename MT // Type of the target dense matrix
1085 , bool SO2 > // Storage order of the target dense matrix
1086 friend inline auto smpSubAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
1087 -> EnableIf_t< UseSMPAssign_v<MT> >
1088 {
1090
1094
1095 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
1096 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
1097
1098 const ResultType tmp( rhs );
1099 smpSubAssign( *lhs, tmp );
1100 }
1102 //**********************************************************************************************
1103
1104 //**SMP subtraction assignment to sparse matrices***********************************************
1105 // No special implementation for the SMP subtraction assignment to sparse matrices.
1106 //**********************************************************************************************
1107
1108 //**SMP Schur product assignment to dense matrices**********************************************
1123 template< typename MT // Type of the target dense matrix
1124 , bool SO2 > // Storage order of the target dense matrix
1125 friend inline auto smpSchurAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
1126 -> EnableIf_t< UseSMPAssign_v<MT> && !IsCommutative_v<MT1,MT2> >
1127 {
1129
1133
1134 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
1135 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
1136
1137 const ResultType tmp( rhs );
1138 smpSchurAssign( *lhs, tmp );
1139 }
1141 //**********************************************************************************************
1142
1143 //**SMP Schur product assignment to dense matrices**********************************************
1158 template< typename MT // Type of the target dense matrix
1159 , bool SO2 > // Storage order of the target dense matrix
1160 friend inline auto smpSchurAssign( DenseMatrix<MT,SO2>& lhs, const DMatDMatSchurExpr& rhs )
1161 -> EnableIf_t< UseSMPAssign_v<MT> && IsCommutative_v<MT1,MT2> >
1162 {
1164
1165 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == rhs.rows() , "Invalid number of rows" );
1166 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == rhs.columns(), "Invalid number of columns" );
1167
1168 if( !RequiresEvaluation_v<MT2> ) {
1169 smpSchurAssign( *lhs, rhs.rhs_ );
1170 smpSchurAssign( *lhs, rhs.lhs_ );
1171 }
1172 else {
1173 smpSchurAssign( *lhs, rhs.lhs_ );
1174 smpSchurAssign( *lhs, rhs.rhs_ );
1175 }
1176 }
1178 //**********************************************************************************************
1179
1180 //**SMP Schur product assignment to sparse matrices*********************************************
1181 // No special implementation for the SMP Schur product assignment to sparse matrices.
1182 //**********************************************************************************************
1183
1184 //**SMP multiplication assignment to dense matrices*********************************************
1185 // No special implementation for the SMP multiplication assignment to dense matrices.
1186 //**********************************************************************************************
1187
1188 //**SMP multiplication assignment to sparse matrices********************************************
1189 // No special implementation for the SMP multiplication assignment to sparse matrices.
1190 //**********************************************************************************************
1191
1192 //**Compile time checks*************************************************************************
1198 //**********************************************************************************************
1199};
1200//*************************************************************************************************
1201
1202
1203
1204
1205//=================================================================================================
1206//
1207// GLOBAL BINARY ARITHMETIC OPERATORS
1208//
1209//=================================================================================================
1210
1211//*************************************************************************************************
1224template< typename MT1 // Type of the left-hand side dense matrix
1225 , typename MT2 // Type of the right-hand side dense matrix
1226 , bool SO // Storage order
1227 , DisableIf_t< ( IsUniLower_v<MT1> && IsUniUpper_v<MT2> ) ||
1228 ( IsUniUpper_v<MT1> && IsUniLower_v<MT2> ) ||
1229 ( IsStrictlyLower_v<MT1> && IsUpper_v<MT2> ) ||
1230 ( IsStrictlyUpper_v<MT1> && IsLower_v<MT2> ) ||
1231 ( IsLower_v<MT1> && IsStrictlyUpper_v<MT2> ) ||
1232 ( IsUpper_v<MT1> && IsStrictlyLower_v<MT2> ) >* = nullptr >
1233inline const DMatDMatSchurExpr<MT1,MT2,SO>
1234 dmatdmatschur( const DenseMatrix<MT1,SO>& lhs, const DenseMatrix<MT2,SO>& rhs )
1235{
1237
1238 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == (*rhs).rows() , "Invalid number of rows" );
1239 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == (*rhs).columns(), "Invalid number of columns" );
1240
1241 return DMatDMatSchurExpr<MT1,MT2,SO>( *lhs, *rhs );
1242}
1244//*************************************************************************************************
1245
1246
1247//*************************************************************************************************
1260template< typename MT1 // Type of the left-hand side dense matrix
1261 , typename MT2 // Type of the right-hand side dense matrix
1262 , bool SO // Storage order
1263 , EnableIf_t< ( IsUniLower_v<MT1> && IsUniUpper_v<MT2> ) ||
1264 ( IsUniUpper_v<MT1> && IsUniLower_v<MT2> ) >* = nullptr >
1265inline decltype(auto)
1266 dmatdmatschur( const DenseMatrix<MT1,SO>& lhs, const DenseMatrix<MT2,SO>& rhs )
1267{
1269
1270 MAYBE_UNUSED( rhs );
1271
1272 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == (*rhs).rows() , "Invalid number of rows" );
1273 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == (*rhs).columns(), "Invalid number of columns" );
1274
1275 using ReturnType = const SchurTrait_t< ResultType_t<MT1>, ResultType_t<MT2> >;
1276
1279
1280 return ReturnType( (*lhs).rows() );
1281}
1283//*************************************************************************************************
1284
1285
1286//*************************************************************************************************
1299template< typename MT1 // Type of the left-hand side dense matrix
1300 , typename MT2 // Type of the right-hand side dense matrix
1301 , bool SO // Storage order
1302 , EnableIf_t< ( IsStrictlyLower_v<MT1> && IsUpper_v<MT2> ) ||
1303 ( IsStrictlyUpper_v<MT1> && IsLower_v<MT2> ) ||
1304 ( IsLower_v<MT1> && IsStrictlyUpper_v<MT2> ) ||
1305 ( IsUpper_v<MT1> && IsStrictlyLower_v<MT2> ) >* = nullptr >
1306inline decltype(auto)
1307 dmatdmatschur( const DenseMatrix<MT1,SO>& lhs, const DenseMatrix<MT2,SO>& rhs )
1308{
1310
1311 MAYBE_UNUSED( rhs );
1312
1313 BLAZE_INTERNAL_ASSERT( (*lhs).rows() == (*rhs).rows() , "Invalid number of rows" );
1314 BLAZE_INTERNAL_ASSERT( (*lhs).columns() == (*rhs).columns(), "Invalid number of columns" );
1315
1316 using ReturnType = const SchurTrait_t< ResultType_t<MT1>, ResultType_t<MT2> >;
1317
1320
1321 return ReturnType( (*lhs).rows(), (*lhs).columns() );
1322}
1324//*************************************************************************************************
1325
1326
1327//*************************************************************************************************
1352template< typename MT1 // Type of the left-hand side dense matrix
1353 , typename MT2 // Type of the right-hand side dense matrix
1354 , bool SO > // Storage order
1355inline decltype(auto)
1356 operator%( const DenseMatrix<MT1,SO>& lhs, const DenseMatrix<MT2,SO>& rhs )
1357{
1359
1360 if( (*lhs).rows() != (*rhs).rows() || (*lhs).columns() != (*rhs).columns() ) {
1361 BLAZE_THROW_INVALID_ARGUMENT( "Matrix sizes do not match" );
1362 }
1363
1364 return dmatdmatschur( *lhs, *rhs );
1365}
1366//*************************************************************************************************
1367
1368
1369
1370
1371//=================================================================================================
1372//
1373// ISALIGNED SPECIALIZATIONS
1374//
1375//=================================================================================================
1376
1377//*************************************************************************************************
1379template< typename MT1, typename MT2, bool SO >
1380struct IsAligned< DMatDMatSchurExpr<MT1,MT2,SO> >
1381 : public BoolConstant< IsAligned_v<MT1> && IsAligned_v<MT2> >
1382{};
1384//*************************************************************************************************
1385
1386
1387
1388
1389//=================================================================================================
1390//
1391// ISPADDED SPECIALIZATIONS
1392//
1393//=================================================================================================
1394
1395//*************************************************************************************************
1397template< typename MT1, typename MT2, bool SO >
1398struct IsPadded< DMatDMatSchurExpr<MT1,MT2,SO> >
1399 : public BoolConstant< IsPadded_v<MT1> && IsPadded_v<MT2> >
1400{};
1402//*************************************************************************************************
1403
1404} // namespace blaze
1405
1406#endif
Header file for auxiliary alias declarations.
typename T::CompositeType CompositeType_t
Alias declaration for nested CompositeType type definitions.
Definition: Aliases.h:110
typename T::ReturnType ReturnType_t
Alias declaration for nested ReturnType type definitions.
Definition: Aliases.h:470
typename T::ResultType ResultType_t
Alias declaration for nested ResultType type definitions.
Definition: Aliases.h:450
typename T::ElementType ElementType_t
Alias declaration for nested ElementType type definitions.
Definition: Aliases.h:190
typename T::OppositeType OppositeType_t
Alias declaration for nested OppositeType type definitions.
Definition: Aliases.h:310
typename T::ConstIterator ConstIterator_t
Alias declaration for nested ConstIterator type definitions.
Definition: Aliases.h:130
typename T::TransposeType TransposeType_t
Alias declaration for nested TransposeType type definitions.
Definition: Aliases.h:550
Header file for run time assertion macros.
Header file for the EnableIf class template.
Header file for the function trace functionality.
Header file for the HasSIMDMult type trait.
Macro for CUDA compatibility.
Constraint on the data type.
Header file for the If class template.
Header file for the IntegralConstant class template.
Header file for the IsAligned type trait.
Header file for the IsCommutative type trait.
Header file for the IsExpression type trait class.
Header file for the IsLower type trait.
Header file for the IsOperation type trait class.
Header file for the IsPadded type trait.
Header file for the IsStrictlyLower type trait.
Header file for the IsStrictlyUpper type trait.
Header file for the IsTemporary type trait class.
Header file for the IsUniLower type trait.
Header file for the IsUniUpper type trait.
Header file for the IsUpper type trait.
Header file for the MAYBE_UNUSED function template.
Header file for all SIMD functionality.
Header file for the Schur product trait.
Constraint on the data type.
Iterator over the elements of the dense matrix.
Definition: DMatDMatSchurExpr.h:196
BLAZE_DEVICE_CALLABLE bool operator<=(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatSchurExpr.h:371
BLAZE_DEVICE_CALLABLE bool operator>(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatSchurExpr.h:360
BLAZE_DEVICE_CALLABLE DifferenceType operator-(const ConstIterator &rhs) const
Calculating the number of elements between two iterators.
Definition: DMatDMatSchurExpr.h:393
ReferenceType reference
Reference return type.
Definition: DMatDMatSchurExpr.h:209
ValueType value_type
Type of the underlying elements.
Definition: DMatDMatSchurExpr.h:207
DifferenceType difference_type
Difference between two iterators.
Definition: DMatDMatSchurExpr.h:210
ptrdiff_t DifferenceType
Difference between two iterators.
Definition: DMatDMatSchurExpr.h:203
BLAZE_DEVICE_CALLABLE ConstIterator & operator++()
Pre-increment operator.
Definition: DMatDMatSchurExpr.h:262
ConstIterator_t< MT1 > LeftIteratorType
ConstIterator type of the left-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:213
BLAZE_DEVICE_CALLABLE ConstIterator(LeftIteratorType left, RightIteratorType right)
Constructor for the ConstIterator class.
Definition: DMatDMatSchurExpr.h:225
BLAZE_DEVICE_CALLABLE ConstIterator & operator-=(size_t dec)
Subtraction assignment operator.
Definition: DMatDMatSchurExpr.h:250
ElementType & ReferenceType
Reference return type.
Definition: DMatDMatSchurExpr.h:202
std::random_access_iterator_tag IteratorCategory
The iterator category.
Definition: DMatDMatSchurExpr.h:199
ReturnType operator*() const
Direct access to the element at the current iterator position.
Definition: DMatDMatSchurExpr.h:306
friend const ConstIterator operator+(const ConstIterator &it, size_t inc)
Addition between a ConstIterator and an integral value.
Definition: DMatDMatSchurExpr.h:405
BLAZE_DEVICE_CALLABLE bool operator!=(const ConstIterator &rhs) const
Inequality comparison between two ConstIterator objects.
Definition: DMatDMatSchurExpr.h:338
BLAZE_DEVICE_CALLABLE bool operator<(const ConstIterator &rhs) const
Less-than comparison between two ConstIterator objects.
Definition: DMatDMatSchurExpr.h:349
PointerType pointer
Pointer return type.
Definition: DMatDMatSchurExpr.h:208
friend const ConstIterator operator+(size_t inc, const ConstIterator &it)
Addition between an integral value and a ConstIterator.
Definition: DMatDMatSchurExpr.h:417
BLAZE_DEVICE_CALLABLE bool operator>=(const ConstIterator &rhs) const
Greater-than comparison between two ConstIterator objects.
Definition: DMatDMatSchurExpr.h:382
LeftIteratorType left_
Iterator to the current left-hand side element.
Definition: DMatDMatSchurExpr.h:436
ElementType * PointerType
Pointer return type.
Definition: DMatDMatSchurExpr.h:201
auto load() const noexcept
Access to the SIMD elements of the matrix.
Definition: DMatDMatSchurExpr.h:316
ElementType ValueType
Type of the underlying elements.
Definition: DMatDMatSchurExpr.h:200
ConstIterator_t< MT2 > RightIteratorType
ConstIterator type of the right-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:216
RightIteratorType right_
Iterator to the current right-hand side element.
Definition: DMatDMatSchurExpr.h:437
BLAZE_DEVICE_CALLABLE ConstIterator & operator--()
Pre-decrement operator.
Definition: DMatDMatSchurExpr.h:284
IteratorCategory iterator_category
The iterator category.
Definition: DMatDMatSchurExpr.h:206
BLAZE_DEVICE_CALLABLE ConstIterator & operator+=(size_t inc)
Addition assignment operator.
Definition: DMatDMatSchurExpr.h:237
friend BLAZE_DEVICE_CALLABLE const ConstIterator operator-(const ConstIterator &it, size_t dec)
Subtraction between a ConstIterator and an integral value.
Definition: DMatDMatSchurExpr.h:429
BLAZE_DEVICE_CALLABLE const ConstIterator operator++(int)
Post-increment operator.
Definition: DMatDMatSchurExpr.h:274
BLAZE_DEVICE_CALLABLE bool operator==(const ConstIterator &rhs) const
Equality comparison between two ConstIterator objects.
Definition: DMatDMatSchurExpr.h:327
BLAZE_DEVICE_CALLABLE const ConstIterator operator--(int)
Post-decrement operator.
Definition: DMatDMatSchurExpr.h:296
Expression object for dense matrix-dense matrix Schur products.
Definition: DMatDMatSchurExpr.h:108
static constexpr size_t SIMDSIZE
The number of elements packed within a single SIMD element.
Definition: DMatDMatSchurExpr.h:453
CompositeType_t< MT2 > CT2
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:116
static constexpr bool returnExpr
Compilation switch for the selection of the subscript operator return type.
Definition: DMatDMatSchurExpr.h:128
ElementType_t< MT1 > ET1
Element type of the left-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:117
static constexpr bool smpAssignable
Compilation switch for the expression template assignment strategy.
Definition: DMatDMatSchurExpr.h:448
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: DMatDMatSchurExpr.h:493
static constexpr bool simdEnabled
Compilation switch for the expression template evaluation strategy.
Definition: DMatDMatSchurExpr.h:444
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatDMatSchurExpr.h:478
LeftOperand lhs_
Left-hand side dense matrix of the Schur product expression.
Definition: DMatDMatSchurExpr.h:630
static constexpr bool useAssign
Compilation switch for the serial evaluation strategy of the Schur product expression.
Definition: DMatDMatSchurExpr.h:142
bool isAligned() const noexcept
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatDMatSchurExpr.h:612
DMatDMatSchurExpr(const MT1 &lhs, const MT2 &rhs) noexcept
Constructor for the DMatDMatSchurExpr class.
Definition: DMatDMatSchurExpr.h:462
ElementType_t< MT2 > ET2
Element type of the right-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:118
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: DMatDMatSchurExpr.h:622
ReturnType_t< MT1 > RN1
Return type of the left-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:113
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: DMatDMatSchurExpr.h:602
OppositeType_t< ResultType > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatDMatSchurExpr.h:175
LeftOperand leftOperand() const noexcept
Returns the left-hand side dense matrix operand.
Definition: DMatDMatSchurExpr.h:567
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: DMatDMatSchurExpr.h:547
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row/column i.
Definition: DMatDMatSchurExpr.h:526
ElementType_t< ResultType > ElementType
Resulting element type.
Definition: DMatDMatSchurExpr.h:177
BLAZE_ALWAYS_INLINE auto load(size_t i, size_t j) const noexcept
Access to the SIMD elements of the matrix.
Definition: DMatDMatSchurExpr.h:511
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: DMatDMatSchurExpr.h:557
If_t< IsExpression_v< MT1 >, const MT1, const MT1 & > LeftOperand
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:186
const If_t< returnExpr, ExprReturnType, ElementType > ReturnType
Return type for expression template evaluations.
Definition: DMatDMatSchurExpr.h:180
RightOperand rightOperand() const noexcept
Returns the right-hand side dense matrix operand.
Definition: DMatDMatSchurExpr.h:577
SchurTrait_t< RT1, RT2 > ResultType
Result type for expression template evaluations.
Definition: DMatDMatSchurExpr.h:174
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: DMatDMatSchurExpr.h:589
ResultType_t< MT2 > RT2
Result type of the right-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:112
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row/column i.
Definition: DMatDMatSchurExpr.h:537
TransposeType_t< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: DMatDMatSchurExpr.h:176
CompositeType_t< MT1 > CT1
Composite type of the left-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:115
If_t< useAssign, const ResultType, const DMatDMatSchurExpr & > CompositeType
Data type for composite expression templates.
Definition: DMatDMatSchurExpr.h:183
If_t< IsExpression_v< MT2 >, const MT2, const MT2 & > RightOperand
Composite type of the right-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:189
ResultType_t< MT1 > RT1
Result type of the left-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:111
decltype(std::declval< RN1 >() *std::declval< RN2 >()) ExprReturnType
Expression return type for the subscript operator.
Definition: DMatDMatSchurExpr.h:131
ReturnType_t< MT2 > RN2
Return type of the right-hand side dense matrix expression.
Definition: DMatDMatSchurExpr.h:114
RightOperand rhs_
Right-hand side dense matrix of the Schur product expression.
Definition: DMatDMatSchurExpr.h:631
Base class for dense matrices.
Definition: DenseMatrix.h:82
SIMD characteristics of data types.
Definition: SIMDTrait.h:297
Constraint on the data type.
Constraint on the data type.
Constraint on the data type.
Header file for the Computation base class.
Header file for the DenseMatrix base class.
Header file for the SchurExpr base class.
decltype(auto) serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:812
#define BLAZE_CONSTRAINT_MATRICES_MUST_HAVE_SAME_STORAGE_ORDER(T1, T2)
Constraint on the data type.
Definition: StorageOrder.h:84
#define BLAZE_CONSTRAINT_MUST_NOT_REQUIRE_EVALUATION(T)
Constraint on the data type.
Definition: RequiresEvaluation.h:81
#define BLAZE_CONSTRAINT_MUST_BE_DENSE_MATRIX_TYPE(T)
Constraint on the data type.
Definition: DenseMatrix.h:61
#define BLAZE_CONSTRAINT_MUST_BE_IDENTITY_MATRIX_TYPE(T)
Constraint on the data type.
Definition: Identity.h:60
#define BLAZE_CONSTRAINT_MUST_FORM_VALID_SCHUREXPR(T1, T2)
Constraint on the data type.
Definition: SchurExpr.h:103
#define BLAZE_CONSTRAINT_MUST_BE_ZERO_TYPE(T)
Constraint on the data type.
Definition: Zero.h:61
#define BLAZE_CONSTRAINT_MUST_BE_MATRIX_WITH_STORAGE_ORDER(T, SO)
Constraint on the data type.
Definition: StorageOrder.h:63
typename SchurTrait< T1, T2 >::Type SchurTrait_t
Auxiliary alias declaration for the SchurTrait class template.
Definition: SchurTrait.h:137
bool isSame(const Matrix< MT1, SO1 > &a, const Matrix< MT2, SO2 > &b) noexcept
Returns whether the two given matrices represent the same observable state.
Definition: Matrix.h:1424
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.
Definition: Assert.h:101
auto smpSubAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs) -> EnableIf_t< IsDenseMatrix_v< MT1 > >
Default implementation of the SMP subtraction assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:162
auto smpAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs) -> EnableIf_t< IsDenseMatrix_v< MT1 > >
Default implementation of the SMP assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:100
auto smpSchurAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs) -> EnableIf_t< IsDenseMatrix_v< MT1 > >
Default implementation of the SMP Schur product assignment of a matrix to dense matrix.
Definition: DenseMatrix.h:194
auto smpAddAssign(Matrix< MT1, SO1 > &lhs, const Matrix< MT2, SO2 > &rhs) -> EnableIf_t< IsDenseMatrix_v< MT1 > >
Default implementation of the SMP addition assignment of a matrix to a dense matrix.
Definition: DenseMatrix.h:131
#define BLAZE_ALWAYS_INLINE
Platform dependent setup of an enforced inline keyword.
Definition: Inline.h:85
#define BLAZE_DEVICE_CALLABLE
Conditional macro that sets host and device attributes when compiled with CUDA.
Definition: HostDevice.h:94
typename EnableIf< Condition, T >::Type EnableIf_t
Auxiliary type for the EnableIf class template.
Definition: EnableIf.h:138
constexpr void MAYBE_UNUSED(const Args &...)
Suppression of unused parameter warnings.
Definition: MaybeUnused.h:81
typename If< Condition >::template Type< T1, T2 > If_t
Auxiliary alias template for the If class template.
Definition: If.h:108
IntegralConstant< bool, B > BoolConstant
Generic wrapper for a compile time constant boolean value.
Definition: IntegralConstant.h:110
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exception.
Definition: Exception.h:331
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.
Definition: Exception.h:235
typename EnableIf<!Condition, T >::Type DisableIf_t
Auxiliary type for the EnableIf class template.
Definition: EnableIf.h:175
#define BLAZE_FUNCTION_TRACE
Function trace macro.
Definition: FunctionTrace.h:94
Header file for the exception macros of the math module.
Constraints on the storage order of matrix types.
Header file for all forward declarations for expression class templates.
Header file for all forward declarations for sparse vectors and matrices.
Header file for the serial shim.
Base class for all compute expression templates.
Definition: Computation.h:68
Base class for all Schur product expression templates.
Definition: SchurExpr.h:68
System settings for the inline keywords.
Header file for the thresholds for matrix/vector and matrix/matrix multiplications.
Header file for the RequiresEvaluation type trait.
Header file for basic type definitions.