Hi!
I would like to incorporate DiffentialEquations for our Accelerator Physicists particle tracking code. The set up is the following: ensembles of particles are moving through the sequence of magnets each having a length, some constant parameters and a Hamiltonian prescribing an ODE. The time like independent variable in this case is space position s, varying from 0 to total length of the accelerator. I wanted to clarify a few aspects of DifferentialEqations package before I start the development.
Are any symplectic methods available to execute on GPU for ensembles DiffEqGPU as it done with GPUTsit5 examples?
With the current set up I can see two possible ways to evolve the particles:
Invoking solve for each accelerator element (magnet) in a loop
Or use one large piecewise Hamiltonian to represent a whole accelerator and use a single solve call.
What is more beneficial, in terms of computational overhead?
Thank you!