I'm getting the somewhat expected result: sage: i=sqrt(-1) sage: for k in range(100): ....: print k; r1=random(); r2=random() ....: try: ....: xx=100*r1-i; yy=100*r2 ....: res=gamma(xx)*gamma(yy)/gamma(xx+yy) ....: except: ....: pass ....: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 *** overflow in t_REAL->double conversion. 17 18 19 20 ... and so on...
Try running your code in sage -gdb and then bt (backtrace) the stack. Since you apparently are running a binary distribution, I'd bet that there is a slight mismatch between the overflow signal handler on your distribution and on the distribution that it was compiled for. So when the overflow occurs it longjmp's into the abyss.... -- To post to this group, send email to [email protected] To unsubscribe from this group, send email to [email protected] For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
