Hi team, Thanks for the feedback on deprecating these functions. There was support for deprecating them, but also some concerns, so the plan has changed to the following: - Deprecate `scipy.integrate.quadrature` (so we can reclaim the name for a more general-purpose function in the future). - Introduce `scipy.integrate.gauss_quad` as an alias for `scipy.integrate.quadrature`, and point users to `gauss_quad` in the `quadrature` deprecation messages. - Mark `scipy.integrate.gauss_quad` and `scipy.integrate.romberg` as legacy functions. We welcome additional feedback in gh-19510 ( https://github.com/scipy/scipy/pull/19510). Thanks, Matt On Sat, Nov 11, 2023 at 4:10 PM Matt Haberland <matt.haberland@gmail.com> wrote:
Hi Team,
`scipy.integrate.romberg` and `scipy.integrate.quadrature` do not reliably hit their target accuracies, and when they do, they typically require more function evaluations and take more wall clock time than `scipy.integrate.quad` (all with default settings). This is documented in https://github.com/scipy/scipy/issues/18574#issuecomment-1703320598, which compares the accuracies, function call counts, and wall clock times on a benchmark set of ~800 integrals with known analytical solutions.
gh-19510 (https://github.com/scipy/scipy/pull/19510) proposes to deprecate these functions, suggesting `scipy.integrate.quad` as a drop-in replacement. Please let us know your thoughts either here or in gh-19510!
Thanks, Matt