[Numpy-discussion] `allclose` vs `assert_allclose`
Chris Barker
chris.barker at noaa.gov
Fri Jul 18 16:32:50 EDT 2014
On Fri, Jul 18, 2014 at 12:43 PM, Pauli Virtanen <pav at iki.fi> wrote:
> 18.07.2014 22:13, Chris Barker kirjoitti:
> [clip]
> > but an appropriate rtol would work there too. If only zero testing is
> > needed, then atol=0 makes sense as a default. (or maybe atol=eps)
>
> There's plenty of room below eps, but finfo(float).tiny ~ 3e-308 (or
> some big multiple) is also reasonable in the scale-freeness sense.
right! brain blip -- eps is the difference between 1 and then next larger
representable number, yes? So a long way away from smallest representable
number. So yes, zero or [something]e-308 -- making zero seem like a good
idea again....
is it totally ridiculous to have the default be dependent on dtype? float32
vs float64?
-Chris
--
Christopher Barker, Ph.D.
Oceanographer
Emergency Response Division
NOAA/NOS/OR&R (206) 526-6959 voice
7600 Sand Point Way NE (206) 526-6329 fax
Seattle, WA 98115 (206) 526-6317 main reception
Chris.Barker at noaa.gov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/numpy-discussion/attachments/20140718/c900403b/attachment.html>
More information about the NumPy-Discussion
mailing list