cla_porpvgrw - compute the reciprocal pivot growth factor norm(A)/norm(U) for a symmetric or Hermitian positive-definite matrix
REAL FUNCTION CLA_PORPVGRW(UPLO, NCOLS, A, LDA, AF, LDAF, WORK) CHARACTER*1 UPLO INTEGER NCOLS, LDA, LDAF COMPLEX A(LDA,*), AF(LDAF,*) REAL WORK(*) REAL FUNCTION CLA_PORPVGRW_64(UPLO, NCOLS, A, LDA, AF, LDAF, WORK) CHARACTER*1 UPLO INTEGER*8 NCOLS, LDA, LDAF COMPLEX A(LDA,*), AF(LDAF,*) REAL WORK(*) F95 INTERFACE REAL FUNCTION LA_PORPVGRW(UPLO, NCOLS, A, LDA, AF, LDAF, WORK) INTEGER :: NCOLS, LDA, LDAF CHARACTER(LEN=1) :: UPLO COMPLEX, DIMENSION(:) :: WORK COMPLEX, DIMENSION(:,:) :: A, AF REAL FUNCTION LA_PORPVGRW_64(UPLO, NCOLS, A, LDA, AF, LDAF, WORK) INTEGER(8) :: NCOLS, LDA, LDAF CHARACTER(LEN=1) :: UPLO COMPLEX, DIMENSION(:) :: WORK COMPLEX, DIMENSION(:,:) :: A, AF C INTERFACE #include <sunperf.h> float cla_porpvgrw (char uplo, int ncols, floatcomplex *a, int info, int lda, floatcomplex *af, int ldaf); float cla_porpvgrw_64 (char uplo, long ncols, floatcomplex *a, long info, long lda, floatcomplex * af, long ldaf);
Oracle Solaris Studio Performance Library cla_porpvgrw(3P)
NAME
cla_porpvgrw - compute the reciprocal pivot growth factor
norm(A)/norm(U) for a symmetric or Hermitian positive-definite matrix
SYNOPSIS
REAL FUNCTION CLA_PORPVGRW(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
CHARACTER*1 UPLO
INTEGER NCOLS, LDA, LDAF
COMPLEX A(LDA,*), AF(LDAF,*)
REAL WORK(*)
REAL FUNCTION CLA_PORPVGRW_64(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
CHARACTER*1 UPLO
INTEGER*8 NCOLS, LDA, LDAF
COMPLEX A(LDA,*), AF(LDAF,*)
REAL WORK(*)
F95 INTERFACE
REAL FUNCTION LA_PORPVGRW(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
INTEGER :: NCOLS, LDA, LDAF
CHARACTER(LEN=1) :: UPLO
COMPLEX, DIMENSION(:) :: WORK
COMPLEX, DIMENSION(:,:) :: A, AF
REAL FUNCTION LA_PORPVGRW_64(UPLO, NCOLS, A, LDA, AF, LDAF, WORK)
INTEGER(8) :: NCOLS, LDA, LDAF
CHARACTER(LEN=1) :: UPLO
COMPLEX, DIMENSION(:) :: WORK
COMPLEX, DIMENSION(:,:) :: A, AF
C INTERFACE
#include <sunperf.h>
float cla_porpvgrw (char uplo, int ncols, floatcomplex *a, int info,
int lda, floatcomplex *af, int ldaf);
float cla_porpvgrw_64 (char uplo, long ncols, floatcomplex *a, long
info, long lda, floatcomplex * af, long ldaf);
PURPOSE
cla_porpvgrw 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.
NCOLS (input)
NCOLS is INTEGER
The number of columns of the matrix A. NCOLS >= 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 triangular factor U or L from the Cholesky factorization
A=U**T*U or A=L*L**T, as computed by CPOTRF.
LDAF (input)
LDAF is INTEGER
The leading dimension of the array AF.
LDAF >= max(1,N).
WORK (input)
WORK is COMPLEX array, dimension (2*N)
7 Nov 2015 cla_porpvgrw(3P)