Merge FastMath branch

  * Use `-ffast-math` to speed up ~all floating-point operations
  * Enable trapping of floating-point exceptions
  * Make MinStepsPerOrbit function parameter again
  * Update docs