optStiefel function

Optimize a function on the Stiefel manifold