y bgcolor="#FFFFFF" text="#000000" link="#0000FF" vlink="#660099">

Driver Routines for Solving Systems of Linear Equations

?gesv computes the solution to the system of linear equations with a square matrix A and multiple right-hand sides.
Fortran 77:
call sgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
call dgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
call cgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
call zgesv(n, nrhs, a, lda, ipiv, b, ldb, info)
Fortran 95:
call gesv(a, b [,ipiv] [,info])

?gesvx computes the solution to the system of linear equations with a square matrix A and multiple right-hand sides, and provides error bounds on the solution.
Fortran 77:
call sgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zgesvx(fact, trans, n, nrhs, a, lda, af, ldaf, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call gesvx(a, b, x [,af] [,ipiv] [,fact] [,trans] [,equed] [,r] [,c] [,ferr] [,berr] [,rcond] [,rpvgrw] [,info])

?gbsv computes the solution to the system of linear equations with a band matrix A and multiple right-hand sides.
Fortran 77:
call sgbsv(n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
call dgbsv(n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
call cgbsv(n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
call zgbsv(n, kl, ku, nrhs, ab, ldab, ipiv, b, ldb, info)
Fortran 95:
call gbsv(a, b [,kl] [,ipiv] [,info])

?gbsvx computes the solution to the real or complex system of linear equations with a band matrix A and multiple right-hand sides, and provides error bounds on the solution.
Fortran 77:
call sgbsvx(fact, trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dgbsvx(fact, trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cgbsvx(fact, trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zgbsvx(fact, trans, n, kl, ku, nrhs, ab, ldab, afb, ldafb, ipiv, equed, r, c, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call gbsvx(a, b, x [,kl] [,af] [,ipiv] [,fact] [,trans] [,equed] [,r] [,c] [,ferr] [,berr] [,rcond] [,rpvgrw] [,info])

?gtsv computes the solution to the system of linear equations with a tridiagonal matrix A and multiple right-hand sides.
Fortran 77:
call sgtsv(n, nrhs, dl, d, du, b, ldb, info)
call dgtsv(n, nrhs, dl, d, du, b, ldb, info)
call cgtsv(n, nrhs, dl, d, du, b, ldb, info)
call zgtsv(n, nrhs, dl, d, du, b, ldb, info)
Fortran 95:
call gtsv(dl, d, du, b [,info])

?gtsvx computes the solution to the real or complex system of linear equations with a tridiagonal matrix A and multiple right-hand sides, and provides error bounds on the solution.
Fortran 77:
call sgtsvx(fact, trans, n, nrhs, dl, d, du, dlf, df, duf, du2, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dgtsvx(fact, trans, n, nrhs, dl, d, du, dlf, df, duf, du2, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cgtsvx(fact, trans, n, nrhs, dl, d, du, dlf, df, duf, du2, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zgtsvx(fact, trans, n, nrhs, dl, d, du, dlf, df, duf, du2, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call gtsvx(dl, d, du, b, x [,dlf] [,df] [,duf] [,du2] [,ipiv] [,fact] [,trans] [,ferr] [,berr] [,rcond] [,info])

?posv computes the solution to the system of linear equations with a symmetric or Hermitian positive definite matrix A and multiple right-hand sides.
Fortran 77:
call sposv(uplo, n, nrhs, a, lda, b, ldb, info)
call dposv(uplo, n, nrhs, a, lda, b, ldb, info)
call cposv(uplo, n, nrhs, a, lda, b, ldb, info)
call zposv(uplo, n, nrhs, a, lda, b, ldb, info)
Fortran 95:
call posv(a, b [,uplo] [,info])

?posvx uses the Cholesky factorization to compute the solution to the system of linear equations with a symmetric or Hermitian positive-definite matrix A, and provides error bounds on the solution.
Fortran 77:
call sposvx(fact, uplo, n, nrhs, a, lda, af, ldaf, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dposvx(fact, uplo, n, nrhs, a, lda, af, ldaf, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cposvx(fact, uplo, n, nrhs, a, lda, af, ldaf, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zposvx(fact, uplo, n, nrhs, a, lda, af, ldaf, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call posvx(a, b, x [,uplo] [,af] [,fact] [,equed] [,s] [,ferr] [,berr] [,rcond] [,info])

?ppsv computes the solution to the system of linear equations with a symmetric (Hermitian) positive-definite packed matrix A and multiple right-hand sides.
Fortran 77:
call sppsv(uplo, n, nrhs, ap, b, ldb, info)
call dppsv(uplo, n, nrhs, ap, b, ldb, info)
call cppsv(uplo, n, nrhs, ap, b, ldb, info)
call zppsv(uplo, n, nrhs, ap, b, ldb, info)
Fortran 95:
call ppsv(a, b [,uplo] [,info])

?ppsvx uses the Cholesky factorization to compute the solution to the system of linear equations with a symmetric (Hermitian) positive-definite packed matrix A and provides error bounds on the solution.
Fortran 77:
call sppsvx(fact, uplo, n, nrhs, ap, afp, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dppsvx(fact, uplo, n, nrhs, ap, afp, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cppsvx(fact, uplo, n, nrhs, ap, afp, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zppsvx(fact, uplo, n, nrhs, ap, afp, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call ppsvx(a, b, x [,uplo] [,af] [,fact] [,equed] [,s] [,ferr] [,berr] [,rcond] [,info])

?pbsv computes the solution to the system of linear equations with a symmetric or Hermitian positive-definite band matrix A and multiple right-hand sides.
Fortran 77:
call spbsv(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
call dpbsv(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
call cpbsv(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
call zpbsv(uplo, n, kd, nrhs, ab, ldab, b, ldb, info)
Fortran 95:
call pbsv(a, b [,uplo] [,info])

?pbsvx uses the Cholesky factorization to compute the solution to the system of linear equations with a symmetric (Hermitian) positive-definite band matrix A and provides error bounds on the solution.
Fortran 77:
call spbsvx(fact, uplo, n, kd, nrhs, ab, ldab, afb, ldafb, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dpbsvx(fact, uplo, n, kd, nrhs, ab, ldab, afb, ldafb, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cpbsvx(fact, uplo, n, kd, nrhs, ab, ldab, afb, ldafb, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zpbsvx(fact, uplo, n, kd, nrhs, ab, ldab, afb, ldafb, equed, s, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call pbsvx(a, b, x [,uplo] [,af] [,fact] [,equed] [,s] [,ferr] [,berr] [,rcond] [,info])

?ptsv computes the solution to the system of linear equations with a symmetric or Hermitian positive-definite tridiagonal matrix A and multiple right-hand sides.
Fortran 77:
call sptsv(n, nrhs, d, e, b, ldb, info)
call dptsv(n, nrhs, d, e, b, ldb, info)
call cptsv(n, nrhs, d, e, b, ldb, info)
call zptsv(n, nrhs, d, e, b, ldb, info)
Fortran 95:
call ptsv(d, e, b [,info])

?ptsvx uses the factorization A=LDLH to compute the solution to the system of linear equations with a symmetric (Hermitian) positive-definite tridiagonal matrix A and provides error bounds on the solution.
Fortran 77:
call sptsvx(fact, n, nrhs, d, e, df, ef, b, ldb, x, ldx, rcond, ferr, berr, work, info)
call dptsvx(fact, n, nrhs, d, e, df, ef, b, ldb, x, ldx, rcond, ferr, berr, work, info)
call cptsvx(fact, n, nrhs, d, e, df, ef, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zptsvx(fact, n, nrhs, d, e, df, ef, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call ptsvx(d, e, b, x [,df] [,ef] [,fact] [,ferr] [,berr] [,rcond] [,info])

?sysv computes the solution to the system of linear equations with a real or complex symmetric matrix A and multiple right-hand sides.
Fortran 77:
call ssysv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
call dsysv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
call csysv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
call zsysv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
Fortran 95:
call sysv(a, b [,uplo] [,ipiv] [,info])

?sysvx uses the diagonal pivoting factorization to compute the solution to the system of linear equations with a real or complex symmetric matrix A and provides error bounds on the solution.
Fortran 77:
call ssysvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, iwork, info)
call dsysvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, iwork, info)
call csysvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, rwork, info)
call zsysvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, rwork, info)
Fortran 95:
call sysvx(a, b, x [,uplo] [,af] [,ipiv] [,fact] [,ferr] [,berr] [,rcond] [,info])

?hesv computes the solution to the system of linear equations with a Hermitian matrix A and multiple right-hand sides.
Fortran 77:
call chesv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
call zhesv(uplo, n, nrhs, a, lda, ipiv, b, ldb, work, lwork, info)
Fortran 95:
call hesv(a, b [,uplo] [,ipiv] [,info])

?hesvx uses the diagonal pivoting factorization to compute the solution to the complex system of linear equations with a Hermitian matrix A and provides error bounds on the solution.
Fortran 77:
call chesvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, rwork, info)
call zhesvx(fact, uplo, n, nrhs, a, lda, af, ldaf, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, lwork, rwork, info)
Fortran 95:
call hesvx(a, b, x [,uplo] [,af] [,ipiv] [,fact] [,ferr] [,berr] [,rcond] [,info])

?spsv computes the solution to the system of linear equations with a real or complex symmetric matrix A stored in packed format, and multiple right-hand sides.
Fortran 77:
call sspsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
call dspsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
call cspsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
call zspsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
Fortran 95:
call spsv(a, b [,uplo] [,ipiv] [,info])

?spsvx uses the diagonal pivoting factorization to compute the solution to the system of linear equations with a real or complex symmetric matrix A stored in packed format and provides error bounds on the solution.
Fortran 77:
call sspsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call dspsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, iwork, info)
call cspsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zspsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call spsvx(a, b, x [,uplo] [,af] [,ipiv] [,fact] [,ferr] [,berr] [,rcond] [,info])

?hpsv computes the solution to the system of linear equations with a Hermitian matrix A stored in packed format, and multiple right-hand sides.
Fortran 77:
call chpsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
call zhpsv(uplo, n, nrhs, ap, ipiv, b, ldb, info)
Fortran 95:
call hpsv(a, b [,uplo] [,ipiv] [,info])

?hpsvx uses the diagonal pivoting factorization to compute the solution to the system of linear equations with a Hermitian matrix A stored in packed format and provides error bounds on the solution.
Fortran 77:
call chpsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
call zhpsvx(fact, uplo, n, nrhs, ap, afp, ipiv, b, ldb, x, ldx, rcond, ferr, berr, work, rwork, info)
Fortran 95:
call hpsvx(a, b, x [,uplo] [,af] [,ipiv] [,fact] [,ferr] [,berr] [,rcond] [,info])

* Legal Information © 1999, 2002-2006, Intel Corporation