Hi Team,
gh-21609 is a Python translation of the 2-clause BSD-licensed getLebedevSphere, which computes weights and points for Lebedev quadrature over a spherical surface.
The PR followed the request in gh-11929, which was to add the function as scipy.special.roots_lebedev. This makes sense from the perspective that the interface and usage would be similar to the other scipy.special.roots_ functions; however, all of those functions involve the roots of orthogonal polynomials, and the Lebedev points don’t seem to fit this pattern.
The proposal now is to expose the function as scipy.integrate.lebedev_rule, since the primary purpose is for quadrature. We will need to consider how this fits alongside scipy.integrate.newton_cotes, which returns weights for Newton-Cotes quadrature, and the precedent it sets for other quadrature rule functions in scipy.integrate.
Please join the discussion in gh-21609!
Matt