Major changes in the way surfaces are handled. Now each possible contact is defined by a different surface, making it easier to mix global objects with surfaces that don't necessarily correspond to a physical object. Also, big improvements in the release solver. Now, releases are detected and solved by monitoring the continuous time Lagrange multipliers instead of the discrete multipliers. This is much easier to implement and gives a true one-step mapping for every step of the integrator.