A Method to Compute Averages over the Compact Stiefel Manifold