now diagnostic only prints absolute value of center of mass and momentum; also, some improvements in leapfrog and (possibly) in hts (needs more tests).
added methods to get min/max/mean/harmonic_mean timestep.
first attempt to implement a PBaSS in HTS integrator. Some improvements are still needed.
added argparse option 'use_cl'.
more improvements in PBaSS implementation. It seems to work...
added methods to set/get the state of particles. added the 'n' attribute in Particles class so we don't have to updade it every time we want to know his value. also, psdf format is working again but improvements are still necessary.
sets default snap_freq=0 in leapfrog if kwargs is not given.
a new attempt to use t_i=1/sum(1/t_ij) in the timestep calculation together with HTS integrator. It seems to work now.
some small bug fixes.
refactoring of split/merge methods in HTS so that 'indexing' isn't necessary anymore. Also, include an attempt to substitute the recursive function by an iterative one using heapq module but this is not working yet.
added a method to finalize integrators. 'merge' method in HTS has been renamed to join.
instantiation of classes inheriting from Pbase looks much faster now. This speed up considerably the overall performance of HTS integrator but the recursive calls still slowdown the performance for deep levels of recursion.
timestep calculation is using t_i = min(t_ij) again. Also, adds several miscellaneous improvements in different modulos.
LeapFrog has been splited in a constant timestep and an adaptive timestep methods.
'models' has been renamed to 'ics'.
some files in 'ics' has been renamed.
added two new initial conditions for a few body systems.
removes some old files in ics.
added a class decorator which apply the provided decorator to all methods in a decorated class.
now all classes are decorated with timings, but it is not enabled by default.