dzhang314 / MultiFloats.jl

Fast, SIMD-accelerated extended-precision arithmetic for Julia
MIT License
77 stars 10 forks source link

Better stopgaps for log1p and expm1 #19

Closed lrnv closed 1 year ago

lrnv commented 3 years ago

Hey,

Thanks a lot for what you did a while back about these stopgaps: this improved a lot my workflow and my runtime ! However i saw that the only two that i am calling these days are log1p and exp1m, and those might not need BigFloat at all.

This two functions are easier to remove from the BigFloat stopgaps than others, even if the algorithms are not the expected ones the result should still be OK, and way faster.

This should NOT be the final implementation: these functions exists because of catastrophic cancellation, and they should be re-implmeented properly. But in the mean time, these are way better than the BigFloat stopgaps.

lrnv commented 3 years ago

Thanks @oscardssmith :)