[SciPy-Dev] distributions.ncf.fit -- never converges?

[Yaroslav Halchenko]

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