124 REAL FUNCTION slangb( NORM, N, KL, KU, AB, LDAB,
134 INTEGER kl, ku, ldab, n
137 REAL ab( ldab, * ), work( * )
145 parameter( one = 1.0e+0, zero = 0.0e+0 )
149 REAL scale, sum,
value, temp
159 INTRINSIC abs, max, min, sqrt
165 ELSE IF(
lsame( norm,
'M' ) )
THEN
171 DO 10 i = max( ku+2-
j, 1 ), min( n+ku+1-
j, kl+ku+1 )
172 temp = abs( ab( i,
j ) )
176 ELSE IF( (
lsame( norm,
'O' ) ) .OR. ( norm.EQ.
'1' ) )
THEN
183 DO 30 i = max( ku+2-
j, 1 ), min( n+ku+1-
j, kl+ku+1 )
184 sum = sum + abs( ab( i,
j ) )
188 ELSE IF(
lsame( norm,
'I' ) )
THEN
197 DO 60 i = max( 1,
j-ku ), min( n,
j+kl )
198 work( i ) = work( i ) + abs( ab( k+i,
j ) )
206 ELSE IF( (
lsame( norm,
'F' ) ) .OR. (
lsame( norm,
'E' ) ) )
THEN
215 CALL
slassq( min( n,
j+kl )-l+1, ab( k,
j ), 1, scale, sum )
217 value = scale*sqrt( sum )
subroutine slassq(N, X, INCX, SCALE, SUMSQ)
SLASSQ updates a sum of squares represented in scaled form.
LOGICAL function lsame(CA, CB)
LSAME
input scalars passed by value
REAL function slangb(NORM, N, KL, KU, AB, LDAB, WORK)
SLANGB returns the value of the 1-norm, Frobenius norm, infinity-norm, or the largest absolute value ...
LOGICAL function sisnan(SIN)
SISNAN tests input for NaN.
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j