Physics 410: Linear Systems
Please report all errors/typos. etc to
choptuik@physics.ubc.ca
Last updated October 8 2000
- Most of the source code covered in class can be found on
the phy410 account on sgi1 in
~phys410/linsys/ex1, ... (organized by lecture).
- General References
- Combined handouts (PS)
and slides (PS)
for the following lectures (does NOT include
Solving General Linear Banded Systems from Lecture 4).
- Lectures 1 & 2 (Oct 10 and 12)
- 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 sgi1.
Output on vnfe1.
Output on physics.
- Maple worksheet (PS)
which verifies above test.
- Lectures 2 & 3 (Oct 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 sgi1.
Plots showing level-8
solution and
error for levels 5, 6 and
7. Gnuplot "scripts"
gpsoln8
and
gperr,
for making solution and error plots.
- Lecture 4 (Oct 17)
- 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 sgi1.
Plots showing level-4
solution and
error for levels 4, 5 and
6.