compute the discrepancy of point sets on spheres