Parallel transport in shape analysis: a scalable numerical scheme