On Fri, Jul 18, 2014 at 2:32 PM, Chris Barker <[email protected]> wrote:
> 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? > > Whatever the final decision is, if the defaults change we should start with a FutureWarning. How we can make that work is uncertain, because I don't know of any reliable way to detect if we are using the default value or if a value was passed in. Maybe just warn if `atol == 0` ? Chuck
_______________________________________________ NumPy-Discussion mailing list [email protected] http://mail.scipy.org/mailman/listinfo/numpy-discussion
