On Fri, Jul 18, 2014 at 12:43 PM, Pauli Virtanen <pav@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@noaa.gov