1 SUBROUTINE dlassq( N, X, INCX, SCALE, SUMSQ )
10 DOUBLE PRECISION SCALE, SUMSQ
13 DOUBLE PRECISION X( * )
61 parameter( zero = 0.0d+0 )
65 DOUBLE PRECISION ABSXI
73 DO 10 ix = 1, 1 + ( n-1 )*incx, incx
74 IF( x( ix ).NE.zero )
THEN
75 absxi = abs( x( ix ) )
76 IF( scale.LT.absxi )
THEN
77 sumsq = 1 + sumsq*( scale / absxi )**2
80 sumsq = sumsq + ( absxi / scale )**2
subroutine dlassq(N, X, INCX, SCALE, SUMSQ)