ssbevx - compute selected eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A
SUBROUTINE SSBEVX( JOBZ, RANGE, UPLO, N, NDIAG, A, LDA, Q, LDQ, VL, * VU, IL, IU, ABTOL, NFOUND, W, Z, LDZ, WORK, IWORK2, IFAIL, INFO) CHARACTER * 1 JOBZ, RANGE, UPLO INTEGER N, NDIAG, LDA, LDQ, IL, IU, NFOUND, LDZ, INFO INTEGER IWORK2(*), IFAIL(*) REAL VL, VU, ABTOL REAL A(LDA,*), Q(LDQ,*), W(*), Z(LDZ,*), WORK(*)
SUBROUTINE SSBEVX_64( JOBZ, RANGE, UPLO, N, NDIAG, A, LDA, Q, LDQ, * VL, VU, IL, IU, ABTOL, NFOUND, W, Z, LDZ, WORK, IWORK2, IFAIL, * INFO) CHARACTER * 1 JOBZ, RANGE, UPLO INTEGER*8 N, NDIAG, LDA, LDQ, IL, IU, NFOUND, LDZ, INFO INTEGER*8 IWORK2(*), IFAIL(*) REAL VL, VU, ABTOL REAL A(LDA,*), Q(LDQ,*), W(*), Z(LDZ,*), WORK(*)
SUBROUTINE SBEVX( JOBZ, RANGE, UPLO, [N], NDIAG, A, [LDA], Q, [LDQ], * VL, VU, IL, IU, ABTOL, NFOUND, W, Z, [LDZ], [WORK], [IWORK2], * IFAIL, [INFO]) CHARACTER(LEN=1) :: JOBZ, RANGE, UPLO INTEGER :: N, NDIAG, LDA, LDQ, IL, IU, NFOUND, LDZ, INFO INTEGER, DIMENSION(:) :: IWORK2, IFAIL REAL :: VL, VU, ABTOL REAL, DIMENSION(:) :: W, WORK REAL, DIMENSION(:,:) :: A, Q, Z
SUBROUTINE SBEVX_64( JOBZ, RANGE, UPLO, [N], NDIAG, A, [LDA], Q, * [LDQ], VL, VU, IL, IU, ABTOL, NFOUND, W, Z, [LDZ], [WORK], * [IWORK2], IFAIL, [INFO]) CHARACTER(LEN=1) :: JOBZ, RANGE, UPLO INTEGER(8) :: N, NDIAG, LDA, LDQ, IL, IU, NFOUND, LDZ, INFO INTEGER(8), DIMENSION(:) :: IWORK2, IFAIL REAL :: VL, VU, ABTOL REAL, DIMENSION(:) :: W, WORK REAL, DIMENSION(:,:) :: A, Q, Z
#include <sunperf.h>
void ssbevx(char jobz, char range, char uplo, int n, int ndiag, float *a, int lda, float *q, int ldq, float vl, float vu, int il, int iu, float abtol, int *nfound, float *w, float *z, int ldz, int *ifail, int *info);
void ssbevx_64(char jobz, char range, char uplo, long n, long ndiag, float *a, long lda, float *q, long ldq, float vl, float vu, long il, long iu, float abtol, long *nfound, float *w, float *z, long ldz, long *ifail, long *info);
ssbevx computes selected eigenvalues and, optionally, eigenvectors of a real symmetric band matrix A. Eigenvalues and eigenvectors can be selected by specifying either a range of values or a range of indices for the desired eigenvalues.
= 'N': Compute eigenvalues only;
= 'V': Compute eigenvalues and eigenvectors.
= 'A': all eigenvalues will be found;
= 'V': all eigenvalues in the half-open interval (VL,VU] will be found; = 'I': the IL-th through IU-th eigenvalues will be found.
= 'U': Upper triangle of A is stored;
= 'L': Lower triangle of A is stored.
A(kd+1+i-j,j)
= A(i,j)
for max(1,j-kd)
< =i < =j;
if UPLO = 'L', A(1+i-j,j)
= A(i,j)
for j < =i < =min(n,j+kd).
On exit, A is overwritten by values generated during the reduction to tridiagonal form. If UPLO = 'U', the first superdiagonal and the diagonal of the tridiagonal matrix T are returned in rows NDIAG and NDIAG+1 of A, and if UPLO = 'L', the diagonal and first subdiagonal of T are returned in the first two rows of A.
ABTOL + EPS * max( |a|,|b| ) ,
where EPS is the machine precision. If ABTOL is less than or equal to zero, then EPS*|T| will be used in its place, where |T| is the 1-norm of the tridiagonal matrix obtained by reducing A to tridiagonal form.
Eigenvalues will be computed most accurately when ABTOL is set to twice the underflow threshold 2*SLAMCH('S'), not zero. If this routine returns with INFO >0, indicating that some eigenvectors did not converge, try setting ABTOL to 2*SLAMCH('S').
See ``Computing Small Singular Values of Bidiagonal Matrices with Guaranteed High Relative Accuracy,'' by Demmel and Kahan, LAPACK Working Note #3.
max(1,NFOUND)
columns are
supplied in the array Z; if RANGE = 'V', the exact value of NFOUND
is not known in advance and an upper bound must be used.
dimension(7*N)
= 0: successful exit.
< 0: if INFO = -i, the i-th argument had an illegal value.
> 0: if INFO = i, then i eigenvectors failed to converge. Their indices are stored in array IFAIL.