constant matrix Y generation for global solve
Merge branch 'master' of https://github.com/rarietta/ProjectionBasedDynamics
Y generation complete