Add a missing reduction for the number of grid points in orbital sets for +U
At the moment only the MPI root taks prints the number of grid points in the spherical mesh corresponding to orbital sets. This leads to a wrong display in domain parallelization.
Edited by Nicolas Tancogne-Dejean