HP MLIB User's Guide Vol. 2 7th Ed.

Chapter 13 Sparse Linear Equations 907
One-call usage DSLEFS
Notes Calling DSLEFS is equivalent to calling DSLEIN, DSLEIM, DSLEOR,
DSLEVM, DSLECO, and DSLESL in sequence. You can follow the call to
DSLEFS with other calls to subprograms in the package to solve additional
right-hand sides, to input new matrix values, to input a new matrix structure,
or to perform utility functions.
Example Solve the small (order 6) sparse system of linear equations where the matrix A
and right-hand side b are
A supernode relaxation factor of 0 (no extra fill allowed) and a pivoting
tolerance of 0.1 are used. Printed output is to Fortran logical unit 6. Only
printed error messages and some runtime statistics are printed.
A =
4.0 0.0 1.0 1.0 0.0 0.0
b =
1.0
0.0 4.0 1.0 0.0 1.0 0.0 0.0
1.0 1.0 4.0 0.0 1.0 0.0 1.0
1.0 0.0 0.0 4.0 1.0 1.0 1.0
0.0 1.0 1.0 1.0 4.0 0.0 0.0
0.0 0.0 0.0 1.0 0.0 1.0 0.0