dsyev - compute all eigenvalues and, optionally, eigenvectors of a real symmetric matrix A
SUBROUTINE DSYEV( JOBZ, UPLO, N, A, LDA, W, WORK, LDWORK, INFO) CHARACTER * 1 JOBZ, UPLO INTEGER N, LDA, LDWORK, INFO DOUBLE PRECISION A(LDA,*), W(*), WORK(*)
SUBROUTINE DSYEV_64( JOBZ, UPLO, N, A, LDA, W, WORK, LDWORK, INFO) CHARACTER * 1 JOBZ, UPLO INTEGER*8 N, LDA, LDWORK, INFO DOUBLE PRECISION A(LDA,*), W(*), WORK(*)
SUBROUTINE SYEV( JOBZ, UPLO, [N], A, [LDA], W, [WORK], [LDWORK], * [INFO]) CHARACTER(LEN=1) :: JOBZ, UPLO INTEGER :: N, LDA, LDWORK, INFO REAL(8), DIMENSION(:) :: W, WORK REAL(8), DIMENSION(:,:) :: A
SUBROUTINE SYEV_64( JOBZ, UPLO, [N], A, [LDA], W, [WORK], [LDWORK], * [INFO]) CHARACTER(LEN=1) :: JOBZ, UPLO INTEGER(8) :: N, LDA, LDWORK, INFO REAL(8), DIMENSION(:) :: W, WORK REAL(8), DIMENSION(:,:) :: A
#include <sunperf.h>
void dsyev(char jobz, char uplo, int n, double *a, int lda, double *w, int *info);
void dsyev_64(char jobz, char uplo, long n, double *a, long lda, double *w, long *info);
dsyev computes all eigenvalues and, optionally, eigenvectors of a real symmetric matrix A.
= 'N': Compute eigenvalues only;
= 'V': Compute eigenvalues and eigenvectors.
= 'U': Upper triangle of A is stored;
= 'L': Lower triangle of A is stored.
WORK(1)
returns the optimal LDWORK.
If LDWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of the WORK array, returns this value as the first entry of the WORK array, and no error message related to LDWORK is issued by XERBLA.
= 0: successful exit
< 0: if INFO = -i, the i-th argument had an illegal value
> 0: if INFO = i, the algorithm failed to converge; i off-diagonal elements of an intermediate tridiagonal form did not converge to zero.