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