Implement the derivative of the projection onto the set of unit vectors

This is needed for projection-based finite elements.

The code in this patch has not been tested yet.

[[Imported from SVN: r9887]]