184 $ IPIV, IPIV2, B, LDB, WORK, LWORK,
195 INTEGER N, NRHS, LDA, LDB, LTB, LWORK, INFO
198 INTEGER IPIV( * ), IPIV2( * )
199 COMPLEX*16 A( LDA, * ), B( LDB, * ), TB( * ), WORK( * )
205 LOGICAL UPPER, TQUERY, WQUERY
224 upper = lsame( uplo,
'U' )
225 wquery = ( lwork.EQ.-1 )
226 tquery = ( ltb.EQ.-1 )
227 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
229 ELSE IF( n.LT.0 )
THEN
231 ELSE IF( nrhs.LT.0 )
THEN
233 ELSE IF( lda.LT.max( 1, n ) )
THEN
235 ELSE IF( ltb.LT.( 4*n ) .AND. .NOT.tquery )
THEN
237 ELSE IF( ldb.LT.max( 1, n ) )
THEN
239 ELSE IF( lwork.LT.n .AND. .NOT.wquery )
THEN
245 $ ipiv2, work, -1, info )
246 lwkopt = int( work(1) )
250 CALL xerbla(
'ZSYSV_AA_2STAGE', -info )
252 ELSE IF( wquery .OR. tquery )
THEN
260 $ work, lwork, info )
266 $ ipiv2, b, ldb, info )
subroutine xerbla(srname, info)
subroutine zsysv_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, work, lwork, info)
ZSYSV_AA_2STAGE computes the solution to system of linear equations A * X = B for SY matrices
subroutine zsytrf_aa_2stage(uplo, n, a, lda, tb, ltb, ipiv, ipiv2, work, lwork, info)
ZSYTRF_AA_2STAGE
subroutine zsytrs_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, info)
ZSYTRS_AA_2STAGE