Generating Optimal Initial Conditions for Smoothed Particle
  Hydrodynamics Simulations
  We review existing SPH setup methods and outline their advantages,
limitations and drawbacks. We present a new method for constructing initial
conditions for smoothed particle hydrodynamics (SPH) simulations, which may
also be of interest for N-body simulations, and demonstrate this method on a
number of applications. This new method is inspired by adaptive binning
techniques using weighted Voronoi tesselations. Particles are placed and
iteratively moved based on their proximity to neighboring particles and the
desired spatial resolution. This new method can satisfy arbitrarily complex
spatial resolution requirements.
