Physics 329: Linear Systems
Please report all errors/typos. etc to
matt@infeld.ph.utexas.edu
- Most of the source code covered in class can be found on
the phy329 account on einstein in
~phy329/linsys/ex1, ... (organized by lecture).
- General References
- Lectures 1 & 2 (Oct 5 and 7)
- Handout (PS) and
slides (PS).
- dgesv.f:
LAPACK driver routine for solving general linear system via
LU decomposition with partial pivoting.
- tdgesv.f: Test
program illustrating use of dgesv.
Makefile
and output on SGIs.
Output on Cray J90.
- Maple worksheet (PS)
which verifies above test.
- Lectures 3 and 4 (Oct 9 and 12)
- Handout (PS) and
slides (PS).
- dgtsv.f:
LAPACK driver routine for solving general tridiagonal system with
no pivoting.
- bvp1d.f: Test
program illustrating use of dgtsv to solve a 1d
boundary-value problem using second-order finite-difference techniques.
Makefile
and sample output on SGIs.
Plots showing level-8
solution and
error for levels 5, 6 and
7. Gnuplot "scripts"
gpsoln8
and
gperr,
for making solution and error plots.
- Lectures 4 and 5 (Oct 12 and 14)
- Handout: parts 1 (PS) and
2 (PS) and
slides: parts 1 (PS) and
2 (PS).
- dgbsv.f:
LAPACK driver routine for solving general banded system via
LU decomposition with partial pivoting.
- bvp1d4.f: Test
program illustrating use of dgbsv to solve a 1d
boundary-value problem using mixed second and fourth order
finite-difference techniques.
Makefile
and sample output on SGIs.
Plots showing level-4
solution and
error for levels 4, 5 and
6.