xexpy(0, 1000) should be zero
cossio opened this issue · comments
Jorge Fernandez-de-Cossio-Diaz commented
julia> xexpy(0, 1000)
NaN
I think xexpy
can handle this better and return zero?
Jorge Fernandez-de-Cossio-Diaz commented
David Widmann commented
Why should it be 0? At least the current implementation agrees with the docstring.
Jorge Fernandez-de-Cossio-Diaz commented
Mathematically, 0 * exp(1000) = 0
. Why should it not be zero?
Tamas K. Papp commented
I also think that we should make xexpy
finite if y
is finite and x
is zero. The purpose of writing a function is precisely to handle these kind of cases.
David Widmann commented
Fixed by #80.