110 SUBROUTINE cpotf2( UPLO, N, A, LDA, INFO )
129 parameter( one = 1.0e+0, zero = 0.0e+0 )
131 parameter( cone = ( 1.0e+0, 0.0e+0 ) )
147 INTRINSIC max,
REAL, sqrt
154 upper =
lsame( uplo,
'U' )
155 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
157 ELSE IF( n.LT.0 )
THEN
159 ELSE IF( lda.LT.max( 1, n ) )
THEN
163 CALL
xerbla(
'CPOTF2', -info )
180 ajj =
REAL( A( J, J ) ) -
cdotc(
j-1, a( 1,
j ), 1,
182 IF( ajj.LE.zero.OR.
sisnan( ajj ) )
THEN
193 CALL
cgemv(
'Transpose',
j-1, n-
j, -cone, a( 1,
j+1 ),
194 $ lda, a( 1,
j ), 1, cone, a(
j,
j+1 ), lda )
196 CALL
csscal( n-
j, one / ajj, a(
j,
j+1 ), lda )
207 ajj =
REAL( A( J, J ) ) -
cdotc(
j-1, a(
j, 1 ), lda,
209 IF( ajj.LE.zero.OR.
sisnan( ajj ) )
THEN
220 CALL
cgemv(
'No transpose', n-
j,
j-1, -cone, a(
j+1, 1 ),
221 $ lda, a(
j, 1 ), lda, cone, a(
j+1,
j ), 1 )
223 CALL
csscal( n-
j, one / ajj, a(
j+1,
j ), 1 )
LOGICAL function lsame(CA, CB)
LSAME
subroutine xerbla(SRNAME, INFO)
XERBLA
LOGICAL function sisnan(SIN)
SISNAN tests input for NaN.
subroutine clacgv(N, X, INCX)
CLACGV conjugates a complex vector.
subroutine cgemv(TRANS, M, N, ALPHA, A, LDA, X, INCX, BETA, Y, INCY)
CGEMV
subroutine cpotf2(UPLO, N, A, LDA, INFO)
CPOTF2 computes the Cholesky factorization of a symmetric/Hermitian positive definite matrix (unblock...
COMPLEX function cdotc(N, CX, INCX, CY, INCY)
CDOTC
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real j
subroutine csscal(N, SA, CX, INCX)
CSSCAL