A geometric multigrid Poisson solver for domains containing solid inclusions

Abstract A Cartesian grid method for the fast solution of the Poisson equation in three-dimensional domains with embedded solid inclusions is presented and its performance analyzed. The efficiency of the method, which assume Neumann conditions at the immersed boundaries, is comparable to that of a multigrid method for regular domains. The method is light in terms of memory usage, and easily adaptable to parallel architectures. Tests with random and ordered arrays of solid inclusions, including spheres... Title of program: MG-Inc Catalogue Id: AEOE_v1_0 Nature of problem Poisson equation in domains containing inclusions; Neumann boundary conditions at immersed boundaries. Versions of this program held in the CPC repository in Mendeley Data AEOE_v1_0; MG-Inc; 10.1016/j.cpc.2012.11.008 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2018)