21 Apr
2014
21 Apr
'14
9:24 a.m.
On 21 April 2014 09:38, John Lato
I was just wondering, why not simply numerically robust algorithms as defaults for these functions? No crashes, no errors, no loss of precision, everything would just work. They aren't particularly complicated, so the performance should even be reasonable.
I think it's best option. log1p and exp1m come with guarantees about precision. log(1+p) default makes it impossible to depend in such guarantees. They will silenly give wrong answer