last modified: 21-JUL-1992 | catalog | categories | new | search |

NESC9577 CPDES2.

CPDES2, Coupled 2-D Partial Differential Equation Solution

top ]
1. NAME OR DESIGNATION OF PROGRAM:  CPDES2.
top ]
2. COMPUTERS
To submit a request, click below on the link of the version you wish to order. Only liaison officers are authorised to submit online requests. Rules for requesters are available here.
Program name Package id Status Status date
CPDES2 NESC9577/01 Tested 21-JUL-1992

Machines used:

Package ID Orig. computer Test computer
NESC9577/01 CRAY 2 DEC VAX 8810
top ]
3. DESCRIPTION OF PROGRAM OR FUNCTION

CPDES2 solves the linear asymmetric equations arising from coupled partial differential equations in two dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils, permits general coupling between all of the component PDE's, and automatically generates the  matrix structures needed to perform the algorithm.
top ]
4. METHOD OF SOLUTION

The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned biconjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used.
top ]
5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

The discretization of the coupled two-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective.
top ]
6. TYPICAL RUNNING TIME

Execution timeS are problem dependent because ill-conditioned matrices require more iterations than well- conditioned ones. The sample problems converged to relative errors of 1.0E-10 in a minute or less on the Cray2.
top ]
7. UNUSUAL FEATURES OF THE PROGRAM

All inner loops are arranged to vectorize on the Cray2. A switch is included to allow the user to choose between the biconjugate (BCG) or conjugate (CG) gradient algorithm.
top ]
8. RELATED AND AUXILIARY PROGRAMS

The preconditioning technique is the same as that used in ILUCG2 (NESC 9929). CPDES3 (NESC 9576) is a preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled PDE's in three dimensions.
top ]
9. STATUS
Package ID Status date Status
NESC9577/01 21-JUL-1992 Screened
top ]
10. REFERENCES:
NESC9577/01, included references:
- D.V. Anderson, A.E. Koniges, and D.E. Shumaker:
  CPDES2 - A Preconditioned Conjugate Gradient Solver for Linear
  Asymmetric Matrix Equations Arising from Coupled Partial
  Differential Equations in Two Dimensions
  UCRL-96617 Preprint (May 1987).
top ]
11. MACHINE REQUIREMENTS

From 8*(n-1)*m to 32*n*m words of memory depending on the matrix sparsity pattern, where n is the number of partial differential equations and m is the number of linear equations.
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NESC9577/01 FORTRAN
top ]
13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  CTSS.
top ]
14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS

All COMMON and NAMELIST declarations are set in macros, and PARAMETER statements are used. The CPDES2 source is processed by a pre- compiler to obtain standard CFT FORTRAN source; the pre-compiler is  not included.
top ]
15. NAME AND ESTABLISHMENT OF AUTHORS

          D.V. Anderson, A.E. Koniges, D.E. Shumaker
          National MFE Computer Center
          Lawrence Livermore National Laboratory
top ]
16. MATERIAL AVAILABLE
NESC9577/01
File name File description Records
NESC9577_01.001 Information file 38
NESC9577_01.002 CPDES2 LRLTRAN source program 1080
top ]
17. CATEGORIES
  • P. General Mathematical and Computing System Routines

Keywords: algorithms, finite difference method, finite element method, magnetohydrodynamics, matrices, partial differential equations, two-dimensional.