121 SUBROUTINE dgesv( N, NRHS, A, LDA, IPIV, B, LDB, INFO )
128 INTEGER INFO, LDA, LDB, N, NRHS
132 DOUBLE PRECISION A( LDA, * ), B( LDB, * )
150 ELSE IF( nrhs.LT.0 )
THEN
152 ELSE IF( lda.LT.max( 1, n ) )
THEN
154 ELSE IF( ldb.LT.max( 1, n ) )
THEN
158 CALL xerbla(
'DGESV ', -info )
164 CALL dgetrf( n, n, a, lda, ipiv, info )
169 CALL dgetrs(
'No transpose', n, nrhs, a, lda, ipiv, b, ldb,
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine dgetrf(M, N, A, LDA, IPIV, INFO)
DGETRF
subroutine dgetrs(TRANS, N, NRHS, A, LDA, IPIV, B, LDB, INFO)
DGETRS
subroutine dgesv(N, NRHS, A, LDA, IPIV, B, LDB, INFO)
DGESV computes the solution to system of linear equations A * X = B for GE matrices