Rebound

* Added new high order symplectic integrators from Wisdom et al. (1996) and Laskar & Robutel (2001). The implementation of these integrators are discussed in Rein, Tamayo & Brown (2019).
* Implemented new bit-wise comparison functions for simulations. Python syntax is simply sim1==sim2.
* Fixed a bug in IAS15 which prevented a restarted simulation to reproduce the original simulation exactly.

* Improves and fixes various issues related to variational equations and MEGNO.

* Fixes a bug which resulted in duplicate snapshots in Simulationarchives when restarting simulations.

* Syntax change on the python side to create a simulation from a binary file or Simulationarchive:

python
rebound.Simulation.from_file("test.bin") becomes rebound.Simulation("test.bin")
rebound.Simulation.from_archive("test.bin",5) becomes rebound.Simulation("test.bin",5)

* The hybrid integrator MERCURIUS has been completely rewritten. It can now much more easily be used in simulations where physical collisions occur. There are no more hidden particle arrays in the background, meaning adding and removing particles can occur in the same way as for other integrators. It also works reliably with any additional forces.
* The old hybrid integrator HERMES has been removed. MERCURIUS should always be equal or better in performance and accuracy.

* Added getBezierPaths to Simulationarchive to allow for easy plotting of complicated trajectories. To do this, store a lot of snapshots in the Simulationarchive (several per orbit!).
* Added functionality to add, subtract, multiply and divide simulations. This might be useful when developing new algorithms, but is most likely not useful for most users.