hannorein / rebound

💫 An open-source multi-purpose N-body code.
https://rebound.readthedocs.io/
GNU General Public License v3.0
820 stars 217 forks source link

Fixed machine independence bug. #685

Closed PeterBartram closed 1 year ago

PeterBartram commented 1 year ago

Issue was the use of the math.pow function. Replaced it with ias15 version of x^1/7. Tested for bitwise energy violation matching for the longterm integration test is the unit test file (1e6 orbits for inner solar system) using a ubuntu desktop and a macbook M1.

hannorein commented 1 year ago

Nice catch! This shouldn't matter for any practical purpose, but it will help a lot with the unit tests!