124 REAL FUNCTION clansy( NORM, UPLO, N, A, LDA, WORK )
144 parameter( one = 1.0e+0, zero = 0.0e+0 )
148 REAL absa, scale, sum,
value
164 ELSE IF(
lsame( norm,
'M' ) )
THEN
169 IF(
lsame( uplo,
'U' ) )
THEN
172 sum = abs( a( i,
j ) )
179 sum = abs( a( i,
j ) )
184 ELSE IF( (
lsame( norm,
'I' ) ) .OR. (
lsame( norm,
'O' ) ) .OR.
185 $ ( norm.EQ.
'1' ) )
THEN
190 IF(
lsame( uplo,
'U' ) )
THEN
194 absa = abs( a( i,
j ) )
196 work( i ) = work( i ) + absa
198 work(
j ) = sum + abs( a(
j,
j ) )
209 sum = work(
j ) + abs( a(
j,
j ) )
211 absa = abs( a( i,
j ) )
213 work( i ) = work( i ) + absa
218 ELSE IF( (
lsame( norm,
'F' ) ) .OR. (
lsame( norm,
'E' ) ) )
THEN
224 IF(
lsame( uplo,
'U' ) )
THEN
226 CALL
classq(
j-1, a( 1,
j ), 1, scale, sum )
230 CALL
classq( n-
j, a(
j+1,
j ), 1, scale, sum )
234 CALL
classq( n, a, lda+1, scale, sum )
235 value = scale*sqrt( sum )