Sorry for the typo but for example 1) listed in the previous email: gsl_sf_hyperg_U(1, -500, 0.1) gives the correct value: ~1.99e-3.
The problem happens when calculating gsl_sf_hyperg_U(1, -500, 40): gsl: gsl: hyperg.c:165: ERROR: overflow. The value should be ~0.00184 so I suppose this is a bug. Again, thanks very much if you can help me solve this problem since I very much need it in my research. Best Regards, Feng 2014-03-10 19:53 GMT-04:00 Feng Gao <[email protected]>: > Hi, > > I have been using the gsl_sf_hyperg_U function for a while and I found > several bugs when calculating normal values. For example, > > 1) gsl_sf_hyperg_U(1, -500, 0.1) gives incorrect value: 9.3246e+03, which > actually should be around 0.19. > > 2) Many domain error bugs. For example, gsl_sf_hyperg_U(1, -6.67, 1) > should have value ~0.128 but will invoke domain error for gsl. > > Hope that you can further investigate this and come up with solutions. > Thanks! > > Best Regards, > > Feng >
