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
135 DO 10 ix = 1, 1 + ( n-1 )*incx, incx
136 IF( x( ix ).NE.zero )
THEN 137 absxi = abs( x( ix ) )
138 IF( scale.LT.absxi )
THEN 139 sumsq = 1 + sumsq*( scale / absxi )**2
142 sumsq = sumsq + ( absxi / scale )**2