171 SUBROUTINE dsysv( UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK,
181 INTEGER INFO, LDA, LDB, LWORK, N, NRHS
185 DOUBLE PRECISION A( lda, * ), B( ldb, * ), WORK( * )
209 lquery = ( lwork.EQ.-1 )
210 IF( .NOT.lsame( uplo,
'U' ) .AND. .NOT.lsame( uplo,
'L' ) )
THEN
212 ELSE IF( n.LT.0 )
THEN
214 ELSE IF( nrhs.LT.0 )
THEN
216 ELSE IF( lda.LT.max( 1, n ) )
THEN
218 ELSE IF( ldb.LT.max( 1, n ) )
THEN
220 ELSE IF( lwork.LT.1 .AND. .NOT.lquery )
THEN
228 CALL
dsytrf( uplo, n, a, lda, ipiv, work, -1, info )
235 CALL
xerbla(
'DSYSV ', -info )
237 ELSE IF( lquery )
THEN
243 CALL
dsytrf( uplo, n, a, lda, ipiv, work, lwork, info )
248 IF ( lwork.LT.n )
THEN
252 CALL
dsytrs( uplo, n, nrhs, a, lda, ipiv, b, ldb, info )
258 CALL
dsytrs2( uplo,n,nrhs,a,lda,ipiv,b,ldb,work,info )
subroutine dsytrs2(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, INFO)
DSYTRS2
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dsysv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
DSYSV computes the solution to system of linear equations A * X = B for SY matrices ...
subroutine dsytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
DSYTRS
subroutine dsytrf(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)
DSYTRF