Cosmic density and displacement field components as obtained with different methods. Upper panels: log-density slices at redshift z = 0 (log(1 + δ)) corresponding to different gravity solvers through the three-dimensional simulation cubical box of 200 h⁻¹ Mpc side length with a mesh of 256³ cells. Lower panels: Same as the upper panels but for the divergence of the irrotational (∇_q ⋅ Ψ_{∇ × Ψ = 0}(q)) and the absolute value of the solenoidal component of the displacement field, |Ψ_{∇ ⋅ Ψ = 0}(q)| (with power spectra within percentage agreement up to k ≃ 0.4 h Mpc⁻¹). The resulting maps have been saturated to the same colour scale after averaging over five slices corresponding to one simulation.