  The effect on the convergence of warm starting the projected Gauss-Seidel
solver for nonsmooth discrete element simulation of granular matter are
investigated. It is found that the computational performance can be increased
by a factor 2 to 5.
