Add exponential-midpoint propagator for multi-system propagation
This merge request adds only the algorithmic steps for the exponential-midpoint propagator. The system specific implementations for the celestial particles and the Maxwell case together with their regression tests will follow in separate merge requests.
- I have checked that my code follows the Octopus coding standards
- I have added tests for all the new features added in this request.