104 SUBROUTINE dlassq( N, X, INCX, SCALE, SUMSQ )
113 DOUBLE PRECISION SCALE, SUMSQ
116 DOUBLE PRECISION X( * )
122 DOUBLE PRECISION ZERO
123 parameter ( zero = 0.0d+0 )
127 DOUBLE PRECISION ABSXI
139 DO 10 ix = 1, 1 + ( n-1 )*incx, incx
140 absxi = abs( x( ix ) )
141 IF( absxi.GT.zero.OR.disnan( absxi ) )
THEN
142 IF( scale.LT.absxi )
THEN
143 sumsq = 1 + sumsq*( scale / absxi )**2
146 sumsq = sumsq + ( absxi / scale )**2
subroutine dlassq(N, X, INCX, SCALE, SUMSQ)
DLASSQ updates a sum of squares represented in scaled form.