Beyond R-barycenters: an effective averaging method on Stiefel and Grassmann manifolds