186 $ IPIV, IPIV2, B, LDB, WORK, LWORK,
197 INTEGER N, NRHS, LDA, LDB, LTB, LWORK, INFO
200 INTEGER IPIV( * ), IPIV2( * )
201 DOUBLE PRECISION A( LDA, * ), B( LDB, * ), TB( * ), WORK( * )
207 LOGICAL UPPER, TQUERY, WQUERY
226 upper = lsame( uplo,
'U' )
227 wquery = ( lwork.EQ.-1 )
228 tquery = ( ltb.EQ.-1 )
229 IF( .NOT.upper .AND. .NOT.lsame( uplo,
'L' ) )
THEN
231 ELSE IF( n.LT.0 )
THEN
233 ELSE IF( nrhs.LT.0 )
THEN
235 ELSE IF( lda.LT.max( 1, n ) )
THEN
237 ELSE IF( ltb.LT.( 4*n ) .AND. .NOT.tquery )
THEN
239 ELSE IF( ldb.LT.max( 1, n ) )
THEN
241 ELSE IF( lwork.LT.n .AND. .NOT.wquery )
THEN
247 $ ipiv2, work, -1, info )
248 lwkopt = int( work(1) )
252 CALL xerbla(
'DSYSV_AA_2STAGE', -info )
254 ELSE IF( wquery .OR. tquery )
THEN
262 $ work, lwork, info )
268 $ ipiv2, b, ldb, info )
subroutine xerbla(srname, info)
subroutine dsysv_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, work, lwork, info)
DSYSV_AA_2STAGE computes the solution to system of linear equations A * X = B for SY matrices
subroutine dsytrf_aa_2stage(uplo, n, a, lda, tb, ltb, ipiv, ipiv2, work, lwork, info)
DSYTRF_AA_2STAGE
subroutine dsytrs_aa_2stage(uplo, n, nrhs, a, lda, tb, ltb, ipiv, ipiv2, b, ldb, info)
DSYTRS_AA_2STAGE