Coastal Ocean Simulations with Unstructured and Block-Structured Grids on CPUs, GPUs, and FPGAs

The shallow water equations model fluid flow in regions where the horizontal scale is much larger than the vertical depth. Applications of the shallow water equations include the modelling of tides, tsunamis or atmospheric flows. In this talk, we will discuss a SYCL implementation of a discontinuous Galerkin discretization for the two dimensional shallow water equations with application to coastal ocean simulations. We will first have a brief discussion of the discontinuous Galerkin discretization and important performance characteristics of the resulting scheme on unstructured grids. We will then show how this scheme can be implemented in SYCL, what abstractions are necessary for different architectures, and highlight differences betweencompilers and backends. At last, we briefly discuss energy efficiency for our code across the tested architectures and give a brief outlook on how block-structured grids can be used to improve simulation performance.