have also been combined with a lattice for deformation during hair styling [Gupta et al. 2006] or with an Eulerian fluid solver to give volume and provide a better initial position for the particle contacts [McAdams et al. 2009]. More recently, hair-body and hair-hair contacts have been more accurately modeled by using a non-smooth Newton solver for the Coulomb friction law [Bertails -Descoubes et al. 2011]. Daviet et al. [2011] globally solves for Coulomb friction with a hybrid Gauss-Seidel algorithm, using an analytic solver to ensure convergence.