Further clarification of Facet code
Job may not be completely complete: I do not detail how contact kinematics are computed between a facet and a sphere when sphere center lies, after projection, outside of the facet. (I'm not planning on doing it in a near future, though. It took me long enough to see a connection between the code and these trilinear and barycentric coordinates)