228 SUBROUTINE dsysv_rk( UPLO, N, NRHS, A, LDA, E, IPIV, B, LDB,
229 $ work, lwork, info )
238 INTEGER INFO, LDA, LDB, LWORK, N, NRHS
242 DOUBLE PRECISION A( lda, * ), B( ldb, * ), E( * ), WORK( * )
266 lquery = ( lwork.EQ.-1 )
267 IF( .NOT.lsame( uplo,
'U' ) .AND. .NOT.lsame( uplo,
'L' ) )
THEN
269 ELSE IF( n.LT.0 )
THEN
271 ELSE IF( nrhs.LT.0 )
THEN
273 ELSE IF( lda.LT.max( 1, n ) )
THEN
275 ELSE IF( ldb.LT.max( 1, n ) )
THEN
277 ELSE IF( lwork.LT.1 .AND. .NOT.lquery )
THEN
285 CALL dsytrf_rk( uplo, n, a, lda, e, ipiv, work, -1, info )
292 CALL xerbla(
'DSYSV_RK ', -info )
294 ELSE IF( lquery )
THEN
301 CALL dsytrf_rk( uplo, n, a, lda, e, ipiv, work, lwork, info )
307 CALL dsytrs_3( uplo, n, nrhs, a, lda, e, ipiv, b, ldb, info )
subroutine dsytrf_rk(UPLO, N, A, LDA, E, IPIV, WORK, LWORK, INFO)
DSYTRF_RK computes the factorization of a real symmetric indefinite matrix using the bounded Bunch-Ka...
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dsytrs_3(UPLO, N, NRHS, A, LDA, E, IPIV, B, LDB, INFO)
DSYTRS_3
subroutine dsysv_rk(UPLO, N, NRHS, A, LDA, E, IPIV, B, LDB, WORK, LWORK, INFO)
DSYSV_RK computes the solution to system of linear equations A * X = B for SY matrices ...