On Fri, Jul 18, 2014 at 12:43 PM, Pauli Virtanen <[email protected]> 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 [email protected]
_______________________________________________ NumPy-Discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
