1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
|
.TH DLAED5 l "15 June 2000" "LAPACK version 3.0" ")"
.SH NAME
DLAED5 - subroutine computes the I-th eigenvalue of a symmetric rank-one modification of a 2-by-2 diagonal matrix diag( D ) + RHO * Z * transpose(Z)
.SH SYNOPSIS
.TP 19
SUBROUTINE DLAED5(
I, D, Z, DELTA, RHO, DLAM )
.TP 19
.ti +4
INTEGER
I
.TP 19
.ti +4
DOUBLE
PRECISION DLAM, RHO
.TP 19
.ti +4
DOUBLE
PRECISION D( 2 ), DELTA( 2 ), Z( 2 )
.SH PURPOSE
This subroutine computes the I-th eigenvalue of a symmetric rank-one modification of a 2-by-2 diagonal matrix diag( D ) + RHO * Z * transpose(Z) .
The diagonal elements in the array D are assumed to satisfy
D(i) < D(j) for i < j .
.br
We also assume RHO > 0 and that the Euclidean norm of the vector
Z is one.
.br
.SH ARGUMENTS
.TP 7
I (input) INTEGER
The index of the eigenvalue to be computed. I = 1 or I = 2.
.TP 7
D (input) DOUBLE PRECISION array, dimension (2)
The original eigenvalues. We assume D(1) < D(2).
.TP 7
Z (input) DOUBLE PRECISION array, dimension (2)
The components of the updating vector.
.TP 7
DELTA (output) DOUBLE PRECISION array, dimension (2)
The vector DELTA contains the information necessary
to construct the eigenvectors.
.TP 7
RHO (input) DOUBLE PRECISION
The scalar in the symmetric updating formula.
.TP 7
DLAM (output) DOUBLE PRECISION
The computed lambda_I, the I-th updated eigenvalue.
.SH FURTHER DETAILS
Based on contributions by
.br
Ren-Cang Li, Computer Science Division, University of California
at Berkeley, USA
.br
|
