high performance computing on graphics processing units: hgpu.org

Field equations can be numerically simulated by approximating a continuous space field by a discrete lattice. There are a number of different lattice geometries that can be used to approximate continuous space which may cause structural artefacts in the simulation. These different lattice structures require the use of different stencil operators to approximate the spatial terms of the field equations. We show how different lattice geometries and associated stencil operators can be implemented in a stencil library which is used in conjunction with a code generator to produce code for field equations simulations that can run on a CPU or GPU.