1DLAS2(1) LAPACK auxiliary routine (version 3.1) DLAS2(1)
2
3
4
6 DLAS2 - the singular values of the 2-by-2 matrix [ F G ] [ 0 H ]
7
9 SUBROUTINE DLAS2( F, G, H, SSMIN, SSMAX )
10
11 DOUBLE PRECISION F, G, H, SSMAX, SSMIN
12
14 DLAS2 computes the singular values of the 2-by-2 matrix
15 [ F G ]
16 [ 0 H ]. On return, SSMIN is the smaller singular value and
17 SSMAX is the larger singular value.
18
19
21 F (input) DOUBLE PRECISION
22 The (1,1) element of the 2-by-2 matrix.
23
24 G (input) DOUBLE PRECISION
25 The (1,2) element of the 2-by-2 matrix.
26
27 H (input) DOUBLE PRECISION
28 The (2,2) element of the 2-by-2 matrix.
29
30 SSMIN (output) DOUBLE PRECISION
31 The smaller singular value.
32
33 SSMAX (output) DOUBLE PRECISION
34 The larger singular value.
35
37 Barring over/underflow, all output quantities are correct to within a
38 few units in the last place (ulps), even in the absence of a guard
39 digit in addition/subtraction.
40
41 In IEEE arithmetic, the code works correctly if one matrix element is
42 infinite.
43
44 Overflow will not occur unless the largest singular value itself over‐
45 flows, or is within a few ulps of overflow. (On machines with partial
46 overflow, like the Cray, overflow may occur if the largest singular
47 value is within a factor of 2 of overflow.)
48
49 Underflow is harmless if underflow is gradual. Otherwise, results may
50 correspond to a matrix modified by perturbations of size near the
51 underflow threshold.
52
53
54
55
56 LAPACK auxiliary routine (versionNo3v.e1m)ber 2006 DLAS2(1)