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 )
subroutine zlassq(N, X, INCX, SCALE, SUMSQ)
ZLASSQ updates a sum of squares represented in scaled form.
input scalars passed by value
logical function lsame(CA, CB)
LSAME
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 ...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
logical function disnan(DIN)
DISNAN tests input for NaN.