In order to extract the maximum amount of cosmological information from large scale structure surveys we need to make sure that we are able to model cosmological observables and their statistical errors at this unprecedented level of accuracy on the smallest scales that will be accessible to these surveys. This is where cosmological simulations are indispensable, since these scales are in the non-linear regime of gravitational collapse of cosmic structures, which cannot be modelled with analytical techniques.