Random orthonormal matrix generation and optimization on the Stiefel manifold
