1SSTEQR(1) LAPACK routine (version 3.2) SSTEQR(1)
2
3
4
6 SSTEQR - computes all eigenvalues and, optionally, eigenvectors of a
7 symmetric tridiagonal matrix using the implicit QL or QR method
8
10 SUBROUTINE SSTEQR( COMPZ, N, D, E, Z, LDZ, WORK, INFO )
11
12 CHARACTER COMPZ
13
14 INTEGER INFO, LDZ, N
15
16 REAL D( * ), E( * ), WORK( * ), Z( LDZ, * )
17
19 SSTEQR computes all eigenvalues and, optionally, eigenvectors of a sym‐
20 metric tridiagonal matrix using the implicit QL or QR method. The
21 eigenvectors of a full or band symmetric matrix can also be found if
22 SSYTRD or SSPTRD or SSBTRD has been used to reduce this matrix to
23 tridiagonal form.
24
26 COMPZ (input) CHARACTER*1
27 = 'N': Compute eigenvalues only.
28 = 'V': Compute eigenvalues and eigenvectors of the original
29 symmetric matrix. On entry, Z must contain the orthogonal
30 matrix used to reduce the original matrix to tridiagonal form.
31 = 'I': Compute eigenvalues and eigenvectors of the tridiagonal
32 matrix. Z is initialized to the identity matrix.
33
34 N (input) INTEGER
35 The order of the matrix. N >= 0.
36
37 D (input/output) REAL array, dimension (N)
38 On entry, the diagonal elements of the tridiagonal matrix. On
39 exit, if INFO = 0, the eigenvalues in ascending order.
40
41 E (input/output) REAL array, dimension (N-1)
42 On entry, the (n-1) subdiagonal elements of the tridiagonal
43 matrix. On exit, E has been destroyed.
44
45 Z (input/output) REAL array, dimension (LDZ, N)
46 On entry, if COMPZ = 'V', then Z contains the orthogonal
47 matrix used in the reduction to tridiagonal form. On exit, if
48 INFO = 0, then if COMPZ = 'V', Z contains the orthonormal
49 eigenvectors of the original symmetric matrix, and if COMPZ =
50 'I', Z contains the orthonormal eigenvectors of the symmetric
51 tridiagonal matrix. If COMPZ = 'N', then Z is not referenced.
52
53 LDZ (input) INTEGER
54 The leading dimension of the array Z. LDZ >= 1, and if eigen‐
55 vectors are desired, then LDZ >= max(1,N).
56
57 WORK (workspace) REAL array, dimension (max(1,2*N-2))
58 If COMPZ = 'N', then WORK is not referenced.
59
60 INFO (output) INTEGER
61 = 0: successful exit
62 < 0: if INFO = -i, the i-th argument had an illegal value
63 > 0: the algorithm has failed to find all the eigenvalues in a
64 total of 30*N iterations; if INFO = i, then i elements of E
65 have not converged to zero; on exit, D and E contain the ele‐
66 ments of a symmetric tridiagonal matrix which is orthogonally
67 similar to the original matrix.
68
69
70
71 LAPACK routine (version 3.2) November 2008 SSTEQR(1)