Software Composability for Implicit Multiphysics with High Throughput
Thursday, February 23, 3:30PM – 5PM
Towards Algorithmic and Software Composability for Implicit >> Multiphysics with High Throughput
Implicit and IMEX solution methods for multiphysics problems are either monolithic (treating the coupled problem directly) or split (solving reduced systems independently). Software for monolithic multigrid is typically more intrusive, requiring a monolithic hierarchy instead of independent hierarchies, but the required number of iterations may be smaller due to coupling on all levels. We describe generic software support for "multigrid inside splitting" and "splitting inside multigrid" and provide comparisons for multiphysics applications in glaciology and geodynamics. The software for each "physics" component has no knowledge of its coupled environment, and indeed, can be used directly for single-physics simulation. In this process, we prefer high-order finite element methods while minimizing the use of assembled matrices, which exploits the nearly "free flops" available on recent hardware to deliver more accurate results with memory bandwidth and communication requirements.
Hosted by Robert Moser, Karl Schultz