Warm starting the projected Gauss-Seidel algorithm for granular matter simulation

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.