Saved in:
| Main Authors: | , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2508.09568 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866908487726923776 |
|---|---|
| author | Haahr, Michael Haugbølle, Troels Nordlund, Åke Karlsson, Sven Richard, Eloi Martaillé |
| author_facet | Haahr, Michael Haugbølle, Troels Nordlund, Åke Karlsson, Sven Richard, Eloi Martaillé |
| contents | We present a GPU-accelerated implementation of a magnetohydrodynamic (MHD) solver using directive-based programming with OpenMP target offloading. The solver is integrated into the DISPATCH framework, which organises the computational domain into a collection of asynchronously updated patches. To reduce GPU kernel launch overhead, patches are grouped into "bunches" that are updated collectively. While porting the particular solver required a complete code refactoring, it yielded performance gains on both GPU and CPU. A stand-alone mini-app achieved a 7.3x speed-up compared to a single NVIDIA A100 GPU to seven AMD 7F72 Rome CPU cores. Within the full DISPATCH framework, the GPU-accelerated MHD Bunch solver showed excellent agreement with the CPU-based reference implementation on standard test problems such as the Sod shock tube and Orszag-Tang vortex. In large-scale 3D tests, the GPU implementation achieved a 9.8x overall speedup, comparing one GPU to 12 CPU cores, with the core MHD update routine being two orders of magnitude faster on the GPU than on a single CPU core. These results demonstrate that OpenMP offloading can provide substantial performance improvements for astrophysical codes while maintaining portability and accuracy. The work also demonstrates how new codes should be structured to allow simple and efficient directive-based GPU offloading. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_09568 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | GPU accelerated MHD in the DISPATCH framework using directive-based programming Haahr, Michael Haugbølle, Troels Nordlund, Åke Karlsson, Sven Richard, Eloi Martaillé Instrumentation and Methods for Astrophysics Computational Physics We present a GPU-accelerated implementation of a magnetohydrodynamic (MHD) solver using directive-based programming with OpenMP target offloading. The solver is integrated into the DISPATCH framework, which organises the computational domain into a collection of asynchronously updated patches. To reduce GPU kernel launch overhead, patches are grouped into "bunches" that are updated collectively. While porting the particular solver required a complete code refactoring, it yielded performance gains on both GPU and CPU. A stand-alone mini-app achieved a 7.3x speed-up compared to a single NVIDIA A100 GPU to seven AMD 7F72 Rome CPU cores. Within the full DISPATCH framework, the GPU-accelerated MHD Bunch solver showed excellent agreement with the CPU-based reference implementation on standard test problems such as the Sod shock tube and Orszag-Tang vortex. In large-scale 3D tests, the GPU implementation achieved a 9.8x overall speedup, comparing one GPU to 12 CPU cores, with the core MHD update routine being two orders of magnitude faster on the GPU than on a single CPU core. These results demonstrate that OpenMP offloading can provide substantial performance improvements for astrophysical codes while maintaining portability and accuracy. The work also demonstrates how new codes should be structured to allow simple and efficient directive-based GPU offloading. |
| title | GPU accelerated MHD in the DISPATCH framework using directive-based programming |
| topic | Instrumentation and Methods for Astrophysics Computational Physics |
| url | https://arxiv.org/abs/2508.09568 |