[SciPy-User] Default tolerances for scipy.integrate.quadrature
Wed Apr 24 08:35:35 CDT 2013
On Wed, Apr 24, 2013 at 9:21 AM, Andrew Nelson <firstname.lastname@example.org> wrote:
> Dear list,
> I am using scipy.integrate.quadrature to do some adaptive gaussian
> quadrature. I am wondering about the default options for the
> tol=1.49e-8, rtol=1.49e-8
> I was wondering where these numbers come from? Are they related to
> some sort of hardware feature (e.g. related to machine precision)?
Yes. Those values are the square root of numpy.finfo(numpy.float64).eps.
eps is the smallest representable positive number such that 1.0 + eps !=
> Whilst I'm on this topic, I've noticed that a minimum number of
> adaption iterations for my problem using this function is
> approximately 10. I was wondering if it would be possible to add a
> 'startingOrder' keyword, so that the iteration could start from a 10th
> order integration? I don't mind trying to add this myself, if it
> would prove useful for scipy.
> Dr. Andrew Nelson
> SciPy-User mailing list
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the SciPy-User