115 DOUBLE PRECISION FUNCTION dlange( NORM, M, N, A, LDA, WORK )
127 DOUBLE PRECISION A( lda, * ), WORK( * )
133 DOUBLE PRECISION ONE, ZERO
134 parameter ( one = 1.0d+0, zero = 0.0d+0 )
138 DOUBLE PRECISION SCALE, SUM,
VALUE, TEMP
144 LOGICAL LSAME, DISNAN
145 EXTERNAL lsame, disnan
148 INTRINSIC abs, min, sqrt
152 IF( min( m, n ).EQ.0 )
THEN
154 ELSE IF( lsame( norm,
'M' ) )
THEN
161 temp = abs( a( i, j ) )
162 IF(
VALUE.LT.temp .OR. disnan( temp ) )
VALUE = temp
165 ELSE IF( ( lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
173 sum = sum + abs( a( i, j ) )
175 IF(
VALUE.LT.sum .OR. disnan( sum ) )
VALUE = sum
177 ELSE IF( lsame( norm,
'I' ) )
THEN
186 work( i ) = work( i ) + abs( a( i, j ) )
192 IF(
VALUE.LT.temp .OR. disnan( temp ) )
VALUE = temp
194 ELSE IF( ( lsame( norm,
'F' ) ) .OR. ( lsame( norm,
'E' ) ) )
THEN
201 CALL dlassq( m, a( 1, j ), 1, scale, sum )
203 VALUE = scale*sqrt( sum )
double precision function dlange(NORM, M, N, A, LDA, WORK)
DLANGE returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
subroutine dlassq(N, X, INCX, SCALE, SUMSQ)
DLASSQ updates a sum of squares represented in scaled form.