JuliaMath / JuliaMath/DoubleFloats.jl

replacing `expm1` implementation

Open
#131 16 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
172
Forks
37
PR merge metrics
No merged PRs in 30d

Description

`expm1(x::Double64)` loses all digits of accuracy below `abs(x) < 1e-163` and returns `0.0` instead of `x`.

This regresses even beyond Float64, which correctly returns `x`.

```
using DoubleFloats
x = 1e-200
@assert x == expm1(x)
xx = Double64(x)
@assert xx == expm1(xx)
```

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.