104 SUBROUTINE slassq( N, X, INCX, SCALE, SUMSQ )
123 parameter ( zero = 0.0e+0 )
139 DO 10 ix = 1, 1 + ( n-1 )*incx, incx
140 absxi = abs( x( ix ) )
141 IF( absxi.GT.zero.OR.sisnan( absxi ) )
THEN
142 IF( scale.LT.absxi )
THEN
143 sumsq = 1 + sumsq*( scale / absxi )**2
146 sumsq = sumsq + ( absxi / scale )**2
subroutine slassq(N, X, INCX, SCALE, SUMSQ)
SLASSQ updates a sum of squares represented in scaled form.