171 SUBROUTINE ssysv( UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK,
181 INTEGER info, lda, ldb, lwork, n, nrhs
185 REAL 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
ssytrf( uplo, n, a, lda, ipiv, work, -1, info )
235 CALL
xerbla(
'SSYSV ', -info )
237 ELSE IF( lquery )
THEN
243 CALL
ssytrf( uplo, n, a, lda, ipiv, work, lwork, info )
248 IF ( lwork.LT.n )
THEN
252 CALL
ssytrs( uplo, n, nrhs, a, lda, ipiv,
b, ldb, info )
258 CALL
ssytrs2( uplo,n,nrhs,a,lda,ipiv,
b,ldb,work,info )
subroutine ssysv(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, LWORK, INFO)
SSYSV computes the solution to system of linear equations A * X = B for SY matrices ...
LOGICAL function lsame(CA, CB)
LSAME
subroutine xerbla(SRNAME, INFO)
XERBLA
set ue cd $ADTTMP cat<< EOF > tmp f Program LinearEquations Implicit none Real b(3) integer i
subroutine ssytrf(UPLO, N, A, LDA, IPIV, WORK, LWORK, INFO)
SSYTRF
subroutine ssytrs(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
SSYTRS
subroutine ssytrs2(UPLO, N, NRHS, A, LDA, IPIV, B, LDB, WORK, INFO)
SSYTRS2