Blaze 3.9
sytri.h
Go to the documentation of this file.
1//=================================================================================================
33//=================================================================================================
34
35#ifndef _BLAZE_MATH_LAPACK_SYTRI_H_
36#define _BLAZE_MATH_LAPACK_SYTRI_H_
37
38
39//*************************************************************************************************
40// Includes
41//*************************************************************************************************
42
43#include <memory>
44#include <blaze/math/Aliases.h>
54#include <blaze/util/Assert.h>
56
57
58namespace blaze {
59
60//=================================================================================================
61//
62// LAPACK LDLT-BASED INVERSION FUNCTIONS (SYTRI)
63//
64//=================================================================================================
65
66//*************************************************************************************************
69template< typename MT, bool SO >
70void sytri( DenseMatrix<MT,SO>& A, char uplo, const blas_int_t* ipiv );
72//*************************************************************************************************
73
74
75//*************************************************************************************************
113template< typename MT // Type of the dense matrix
114 , bool SO > // Storage order of the dense matrix
115inline void sytri( DenseMatrix<MT,SO>& A, char uplo, const blas_int_t* ipiv )
116{
122
123 using ET = ElementType_t<MT>;
124
125 if( !isSquare( *A ) ) {
126 BLAZE_THROW_INVALID_ARGUMENT( "Invalid non-square matrix provided" );
127 }
128
129 if( uplo != 'L' && uplo != 'U' ) {
130 BLAZE_THROW_INVALID_ARGUMENT( "Invalid uplo argument provided" );
131 }
132
133 blas_int_t n ( numeric_cast<blas_int_t>( (*A).columns() ) );
134 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
135 blas_int_t info( 0 );
136
137 if( n == 0 ) {
138 return;
139 }
140
141 if( IsRowMajorMatrix_v<MT> ) {
142 ( uplo == 'L' )?( uplo = 'U' ):( uplo = 'L' );
143 }
144
145 const std::unique_ptr<ET[]> work( new ET[n] );
146
147 sytri( uplo, n, (*A).data(), lda, ipiv, work.get(), &info );
148
149 BLAZE_INTERNAL_ASSERT( info >= 0, "Invalid argument for matrix inversion" );
150
151 if( info > 0 ) {
152 BLAZE_THROW_LAPACK_ERROR( "Inversion of singular matrix failed" );
153 }
154}
155//*************************************************************************************************
156
157} // namespace blaze
158
159#endif
Constraint on the data type.
Header file for auxiliary alias declarations.
typename T::ElementType ElementType_t
Alias declaration for nested ElementType type definitions.
Definition: Aliases.h:190
Header file for run time assertion macros.
Constraint on the data type.
Constraint on the data type.
Header file for the IsRowMajorMatrix type trait.
Constraint on the data type.
Cast operators for numeric types.
Header file for the CLAPACK sytri wrapper functions.
Base class for dense matrices.
Definition: DenseMatrix.h:82
Constraint on the data type.
Header file for the DenseMatrix base class.
void sytri(DenseMatrix< MT, SO > &A, char uplo, const blas_int_t *ipiv)
LAPACK kernel for the inversion of the given dense symmetric indefinite matrix.
Definition: sytri.h:115
#define BLAZE_CONSTRAINT_MUST_BE_BLAS_COMPATIBLE_TYPE(T)
Constraint on the data type.
Definition: BLASCompatible.h:61
#define BLAZE_CONSTRAINT_MUST_BE_CONTIGUOUS_TYPE(T)
Constraint on the data type.
Definition: Contiguous.h:61
#define BLAZE_CONSTRAINT_MUST_NOT_BE_COMPUTATION_TYPE(T)
Constraint on the data type.
Definition: Computation.h:81
#define BLAZE_CONSTRAINT_MUST_NOT_BE_ADAPTOR_TYPE(T)
Constraint on the data type.
Definition: Adaptor.h:81
#define BLAZE_CONSTRAINT_MUST_HAVE_MUTABLE_DATA_ACCESS(T)
Constraint on the data type.
Definition: MutableDataAccess.h:61
int32_t blas_int_t
Signed integer type used in the BLAS/LAPACK wrapper functions.
Definition: Types.h:64
#define BLAZE_THROW_LAPACK_ERROR(MESSAGE)
Macro for the emission of an exception on detection of a LAPACK error.
Definition: Exception.h:146
bool isSquare(const Matrix< MT, SO > &matrix) noexcept
Checks if the given matrix is a square matrix.
Definition: Matrix.h:1383
#define BLAZE_INTERNAL_ASSERT(expr, msg)
Run time assertion macro for internal checks.
Definition: Assert.h:101
#define BLAZE_THROW_INVALID_ARGUMENT(MESSAGE)
Macro for the emission of a std::invalid_argument exception.
Definition: Exception.h:235
Header file for the exception macros of the math module.