35 #ifndef _BLAZE_MATH_EXPRESSIONS_DMATINVEXPR_H_
36 #define _BLAZE_MATH_EXPRESSIONS_DMATINVEXPR_H_
121 enum :
bool { simdEnabled =
false };
124 enum :
bool { smpAssignable =
false };
142 inline size_t rows() const noexcept {
143 return dm_.columns();
173 template<
typename T >
174 inline bool canAlias(
const T* alias )
const noexcept {
175 return dm_.isAliased( alias );
185 template<
typename T >
186 inline bool isAliased(
const T* alias )
const noexcept {
187 return dm_.isAliased( alias );
227 template<
typename MT2
237 assign( ~lhs, rhs.
dm_ );
240 invert< DMatInvExpr::getInversionFlag() >( ~lhs );
257 template<
typename MT2
259 friend inline void assign( SparseMatrix<MT2,SO2>& lhs,
const DMatInvExpr& rhs )
263 typedef IfTrue_< SO == SO2, ResultType, OppositeType > TmpType;
275 const TmpType tmp(
serial( rhs ) );
293 template<
typename MT2
295 friend inline void addAssign( DenseMatrix<MT2,SO2>& lhs,
const DMatInvExpr& rhs )
302 const ResultType tmp(
serial( rhs ) );
303 addAssign( ~lhs, tmp );
324 template<
typename MT2
326 friend inline void subAssign( DenseMatrix<MT2,SO2>& lhs,
const DMatInvExpr& rhs )
333 const ResultType tmp(
serial( rhs ) );
334 subAssign( ~lhs, tmp );
399 template<
typename MT
442 template<
typename MT
444 inline typename DMatInvExpr<MT,SO>::Operand
inv(
const DMatInvExpr<MT,SO>& dm )
464 template<
typename MT,
bool SO >
465 struct Rows< DMatInvExpr<MT,SO> > :
public Rows<MT>
481 template<
typename MT,
bool SO >
482 struct Columns< DMatInvExpr<MT,SO> > :
public Columns<MT>
498 template<
typename MT,
bool SO >
499 struct IsSymmetric< DMatInvExpr<MT,SO> >
516 template<
typename MT,
bool SO >
517 struct IsHermitian< DMatInvExpr<MT,SO> >
534 template<
typename MT,
bool SO >
535 struct IsLower< DMatInvExpr<MT,SO> >
552 template<
typename MT,
bool SO >
553 struct IsUniLower< DMatInvExpr<MT,SO> >
570 template<
typename MT,
bool SO >
571 struct IsUpper< DMatInvExpr<MT,SO> >
588 template<
typename MT,
bool SO >
589 struct IsUniUpper< DMatInvExpr<MT,SO> >
606 template<
typename MT >
607 struct DMatInvExprTrait< DMatInvExpr<MT,false> >
611 using Type = If_< And< IsDenseMatrix<MT>, IsRowMajorMatrix<MT> >
612 , Operand_< DMatInvExpr<MT,false> >
622 template<
typename MT >
623 struct TDMatInvExprTrait< DMatInvExpr<MT,true> >
627 using Type = If_< And< IsDenseMatrix<MT>, IsColumnMajorMatrix<MT> >
628 , Operand_< DMatInvExpr<MT,true> >
#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.
TransposeType_< MT > TransposeType
Transpose type for expression template evaluations.
Definition: DMatInvExpr.h:108
Operand dm_
Dense matrix of the inversion expression.
Definition: DMatInvExpr.h:193
Header file for the Rows type trait.
Header file for the IsUniUpper type trait.
BLAZE_ALWAYS_INLINE 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:653
Header file for basic type definitions.
ResultType_< MT > ResultType
Result type for expression template evaluations.
Definition: DMatInvExpr.h:106
Header file for the serial shim.
#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.
Header file for the TDMatInvExprTrait class template.
Flag for the inversion of a Hermitian matrix (same as byLDLH).
Definition: InversionFlag.h:110
#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
bool isAliased(const T *alias) const noexcept
Returns whether the expression is aliased with the given address alias.
Definition: DMatInvExpr.h:186
Flag for the inversion of a diagonal matrix.
Definition: InversionFlag.h:115
Header file for the IsColumnMajorMatrix type trait.
Flag for the inversion of a general matrix (same as byLU).
Definition: InversionFlag.h:108
size_t columns() const noexcept
Returns the current number of columns of the matrix.
Definition: DMatInvExpr.h:152
Flag for the inversion of a upper unitriangular matrix.
Definition: InversionFlag.h:114
Header file for the And class template.
Compile time check for lower triangular matrices.This type trait tests whether or not the given templ...
Definition: IsLower.h:88
const DMatSerialExpr< MT, SO > serial(const DenseMatrix< MT, SO > &dm)
Forces the serial evaluation of the given dense matrix expression dm.
Definition: DMatSerialExpr.h:723
Header file for the Computation base class.
Compile time check for upper triangular matrices.This type trait tests whether or not the given templ...
Definition: IsUpper.h:88
Flag for the inversion of a upper triangular matrix.
Definition: InversionFlag.h:113
Header file for the MatInvExpr base class.
Header file for the IsUniLower type trait.
typename T::ResultType ResultType_
Alias declaration for nested ResultType type definitions.The ResultType_ alias declaration provides a...
Definition: Aliases.h:323
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.
typename T::ReturnType ReturnType_
Alias declaration for nested ReturnType type definitions.The ReturnType_ alias declaration provides a...
Definition: Aliases.h:343
If_< IsExpression< MT >, const MT, const MT & > Operand
Composite data type of the dense matrix expression.
Definition: DMatInvExpr.h:116
DMatInvExpr(const MT &dm) noexcept
Constructor for the DMatInvExpr class.
Definition: DMatInvExpr.h:132
Constraint on the data type.
typename T::CompositeType CompositeType_
Alias declaration for nested CompositeType type definitions.The CompositeType_ alias declaration prov...
Definition: Aliases.h:83
Compile time check for upper unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniUpper.h:86
Flag for the inversion of a symmetric matrix (same as byLDLT).
Definition: InversionFlag.h:109
Header file for the IsSymmetric type trait.
Namespace of the Blaze C++ math library.
Definition: Blaze.h:57
Header file for the If class template.
OppositeType_< MT > OppositeType
Result type with opposite storage order for expression template evaluations.
Definition: DMatInvExpr.h:107
Header file for the DenseMatrix base class.
Header file for the Columns type trait.
typename T::ElementType ElementType_
Alias declaration for nested ElementType type definitions.The ElementType_ alias declaration provides...
Definition: Aliases.h:163
Header file for the IsLower type trait.
Compile time check for diagonal matrices.This type trait tests whether or not the given template para...
Definition: IsDiagonal.h:90
ReturnType_< MT > ReturnType
Return type for expression template evaluations.
Definition: DMatInvExpr.h:110
#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:60
Operand operand() const noexcept
Returns the dense matrix operand.
Definition: DMatInvExpr.h:162
Constraints on the storage order of matrix types.
Compile time check for symmetric matrices.This type trait tests whether or not the given template par...
Definition: IsSymmetric.h:85
Header file for the exception macros of the math module.
Header file for all forward declarations for expression class templates.
Header file for the IsDenseMatrix type trait.
InversionFlag
Inversion flag.The InversionFlag type enumeration represents the different types of matrix inversion ...
Definition: InversionFlag.h:101
Compile time check for lower unitriangular matrices.This type trait tests whether or not the given te...
Definition: IsUniLower.h:86
Flag for the inversion of a lower triangular matrix.
Definition: InversionFlag.h:111
Header file for the DMatInvExprTrait class template.
Header file for run time assertion macros.
Utility type for generic codes.
typename If< T1, T2, T3 >::Type If_
Auxiliary alias declaration for the If class template.The If_ alias declaration provides a convenient...
Definition: If.h:160
Compile time check for Hermitian matrices.This type trait tests whether or not the given template par...
Definition: IsHermitian.h:85
Base class for all matrix inversion expression templates.The MatInvExpr class serves as a tag for all...
Definition: MatInvExpr.h:65
DMatInvExpr< MT, SO > This
Type of this DMatInvExpr instance.
Definition: DMatInvExpr.h:105
IntegralConstant< bool, B > BoolConstant
Generic wrapper for a compile time constant boolean value.The BoolConstant class template represents ...
Definition: IntegralConstant.h:100
bool canAlias(const T *alias) const noexcept
Returns whether the expression can alias with the given address alias.
Definition: DMatInvExpr.h:174
typename T::OppositeType OppositeType_
Alias declaration for nested OppositeType type definitions.The OppositeType_ alias declaration provid...
Definition: Aliases.h:243
#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:84
Expression object for dense matrix inversions.The DMatInvExpr class represents the compile time expre...
Definition: DMatInvExpr.h:93
ElementType_< MT > ElementType
Resulting element type.
Definition: DMatInvExpr.h:109
Header file for the IsRowMajorMatrix type trait.
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
const ResultType CompositeType
Data type for composite expression templates.
Definition: DMatInvExpr.h:113
Header file for the IntegralConstant class template.
typename T::TransposeType TransposeType_
Alias declaration for nested TransposeType type definitions.The TransposeType_ alias declaration prov...
Definition: Aliases.h:403
Header file for the IsUpper type trait.
Flag for the inversion of a lower unitriangular matrix.
Definition: InversionFlag.h:112
Header file for the IsHermitian type trait.
const DMatInvExpr< MT, SO > inv(const DenseMatrix< MT, SO > &dm)
Calculation of the inverse of the given dense matrix.
Definition: DMatInvExpr.h:401
BLAZE_ALWAYS_INLINE bool isSquare(const Matrix< MT, SO > &matrix) noexcept
Checks if the given matrix is a square matrix.
Definition: Matrix.h:609
ResultType_< MT > RT
Result type of the dense matrix expression.
Definition: DMatInvExpr.h:99
#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
Header file for the IsExpression type trait class.
CompositeType_< MT > CT
Composite type of the dense matrix expression.
Definition: DMatInvExpr.h:100
size_t rows() const noexcept
Returns the current number of rows of the matrix.
Definition: DMatInvExpr.h:142
Header file for the FunctionTrace class.