We introduce a CUDA-Accelerated Timestepper for Alpha Particles Using Local Tricubics (CATAPULT) for use in Monte Carlo calculations of alpha particle confinement in stellarators. Our GPU implementation is significantly faster than existing parallelized CPU implementations, and handles both equilibrium magnetic fields and Shear Alfven Waves. We test our implementation on several example stellarators to exhibit both the speed and correctness of our code. The source code is included in the firm3d Python package.
@article{arxiv.2604.07617,
title = {CATAPULT: A CUDA-Accelerated Timestepper for Alpha Particles Using Local Tricubics},
author = {Michael Czekanski and Alexey R. Knyazev and David Bindel and Elizabeth J. Paul},
journal= {arXiv preprint arXiv:2604.07617},
year = {2026}
}