JuliaLang/julia — comment on commit e8796226c4ba144bc371b9b120bb80d782d79c85
File:  line 0
By JeffBezanson:

Here are 3 cases where amos gets closer than fdlibm:

```
julia> lgamma(11)
15.104412573075516

julia> ccall(dlsym(_jl_libfdm,:lgamma),Float64,(Float64,),11.0)
15.104412573075518

julia> lgamma(pi)
0.8276945923234371

julia> ccall(dlsym(_jl_libfdm,:lgamma),Float64,(Float64,),pi)
0.827694592323437

julia> lgamma(50)
144.5657439463449

julia> ccall(dlsym(_jl_libfdm,:lgamma),Float64,(Float64,),50.0)
144.56574394634487
```

So far I've only tried `lgamma`.
libm doesn't have the airy function, but on that amos did better than NUMAL and cephes.
I believe the msun code for lgamma is identical to what we have in libfdm, but diff them to be sure.