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

NESC9576 CPDES3.

CPDES3, Coupled 3-D Partial Differential Equation Solution

top ]
1. NAME OR DESIGNATION OF PROGRAM:  CPDES3.
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
CPDES3 NESC9576/01 Tested 21-JUL-1992

Machines used:

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

CPDES3 solves the linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils, permits general couplings 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 three-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-conditionel matrices require more iterations than well- conditioned ones. The sample problems converged to relative errors of 1.0E-10 in 10 to 90 seconds 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 ILUCG3 (NESC 9927). CPDES2 (NESC 9577) is a preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled PDE's in two dimensions.
top ]
9. STATUS
Package ID Status date Status
NESC9576/01 21-JUL-1992 Screened
top ]
10. REFERENCES:
NESC9576/01, included references:
- D.V. Anderson, A.E. Koniges, and D.E. Shumaker:
  CPDES3 - A Preconditioned Conjugate Gradient Solver for Linear
  Asymmetric Matrix Equations Arising from Coupled Partial
  Differential Equations in Three Dimensions
  UCRL-96618 Preprint (May 1987).
top ]
11. MACHINE REQUIREMENTS

From 24*n*m to (256*n-76)*m words of the memory depending on the matrix sparsity pattern, where n is the number of coupled partial differential equations and m is the number of linear equations.
top ]
12. PROGRAMMING LANGUAGE(S) USED
Package ID Computer language
NESC9576/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 CPDES3 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
NESC9576/01
File name File description Records
NESC9576_01.001 Information file 38
NESC9576_01.002 CPDES3 LRLTRAN source program 1137
top ]
17. CATEGORIES
  • P. General Mathematical and Computing System Routines

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