Radial density profiles of the 3D Sedov blast wave test at t = 0.06. The solid black line indicates the analytic solution, while the simulation profiles are obtained by averaging for each radial bin SPH densities calculated from the particle distributions over a set of grid points located at the surface of a spherical shell and uniformly spaced in angular coordinates. All of the runs were performed using N = 524   288 particles.