1DLASQ5(1)LAPACK routine (version 3.2) DLASQ5(1)
2
3
4
6 DLASQ5 - computes one dqds transform in ping-pong form, one version for
7 IEEE machines another for non IEEE machines
8
10 SUBROUTINE DLASQ5( I0, N0, Z, PP, TAU, DMIN, DMIN1, DMIN2, DN, DNM1,
11 DNM2, IEEE )
12
13 LOGICAL IEEE
14
15 INTEGER I0, N0, PP
16
17 DOUBLE PRECISION DMIN, DMIN1, DMIN2, DN, DNM1, DNM2, TAU
18
19 DOUBLE PRECISION Z( * )
20
22 DLASQ5 computes one dqds transform in ping-pong form, one version for
23 IEEE machines another for non IEEE machines.
24
26 I0 (input) INTEGER
27 First index.
28
29 N0 (input) INTEGER
30 Last index.
31
32 Z (input) DOUBLE PRECISION array, dimension ( 4*N )
33 Z holds the qd array. EMIN is stored in Z(4*N0) to avoid an extra
34 argument.
35
36 PP (input) INTEGER
37 PP=0 for ping, PP=1 for pong.
38
39 TAU (input) DOUBLE PRECISION
40 This is the shift.
41
42 DMIN (output) DOUBLE PRECISION
43 Minimum value of d. DMIN1 (output) DOUBLE PRECISION Minimum
44 value of d, excluding D( N0 ). DMIN2 (output) DOUBLE PRECISION
45 Minimum value of d, excluding D( N0 ) and D( N0-1 ).
46
47 DN (output) DOUBLE PRECISION
48 d(N0), the last value of d.
49
50 DNM1 (output) DOUBLE PRECISION
51 d(N0-1).
52
53 DNM2 (output) DOUBLE PRECISION
54 d(N0-2).
55
56 IEEE (input) LOGICAL
57 Flag for IEEE or non IEEE arithmetic.
58
59
60
61 LAPACK routine (version 3.2) November 2008 DLASQ5(1)