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
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 ) )
165 ELSE IF( (
lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
173 sum = sum + abs( a( i,
j ) )
177 ELSE IF(
lsame( norm,
'I' ) )
THEN
186 work( i ) = work( i ) + abs( a( i,
j ) )
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 )
input scalars passed by value
logical function lsame(CA, CB)
LSAME
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.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
logical function disnan(DIN)
DISNAN tests input for NaN.