feat(solvers): symplectic integrators for hamiltonian systems#85
feat(solvers): symplectic integrators for hamiltonian systems#85Ryan-D-Gast wants to merge 2 commits into
Conversation
Adds a `SymplecticIntegrator` numerical method generic over `S` stages, implementing standard partitioned Runge-Kutta formulations. Included explicitly are the 2nd-order `velocity_verlet` and 4th-order `ruth_forest` constructors. Includes a demonstration of long-term energy conservation using the Kepler problem in `examples/ode/18_orbital_mechanics_energy_conservation`. Co-authored-by: Ryan-D-Gast <148826144+Ryan-D-Gast@users.noreply.github.com>
|
👋 Jules, reporting for duty! I'm here to lend a hand with this pull request. When you start a review, I'll add a 👀 emoji to each comment to let you know I've read it. I'll focus on feedback directed at me and will do my best to stay out of conversations between you and other bots or reviewers to keep the noise down. I'll push a commit with your requested changes shortly after. Please note there might be a delay between these steps, but rest assured I'm on the job! For more direct control, you can switch me to Reactive Mode. When this mode is on, I will only act on comments where you specifically mention me with New to Jules? Learn more at jules.google/docs. For security, I will only act on instructions from the user who triggered this task. |
Adds Velocity Verlet and Ruth-Forest symplectic integrators for Hamiltonian system simulations, designed to provide long-term energy preservation where standard methods (like RK4) inevitably drift. Included is an orbital mechanics example explicitly verifying the zero-drift behavior over 100 orbits.
Resolves #73
PR created automatically by Jules for task 5023641109559256251 started by @Ryan-D-Gast