just in case it would come useful in troubleshooting -- gets stuck in cumfnc_ (not modified since 2007) and blunt guess could be that reason is could be related to cdffnc1_wrap (dfn=-nan(0x8000000000000), dfd=-nan(0x8000000000000), ... ? *(gdb) bt full 10 #0 0x00007f69103a44eb in cumfnc_ () from /home/yoh/proj/scipy/install/lib/python2.6/site-packages/scipy/special/_cephes.so No symbol table info available. #1 0x00007f691039fab0 in cdffnc_ () from /home/yoh/proj/scipy/install/lib/python2.6/site-packages/scipy/special/_cephes.so No symbol table info available. #2 0x00007f6910363b84 in cdffnc1_wrap (dfn=-nan(0x8000000000000), dfd=-nan(0x8000000000000), nc=1, f=0.5) at scipy/special/cdf_wrappers.c:254 which = 1 q = 8.4886169591766247e-317 p = <value optimized out> bound = <value optimized out> status = 6 #3 0x00007f69103668fe in PyUFunc_dddd_d (args=0x7f691067a750, dimensions=<value optimized out>, steps=0x7ff80000, func=<value optimized out>) at scipy/special/ufunc_extras.c:510 i = 1 ip1 = 0xf6e2c0 "" ip2 = 0x19df1e0 "" ip3 = 0x19ee0c0 "" ip4 = 0x190c7b0 "" op = 0x17dedf0 "\b\357\374\023i\177" n = 10 #4 0x00007f6913257209 in ?? () from /usr/lib/pymodules/python2.6/numpy/core/umath.so No symbol table info available. #5 0x00007f6913257d13 in ?? () from /usr/lib/pymodules/python2.6/numpy/core/umath.so No symbol table info available. #6 0x000000000041ef47 in PyObject_Call (func=<numpy.ufunc at remote 0x14dba50>, arg=<unknown at remote 0x43680000>, kw=<unknown at remote 0x7ff80000>) at ../Objects/abstract.c:2492 result = < at remote 0xedd180> call = 0x7f6913257cb0 #7 0x00000000004a72b8 in do_call (f= Frame 0x19e4b60, for file /home/yoh/proj/scipy/install/lib/python2.6/site-packages/scipy/stats/distributions.py, line 4241, in _cdf (self=<ncf_gen(a=<float at remote 0x15f1a88>, moment_type=1, b=<float at remote 0xefa5d0>, vecfunc=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890640>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1891f50>, name='ncf', xb=<float at remote 0x1589d58>, xa=<float at remote 0x1589d70>, m=<float at remote 0x1589d28>, vecentropy=<vectorize(nin=4, nin_wo_defaults=None, thefunc=<instancemethod at remote 0x18906e0>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1891fd0>, numargs=3, generic_moment=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890fa0>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1893090>, veccdf=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890e10>, otypes='d', lastcallargs=0, ufunc=None, __doc_...(truncated), throwflag=<value optimized out>) at ../Python/ceval.c:3968 callargs = (<numpy.float64 at remote 0x1643710>, <numpy.float64 at remote 0x1553990>, <numpy.float64 at remote 0x18d1b80>, <numpy.ndarray at remote 0x19f4590>) kwdict = 0x0 #8 call_function (f= Frame 0x19e4b60, for file /home/yoh/proj/scipy/install/lib/python2.6/site-packages/scipy/stats/distributions.py, line 4241, in _cdf (self=<ncf_gen(a=<float at remote 0x15f1a88>, moment_type=1, b=<float at remote 0xefa5d0>, vecfunc=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890640>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1891f50>, name='ncf', xb=<float at remote 0x1589d58>, xa=<float at remote 0x1589d70>, m=<float at remote 0x1589d28>, vecentropy=<vectorize(nin=4, nin_wo_defaults=None, thefunc=<instancemethod at remote 0x18906e0>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1891fd0>, numargs=3, generic_moment=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890fa0>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1893090>, veccdf=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890e10>, otypes='d', lastcallargs=0, ufunc=None, __doc_...(truncated), throwflag=<value optimized out>) at ../Python/ceval.c:3773 func = <numpy.ufunc at remote 0x14dba50> w = <value optimized out> nk = <value optimized out> n = <value optimized out> pfunc = 0x19e4d00 #9 PyEval_EvalFrameEx (f= Frame 0x19e4b60, for file /home/yoh/proj/scipy/install/lib/python2.6/site-packages/scipy/stats/distributions.py, line 4241, in _cdf (self=<ncf_gen(a=<float at remote 0x15f1a88>, moment_type=1, b=<float at remote 0xefa5d0>, vecfunc=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890640>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1891f50>, name='ncf', xb=<float at remote 0x1589d58>, xa=<float at remote 0x1589d70>, m=<float at remote 0x1589d28>, vecentropy=<vectorize(nin=4, nin_wo_defaults=None, thefunc=<instancemethod at remote 0x18906e0>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1891fd0>, numargs=3, generic_moment=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890fa0>, otypes='d', lastcallargs=0, ufunc=None, __doc__=None, nout=None) at remote 0x1893090>, veccdf=<vectorize(nin=4, nin_wo_defaults=1, thefunc=<instancemethod at remote 0x1890e10>, otypes='d', lastcallargs=0, ufunc=None, __doc_...(truncated), throwflag=<value optimized out>) at ../Python/ceval.c:2412 sp = 0x19e4d08 stack_pointer = <value optimized out> next_instr = 0x1615229 "S" opcode = <value optimized out> oparg = <value optimized out> why = <value optimized out> err = <value optimized out> x = <value optimized out> v = <value optimized out> w = <value optimized out> u = <value optimized out> t = <value optimized out> stream = <value optimized out> freevars = 0x19e4d00 retval = 0x0 tstate = 0xedd180 co = 0x1614828 instr_ub = -1 instr_lb = 0 instr_prev = -1 first_instr = 0x1615214 "t" names = ('special', 'ncfdtr') consts = (None,) (More stack frames follow...) -- =------------------------------------------------------------------= Keep in touch www.onerussian.com Yaroslav Halchenko www.ohloh.net/accounts/yarikoptic