Yaravi is a python code to solve the astrophysical N-body problem in arbitrary precision arithmetic, with a default target of the python float( 64-bit floating point) precision. It uses, like the Brutus code, an implementation of the Bulirsch-Stoer integrator with adaptive timestepping. It has crude support for parallelization.

It can be used through the AMUSE interface (included) or seperately if desired.

references: Boekholt & Portegies Zwart ,Computational Astrophysics and Cosmology, Volume 2, article id.2, 21 pp.

Note that (re)commit_particles calls generate a checkpoint in the code: the solution found is guaranteed to floating point precision from the last check point.