TESTING/LIN/sppt03.f(3) | Library Functions Manual | TESTING/LIN/sppt03.f(3) |
NAME
TESTING/LIN/sppt03.f
SYNOPSIS
Functions/Subroutines
subroutine sppt03 (uplo, n, a, ainv, work, ldwork, rwork,
rcond, resid)
SPPT03
Function/Subroutine Documentation
subroutine sppt03 (character uplo, integer n, real, dimension( * ) a, real, dimension( * ) ainv, real, dimension( ldwork, * ) work, integer ldwork, real, dimension( * ) rwork, real rcond, real resid)
SPPT03
Purpose:
SPPT03 computes the residual for a symmetric packed matrix times its inverse: norm( I - A*AINV ) / ( N * norm(A) * norm(AINV) * EPS ), where EPS is the machine epsilon.
Parameters
UPLO
UPLO is CHARACTER*1 Specifies whether the upper or lower triangular part of the symmetric matrix A is stored: = 'U': Upper triangular = 'L': Lower triangular
N
N is INTEGER The number of rows and columns of the matrix A. N >= 0.
A
A is REAL array, dimension (N*(N+1)/2) The original symmetric matrix A, stored as a packed triangular matrix.
AINV
AINV is REAL array, dimension (N*(N+1)/2) The (symmetric) inverse of the matrix A, stored as a packed triangular matrix.
WORK
WORK is REAL array, dimension (LDWORK,N)
LDWORK
LDWORK is INTEGER The leading dimension of the array WORK. LDWORK >= max(1,N).
RWORK
RWORK is REAL array, dimension (N)
RCOND
RCOND is REAL The reciprocal of the condition number of A, computed as ( 1/norm(A) ) / norm(AINV).
RESID
RESID is REAL norm(I - A*AINV) / ( N * norm(A) * norm(AINV) * EPS )
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.
Definition at line 108 of file sppt03.f.
Author
Generated automatically by Doxygen for LAPACK from the source code.
Version 3.12.0 | LAPACK |