1 SUBROUTINE dlasq5( I0, N0, Z, PP, TAU, DMIN, DMIN1, DMIN2, DN,
12 DOUBLE PRECISION DMIN, DMIN1, DMIN2, DN, DNM1, DNM2, TAU
15 DOUBLE PRECISION Z( * )
68 parameter( zero = 0.0d0 )
72 DOUBLE PRECISION D, EMIN, TEMP
79 IF( ( n0-i0-1 ).LE.0 )
93 DO 10 j4 = 4*i0, 4*( n0-3 ), 4
94 z( j4-2 ) = d + z( j4-1 )
95 temp = z( j4+1 ) / z( j4-2 )
98 z( j4 ) = z( j4-1 )*temp
99 emin = min( z( j4 ), emin )
102 DO 20 j4 = 4*i0, 4*( n0-3 ), 4
103 z( j4-3 ) = d + z( j4 )
104 temp = z( j4+2 ) / z( j4-3 )
106 dmin = min( dmin, d )
107 z( j4-1 ) = z( j4 )*temp
108 emin = min( z( j4-1 ), emin )
118 z( j4-2 ) = dnm2 + z( j4p2 )
119 z( j4 ) = z( j4p2+2 )*( z( j4p2 ) / z( j4-2 ) )
120 dnm1 = z( j4p2+2 )*( dnm2 / z( j4-2 ) ) - tau
121 dmin = min( dmin, dnm1 )
126 z( j4-2 ) = dnm1 + z( j4p2 )
127 z( j4 ) = z( j4p2+2 )*( z( j4p2 ) / z( j4-2 ) )
128 dn = z( j4p2+2 )*( dnm1 / z( j4-2 ) ) - tau
129 dmin = min( dmin, dn )
136 DO 30 j4 = 4*i0, 4*( n0-3 ), 4
137 z( j4-2 ) = d + z( j4-1 )
141 z( j4 ) = z( j4+1 )*( z( j4-1 ) / z( j4-2 ) )
142 d = z( j4+1 )*( d / z( j4-2 ) ) - tau
144 dmin = min( dmin, d )
145 emin = min( emin, z( j4 ) )
148 DO 40 j4 = 4*i0, 4*( n0-3 ), 4
149 z( j4-3 ) = d + z( j4 )
153 z( j4-1 ) = z( j4+2 )*( z( j4 ) / z( j4-3 ) )
154 d = z( j4+2 )*( d / z( j4-3 ) ) - tau
156 dmin = min( dmin, d )
157 emin = min( emin, z( j4-1 ) )
167 z( j4-2 ) = dnm2 + z( j4p2 )
168 IF( dnm2.LT.zero )
THEN
171 z( j4 ) = z( j4p2+2 )*( z( j4p2 ) / z( j4-2 ) )
172 dnm1 = z( j4p2+2 )*( dnm2 / z( j4-2 ) ) - tau
174 dmin = min( dmin, dnm1 )
179 z( j4-2 ) = dnm1 + z( j4p2 )
180 IF( dnm1.LT.zero )
THEN
183 z( j4 ) = z( j4p2+2 )*( z( j4p2 ) / z( j4-2 ) )
184 dn = z( j4p2+2 )*( dnm1 / z( j4-2 ) ) - tau
186 dmin = min( dmin, dn )
subroutine dlasq5(I0, N0, Z, PP, TAU, DMIN, DMIN1, DMIN2, DN, DNM1, DNM2, IEEE)