Energy conserving (or symplectic) solvers for differential equations?

Specific Domains Modelling & Simulations
1

Hello there,
I am relatively new to Julia and am currently trying to calculate trajectories of a physical system on GPUs, which in generally works like a charm, as long as I do not include an minimum coupling term in my problem. In this case the energy conservation of the system breaks and the results are all wrong. So my question is, is there any symplectic / energy conserving integrator available on GPUs that I have missed, or are there other known ways to avoid this kind of problems ?

Already thanks in advance for any kind of input!

2

I don’t know if you are simulating atomic systems, but there might be something useful in Molly.jl.

3

Thanks for taking the time to answer. I found a solution in enforcing the energy conservation inside my definition of the ODE. That should work for now, but thanks for your answer!

1 Like