110 SUBROUTINE spotf2( UPLO, N, A, LDA, INFO )
129 parameter( one = 1.0e+0, zero = 0.0e+0 )
152 upper =
lsame( uplo,
'U' )
153 IF( .NOT.upper .AND. .NOT.
lsame( uplo,
'L' ) )
THEN
155 ELSE IF( n.LT.0 )
THEN
157 ELSE IF( lda.LT.max( 1, n ) )
THEN
161 CALL
xerbla(
'SPOTF2', -info )
178 ajj = a( j, j ) -
sdot( j-1, a( 1, j ), 1, a( 1, j ), 1 )
179 IF( ajj.LE.zero.OR.
sisnan( ajj ) )
THEN
189 CALL
sgemv(
'Transpose', j-1, n-j, -one, a( 1, j+1 ),
190 $ lda, a( 1, j ), 1, one, a( j, j+1 ), lda )
191 CALL
sscal( n-j, one / ajj, a( j, j+1 ), lda )
202 ajj = a( j, j ) -
sdot( j-1, a( j, 1 ), lda, a( j, 1 ),
204 IF( ajj.LE.zero.OR.
sisnan( ajj ) )
THEN
214 CALL
sgemv(
'No transpose', n-j, j-1, -one, a( j+1, 1 ),
215 $ lda, a( j, 1 ), lda, one, a( j+1, j ), 1 )
216 CALL
sscal( n-j, one / ajj, a( j+1, j ), 1 )