[SciPy-Dev] severe bug in scipy.integrate.quadrature?

Pauli Virtanen pav@iki...
Wed Jul 28 17:51:04 CDT 2010


Wed, 28 Jul 2010 17:32:28 -0400, Thomas Robitaille wrote:
[clip: scipy.integrate.quadrature]
> The error inside the loop is the absolute error, so it looks like the
> tolerance to specify is the absolute tolerance, not relative. However,
> the original err is set to 100, so if I compute an integral that has a
> result of 1e20, I might want a tolerance of say 1e15, in which case the
> loop will be skipped, because err is not greater than tol. In its
> present form, I think that err should be initially set to infinity.

Yep, it's wrong. Fixed.

-- 
Pauli Virtanen



More information about the SciPy-Dev mailing list