[metapost] bug in arctime calculation

[Hans Hagen]

On 8/16/2020 2:19 PM, luigi scarso wrote:

> Both binary and decimal uses and epsilon = 1E-52
> buth double has epsilon = 2^-52 and this matches with the binary search.
> Locally I have set  epsilon = 2^-52 for both binary and decimal
> and the results  are now:
I had been staring at that but wondered if

@d epsilonf                pow(2.0,-52.0)

was to be uised instead. There are a few more 'looking scaled related' 
constants, like the

fraction_threshold      0.04096

shouldn't these be more accurate?

Hans

-----------------------------------------------------------------
                                           Hans Hagen | PRAGMA ADE
               Ridderstraat 27 | 8061 GH Hasselt | The Netherlands
        tel: 038 477 53 69 | www.pragma-ade.nl | www.pragma-pod.nl
-----------------------------------------------------------------