Hi all, I have noticed that many of the constants related to pi in the file sequences.jl of pkg DoubleFloats are flawed. For example, the DD values for (15/32)*%pi read:
"Double64(1.4726215563702156, -5.3616983752483455e-17),"
whereas my PY/MX/JL-scripts as well as the XL-addin XNumbers compute:
DDhi: +1.4726215563702156341e+00... , 3F F7 8F DB 9E FF EA 47
DDlo: -5.3616983752483473226e-17 ; BC 8E E8 78 60 39 D3 73 .
So in this example, DoubleFloats-DDlo is off by 3 ULPs. Any ideas??
Hello. You earned this
I am regenerating every constant in the repository.
(The current implementation does meet documented accuracies.)
The original values were transcribed from constants given by
one of the original authors. There was a note that many people
were writing to explain that some of the constants were inexact.
The note stated the values as written are as intended.
Any very low order mismatches were the result of subtle
misapplication of representation interconversion, and maybe
an incorrect assumption about the rounding techniques used
with each part.
I will replace the constants and rerun the accuracy tests.
3 Likes