1 SUBROUTINE slacn2( N, V, X, ISGN, EST, KASE, ISAVE )
12 INTEGER isgn( * ), isave( 3 )
77 parameter( itmax = 5 )
79 parameter( zero = 0.0e+0, one = 1.0e+0, two = 2.0e+0 )
83 REAL altsgn, estold, temp
88 EXTERNAL isamax, sasum
94 INTRINSIC abs, nint,
REAL, sign
100 x( i ) = one /
REAL( n )
107 go to( 20, 40, 70, 110, 140 )isave( 1 )
119 est = sasum( n, x, 1 )
122 x( i ) = sign( one, x( i ) )
123 isgn( i ) = nint( x( i ) )
133 isave( 2 ) = isamax( n, x, 1 )
142 x( isave( 2 ) ) = one
151 CALL scopy( n, x, 1, v, 1 )
153 est = sasum( n, v, 1 )
155 IF( nint( sign( one, x( i ) ) ).NE.isgn( i ) )
167 x( i ) = sign( one, x( i ) )
168 isgn( i ) = nint( x( i ) )
179 isave( 2 ) = isamax( n, x, 1 )
180 IF( ( x( jlast ).NE.abs( x( isave( 2 ) ) ) ) .AND.
181 $ ( isave( 3 ).LT.itmax ) )
THEN
182 isave( 3 ) = isave( 3 ) + 1
191 x( i ) = altsgn*( one+
REAL( I-1 ) /
REAL( N-1 ) )
202 temp = two*( sasum( n, x, 1 ) /
REAL( 3*N ) )
203 IF( temp.GT.est )
THEN
204 CALL scopy( n, x, 1, v, 1 )