NAG Library Routine Document

F07AWF (ZGETRI)

+− Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Further Comments

+− 9 Example

9.1 Program Text

9.2 Program Data

9.3 Program Results

1 Purpose

F07AWF (ZGETRI) computes the inverse of a complex matrix

A

, where

A

has been factorized by F07ARF (ZGETRF).

2 Specification

SUBROUTINE F07AWF (

N, A, LDA, IPIV, WORK, LWORK, INFO)

INTEGER	N, LDA, IPIV(*), LWORK, INFO
COMPLEX (KIND=nag_wp)	A(LDA,*), WORK(max(1,LWORK))

The routine may be called by its LAPACK name zgetri.

3 Description

F07AWF (ZGETRI) is used to compute the inverse of a complex matrix

A

, the routine must be preceded by a call to F07ARF (ZGETRF), which computes the

L U

factorization of

A

A = P L U

. The inverse of

A

is computed by forming

U^{- 1}

and then solving the equation

X P L = U^{- 1}

for

X

4 References

Du Croz J J and Higham N J (1992) Stability of methods for matrix inversion IMA J. Numer. Anal. 12 1–19

5 Parameters

1: N – INTEGERInput: On entry: $n$ , the order of the matrix $A$ .
Constraint: $N \geq 0$ .
2: A(LDA, $*$ ) – COMPLEX (KIND=nag_wp) arrayInput/Output: Note: the second dimension of the array A must be at least $\max (1, N)$ .
On entry: the $L U$ factorization of $A$ , as returned by F07ARF (ZGETRF).

On exit: the factorization is overwritten by the $n$ by $n$ matrix $A^{- 1}$ .
3: LDA – INTEGERInput: On entry: the first dimension of the array A as declared in the (sub)program from which F07AWF (ZGETRI) is called.
Constraint: $LDA \geq \max (1, N)$ .
4: IPIV( $*$ ) – INTEGER arrayInput: Note: the dimension of the array IPIV must be at least $\max (1, N)$ .
On entry: the pivot indices, as returned by F07ARF (ZGETRF).
5: WORK( $\max (1, LWORK)$ ) – COMPLEX (KIND=nag_wp) arrayWorkspace: On exit: if $INFO = 0$ , $WORK (1)$ contains the minimum value of LWORK required for optimum performance.
6: LWORK – INTEGERInput: On entry: the dimension of the array WORK as declared in the (sub)program from which F07AWF (ZGETRI) is called, unless $LWORK = - 1$ , in which case a workspace query is assumed and the routine only calculates the optimal dimension of WORK (using the formula given below).

Suggested value: for optimum performance LWORK should be at least $N \times nb$ , where $nb$ is the block size.
Constraint: $LWORK \geq \max (1, N)$ or $LWORK = - 1$ .
7: INFO – INTEGEROutput: On exit: $INFO = 0$ unless the routine detects an error (see Section 6).

6 Error Indicators and Warnings

Errors or warnings detected by the routine:

$INFO < 0$: If $INFO = - i$ , the $i$ th parameter had an illegal value. An explanatory message is output, and execution of the program is terminated.

$INFO > 0$: If $INFO = i$ , the $i$ th diagonal element of the factor $U$ is zero, $U$ is singular, and the inverse of $A$ cannot be computed.

7 Accuracy

The computed inverse

X

satisfies a bound of the form:

|X A - I| \leq c (n) ε |X| P |L| |U|,

where

c (n)

is a modest linear function of

n

, and

ε

is the machine precision.

Note that a similar bound for

|A X - I|

cannot be guaranteed, although it is almost always satisfied. See Du Croz and Higham (1992).

8 Further Comments

The total number of real floating point operations is approximately

\frac{16}{3} n^{3}

The real analogue of this routine is F07AJF (DGETRI).

9 Example

This example computes the inverse of the matrix

A

, where

A = (\begin{array}{r} - 1.34 + 2.55 i & 0.28 + 3.17 i & - 6.39 - 2.20 i & 0.72 - 0.92 i \\ - 0.17 - 1.41 i & 3.31 - 0.15 i & - 0.15 + 1.34 i & 1.29 + 1.38 i \\ - 3.29 - 2.39 i & - 1.91 + 4.42 i & - 0.14 - 1.35 i & 1.72 + 1.35 i \\ 2.41 + 0.39 i & - 0.56 + 1.47 i & - 0.83 - 0.69 i & - 1.96 + 0.67 i \end{array}) .

Here

A

is nonsymmetric and must first be factorized by F07ARF (ZGETRF).

NAG Library Routine DocumentF07AWF (ZGETRI)