Domain Decomposition for Elliptic Partial Differential Equations with Neumann Boundary Conditions

Discretization of a self-adjoint elliptic partial differential equation by finite differences or finite elements yields a large, sparse, symmetric system of equations, <em>Ax=b</em>. We use the preconditioned conjugate gradient method with domain decomposition to develop an effective, vectorizable preconditioner which is suitable for solving large two-dimensional problems on vector and parallel machines