- Fixed constraint sign convention-related bugs in position projection.
- Preallocate some temps in the PLUS solver to avoid heap activity.
- Added code to get the velocity-level solver to produce acceleration-level values such as reaction forces and body accelerations. Issue: reaction forces are somewhat incorrect since they are saved mid-step but then finally evaluated after position changes.
- Minor comment tweaking here and there.