TOMS565
PDETWO: Approximation of a PDE in 1 Time and 2 Space Dimensions

TOMS565 is a FORTRAN77 library which approximates a time-dependent 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:

1. David Melgaard, Richard Sincovec,
   Algorithm 565: PDETWO/PSETM/GEARB: Solution of systems of two-dimensional nonlinear partial differential equations,
   Volume 7, Number 1, March 1981, pages 126-135.

Source Code:

• toms565.f, the source code;
• toms565.sh, commands to compile the source code;

Examples and Tests:

Problem 1 is a simple elliptic equation.

• toms565_prb1.f, a sample calling program;
• toms565_prb1.sh, commands to compile, link and run the sample calling program;
• toms565_prb1_output.txt, the output file.

Problem 2 is the Burgers equation.

• toms565_prb2.f, a sample calling program;
• toms565_prb2.sh, commands to compile, link and run the sample calling program;
• toms565_prb2_output.txt, the output file.

Problem 3 is a coupled system of PDE's. This problem currently causes the time integrator to fail.

• toms565_prb3.f, a sample calling program;
• toms565_prb3.sh, commands to compile, link and run the sample calling program;
• toms565_prb3_output.txt, the output file.

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.