StagBL: A Scalable, Portable, High-Performance Discretization and Solver Layer for Geodynamic Simulation

StagBL is an open-source parallel solver and discretization library for geodynamic simulation, encapsulating and optimizing operations essential to staggered-grid finite volume Stokes flow solvers. It provides a parallel staggered-grid abstraction with a high-level interface in C and Fortran. On top of this abstraction, tools are available to define boundary conditions and interact with particle systems.