Cosmological N-body simulations are done on massively parallel computers. This necessitates the use of simple time integrators, and, additionally, of mesh-grid approximations of the potentials. Recently, Adamek et al. (2015); Barrera-Hinojosa et al. (2019) have developed general relativistic N-body simulations to capture relativistic effects mainly for cosmological purposes. We therefore ask whether, with the available technology, relativistic effects like perihelion advance can be detected numerically to a relevant precision. We first study the spurious perihelion shift in the Kepler problem, as a function of the integration method used, and then as a function of an additional interpolation of forces on a 2-dimensional lattice. This is done for several choices of eccentricities and semi-major axes. Using these results, we can predict which precisions and lattice constants allow for a detection of the relativistic perihelion advance in N-body simulation. We find that there are only small windows of parameters---such as eccentricity, distance from the central object and the Schwarzschild radius---for which the corrections can be detected in the numerics.

Estimates of the fractal dimension $D$ of the set of galaxies in the universe, based on ever improving data sets, tend to settle on $D\approx 2$. This result raised a raging debate due to its glaring contradiction with astrophysical models that expect a homogeneous universe. A recent mathematical result indicates that there is no contradiction, since measurements of the dimension of the {\em visible} subset of galaxies is bounded from above by D=2 even if the true dimension is anything between D=2 and D=3. We demonstrate this result in the context of a simple fractal model, and explain how to proceed in order to find a better estimate of the true dimension of the set of galaxies.

We introduce a concept of porosity for measures and study relations between dimensions and porosities for two classes of measures: measures on $R^n$ which satisfy the doubling condition and strongly porous measures on $R$.

We discuss various questions which arise when one considers the central projection of three dimensional fractal sets (galaxy catalogs) onto the celestial globe. The issues are related to how fractal such projections look. First we show that the lacunarity in the projection can be arbitrarily small. Further characteristics of the projected set---in particular scaling---depend sensitively on how the apparent sizes of galaxies are taken into account. Finally, we discuss the influence of opacity of galaxies. Combining these ideas, seemingly contradictory statements about lacunarity and apparent projections can be reconciled.