116 SUBROUTINE slacon( N, V, X, ISGN, EST, KASE )
136 parameter( itmax = 5 )
138 parameter( zero = 0.0e+0, one = 1.0e+0, two = 2.0e+0 )
141 INTEGER i, iter,
j, jlast, jump
142 REAL altsgn, estold, temp
153 INTRINSIC abs, nint,
REAL, sign
162 x( i ) = one /
REAL( n )
169 go to( 20, 40, 70, 110, 140 )jump
184 x( i ) = sign( one,
x( i ) )
185 isgn( i ) = nint(
x( i ) )
213 CALL
scopy( n,
x, 1, v, 1 )
215 est =
sasum( n, v, 1 )
217 IF( nint( sign( one,
x( i ) ) ).NE.isgn( i ) )
229 x( i ) = sign( one,
x( i ) )
230 isgn( i ) = nint(
x( i ) )
242 IF( (
x( jlast ).NE.abs(
x(
j ) ) ) .AND. ( iter.LT.itmax ) )
THEN
252 x( i ) = altsgn*( one+
REAL( I-1 ) /
REAL( N-1 ) )
263 temp = two*(
sasum( n,
x, 1 ) /
REAL( 3*N ) )
264 IF( temp.GT.est )
THEN
265 CALL
scopy( n,
x, 1, v, 1 )
INTEGER function isamax(N, SX, INCX)
ISAMAX
subroutine scopy(N, SX, INCX, SY, INCY)
SCOPY
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
REAL function sasum(N, SX, INCX)
SASUM
subroutine slacon(N, V, X, ISGN, EST, KASE)
SLACON estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...