Iterative Poisson Solvers for Self-gravity with the GPU Code Astaroth
Fuente:
arXiv
Saved in:
| Main Authors: | , , , , , , , |
|---|---|
| Format: | Preprint |
| Published: |
2026
|
| Subjects: | |
| Online Access: | |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866915980809076736 |
|---|---|
| author | Krasnopolsky, Ruben Puro, Touko Li, Wei-Wen Shang, Hsien Väisälä, Miikka S. Mac Low, Mordecai-Mark Rheinhardt, Matthias Korpi-Lagg, Maarit |
| author_facet | Krasnopolsky, Ruben Puro, Touko Li, Wei-Wen Shang, Hsien Väisälä, Miikka S. Mac Low, Mordecai-Mark Rheinhardt, Matthias Korpi-Lagg, Maarit |
| contents | We present the development and benchmarking of Poisson solvers for graphics processing units (GPUs). Implemented in the Astaroth platform, the solvers feature high computational efficiency. We present novel combinations of discretizations and smoothers and document practical and performance-focused implementations aimed at reducing time-to-solution for self-gravitating systems. We describe the solver architectures and validate their accuracy against known analytic solutions. We measure convergence and timing per iteration for various solver algorithms, including conjugate gradient, successive overrelaxation, and multigrid in Cartesian coordinates, along with biconjugate gradient stabilized in spherical coordinates. We also couple the solvers to the Astaroth hydrodynamics to simulate a classic time-dependent problem in star formation, measuring accuracy and time-to-solution, for self-gravity on three-dimensional structured grids. Our results demonstrate that the solvers achieve performance similar to other algorithms implemented in Astaroth, and provide a solid foundation for integration into production-scale astrophysical simulations. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2605_03563 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Iterative Poisson Solvers for Self-gravity with the GPU Code Astaroth Krasnopolsky, Ruben Puro, Touko Li, Wei-Wen Shang, Hsien Väisälä, Miikka S. Mac Low, Mordecai-Mark Rheinhardt, Matthias Korpi-Lagg, Maarit Instrumentation and Methods for Astrophysics We present the development and benchmarking of Poisson solvers for graphics processing units (GPUs). Implemented in the Astaroth platform, the solvers feature high computational efficiency. We present novel combinations of discretizations and smoothers and document practical and performance-focused implementations aimed at reducing time-to-solution for self-gravitating systems. We describe the solver architectures and validate their accuracy against known analytic solutions. We measure convergence and timing per iteration for various solver algorithms, including conjugate gradient, successive overrelaxation, and multigrid in Cartesian coordinates, along with biconjugate gradient stabilized in spherical coordinates. We also couple the solvers to the Astaroth hydrodynamics to simulate a classic time-dependent problem in star formation, measuring accuracy and time-to-solution, for self-gravity on three-dimensional structured grids. Our results demonstrate that the solvers achieve performance similar to other algorithms implemented in Astaroth, and provide a solid foundation for integration into production-scale astrophysical simulations. |
| title | Iterative Poisson Solvers for Self-gravity with the GPU Code Astaroth |
| topic | Instrumentation and Methods for Astrophysics |
| url | https://arxiv.org/abs/2605.03563 |