116 DOUBLE PRECISION FUNCTION zlange( NORM, M, N, A, LDA, WORK )
128 DOUBLE PRECISION work( * )
129 COMPLEX*16 a( lda, * )
135 DOUBLE PRECISION one, zero
136 parameter( one = 1.0d+0, zero = 0.0d+0 )
140 DOUBLE PRECISION scale, sum,
value, temp
150 INTRINSIC abs, min, sqrt
154 IF( min( m, n ).EQ.0 )
THEN
156 ELSE IF(
lsame( norm,
'M' ) )
THEN
163 temp = abs( a( i,
j ) )
167 ELSE IF( (
lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
175 sum = sum + abs( a( i,
j ) )
179 ELSE IF(
lsame( norm,
'I' ) )
THEN
188 work( i ) = work( i ) + abs( a( i,
j ) )
196 ELSE IF( (
lsame( norm,
'F' ) ) .OR. (
lsame( norm,
'E' ) ) )
THEN
203 CALL
zlassq( m, a( 1,
j ), 1, scale, sum )
205 value = scale*sqrt( sum )
LOGICAL function lsame(CA, CB)
LSAME
subroutine zlassq(N, X, INCX, SCALE, SUMSQ)
ZLASSQ updates a sum of squares represented in scaled form.
input scalars passed by value
DOUBLE PRECISION function zlange(NORM, M, N, A, LDA, WORK)
ZLANGE returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
LOGICAL function disnan(DIN)
DISNAN tests input for NaN.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j