1DLASQ5(1) LAPACK auxiliary routine (version 3.1) DLASQ5(1)
2
3
4
6 DLASQ5 - one dqds transform in ping-pong form, one version for IEEE
7 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
25
27 I0 (input) INTEGER
28 First index.
29
30 N0 (input) INTEGER
31 Last index.
32
33 Z (input) DOUBLE PRECISION array, dimension ( 4*N )
34 Z holds the qd array. EMIN is stored in Z(4*N0) to avoid an extra
35 argument.
36
37 PP (input) INTEGER
38 PP=0 for ping, PP=1 for pong.
39
40 TAU (input) DOUBLE PRECISION
41 This is the shift.
42
43 DMIN (output) DOUBLE PRECISION
44 Minimum value of d.
45
46 DMIN1 (output) DOUBLE PRECISION Minimum value of d, excluding D(
47 N0 ).
48
49 DMIN2 (output) DOUBLE PRECISION Minimum value of d, excluding D(
50 N0 ) and D( N0-1 ).
51
52 DN (output) DOUBLE PRECISION
53 d(N0), the last value of d.
54
55 DNM1 (output) DOUBLE PRECISION
56 d(N0-1).
57
58 DNM2 (output) DOUBLE PRECISION
59 d(N0-2).
60
61 IEEE (input) LOGICAL
62 Flag for IEEE or non IEEE arithmetic.
63
64
65
66 LAPACK auxiliary routine (versionNo3v.e1m)ber 2006 DLASQ5(1)