The time-symmetric block time--step (TSBTS) algorithm is a newly developed efficient scheme for $N$--body integrations. It is constructed on an era-based iteration. In this work, we re-designed the TSBTS integration scheme with dynamically changing era size. A number of numerical tests were performed to show the importance of choosing the size of the era, especially for long time integrations. Our second aim was to show that the TSBTS scheme is as suitable as previously known schemes for developing parallel $N$--body codes. In this work, we relied on a parallel scheme using the copy algorithm for the time-symmetric scheme. We implemented a hybrid of data and task parallelization for force calculation to handle load balancing problems that can appear in practice. Using the Plummer model initial conditions for different numbers of particles, we obtained the expected efficiency and speedup for a small number of particles. Although parallelization of the direct $N$--body codes is negatively affected by the communication/calculation ratios, we obtained good load balance results. Moreover, we were able to conserve the advantages of the algorithm (e.g., energy conservation for long--term simulations).

The method of choice for integrating the equations of motion of the general N-body problem has been to use an individual time step scheme. For the sake of efficiency, block time steps have been the most popular, where all time step sizes are smaller than a maximum time step size by an integer power of two. We present the first successful attempt to construct a time-symmetric integration scheme, based on block time steps. We demonstrate how our scheme shows a vastly better long-time behavior of energy errors, in the form of a random walk rather than a linear drift. Increasing the number of particles makes the improvement even more pronounced.

Time-symmetric integration schemes share with symplectic schemes the property that their energy errors show a much better behavior than is the case for generic integration schemes. Allowing adaptive time steps typically leads to a loss of symplecticity. In contrast, time symmetry can be easily maintained, at least for a continuous choice of time step size. In large-scale N-body simulations, however, one often uses block time steps, where all time steps are forced to take on values as powers of two. This greatly facilitates parallelization, and hence code efficiency. Straightforward implementation of time-symmetry, translated to block time steps, faces significant hurdles. For example, iteration can lead to oscillatory behavior, and even when such behavior is suppressed, energy errors show a linear drift in time. We present an approach that circumvents these problems.