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/
_______________________________________________
Bug-gsl mailing list
[email protected]
http://lists.gnu.org/mailman/listinfo/bug-gsl
