35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATDECLDIAGEXPR_H_ 36 #define _BLAZE_MATH_EXPRESSIONS_DMATDECLDIAGEXPR_H_ 95 :
public DeclDiagExpr< DenseMatrix< DMatDeclDiagExpr<MT,SO>, SO > >
114 static constexpr
bool useAssign = RequiresEvaluation_v<MT>;
117 template<
typename MT2 >
119 static constexpr
bool UseAssign_v =
useAssign;
131 template<
typename MT2 >
204 if( i >=
dm_.rows() ) {
207 if( j >=
dm_.columns() ) {
226 return dm_.load(i,j);
267 inline size_t rows() const noexcept {
278 return dm_.columns();
298 template<
typename T >
299 inline bool canAlias(
const T* alias )
const noexcept {
300 return dm_.canAlias( alias );
310 template<
typename T >
311 inline bool isAliased(
const T* alias )
const noexcept {
312 return dm_.isAliased( alias );
322 return dm_.isAligned();
332 return dm_.canSMPAssign();
353 template<
typename MT2
363 assign( ~lhs, rhs.dm_ );
380 template<
typename MT2
390 assign( ~lhs, rhs.dm_ );
407 template<
typename MT2
409 friend inline auto addAssign( DenseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
410 -> EnableIf_t< UseAssign_v<MT2> >
417 addAssign( ~lhs, rhs.dm_ );
434 template<
typename MT2
436 friend inline auto addAssign( SparseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
437 -> EnableIf_t< UseAssign_v<MT2> >
444 addAssign( ~lhs, rhs.dm_ );
461 template<
typename MT2
463 friend inline auto subAssign( DenseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
464 -> EnableIf_t< UseAssign_v<MT2> >
471 subAssign( ~lhs, rhs.dm_ );
488 template<
typename MT2
490 friend inline auto subAssign( SparseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
491 -> EnableIf_t< UseAssign_v<MT2> >
498 subAssign( ~lhs, rhs.dm_ );
515 template<
typename MT2
517 friend inline auto schurAssign( DenseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
518 -> EnableIf_t< UseAssign_v<MT2> >
525 schurAssign( ~lhs, rhs.dm_ );
542 template<
typename MT2
544 friend inline auto schurAssign( SparseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
545 -> EnableIf_t< UseAssign_v<MT2> >
552 schurAssign( ~lhs, rhs.dm_ );
569 template<
typename MT2
571 friend inline auto multAssign( DenseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
572 -> EnableIf_t< UseAssign_v<MT2> >
579 multAssign( ~lhs, rhs.dm_ );
596 template<
typename MT2
598 friend inline auto multAssign( SparseMatrix<MT2,SO2>& lhs,
const DMatDeclDiagExpr& rhs )
599 -> EnableIf_t< UseAssign_v<MT2> >
606 multAssign( ~lhs, rhs.dm_ );
623 template<
typename MT2
626 -> EnableIf_t< UseSMPAssign_v<MT2> >
650 template<
typename MT2
653 -> EnableIf_t< UseSMPAssign_v<MT2> >
677 template<
typename MT2
680 -> EnableIf_t< UseSMPAssign_v<MT2> >
704 template<
typename MT2
707 -> EnableIf_t< UseSMPAssign_v<MT2> >
731 template<
typename MT2
734 -> EnableIf_t< UseSMPAssign_v<MT2> >
758 template<
typename MT2
761 -> EnableIf_t< UseSMPAssign_v<MT2> >
785 template<
typename MT2
788 -> EnableIf_t< UseSMPAssign_v<MT2> >
812 template<
typename MT2
815 -> EnableIf_t< UseSMPAssign_v<MT2> >
840 template<
typename MT2
843 -> EnableIf_t< UseSMPAssign_v<MT2> >
868 template<
typename MT2
871 -> EnableIf_t< UseSMPAssign_v<MT2> >
914 template<
typename MT
916 , DisableIf_t< IsDiagonal_v<MT> >* =
nullptr >
917 inline const DMatDeclDiagExpr<MT,SO> decldiag_backend(
const DenseMatrix<MT,SO>& dm )
923 return DMatDeclDiagExpr<MT,SO>( ~dm );
940 template<
typename MT
942 , EnableIf_t< IsDiagonal_v<MT> >* =
nullptr >
943 inline const MT& decldiag_backend(
const DenseMatrix<MT,SO>& dm )
973 template<
typename MT
983 return decldiag_backend( ~dm );
1010 template<
typename MT
1013 , DisableIf_t< IsDiagonal_v<MT> >* =
nullptr >
1014 inline decltype(
auto)
decldiag( const DMatScalarMultExpr<MT,ST,SO>& dm )
1022 return decldiag( dm.leftOperand() ) * dm.rightOperand();
1038 template<
typename MT,
bool SO >
1039 struct HasConstDataAccess< DMatDeclDiagExpr<MT,SO> >
1040 :
public HasConstDataAccess<MT>
1056 template<
typename MT,
bool SO >
1057 struct IsAligned< DMatDeclDiagExpr<MT,SO> >
1058 :
public IsAligned<MT>
1074 template<
typename MT,
bool SO >
1075 struct IsSymmetric< DMatDeclDiagExpr<MT,SO> >
1092 template<
typename MT,
bool SO >
1093 struct IsLower< DMatDeclDiagExpr<MT,SO> >
1110 template<
typename MT,
bool SO >
1111 struct IsUpper< DMatDeclDiagExpr<MT,SO> >
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.This macro encapsulates the default way o...
Definition: Exception.h:235
Header file for auxiliary alias declarations.
Header file for the decldiag trait.
decltype(auto) decldiag(const DenseMatrix< MT, SO > &dm)
Declares the given dense matrix expression dm as diagonal.
Definition: DMatDeclDiagExpr.h:975
Header file for basic type definitions.
typename If< Condition, T1, T2 >::Type If_t
Auxiliary alias declaration for the If class template.The If_t alias declaration provides a convenien...
Definition: If.h:109
typename T::ResultType ResultType_t
Alias declaration for nested ResultType type definitions.The ResultType_t alias declaration provides ...
Definition: Aliases.h:390
BLAZE_ALWAYS_INLINE auto load(size_t i, size_t j) const noexcept
Access to the SIMD elements of the matrix.
Definition: DMatDeclDiagExpr.h:221
#define BLAZE_CONSTRAINT_MUST_BE_MATRIX_WITH_STORAGE_ORDER(T, SO)
Constraint on the data type.In case the given data type T is not a dense or sparse matrix type and in...
Definition: StorageOrder.h:63
Header file for the IsDiagonal type trait.
#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:61
Base class for all declaration expression templates.The Declaration class serves as a tag for all dec...
Definition: Declaration.h:70
If_t< IsExpression_v< MT >, const MT, const MT &> Operand
Composite data type of the dense matrix expression.
Definition: DMatDeclDiagExpr.h:153
static constexpr bool smpAssignable
Compilation flag for SMP assignments.
Definition: CompressedMatrix.h:3113
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: DMatDeclDiagExpr.h:277
BLAZE_CREATE_GET_TYPE_MEMBER_TYPE_TRAIT(GetConstIterator, ConstIterator, INVALID_TYPE)
Definition of the GetConstIterator type trait.
static constexpr bool simdEnabled
Compilation switch for the expression template evaluation strategy.
Definition: DMatDeclDiagExpr.h:158
BoolConstant< true > TrueType
Type traits base class.The TrueType class is used as base class for type traits and value traits that...
Definition: TrueType.h:61
Header file for the RequiresEvaluation type trait.
Expression object for the explicit diagonal declaration of dense matrices.The DMatDeclDiagExpr class ...
Definition: DMatDeclDiagExpr.h:94
typename T::ReturnType ReturnType_t
Alias declaration for nested ReturnType type definitions.The ReturnType_t alias declaration provides ...
Definition: Aliases.h:410
Header file for the SIMD trait.
Base class for dense matrices.The DenseMatrix class is a base class for all dense matrix classes...
Definition: DenseMatrix.h:80
Base class for sparse matrices.The SparseMatrix class is a base class for all sparse matrix classes...
Definition: Forward.h:137
typename T::ElementType ElementType_t
Alias declaration for nested ElementType type definitions.The ElementType_t alias declaration provide...
Definition: Aliases.h:170
Header file for all forward declarations for sparse vectors and matrices.
Base class for all decldiag expression templates.The DeclDiagExpr class serves as a tag for all expre...
Definition: DeclDiagExpr.h:66
Constraint on the data type.
typename EnableIf< Condition, T >::Type EnableIf_t
Auxiliary type for the EnableIf class template.The EnableIf_t alias declaration provides a convenient...
Definition: EnableIf.h:138
TransposeType_t< ResultType > TransposeType
Transpose type for expression template evaluations.
Definition: DMatDeclDiagExpr.h:142
static constexpr size_t SIMDSIZE
The number of elements packed within a single SIMD element.
Definition: DMatDeclDiagExpr.h:166
Header file for the DisableIf class template.
Header file for the IsSymmetric type trait.
ReturnType operator()(size_t i, size_t j) const
2D-access to the matrix elements.
Definition: DMatDeclDiagExpr.h:188
Namespace of the Blaze C++ math library.
Definition: Blaze.h:58
#define BLAZE_ALWAYS_INLINE
Platform dependent setup of an enforced inline keyword.
Definition: Inline.h:85
Header file for the If class template.
Header file for the DeclDiagExpr base class.
bool canSMPAssign() const noexcept
Returns whether the expression can be used in SMP assignments.
Definition: DMatDeclDiagExpr.h:331
#define BLAZE_THROW_OUT_OF_RANGE(MESSAGE)
Macro for the emission of a std::out_of_range exception.This macro encapsulates the default way of Bl...
Definition: Exception.h:331
Header file for the DenseMatrix base class.
ReturnType_t< MT > ReturnType
Return type for expression template evaluations.
Definition: DMatDeclDiagExpr.h:144
GetConstIterator_t< MT > ConstIterator
Iterator over the elements of the dense matrix.
Definition: DMatDeclDiagExpr.h:150
#define BLAZE_CONSTRAINT_MUST_NOT_BE_DIAGONAL_MATRIX_TYPE(T)
Constraint on the data type.In case the given data type T is a diagonal matrix type, a compilation error is created.
Definition: Diagonal.h:79
Header file for the GetMemberType type trait.
Header file for the IsLower type trait.
Header file for the IsAligned type trait.
OppositeType_t< ResultType > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatDeclDiagExpr.h:141
#define BLAZE_CONSTRAINT_MUST_NOT_BE_MATMATMULTEXPR_TYPE(T)
Constraint on the data type.In case the given data type T is a matrix/matrix multiplication expressio...
Definition: MatMatMultExpr.h:83
If_t< RequiresEvaluation_v< MT >, const ResultType, const DMatDeclDiagExpr &> CompositeType
Data type for composite expression templates.
Definition: DMatDeclDiagExpr.h:147
bool isAligned() const noexcept
Returns whether the operands of the expression are properly aligned in memory.
Definition: DMatDeclDiagExpr.h:321
Constraints on the storage order of matrix types.
Operand operand() const noexcept
Returns the dense matrix operand.
Definition: DMatDeclDiagExpr.h:287
const ElementType * data() const noexcept
Low-level data access to the matrix elements.
Definition: DMatDeclDiagExpr.h:235
DMatDeclDiagExpr(const MT &dm) noexcept
Constructor for the DMatDeclDiagExpr class.
Definition: DMatDeclDiagExpr.h:174
Header file for the exception macros of the math module.
Header file for all forward declarations for expression class templates.
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: DMatDeclDiagExpr.h:267
Header file for the EnableIf class template.
ElementType_t< MT > ElementType
Resulting element type.
Definition: DMatDeclDiagExpr.h:143
typename T::OppositeType OppositeType_t
Alias declaration for nested OppositeType type definitions.The OppositeType_t alias declaration provi...
Definition: Aliases.h:270
Header file for the HasConstDataAccess type trait.
Header file for the Declaration base class.
typename T::TransposeType TransposeType_t
Alias declaration for nested TransposeType type definitions.The TransposeType_t alias declaration pro...
Definition: Aliases.h:470
Header file for run time assertion macros.
Utility type for generic codes.
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
SIMD characteristics of data types.The SIMDTrait class template provides the SIMD characteristics of ...
Definition: SIMDTrait.h:295
#define BLAZE_FUNCTION_TRACE
Function trace macro.This macro can be used to reliably trace function calls. In case function tracin...
Definition: FunctionTrace.h:94
Constraint on the data type.
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
static constexpr bool useAssign
Compilation switch for the serial evaluation strategy of the diagonal declaration expression...
Definition: DMatDeclDiagExpr.h:114
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
ResultType_t< MT > RT
Result type of the dense matrix expression.
Definition: DMatDeclDiagExpr.h:100
static constexpr bool smpAssignable
Compilation switch for the expression template assignment strategy.
Definition: DMatDeclDiagExpr.h:161
ConstIterator begin(size_t i) const
Returns an iterator to the first non-zero element of row/column i.
Definition: DMatDeclDiagExpr.h:246
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: DMatDeclDiagExpr.h:311
Header file for the implementation of the base template of the DiagonalMatrix.
ReturnType at(size_t i, size_t j) const
Checked access to the matrix elements.
Definition: DMatDeclDiagExpr.h:203
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
Operand dm_
Dense matrix of the decldiag expression.
Definition: DMatDeclDiagExpr.h:338
DeclDiagTrait_t< RT > ResultType
Result type for expression template evaluations.
Definition: DMatDeclDiagExpr.h:140
ConstIterator end(size_t i) const
Returns an iterator just past the last non-zero element of row/column i.
Definition: DMatDeclDiagExpr.h:257
Header file for the IsUpper type trait.
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: DMatDeclDiagExpr.h:299
Constraint on the data type.
typename DeclDiagTrait< MT >::Type DeclDiagTrait_t
Auxiliary alias declaration for the DeclDiagTrait type trait.The DeclDiagTrait_t alias declaration pr...
Definition: DeclDiagTrait.h:160
bool isSquare(const Matrix< MT, SO > &matrix) noexcept
Checks if the given matrix is a square matrix.
Definition: Matrix.h:951
#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
auto smpMultAssign(Vector< VT1, TF1 > &lhs, const Vector< VT2, TF2 > &rhs) -> EnableIf_t< IsDenseVector_v< VT1 > >
Default implementation of the SMP multiplication assignment of a vector to a dense vector...
Definition: DenseVector.h:191
Header file for the TrueType type/value trait base class.
Header file for the IsExpression type trait class.
Header file for the function trace functionality.