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 )