Numerical Computing with Modern Fortran

Richard J. Hanson and Tim Hopkins

Chapter 12: Case Study: Two Sparse Least-Squares System Examples

Source Code:

  1. in zip format
  2. in gzipped tar format

Two application routines are constructed that lead to the solution of sparse linear systems of algebraic equations with size (M+N) by (M+N). The operations to build the matrices and solve the equations use the defined operations for sparse matrices and the provided utilities to build the data structures.

There are two executables which may be generated using the makefile.