[SciPy-Dev] severe bug in scipy.integrate.quadrature?
Wed Jul 28 17:51:04 CDT 2010
Wed, 28 Jul 2010 17:32:28 -0400, Thomas Robitaille wrote:
> 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.
More information about the SciPy-Dev