CPDES2: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in two dimensions
Description
Abstract Many physical problems require the solution of coupled partial differential equations on two-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. CPDES2 allows each spatial operator to have 5 or 9 point stencils and allows for general couplings b... Title of program: CPDES2 Catalogue Id: ABFD_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 two 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 ABFD_v1_0; CPDES2; 10.1016/0010-4655(88)90152-X This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)