A parallel algorithm for scattered data fitting on the sphere