1SGETRI(1)                LAPACK routine (version 3.1)                SGETRI(1)
2
3
4

NAME

6       SGETRI - the inverse of a matrix using the LU factorization computed by
7       SGETRF
8

SYNOPSIS

10       SUBROUTINE SGETRI( N, A, LDA, IPIV, WORK, LWORK, INFO )
11
12           INTEGER        INFO, LDA, LWORK, N
13
14           INTEGER        IPIV( * )
15
16           REAL           A( LDA, * ), WORK( * )
17

PURPOSE

19       SGETRI computes the inverse of a matrix using the LU factorization com‐
20       puted by SGETRF.
21
22       This  method  inverts  U and then computes inv(A) by solving the system
23       inv(A)*L = inv(U) for inv(A).
24
25

ARGUMENTS

27       N       (input) INTEGER
28               The order of the matrix A.  N >= 0.
29
30       A       (input/output) REAL array, dimension (LDA,N)
31               On entry, the factors L and U from the factorization A =  P*L*U
32               as  computed  by  SGETRF.  On exit, if INFO = 0, the inverse of
33               the original matrix A.
34
35       LDA     (input) INTEGER
36               The leading dimension of the array A.  LDA >= max(1,N).
37
38       IPIV    (input) INTEGER array, dimension (N)
39               The pivot indices from SGETRF; for 1<=i<=N, row i of the matrix
40               was interchanged with row IPIV(i).
41
42       WORK    (workspace/output) REAL array, dimension (MAX(1,LWORK))
43               On exit, if INFO=0, then WORK(1) returns the optimal LWORK.
44
45       LWORK   (input) INTEGER
46               The dimension of the array WORK.  LWORK >= max(1,N).  For opti‐
47               mal performance LWORK >= N*NB, where NB is the  optimal  block‐
48               size returned by ILAENV.
49
50               If  LWORK  = -1, then a workspace query is assumed; the routine
51               only calculates the optimal size of  the  WORK  array,  returns
52               this  value  as the first entry of the WORK array, and no error
53               message related to LWORK is issued by XERBLA.
54
55       INFO    (output) INTEGER
56               = 0:  successful exit
57               < 0:  if INFO = -i, the i-th argument had an illegal value
58               > 0:  if INFO = i, U(i,i) is exactly zero; the matrix is singu‐
59               lar and its inverse could not be computed.
60
61
62
63 LAPACK routine (version 3.1)    November 2006                       SGETRI(1)
Impressum