108 REAL function
clanhs( norm, n, a, lda, work )
127 parameter( one = 1.0e+0, zero = 0.0e+0 )
131 REAL scale, sum, value
141 INTRINSIC abs, min, sqrt
147 ELSE IF(
lsame( norm,
'M' ) )
THEN
153 DO 10 i = 1, min( n, j+1 )
154 sum = abs( a( i, j ) )
155 IF(
VALUE .LT. sum .OR.
sisnan( sum ) )
VALUE = sum
158 ELSE IF( (
lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
165 DO 30 i = 1, min( n, j+1 )
166 sum = sum + abs( a( i, j ) )
168 IF(
VALUE .LT. sum .OR.
sisnan( sum ) )
VALUE = sum
170 ELSE IF(
lsame( norm,
'I' ) )
THEN
178 DO 60 i = 1, min( n, j+1 )
179 work( i ) = work( i ) + abs( a( i, j ) )
185 IF(
VALUE .LT. sum .OR.
sisnan( sum ) )
VALUE = sum
187 ELSE IF( (
lsame( norm,
'F' ) ) .OR. (
lsame( norm,
'E' ) ) )
THEN
194 CALL classq( min( n, j+1 ), a( 1, j ), 1, scale, sum )
196 VALUE = scale*sqrt( sum )
subroutine classq(n, x, incx, scl, sumsq)
CLASSQ 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 clanhs(NORM, N, A, LDA, WORK)
CLANHS returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...