Maybe you missed it, but this point was already challenged and expanded upon in the previous thread. The library function that Fortran is using does indeed have on average better accuracy bounds, but it also has a significantly longer tail of low accuracy results than the algorithm that Bessels.jl is using.
This trade off is much more subtle, and it’s hard to say what’s ‘better’ in general without discussing a use-case.
