cla_herpvgrw - compute the reciprocal pivot growth factor using the "max absolute element" norm
REAL FUNCTION CLA_HERPVGRW(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) CHARACTER*1 UPLO INTEGER N, INFO, LDA, LDAF INTEGER IPIV(*) COMPLEX A(LDA,*), AF(LDAF,*) REAL WORK(*) REAL FUNCTION CLA_HERPVGRW_64(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) CHARACTER*1 UPLO INTEGER*8 N, INFO, LDA, LDAF INTEGER*8 IPIV(*) COMPLEX A(LDA,*), AF(LDAF,*) REAL WORK(*) F95 INTERFACE REAL FUNCTION LA_HERPVGRW(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) INTEGER :: N, INFO, LDA, LDAF CHARACTER(LEN=1) :: UPLO INTEGER, DIMENSION(:) :: IPIV COMPLEX, DIMENSION(:) :: WORK COMPLEX, DIMENSION(:,:) :: A, AF REAL FUNCTION LA_HERPVGRW_64(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) INTEGER(8) :: N, INFO, LDA, LDAF CHARACTER(LEN=1) :: UPLO INTEGER(8), DIMENSION(:) :: IPIV COMPLEX, DIMENSION(:) :: WORK COMPLEX, DIMENSION(:,:) :: A, AF C INTERFACE #include <sunperf.h> float cla_herpvgrw (char uplo, int n, int info, floatcomplex *a, int lda, floatcomplex *af, int ldaf, int *ipiv); float cla_herpvgrw_64 (char uplo, long n, long info, floatcomplex *a, long lda, floatcomplex *af, long ldaf, long *ipiv);
Oracle Solaris Studio Performance Library cla_herpvgrw(3P) NAME cla_herpvgrw - compute the reciprocal pivot growth factor using the "max absolute element" norm SYNOPSIS REAL FUNCTION CLA_HERPVGRW(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) CHARACTER*1 UPLO INTEGER N, INFO, LDA, LDAF INTEGER IPIV(*) COMPLEX A(LDA,*), AF(LDAF,*) REAL WORK(*) REAL FUNCTION CLA_HERPVGRW_64(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) CHARACTER*1 UPLO INTEGER*8 N, INFO, LDA, LDAF INTEGER*8 IPIV(*) COMPLEX A(LDA,*), AF(LDAF,*) REAL WORK(*) F95 INTERFACE REAL FUNCTION LA_HERPVGRW(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) INTEGER :: N, INFO, LDA, LDAF CHARACTER(LEN=1) :: UPLO INTEGER, DIMENSION(:) :: IPIV COMPLEX, DIMENSION(:) :: WORK COMPLEX, DIMENSION(:,:) :: A, AF REAL FUNCTION LA_HERPVGRW_64(UPLO, N, INFO, A, LDA, AF, LDAF, IPIV, WORK) INTEGER(8) :: N, INFO, LDA, LDAF CHARACTER(LEN=1) :: UPLO INTEGER(8), DIMENSION(:) :: IPIV COMPLEX, DIMENSION(:) :: WORK COMPLEX, DIMENSION(:,:) :: A, AF C INTERFACE #include <sunperf.h> float cla_herpvgrw (char uplo, int n, int info, floatcomplex *a, int lda, floatcomplex *af, int ldaf, int *ipiv); float cla_herpvgrw_64 (char uplo, long n, long info, floatcomplex *a, long lda, floatcomplex *af, long ldaf, long *ipiv); PURPOSE cla_herpvgrw computes the reciprocal pivot growth factor norm(A)/norm(U). The "max absolute element" norm is used. If this is much less than 1, the stability of the LU factorization of the (equili- brated) matrix A could be poor. This also means that the solution X, estimated condition numbers, and error bounds could be unreliable. ARGUMENTS UPLO (input) UPLO is CHARACTER*1 = 'U': Upper triangle of A is stored; = 'L': Lower triangle of A is stored. N (input) N is INTEGER The number of linear equations, i.e., the order of the matrix A. N >= 0. INFO (input) INFO is INTEGER The value of INFO returned from SSYTRF, i.e., the pivot in column INFO is exactly 0. A (input) A is COMPLEX array, dimension (LDA,N) On entry, the N-by-N matrix A. LDA (input) LDA is INTEGER The leading dimension of the array A. LDA >= max(1,N). AF (input) AF is COMPLEX array, dimension (LDAF,N) The block diagonal matrix D and the multipliers used to obtain the factor U or L as computed by CHETRF. LDAF (input) LDAF is INTEGER The leading dimension of the array AF. LDAF >= max(1,N). IPIV (input) IPIV is INTEGER array, dimension (N) Details of the interchanges and the block structure of D as determined by CHETRF. WORK (input) WORK is COMPLEX array, dimension (2*N) 7 Nov 2015 cla_herpvgrw(3P)