@DanielSank 1. I can guess what it does but it would be nice if the readme/documentation told me explicitly ;) 2. Isn't the circle/area constraint dangerous because floating point equality is a horrible minefield?
In a non-orthogonal coordinate system, what is the physically significant difference between the components of a vector on the skew axes and its projection onto each axis? Why would one want to find one over the other?