1 DOUBLE PRECISION FUNCTION dlanhs( NORM, N, A, LDA, WORK )
13 DOUBLE PRECISION a( lda, * ), work( * )
66 DOUBLE PRECISION one, zero
67 parameter( one = 1.0d+0, zero = 0.0d+0 )
71 DOUBLE PRECISION scale, sum, value
81 INTRINSIC abs, max, min, sqrt
87 ELSE IF(
lsame( norm,
'M' ) )
THEN
93 DO 10 i = 1, min( n, j+1 )
94 VALUE = max(
VALUE, abs( a( i, j ) ) )
97 ELSE IF( (
lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
104 DO 30 i = 1, min( n, j+1 )
105 sum = sum + abs( a( i, j ) )
107 VALUE = max(
VALUE, sum )
109 ELSE IF(
lsame( norm,
'I' ) )
THEN
117 DO 60 i = 1, min( n, j+1 )
118 work( i ) = work( i ) + abs( a( i, j ) )
123 VALUE = max(
VALUE, work( i ) )
125 ELSE IF( (
lsame( norm,
'F' ) ) .OR. (
lsame( norm,
'E' ) ) )
THEN
132 CALL dlassq( min( n, j+1 ), a( 1, j ), 1,
scale, sum )
134 VALUE =
scale*sqrt( sum )
subroutine scale(xyzl, nl)
double precision function dlanhs(NORM, N, A, LDA, WORK)
subroutine dlassq(N, X, INCX, SCALE, SUMSQ)
logical function lsame(CA, CB)