183 $ IPIV, IPIV2, B, LDB, WORK, LWORK,
194 INTEGER N, NRHS, LDA, LDB, LTB, LWORK, INFO
197 INTEGER IPIV( * ), IPIV2( * )
198 COMPLEX A( LDA, * ), B( LDB, * ), TB( * ), WORK( * )
204 LOGICAL UPPER, TQUERY, WQUERY
205 INTEGER LWKMIN, LWKOPT
210 EXTERNAL lsame, sroundup_lwork
224 upper = lsame( uplo,
'U' )
225 wquery = ( lwork.EQ.-1 )
226 tquery = ( ltb.EQ.-1 )
228 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
230 ELSE IF( n.LT.0 )
THEN
232 ELSE IF( nrhs.LT.0 )
THEN
234 ELSE IF( lda.LT.max( 1, n ) )
THEN
236 ELSE IF( ltb.LT.max( 1, 4*n ) .AND. .NOT.tquery )
THEN
238 ELSE IF( ldb.LT.max( 1, n ) )
THEN
240 ELSE IF( lwork.LT.lwkmin .AND. .NOT.wquery )
THEN
246 $ ipiv2, work, -1, info )
247 lwkopt = max( lwkmin, int( work( 1 ) ) )
248 work( 1 ) = sroundup_lwork( lwkopt )
252 CALL xerbla(
'CHESV_AA_2STAGE', -info )
254 ELSE IF( wquery .OR. tquery )
THEN
261 $ work, lwork, info )
267 $ ipiv2, b, ldb, info )
271 work( 1 ) = sroundup_lwork( lwkopt )
subroutine chesv_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, work, lwork, info)
CHESV_AA_2STAGE computes the solution to system of linear equations A * X = B for HE matrices
subroutine chetrf_aa_2stage(uplo, n, a, lda, tb, ltb, ipiv, ipiv2, work, lwork, info)
CHETRF_AA_2STAGE
subroutine chetrs_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, info)
CHETRS_AA_2STAGE