112 SUBROUTINE cgetc2( N, A, LDA, IPIV, JPIV, INFO )
123 INTEGER ipiv( * ), jpiv( * )
131 parameter( zero = 0.0e+0, one = 1.0e+0 )
134 INTEGER i, ip, ipv,
j, jp, jpv
135 REAL bignum, eps, smin, smlnum, xmax
145 INTRINSIC abs, cmplx, max
153 smlnum =
slamch(
'S' ) / eps
154 bignum = one / smlnum
155 CALL
slabad( smlnum, bignum )
167 IF( abs( a( ip, jp ) ).GE.xmax )
THEN
168 xmax = abs( a( ip, jp ) )
175 $ smin = max( eps*xmax, smlnum )
180 $ CALL
cswap( n, a( ipv, 1 ), lda, a( i, 1 ), lda )
186 $ CALL
cswap( n, a( 1, jpv ), 1, a( 1, i ), 1 )
191 IF( abs( a( i, i ) ).LT.smin )
THEN
193 a( i, i ) = cmplx( smin, zero )
196 a(
j, i ) = a(
j, i ) / a( i, i )
198 CALL
cgeru( n-i, n-i, -cmplx( one ), a( i+1, i ), 1,
199 $ a( i, i+1 ), lda, a( i+1, i+1 ), lda )
202 IF( abs( a( n, n ) ).LT.smin )
THEN
204 a( n, n ) = cmplx( smin, zero )
REAL function slamch(CMACH)
SLAMCH
subroutine cgetc2(N, A, LDA, IPIV, JPIV, INFO)
CGETC2 computes the LU factorization with complete pivoting of the general n-by-n matrix...
subroutine cswap(N, CX, INCX, CY, INCY)
CSWAP
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine slabad(SMALL, LARGE)
SLABAD
subroutine cgeru(M, N, ALPHA, X, INCX, Y, INCY, A, LDA)
CGERU