127 COMPLEX v( n ),
x( n )
134 parameter( itmax = 5 )
136 parameter( one = 1.0e0, two = 2.0e0 )
138 parameter( czero = ( 0.0e0, 0.0e0 ),
139 $ cone = ( 1.0e0, 0.0e0 ) )
142 INTEGER i, iter,
j, jlast, jump
143 REAL absxi, altsgn, estold, safmin, temp
154 INTRINSIC abs, aimag, cmplx, real
161 safmin =
slamch(
'Safe minimum' )
164 x( i ) = cmplx( one /
REAL( N ) )
171 go to( 20, 40, 70, 90, 120 )jump
186 absxi = abs(
x( i ) )
187 IF( absxi.GT.safmin )
THEN
188 x( i ) = cmplx(
REAL( X( I ) ) / absxi,
189 $ aimag(
x( i ) ) / absxi )
220 CALL
ccopy( n,
x, 1, v, 1 )
229 absxi = abs(
x( i ) )
230 IF( absxi.GT.safmin )
THEN
231 x( i ) = cmplx(
REAL( X( I ) ) / absxi,
232 $ aimag(
x( i ) ) / absxi )
247 IF( ( abs(
x( jlast ) ).NE.abs(
x(
j ) ) ) .AND.
248 $ ( iter.LT.itmax ) )
THEN
258 x( i ) = cmplx( altsgn*( one+
REAL( I-1 ) /
REAL( N-1 ) ) )
269 temp = two*(
scsum1( n,
x, 1 ) /
REAL( 3*N ) )
270 IF( temp.GT.est )
THEN
271 CALL
ccopy( n,
x, 1, v, 1 )
subroutine clacon(N, V, X, EST, KASE)
CLACON estimates the 1-norm of a square matrix, using reverse communication for evaluating matrix-vec...
REAL function slamch(CMACH)
SLAMCH
REAL function scsum1(N, CX, INCX)
SCSUM1 forms the 1-norm of the complex vector using the true absolute value.
INTEGER function icmax1(N, CX, INCX)
ICMAX1 finds the index of the vector element whose real part has maximum absolute value...
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine ccopy(N, CX, INCX, CY, INCY)
CCOPY