259 SUBROUTINE csytrf_rk( UPLO, N, A, LDA, E, IPIV, WORK, LWORK,
269 INTEGER INFO, LDA, LWORK, N
273 COMPLEX A( lda, * ), E( * ), WORK( * )
279 LOGICAL LQUERY, UPPER
280 INTEGER I, IINFO, IP, IWS, K, KB, LDWORK, LWKOPT,
286 EXTERNAL lsame, ilaenv
299 upper = lsame( uplo,
'U' )
300 lquery = ( lwork.EQ.-1 )
301 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
303 ELSE IF( n.LT.0 )
THEN
305 ELSE IF( lda.LT.max( 1, n ) )
THEN
307 ELSE IF( lwork.LT.1 .AND. .NOT.lquery )
THEN
315 nb = ilaenv( 1,
'CSYTRF_RK', uplo, n, -1, -1, -1 )
321 CALL xerbla(
'CSYTRF_RK', -info )
323 ELSE IF( lquery )
THEN
329 IF( nb.GT.1 .AND. nb.LT.n )
THEN
331 IF( lwork.LT.iws )
THEN
332 nb = max( lwork / ldwork, 1 )
333 nbmin = max( 2, ilaenv( 2,
'CSYTRF_RK',
334 $ uplo, n, -1, -1, -1 ) )
363 CALL clasyf_rk( uplo, k, nb, kb, a, lda, e,
364 $ ipiv, work, ldwork, iinfo )
369 CALL csytf2_rk( uplo, k, a, lda, e, ipiv, iinfo )
375 IF( info.EQ.0 .AND. iinfo.GT.0 )
391 DO i = k, ( k - kb + 1 ), -1
392 ip = abs( ipiv( i ) )
394 CALL cswap( n-k, a( i, k+1 ), lda,
395 $ a( ip, k+1 ), lda )
431 CALL clasyf_rk( uplo, n-k+1, nb, kb, a( k, k ), lda, e( k ),
432 $ ipiv( k ), work, ldwork, iinfo )
439 CALL csytf2_rk( uplo, n-k+1, a( k, k ), lda, e( k ),
447 IF( info.EQ.0 .AND. iinfo.GT.0 )
448 $ info = iinfo + k - 1
453 IF( ipiv( i ).GT.0 )
THEN
454 ipiv( i ) = ipiv( i ) + k - 1
456 ipiv( i ) = ipiv( i ) - k + 1
470 DO i = k, ( k + kb - 1 ), 1
471 ip = abs( ipiv( i ) )
473 CALL cswap( k-1, a( i, 1 ), lda,
subroutine clasyf_rk(UPLO, N, NB, KB, A, LDA, E, IPIV, W, LDW, INFO)
CLASYF_RK computes a partial factorization of a complex symmetric indefinite matrix using bounded Bun...
subroutine csytrf_rk(UPLO, N, A, LDA, E, IPIV, WORK, LWORK, INFO)
CSYTRF_RK computes the factorization of a complex symmetric indefinite matrix using the bounded Bunch...
subroutine csytf2_rk(UPLO, N, A, LDA, E, IPIV, INFO)
CSYTF2_RK computes the factorization of a complex symmetric indefinite matrix using the bounded Bunch...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine cswap(N, CX, INCX, CY, INCY)
CSWAP