https://gcc.gnu.org/bugzilla/show_bug.cgi?id=105384
--- Comment #5 from Richard Biener <rguenth at gcc dot gnu.org> --- The bessel functions are known to take a lot of time, there are (fixed?) older bugreports about those and we do take some measures to speed up the computation. We might want to give up the idea on producing correctly rounded results and simply invoke the host library routines here ...