Combining CUDA and SYCL for Accelerating Earthquake Simulations in SeisSol

SeisSol is a high-performance simulation package designed for large-scale simulations of earthquakes and seismic phenomena. SeisSol is based on the discontinuous Galerkin method, the ADER time integration and Local Time Stepping. SeisSol uses a custom Domain Specific Language, called YATeTo, to tackle the demand of modeling various wave propagation models as well as providing portability between various CPU and GPU architectures from different vendors which have different capabilities and limitations. In this talk, we are going to report on 1) our recent results obtained with the SYCL compiler backend of YATeTo, 2) discuss benefits of using SYCL instead of OpenMP in the implementation of the dynamic rupture model in SeisSol and thus 3) coupling CUDA/HIP GPU programming models with SYCL.     «

SeisSol is a high-performance simulation package designed for large-scale simulations of earthquakes and seismic phenomena. SeisSol is based on the discontinuous Galerkin method, the ADER time integration and Local Time Stepping. SeisSol uses a custom Domain Specific Language, called YATeTo, to tackle the demand of modeling various wave propagation models as well as providing portability between various CPU and GPU architectures from different vendors which have different capabilities and limita...     »