Blaze 3.9
gesv.h
Go to the documentation of this file.
1//=================================================================================================
33//=================================================================================================
34
35#ifndef _BLAZE_MATH_LAPACK_GESV_H_
36#define _BLAZE_MATH_LAPACK_GESV_H_
37
38
39//*************************************************************************************************
40// Includes
41//*************************************************************************************************
42
43#include <blaze/math/Aliases.h>
53#include <blaze/util/Assert.h>
56
57
58namespace blaze {
59
60//=================================================================================================
61//
62// LAPACK GENERAL LINEAR SYSTEM FUNCTIONS (GESV)
63//
64//=================================================================================================
65
66//*************************************************************************************************
69template< typename MT, bool SO, typename VT, bool TF >
70void gesv( DenseMatrix<MT,SO>& A, DenseVector<VT,TF>& b, blas_int_t* ipiv );
71
72template< typename MT1, bool SO1, typename MT2, bool SO2 >
73void gesv( DenseMatrix<MT1,SO1>& A, DenseMatrix<MT2,SO2>& B, blas_int_t* ipiv );
75//*************************************************************************************************
76
77
78//*************************************************************************************************
132template< typename MT // Type of the system matrix
133 , bool SO // Storage order of the system matrix
134 , typename VT // Type of the right-hand side vector
135 , bool TF > // Transpose flag of the right-hand side vector
137{
143
148
149 if( !isSquare( *A ) ) {
150 BLAZE_THROW_INVALID_ARGUMENT( "Invalid non-square matrix provided" );
151 }
152
153 if( (*b).size() != (*A).rows() ) {
154 BLAZE_THROW_INVALID_ARGUMENT( "Invalid right-hand side vector provided" );
155 }
156
157 blas_int_t n ( numeric_cast<blas_int_t>( (*A).rows() ) );
158 blas_int_t nrhs( 1 );
159 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
160 blas_int_t ldb ( numeric_cast<blas_int_t>( (*b).size() ) );
161 blas_int_t info( 0 );
162
163 if( n == 0 ) {
164 return;
165 }
166
167 gesv( n, nrhs, (*A).data(), lda, ipiv, (*b).data(), ldb, &info );
168
169 BLAZE_INTERNAL_ASSERT( info >= 0, "Invalid function argument" );
170
171 if( info > 0 ) {
172 BLAZE_THROW_LAPACK_ERROR( "Inversion of singular matrix failed" );
173 }
174}
175//*************************************************************************************************
176
177
178//*************************************************************************************************
235template< typename MT1 // Type of the system matrix
236 , bool SO1 // Storage order of the system matrix
237 , typename MT2 // Type of the right-hand side matrix
238 , bool SO2 > // Storage order of the right-hand side matrix
240{
246
252
253 if( !isSquare( *A ) ) {
254 BLAZE_THROW_INVALID_ARGUMENT( "Invalid non-square matrix provided" );
255 }
256
257 blas_int_t n ( numeric_cast<blas_int_t>( (*A).rows() ) );
258 blas_int_t mrhs( numeric_cast<blas_int_t>( SO2 ? (*B).rows() : (*B).columns() ) );
259 blas_int_t nrhs( numeric_cast<blas_int_t>( SO2 ? (*B).columns() : (*B).rows() ) );
260 blas_int_t lda ( numeric_cast<blas_int_t>( (*A).spacing() ) );
261 blas_int_t ldb ( numeric_cast<blas_int_t>( (*B).spacing() ) );
262 blas_int_t info( 0 );
263
264 if( n != mrhs ) {
265 BLAZE_THROW_INVALID_ARGUMENT( "Invalid right-hand side matrix provided" );
266 }
267
268 if( n == 0 ) {
269 return;
270 }
271
272 gesv( n, nrhs, (*A).data(), lda, ipiv, (*B).data(), ldb, &info );
273
274 BLAZE_INTERNAL_ASSERT( info >= 0, "Invalid function argument" );
275
276 if( info > 0 ) {
277 BLAZE_THROW_LAPACK_ERROR( "Inversion of singular matrix failed" );
278 }
279}
280//*************************************************************************************************
281
282} // namespace blaze
283
284#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.
Constraint on the data type.
Cast operators for numeric types.
Data type constraint.
Header file for the CLAPACK gesv wrapper functions.
Base class for dense matrices.
Definition: DenseMatrix.h:82
Base class for N-dimensional dense vectors.
Definition: DenseVector.h:77
Constraint on the data type.
Header file for the DenseMatrix base class.
Header file for the DenseVector base class.
#define BLAZE_CONSTRAINT_MUST_BE_SAME_TYPE(A, B)
Data type constraint.
Definition: SameType.h:71
void gesv(DenseMatrix< MT1, SO1 > &A, DenseMatrix< MT2, SO2 > &B, blas_int_t *ipiv)
LAPACK kernel for solving a general linear system of equations ( ).
Definition: gesv.h:239
#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.