code: def average_over_unit_sphere(self, quad=None):
          
          quad = quad or DEFAULT_QUAD
          weights, points = quad['weights'], quad['points']
          return sum([w * self.project(n) for w, n in zip(weights, points)])
docs: Method for averaging the tensor projection over the unit
          with option for custom quadrature.
  
          Args:
              quad (dict): quadrature for integration, should be
                  dictionary with "points" and "weights" keys defaults
                  to quadpy.sphere.Lebedev(19) as read from file
  
          Returns:
              Average of tensor projected into vectors on the unit sphere