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.