Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] aczl_v1_0.gz(10 Kbytes)
Manuscript Title: SIPSOL: a suite of subprograms for the solution of the linear equations arising from elliptic partial differential equations.
Authors: C.R. Jesshope
Program title: SIPSOL
Catalogue identifier: ACZL_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 17(1979)383
Programming language: Fortran.
Computer: CDC 7600.
Operating system: SCOPE 2.1.4.
RAM: 11K words
Word size: 60
Keywords: General purpose, Stones algorithm, Strongly implicit Procedure, Iterative, Elliptic, Partial, Differential equation, Matrix.
Classification: 4.3, 4.8.

Nature of problem:
SIPSOL is a suite of subprograms for iteratively inverting a set of linear equations obtained from the discretisation of an arbitrary elliptic partial differential equation over a rectangular grid of mesh points.

Solution method:
A strongly implicit procedure using Stone's algorithm.

In this version the mesh must be rectangular, although arbitrary boundary conditions may be specified. A five point difference operator must be used, giving a sparse five diagonal matrix system. The mesh may be non-uniform. A subprogram is supplied for calculating iteration parameters for the method, these are calculated for the general anisotropic diffusion equation; as a function of the mesh, or as a scalar.

Unusual features:
The surface output routine LIST in the test procedure uses A10 format (CDC 7600).

Running time:
For a mesh of 1024 points, a reduction in the L2 residual norm of 6 orders of magnitude can be obtained in 0.2-1.0 s on a CDC 7600.