CPDES3: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions

Abstract Many physical problems require the solution of coupled partial differential equations on three-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils and allows for general... Title of program: CPDES3 Catalogue Id: ABFE_v1_0 Nature of problem Certain coupled elliptic and parabolic partial differential equations that arise in plasma physics and other applications are to be solved in three dimensions. The implicit solution techniques used for these equations give rise to a system of linear equations whose matrix operator is sparse (with a complicated subband structure) and generally asymmetric. We provide a fully vectorized algorithm for their solution. Versions of this program held in the CPC repository in Mendeley Data ABFE_v1_0; CPDES3; 10.1016/0010-4655(88)90153-1 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)