dla_porcond - tive-definite matrix
DOUBLE PRECISION FUNCTION DLA_PORCOND(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) CHARACTER*1 UPLO INTEGER N, LDA, LDAF, INFO, CMODE DOUBLE PRECISION A(LDA,*), AF(LDAF,*), WORK(*), C(*) INTEGER IWORK(*) DOUBLE PRECISION FUNCTION DLA_PORCOND_64(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) CHARACTER*1 UPLO INTEGER*8 N, LDA, LDAF, INFO, CMODE DOUBLE PRECISION A(LDA,*), AF(LDAF,*), WORK(*), C(*) INTEGER*8 IWORK(*) F95 INTERFACE REAL(8) FUNCTION LA_PORCOND(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) INTEGER :: N, LDA, LDAF, CMODE, INFO CHARACTER(LEN=1) :: UPLO INTEGER, DIMENSION(:) :: IWORK REAL(8), DIMENSION(:,:) :: A, AF REAL(8), DIMENSION(:) :: C, WORK REAL(8) FUNCTION LA_PORCOND_64(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) INTEGER(8) :: N, LDA, LDAF, CMODE, INFO CHARACTER(LEN=1) :: UPLO INTEGER(8), DIMENSION(:) :: IWORK REAL(8), DIMENSION(:,:) :: A, AF REAL(8), DIMENSION(:) :: C, WORK C INTERFACE #include <sunperf.h> double dla_porcond (char uplo, int n, double *a, int lda, double *af, int ldaf, int cmode, double *c, int *info); double dla_porcond_64 (char uplo, long n, double *a, long lda, double *af, long ldaf, long cmode, double *c, long *info);
Oracle Solaris Studio Performance Library dla_porcond(3P) NAME dla_porcond - estimate the Skeel condition number for a symmetric posi- tive-definite matrix SYNOPSIS DOUBLE PRECISION FUNCTION DLA_PORCOND(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) CHARACTER*1 UPLO INTEGER N, LDA, LDAF, INFO, CMODE DOUBLE PRECISION A(LDA,*), AF(LDAF,*), WORK(*), C(*) INTEGER IWORK(*) DOUBLE PRECISION FUNCTION DLA_PORCOND_64(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) CHARACTER*1 UPLO INTEGER*8 N, LDA, LDAF, INFO, CMODE DOUBLE PRECISION A(LDA,*), AF(LDAF,*), WORK(*), C(*) INTEGER*8 IWORK(*) F95 INTERFACE REAL(8) FUNCTION LA_PORCOND(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) INTEGER :: N, LDA, LDAF, CMODE, INFO CHARACTER(LEN=1) :: UPLO INTEGER, DIMENSION(:) :: IWORK REAL(8), DIMENSION(:,:) :: A, AF REAL(8), DIMENSION(:) :: C, WORK REAL(8) FUNCTION LA_PORCOND_64(UPLO, N, A, LDA, AF, LDAF, CMODE, C, INFO, WORK, IWORK) INTEGER(8) :: N, LDA, LDAF, CMODE, INFO CHARACTER(LEN=1) :: UPLO INTEGER(8), DIMENSION(:) :: IWORK REAL(8), DIMENSION(:,:) :: A, AF REAL(8), DIMENSION(:) :: C, WORK C INTERFACE #include <sunperf.h> double dla_porcond (char uplo, int n, double *a, int lda, double *af, int ldaf, int cmode, double *c, int *info); double dla_porcond_64 (char uplo, long n, double *a, long lda, double *af, long ldaf, long cmode, double *c, long *info); PURPOSE dla_porcond Estimates the Skeel condition number of op(A)*op2(C) where op2 is determined by CMODE as follows: CMODE = 1 op2(C) = C CMODE = 0 op2(C) = I CMODE = -1 op2(C) = inv(C) The Skeel condition number cond(A)=norminf(|inv(A)||A|) is computed by computing scaling factors R such that diag(R)*A*op2(C) is row equili- brated and computing the standard infinity-norm condition number. 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. A (input) A is DOUBLE PRECISION 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 DOUBLE PRECISION 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 DPOTRF. LDAF (input) LDAF is INTEGER The leading dimension of the array AF. LDAF >= max(1,N). CMODE (input) CMODE is INTEGER Determines op2(C) in the formula op(A)*op2(C) as follows: CMODE = 1 op2(C) = C CMODE = 0 op2(C) = I CMODE = -1 op2(C) = inv(C) C (input) C is DOUBLE PRECISION array, dimension (N) The vector C in the formula op(A)*op2(C). INFO (output) INFO is INTEGER = 0: Successful exit. i > 0: The ith argument is invalid. WORK (input) WORK is DOUBLE PRECISION array, dimension (3*N). Workspace. IWORK (input) IWORK is INTEGER array, dimension (N). Workspace. 7 Nov 2015 dla_porcond(3P)