A c++ code for integrating particle trajectories. Composition schemes of 2nd, 4th, 6th, 8th and 10th order are implemented in it. It should be easily reusable.
The composition schemes are taken from the book "Geometric Numerical Integration", by E. Hairer, C. Lubich and G. Wanner.
By the way, some very interesting stuff related to particle trajectory integration can be found on Professor Hairer's webpage.
I've attached my report on implementing composition schemes for charged particles. There's nothing original in it as far as I know, but it brings together a few pieces of information, and it might be helpful for those interested.
This file is a report I made back in 2007, before I learned about composition schemes and their importance. I think it's a nice example of what can be done with very little programming. And it has chaos and fractals.
This file is an outline of a hierarchy of solvers that I constructed when I couldn't find something better. I implementend some of those a few times, for various tests, but I'm not sure I'll ever use them again. Anyway, I like the nice formulas that I obtained. I can't prove something for it so it's not a "finished product", but I'm posting everything here, maybe someone will find it useful.