107 REAL function
slanhs( norm, n, a, lda, work )
118 REAL a( lda, * ), work( * )
125 parameter( one = 1.0e+0, zero = 0.0e+0 )
129 REAL scale, sum, value
139 INTRINSIC abs, min, sqrt
145 ELSE IF(
lsame( norm,
'M' ) )
THEN
151 DO 10 i = 1, min( n, j+1 )
152 sum = abs( a( i, j ) )
153 IF(
VALUE .LT. sum .OR.
sisnan( sum ) )
VALUE = sum
156 ELSE IF( (
lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
163 DO 30 i = 1, min( n, j+1 )
164 sum = sum + abs( a( i, j ) )
166 IF(
VALUE .LT. sum .OR.
sisnan( sum ) )
VALUE = sum
168 ELSE IF(
lsame( norm,
'I' ) )
THEN
176 DO 60 i = 1, min( n, j+1 )
177 work( i ) = work( i ) + abs( a( i, j ) )
183 IF(
VALUE .LT. sum .OR.
sisnan( sum ) )
VALUE = sum
185 ELSE IF( (
lsame( norm,
'F' ) ) .OR. (
lsame( norm,
'E' ) ) )
THEN
192 CALL slassq( min( n, j+1 ), a( 1, j ), 1, scale, sum )
194 VALUE = scale*sqrt( sum )
subroutine slassq(n, x, incx, scl, sumsq)
SLASSQ updates a sum of squares represented in scaled form.
logical function sisnan(SIN)
SISNAN tests input for NaN.
logical function lsame(CA, CB)
LSAME
real function slanhs(NORM, N, A, LDA, WORK)
SLANHS returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...