http://gcc.gnu.org/bugzilla/show_bug.cgi?id=56216
Bug #: 56216 Summary: TR1 bessel functions bomb at x == 0! Classification: Unclassified Product: gcc Version: unknown Status: UNCONFIRMED Severity: normal Priority: P3 Component: libstdc++ AssignedTo: unassig...@gcc.gnu.org ReportedBy: 3dw...@verizon.net The Bessel functions std::tr1::cyl_bessel_i(double nu, double x) std::tr1::cyl_bessel_j(double nu, double x) return -nan at x=0 instead of 0 (nu==0) or 1 (nu > 0). I have a patch that fixes this and some other things.