Decomposition of stencil update formula into atomic stages

Abstract

In parallel solution of partial differential equations, a complex stencil update formula that accesses multiple layers of neighboring grid points sometimes must be decomposed into atomic stages, ones that access only immediately neighboring grid points. This paper shows that this requirement can be formulated as constraints of an optimization problem, which is equivalent to the dual of a minimum-cost network flow problem. An optimized decomposition of a single stencil on one set of grid points can thereby be computed efficiently.

References in corpus (3)

• The prospect of using LES and DES in engineering design, and the research required to get there
• The swept rule for breaking the latency barrier in time advancing PDEs
• The swept rule for breaking the latency barrier in time advancing two-dimensional PDEs

Cited by in corpus (2)

• Accelerating solutions of one-dimensional unsteady PDEs with GPU-based swept time-space decomposition
• The swept rule for breaking the latency barrier in time advancing two-dimensional PDEs