XLPack 6.1
Excel VBA Numerical Library Reference Manual
Loading...
Searching...
No Matches

◆ Ztpcon()

Sub Ztpcon ( Norm As  String,
Uplo As  String,
Diag As  String,
N As  Long,
Ap() As  Complex,
RCond As  Double,
Info As  Long 
)

Condition number of a complex triangular matrix in packed form

Purpose
This routine estimates the reciprocal of the condition number of a packed triangular matrix A, in either the 1-norm or the infinity-norm.

The norm of A is computed and an estimate is obtained for norm(inv(A)), then the reciprocal of the condition number is computed as
RCond = 1 / (norm(A) * norm(inv(A))).
Parameters
[in]NormSpecifies whether the 1-norm condition number or the infinity-norm condition number is required:
= "1" or "O": 1-norm.
= "I": Infinity-norm.
[in]Uplo= "U": A is upper triangular.
= "L": A is lower triangular.
[in]Diag= "N": A is non-unit triangular.
= "U": A is unit triangular. (Diagonal elements of Ap() are not referenced and are assumed to be 1)
[in]NOrder of the matrix A. (N >= 0) (If N = 0, returns RCond = 1)
[in]Ap()Array Ap(LAp - 1) (LAp >= N(N + 1)/2)
N x N triangular matrix A in packed symmetric matrix form. Upper or lower part is to be stored in accordance with Uplo.
[out]RCondThe reciprocal of the condition number of the matrix A, computed as RCond = 1/(norm(A) * norm(inv(A))).
[out]Info= 0: Successful exit.
= -1: The argument Norm had an illegal value. (Norm <> "1", "O" nor "I")
= -2: The argument Uplo had an illegal value. (Uplo <> "U" nor "L")
= -3: The argument Diag had an illegal value. (Diag <> "N" nor "U")
= -4: The argument N had an illegal value. (N < 0)
= -5: The argument Ap() is invalid.
Reference
LAPACK
Example Program
See example of Ztptrs.