问题
I have been learning about QUADPACK and I used scipy.integrate.quad to calculate an integral from 0 to infinite. It gave a very good result, but now I want to know which integration method (QAGI, QAWF,etc.) the software has applied. Is there any way for printing the technique? Does the software apply some decision tree? Thanks in advance for your time. Regards.
回答1:
As scipy
is open source, you can actually read the code for integrate.quad, which says:
For finite integration limits, the integration is performed using a Clenshaw-Curtis method which uses Chebyshev moments. ... If one of the integration limits is infinite, then a Fourier integral is computed (assuming w neq 0).
The call in _quad
for an infinite bound is to _quadpack._qagie
.
来源:https://stackoverflow.com/questions/21153436/how-to-know-integration-technique-used-by-quadpack