[Bug-gsl] UNEXPECTED BUG in gsl_sf_lngamma_sgn_e !

[Nenad Sakan]

I had some problem with the hypergeometric functions that I use.
when I debug the program an run backtrace the error occures at the point 
before call to lngamma_sgn_sing (gamma.c: leatle before 1240). By all 
possible logic eps SHOULD NOT BE 0 at this point (x=0.0001139406357036744). 

N = -(int) x;
eps = x + N;

I started with casting of N as a double, but with no efect. 
After that I separated the rvalue of eps = ... in several lines but also with 
no efect also.

What is happening?

Please help!

Sincerely Yours,
Nenad Sakan
-------------------------------------
gdb backtrace :

#3  0x400a5bac in gsl_error (reason=0x6 <Address 0x6 out of bounds>, 
file=0x8049150 "@F6@", line=1077297136, gsl_errno=1)
    at error.c:47
#4  0x40166157 in lngamma_sgn_sing (N=1, eps=0, lng=0xbfffef20, 
sgn=0xbfffef00) at gamma.c:769
#5  0x40166543 in gsl_sf_lngamma_sgn_e (x=0.0001139406357036744, 
result_lg=0xbfffef20, sgn=0xbfffef00) at gamma.c:1240

-----------------
Institute of Physics, Belgrade
http://www.phy.bg.ac.yu/