TOMS565
PDETWO: Approximation of a PDE in 1 Time and 2 Space Dimensions
TOMS565
is a FORTRAN77 library which
approximates a timedependent partial differential equation (PDE)
in two spatial dimensions.
The text of many ACM TOMS algorithms is available online
through ACM:
http://www.acm.org/pubs/calgo
or NETLIB:
http://www.netlib.org/toms/index.html.
Languages:
TOMS565 is available in
a FORTRAN77 version.
Related Data and Software:
TOMS494
a FORTRAN77 library which
approximates a 1D PDE as a system of ODE's;
this library is commonly called PDEONE;
this is ACM TOMS algorithm 494.
Reference:

David Melgaard, Richard Sincovec,
Algorithm 565:
PDETWO/PSETM/GEARB: Solution of systems of twodimensional
nonlinear partial differential equations,
Volume 7, Number 1, March 1981, pages 126135.
Source Code:
Examples and Tests:
Problem 1 is a simple elliptic equation.
Problem 2 is the Burgers equation.
Problem 3 is a coupled system of PDE's. This problem currently
causes the time integrator to fail.
List of Routines:

PSETM generates the jacobian matrix.

PDETWO converts the user's PDE system into ODE's for the time integrator.

STRSET sets up storage.

DRIVEP is the driver program for the time integrator.

INTERP interpolates values of the dependent variable.

STIFFP takes one integration step.

COSET sets integration coefficients.

DECBR computes the LU decomposition of a banded matrix.

SOLBR solves a banded linear system.

PDB defines the jacobian if MITER = 1.
You can go up one level to
the FORTRAN77 source codes.
Last revised on 06 February 2011.