question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

15 test failures for full/slow test suite run

See original GitHub issue

When running SCIPY_XSLOW=1 python dev.py -j4 -m full, I see the following test failures (Linux, conda-forge dev env created from environment_meson.yml):

================================= short test summary info =================================
FAILED scipy/linalg/tests/test_lapack.py::test_sgesdd_lwork_bug_workaround - pytest.Pyte...
FAILED scipy/sparse/linalg/tests/test_onenormest.py::TestOnenormest::test_onenormest_table_6_t_1
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_airyai_prime - Assertion...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_airybi - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_ci_complex - AssertionEr...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_igam_fac - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_legenp - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_sinpi - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_struvel - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_riemann_zeta - ValueErro...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_zetac - ValueError: zeta...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_spherical_in - Assertion...
FAILED scipy/stats/tests/test_sampling.py::TestNumericalInverseHermite::test_basic_all_scipy_dists[anglit-shapes1]
FAILED scipy/stats/tests/test_continuous_basic.py::test_moments[vonmises-arg102-True-True-True]
FAILED scipy/spatial/tests/test_hausdorff.py::test_massive_arr_overflow
=== 15 failed, 47544 passed, 2382 skipped, 237 xfailed, 11 xpassed in 633.30s (0:10:33) ===

Full tracebacks:

======================================== FAILURES =========================================
____________________________ test_sgesdd_lwork_bug_workaround _____________________________
[gw1] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/_pytest/runner.py:311: in from_call
    result: Optional[TResult] = func()
        cls        = <class '_pytest.runner.CallInfo'>
        duration   = 5.042588494994561
        excinfo    = <ExceptionInfo PytestUnraisableExceptionWarning('Exception ignored in: <function Popen.__del__ at 0x7f312897e8b0>\n\nT...   _warn("subprocess %s is still running" % self.pid,\nResourceWarning: subprocess 122546 is still running\n') tblen=6>
        func       = <function call_runtest_hook.<locals>.<lambda> at 0x7f30fc4823a0>
        precise_start = 58205.597220761
        precise_stop = 58210.639809256
        reraise    = (<class '_pytest.outcomes.Exit'>, <class 'KeyboardInterrupt'>)
        result     = None
        start      = 1641647516.2764997
        stop       = 1641647521.3190925
        when       = 'call'
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/_pytest/runner.py:255: in <lambda>
    lambda: ihook(item=item, **kwds), when=when, reraise=reraise
        ihook      = <_HookCaller 'pytest_runtest_call'>
        item       = <Function test_sgesdd_lwork_bug_workaround>
        kwds       = {}
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/pluggy/_hooks.py:265: in __call__
    return self._hookexec(self.name, self.get_hookimpls(), kwargs, firstresult)
        argname    = 'item'
        args       = ()
        firstresult = False
        kwargs     = {'item': <Function test_sgesdd_lwork_bug_workaround>}
        self       = <_HookCaller 'pytest_runtest_call'>
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/pluggy/_manager.py:80: in _hookexec
    return self._inner_hookexec(hook_name, methods, kwargs, firstresult)
        firstresult = False
        hook_name  = 'pytest_runtest_call'
        kwargs     = {'item': <Function test_sgesdd_lwork_bug_workaround>}
        methods    = [<HookImpl plugin_name='runner', plugin=<module '_pytest.runner' from '/home/rgommers/anaconda3/envs/scipy-dev/lib/pyt...readexception' from '/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/_pytest/threadexception.py'>>]
        self       = <_pytest.config.PytestPluginManager object at 0x7f31289ac070>
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/_pytest/unraisableexception.py:88: in pytest_runtest_call
    yield from unraisable_exception_runtest_hook()
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/_pytest/unraisableexception.py:78: in unraisable_exception_runtest_hook
    warnings.warn(pytest.PytestUnraisableExceptionWarning(msg))
E   pytest.PytestUnraisableExceptionWarning: Exception ignored in: <function Popen.__del__ at 0x7f312897e8b0>
E   
E   Traceback (most recent call last):
E     File "/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/subprocess.py", line 1052, in __del__
E       _warn("subprocess %s is still running" % self.pid,
E   ResourceWarning: subprocess 122546 is still running
        cm         = <_pytest.unraisableexception.catch_unraisable_exception object at 0x7f30fc46b340>
        err_msg    = 'Exception ignored in'
        msg        = 'Exception ignored in: <function Popen.__del__ at 0x7f312897e8b0>\n\nTraceback (most recent call last):\n  File "/home...__del__\n    _warn("subprocess %s is still running" % self.pid,\nResourceWarning: subprocess 122546 is still running\n'
_______________________ TestOnenormest.test_onenormest_table_6_t_1 ________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/sparse/linalg/tests/test_onenormest.py:183: in test_onenormest_table_6_t_1
    assert_(0.7 < proportion_exact < 0.8)
E   AssertionError
        A          = array([[ 1.0543604 +0.03847635j, -0.33208373-0.22861055j,
         0.40431944+0.0792976j , ..., -0.05681823-0.97139639...     -0.72327521+0.42464523j, ...,  0.86198444+1.35242015j,
        -1.4396941 -0.77624368j, -0.03960976-0.01487513j]])
        A_inv      = array([[0.29119074+0.65104849j, 0.58305885+0.67128938j,
        0.06722355+0.93823077j, ..., 0.6889593 +0.84052724j,
 ...
        0.45835081+0.39487713j, ..., 0.39236538+0.22490301j,
        0.37966209+0.48231849j, 0.52764845+0.90868658j]])
        est        = 286.6717944249122
        expected   = array([ 71.97806709,  70.56821698,  42.24911326, ...,  42.74919002,
        71.59204191, 287.08535285])
        i          = 4999
        itmax      = 5
        max_nresamples = 0
        n          = 100
        nexact     = 2781
        nmult_list = [5, 7, 4, 6, 7, 4, ...]
        nmults     = 4
        nresample_list = [0, 0, 0, 0, 0, 0, ...]
        nresamples = 0
        nsamples   = 5000
        observed   = array([ 71.97806709,  70.56821698,  35.59762277, ...,  42.74919002,
        63.85697738, 286.67179442])
        proportion_exact = 0.5562
        relative_errors = array([1.97433125e-16, 2.01377551e-16, 1.57435032e-01, ...,
       0.00000000e+00, 1.08043636e-01, 1.44054171e-03])
        self       = <scipy.sparse.linalg.tests.test_onenormest.TestOnenormest object at 0x7fa6a4b38d60>
        t          = 1
        underestimation_ratio = array([1.        , 1.        , 0.84256497, ..., 1.        , 0.89195636,
       0.99855946])
        underestimation_ratio_mean = 0.9414421444954721
        v          = array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 0., 0., 0., 0., 0., 0., 0., ...0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
       0., 0., 1., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])
        w          = array([ 1.93101609+4.88613457e-02j,  0.95209743+2.99680743e-01j,
        3.85029705-1.27986728e+00j,  1.89775519+4.326...578313 -9.43259458e-01j,  1.87828784-6.22941868e+00j,
        0.80356713+1.08910034e+00j, -2.96027481+1.61222120e+00j])
____________________________ TestSystematic.test_airyai_prime _____________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:696: in test_airyai_prime
    assert_mpmath_equal(lambda z: sc.airy(z)[1], lambda z:
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a44b9760>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<function TestSystematic.test_airyai_prime.<locals>.<lambda> at 0x7fa6a47144c0>, <function TestSystematic.test_airyai_prime.<locals>.<lambda> at 0x7fa6a4714d30>, [<scipy.special._mptestutils.Arg object at 0x7fa6a44d3e20>])
        d          = <MpmathData: <lambda>>
        kw         = {'rtol': 1e-05}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a4852d30>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a48521f0>
        self       = <MpmathData: <lambda>>
        tb         = <traceback object at 0x7fa6a43ecd80>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 6.71289e-07\nMax |rdiff|: 0.0001235\nBad results (1 out of 4647) for the following poin...628180958 =>           0.005436231808466818 !=           0.005435560519229373  (rdiff           0.000123499542516319)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a4852d30>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a48521f0>
        self       = <MpmathData: <lambda>>
        tb         = <traceback object at 0x7fa6a43ecd80>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 6.71289e-07\nMax |rdiff|: 0.0001235\nBad results (1 out of 4647) for the following poin...628180958 =>           0.005436231808466818 !=           0.005435560519229373  (rdiff           0.000123499542516319)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 6.71289e-07
E   Max |rdiff|: 0.0001235
E   Bad results (1 out of 4647) for the following points (in output 0):
E              -1014845.6628180958 =>           0.005436231808466818 !=           0.005435560519229373  (rdiff           0.000123499542516319)
        atol       = 1e-300
        data       = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for <lambda>>
        func       = <function TestSystematic.test_airyai_prime.<locals>.<lambda> at 0x7fa6a47144c0>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = True
        nparams    = 1
        param_filter = None
        points     = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a4852af0>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a4852af0>
        rtol       = 1e-05
        vectorized = False
_______________________________ TestSystematic.test_airybi ________________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:711: in test_airybi
    assert_mpmath_equal(lambda z: sc.airy(z)[2], lambda z:
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a445efa0>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<function TestSystematic.test_airybi.<locals>.<lambda> at 0x7fa6a47149d0>, <function TestSystematic.test_airybi.<locals>.<lambda> at 0x7fa6a46a8670>, [<scipy.special._mptestutils.Arg object at 0x7fa6a445e640>])
        d          = <MpmathData: <lambda>>
        kw         = {'rtol': 1e-05}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a46a8af0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a46a8040>
        self       = <MpmathData: <lambda>>
        tb         = <traceback object at 0x7fa6a4619400>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 6.66361e-10\nMax |rdiff|: 0.000123499\nBad results (1 out of 4290) for the following po...628180958 =>         -5.396327727623311e-06 !=        -5.3956613664474076e-06  (rdiff              0.000123499443469)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a46a8af0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a46a8040>
        self       = <MpmathData: <lambda>>
        tb         = <traceback object at 0x7fa6a4619400>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 6.66361e-10\nMax |rdiff|: 0.000123499\nBad results (1 out of 4290) for the following po...628180958 =>         -5.396327727623311e-06 !=        -5.3956613664474076e-06  (rdiff              0.000123499443469)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 6.66361e-10
E   Max |rdiff|: 0.000123499
E   Bad results (1 out of 4290) for the following points (in output 0):
E              -1014845.6628180958 =>         -5.396327727623311e-06 !=        -5.3956613664474076e-06  (rdiff              0.000123499443469)
        atol       = 1e-300
        data       = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for <lambda>>
        func       = <function TestSystematic.test_airybi.<locals>.<lambda> at 0x7fa6a47149d0>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = True
        nparams    = 1
        param_filter = None
        points     = array([[-1.00000000e+08],
       [-9.74540742e+07],
       [-9.49729657e+07],
       ...,
       [ 9.49651159e+07],
       [ 9.74500466e+07],
       [ 1.00000000e+08]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a46a8790>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a46a8790>
        rtol       = 1e-05
        vectorized = False
_____________________________ TestSystematic.test_ci_complex ______________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:955: in test_ci_complex
    assert_mpmath_equal(ci,
        ci         = <function TestSystematic.test_ci_complex.<locals>.ci at 0x7fa6a46a80d0>
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a4437760>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<function TestSystematic.test_ci_complex.<locals>.ci at 0x7fa6a46a80d0>, <bound method ci of <mpmath.ctx_mp.MPContext object at 0x7fa6bb4545e0>>, [<scipy.special._mptestutils.ComplexArg object at 0x7fa6a44375b0>])
        d          = <MpmathData: ci (complex)>
        kw         = {'rtol': 1e-08}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[-1.e+08-8.98846567e+307j],
       [-1.e+08-3.84127030e+269j],
       [-1.e+08-1.64158802e+231j],
       ...,
       [ 1.e+08+1.64158802e+231j],
       [ 1.e+08+3.84127030e+269j],
       [ 1.e+08+8.98846567e+307j]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a46a84c0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function mpc2complex at 0x7fa6b5e664c0>
        self       = <MpmathData: ci (complex)>
        tb         = <traceback object at 0x7fa6a4208440>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 1.90917e-16\nMax |rdiff|: 2.04926e-08\nBad results (3 out of 3849) for the following po...577e-09-2.220446049250313e-16j != 9.316390307435769e-09-7.39456652753901e-17j  (rdiff         1.5896601018469327e-08)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-1.e+08-8.98846567e+307j],
       [-1.e+08-3.84127030e+269j],
       [-1.e+08-1.64158802e+231j],
       ...,
       [ 1.e+08+1.64158802e+231j],
       [ 1.e+08+3.84127030e+269j],
       [ 1.e+08+8.98846567e+307j]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a46a84c0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function mpc2complex at 0x7fa6b5e664c0>
        self       = <MpmathData: ci (complex)>
        tb         = <traceback object at 0x7fa6a4208440>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 1.90917e-16\nMax |rdiff|: 2.04926e-08\nBad results (3 out of 3849) for the following po...577e-09-2.220446049250313e-16j != 9.316390307435769e-09-7.39456652753901e-17j  (rdiff         1.5896601018469327e-08)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 1.90917e-16
E   Max |rdiff|: 2.04926e-08
E   Bad results (3 out of 3849) for the following points (in output 0):
E     1.e+08-0.000262352790581753j => 9.316390628054619e-09+9.534595335480844e-13j != 9.316390628054619e-09+9.533509334190527e-13j  (rdiff         1.1656888681346863e-08)
E    1.e+08-1.307750747794633e-07j => 9.316390307435848e-09+6.661338147750939e-16j != 9.316390307435847e-09+4.752171230680926e-16j  (rdiff         2.0492560466753038e-08)
E     1.e+08+2.03491194526928e-08j => 9.31639030743577e-09-2.220446049250313e-16j != 9.316390307435769e-09-7.39456652753901e-17j  (rdiff         1.5896601018469327e-08)
        atol       = 1e-300
        data       = array([[-1.e+08-8.98846567e+307j],
       [-1.e+08-3.84127030e+269j],
       [-1.e+08-1.64158802e+231j],
       ...,
       [ 1.e+08+1.64158802e+231j],
       [ 1.e+08+3.84127030e+269j],
       [ 1.e+08+8.98846567e+307j]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for ci>
        func       = <function TestSystematic.test_ci_complex.<locals>.ci at 0x7fa6a46a80d0>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = True
        nparams    = 1
        param_filter = None
        points     = array([[-1.e+08-8.98846567e+307j],
       [-1.e+08-3.84127030e+269j],
       [-1.e+08-1.64158802e+231j],
       ...,
       [ 1.e+08+1.64158802e+231j],
       [ 1.e+08+3.84127030e+269j],
       [ 1.e+08+8.98846567e+307j]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a46a8f70>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a46a8f70>
        rtol       = 1e-08
        vectorized = False
______________________________ TestSystematic.test_igam_fac _______________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1384: in test_igam_fac
    assert_mpmath_equal(_igam_fac,
        mp_igam_fac = <function TestSystematic.test_igam_fac.<locals>.mp_igam_fac at 0x7fa6a4316040>
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6c4ae4460>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<ufunc '_igam_fac'>, <function TestSystematic.test_igam_fac.<locals>.mp_igam_fac at 0x7fa6a4316040>, [<scipy.special._mptestutils.Arg object at 0x7fa6c4aee100>, <scipy.special._mptestutils.Arg object at 0x7fa6c4aee8e0>])
        d          = <MpmathData: _igam_fac>
        kw         = {'rtol': 1e-10}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[1.00000000e-30, 0.00000000e+00],
       [1.00000000e-30, 1.00000000e-30],
       [1.00000000e-30, 5.39488518e-...   [1.00000000e+14, 2.95520924e+12],
       [1.00000000e+14, 1.71907220e+13],
       [1.00000000e+14, 1.00000000e+14]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a4316670>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a4316790>
        self       = <MpmathData: _igam_fac>
        tb         = <traceback object at 0x7fa6b32eb840>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 0.000206561\nMax |rdiff|: 1.24879e-10\nBad results (1 out of 5041) for the following po...58571e+13 =>             1654082.3391163372 !=             1654082.3389097766  (rdiff         1.2487926997466973e-10)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[1.00000000e-30, 0.00000000e+00],
       [1.00000000e-30, 1.00000000e-30],
       [1.00000000e-30, 5.39488518e-...   [1.00000000e+14, 2.95520924e+12],
       [1.00000000e+14, 1.71907220e+13],
       [1.00000000e+14, 1.00000000e+14]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a4316670>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a4316790>
        self       = <MpmathData: _igam_fac>
        tb         = <traceback object at 0x7fa6b32eb840>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 0.000206561\nMax |rdiff|: 1.24879e-10\nBad results (1 out of 5041) for the following po...58571e+13 =>             1654082.3391163372 !=             1654082.3389097766  (rdiff         1.2487926997466973e-10)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 0.000206561
E   Max |rdiff|: 1.24879e-10
E   Bad results (1 out of 5041) for the following points (in output 0):
E            1.719072201858571e+13           1.719072201858571e+13 =>             1654082.3391163372 !=             1654082.3389097766  (rdiff         1.2487926997466973e-10)
        atol       = 1e-300
        data       = array([[1.00000000e-30, 0.00000000e+00],
       [1.00000000e-30, 1.00000000e-30],
       [1.00000000e-30, 5.39488518e-...   [1.00000000e+14, 2.95520924e+12],
       [1.00000000e+14, 1.71907220e+13],
       [1.00000000e+14, 1.00000000e+14]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for _igam_fac>
        func       = <ufunc '_igam_fac'>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = True
        nparams    = 2
        param_filter = None
        points     = array([[1.00000000e-30, 0.00000000e+00],
       [1.00000000e-30, 1.00000000e-30],
       [1.00000000e-30, 5.39488518e-...   [1.00000000e+14, 2.95520924e+12],
       [1.00000000e+14, 1.71907220e+13],
       [1.00000000e+14, 1.00000000e+14]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a43169d0>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a43169d0>
        rtol       = 1e-10
        vectorized = False
_______________________________ TestSystematic.test_legenp ________________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1547: in test_legenp
    assert_mpmath_equal(lpnm,
        legenp     = <function TestSystematic.test_legenp.<locals>.legenp at 0x7fa6a4852790>
        lpnm       = <function TestSystematic.test_legenp.<locals>.lpnm at 0x7fa6a4852820>
        lpnm_2     = <function TestSystematic.test_legenp.<locals>.lpnm_2 at 0x7fa6a4852670>
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a565f5e0>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<function TestSystematic.test_legenp.<locals>.lpnm at 0x7fa6a4852820>, <function TestSystematic.test_legenp.<locals>....cipy.special._mptestutils.IntArg object at 0x7fa6a56b6bb0>, <scipy.special._mptestutils.Arg object at 0x7fa6a56b6d30>])
        d          = <MpmathData: lpnm>
        kw         = {}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[-1.00000000e+002, -1.00000000e+002, -8.98846567e+307],
       [-1.00000000e+002, -1.00000000e+002, -2.99807700...[ 1.00000000e+001,  1.00000000e+001,  1.00000000e+001],
       [ 1.00000000e+001,  1.00000000e+001,  8.98846567e+307]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a46a8f70>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a46a8af0>
        self       = <MpmathData: lpnm>
        tb         = <traceback object at 0x7fa6a4ae2740>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 0\nMax |rdiff|: 0\nBad results (2 out of 664) for the following points (in output 0):\n...3464e+154 =>                           -inf !=       -2.9980769960613464e+154  (rdiff                             0.)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-1.00000000e+002, -1.00000000e+002, -8.98846567e+307],
       [-1.00000000e+002, -1.00000000e+002, -2.99807700...[ 1.00000000e+001,  1.00000000e+001,  1.00000000e+001],
       [ 1.00000000e+001,  1.00000000e+001,  8.98846567e+307]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a46a8f70>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a46a8af0>
        self       = <MpmathData: lpnm>
        tb         = <traceback object at 0x7fa6a4ae2740>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 0\nMax |rdiff|: 0\nBad results (2 out of 664) for the following points (in output 0):\n...3464e+154 =>                           -inf !=       -2.9980769960613464e+154  (rdiff                             0.)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 0
E   Max |rdiff|: 0
E   Bad results (2 out of 664) for the following points (in output 0):
E                               1.                             -1.        -2.9980769960613464e+154 =>                           -inf !=       -1.4990384980306732e+154  (rdiff                             0.)
E                               1.                              1.        -2.9980769960613464e+154 =>                           -inf !=       -2.9980769960613464e+154  (rdiff                             0.)
        atol       = 1e-300
        data       = array([[-1.00000000e+002, -1.00000000e+002, -8.98846567e+307],
       [-1.00000000e+002, -1.00000000e+002, -2.99807700...[ 1.00000000e+001,  1.00000000e+001,  1.00000000e+001],
       [ 1.00000000e+001,  1.00000000e+001,  8.98846567e+307]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for lpnm>
        func       = <function TestSystematic.test_legenp.<locals>.lpnm at 0x7fa6a4852820>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = True
        nparams    = 3
        param_filter = None
        points     = array([[-1.00000000e+002, -1.00000000e+002, -8.98846567e+307],
       [-1.00000000e+002, -1.00000000e+002, -2.99807700...[ 1.00000000e+001,  1.00000000e+001,  1.00000000e+001],
       [ 1.00000000e+001,  1.00000000e+001,  8.98846567e+307]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a46a8a60>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a46a8a60>
        rtol       = 1e-07
        vectorized = False
________________________________ TestSystematic.test_sinpi ________________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1743: in test_sinpi
    assert_mpmath_equal(_sinpi, mpmath.sinpi,
        eps        = 2.220446049250313e-16
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a4bf1160>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<ufunc '_sinpi'>, <function PythonMPContext._wrap_libmp_function.<locals>.f at 0x7fa6bb478c10>, [<scipy.special._mptestutils.Arg object at 0x7fa6a4beae80>])
        d          = <MpmathData: _sinpi>
        kw         = {'nan_ok': False, 'rtol': 2.220446049250313e-16}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[-8.98846567e+307],
       [-2.90107975e+307],
       [-9.36340418e+306],
       ...,
       [ 9.32952750e+306],
       [ 2.89582696e+307],
       [ 8.98846567e+307]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6b3a295e0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a4305c10>
        self       = <MpmathData: _sinpi>
        tb         = <traceback object at 0x7fa6b3930780>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 1.11022e-16\nMax |rdiff|: 2.27509e-16\nBad results (2 out of 5001) for the following po...     0.92 =>            0.24868988716485463 !=            0.24868988716485468  (rdiff          2.232143488587451e-16)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-8.98846567e+307],
       [-2.90107975e+307],
       [-9.36340418e+306],
       ...,
       [ 9.32952750e+306],
       [ 2.89582696e+307],
       [ 8.98846567e+307]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6b3a295e0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a4305c10>
        self       = <MpmathData: _sinpi>
        tb         = <traceback object at 0x7fa6b3930780>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 1.11022e-16\nMax |rdiff|: 2.27509e-16\nBad results (2 out of 5001) for the following po...     0.92 =>            0.24868988716485463 !=            0.24868988716485468  (rdiff          2.232143488587451e-16)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 1.11022e-16
E   Max |rdiff|: 2.27509e-16
E   Bad results (2 out of 5001) for the following points (in output 0):
E              -2.8377298161470823 =>           -0.48799118192566504 !=           -0.48799118192566515  (rdiff         2.2750882920549882e-16)
E                             0.92 =>            0.24868988716485463 !=            0.24868988716485468  (rdiff          2.232143488587451e-16)
        atol       = 1e-300
        data       = array([[-8.98846567e+307],
       [-2.90107975e+307],
       [-9.36340418e+306],
       ...,
       [ 9.32952750e+306],
       [ 2.89582696e+307],
       [ 8.98846567e+307]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for _sinpi>
        func       = <ufunc '_sinpi'>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = False
        nparams    = 1
        param_filter = None
        points     = array([[-8.98846567e+307],
       [-2.90107975e+307],
       [-9.36340418e+306],
       ...,
       [ 9.32952750e+306],
       [ 2.89582696e+307],
       [ 8.98846567e+307]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a4305ee0>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a4305ee0>
        rtol       = 2.220446049250313e-16
        vectorized = False
_______________________________ TestSystematic.test_struvel _______________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1826: in test_struvel
    assert_mpmath_equal(sc.modstruve,
        mp_struvel = <function TestSystematic.test_struvel.<locals>.mp_struvel at 0x7fa6a40f63a0>
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6b2e9d040>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<ufunc 'modstruve'>, <function exception_to_nan.<locals>.wrap at 0x7fa6a40f6700>, [<scipy.special._mptestutils.Arg object at 0x7fa6a564a100>, <scipy.special._mptestutils.Arg object at 0x7fa6a564a400>])
        d          = <MpmathData: modstruve>
        kw         = {'ignore_inf_sign': True, 'rtol': 5e-10}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[-1.00000000e+04,  0.00000000e+00],
       [-1.00000000e+04,  1.00000000e-30],
       [-1.00000000e+04,  5.3948....00000000e+04,  4.43668733e+03],
       [ 1.00000000e+04,  6.66084629e+03],
       [ 1.00000000e+04,  1.00000000e+04]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a40f6af0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a40f6280>
        self       = <MpmathData: modstruve>
        tb         = <traceback object at 0x7fa6a42f8b80>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 2.93606e+121\nMax |rdiff|: 6.43922e+35\nBad results (1 out of 5036) for the following p...235202888 =>        2.9360586115244534e+121 !=         4.5596491027046297e+85  (rdiff          6.439220530770411e+35)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-1.00000000e+04,  0.00000000e+00],
       [-1.00000000e+04,  1.00000000e-30],
       [-1.00000000e+04,  5.3948....00000000e+04,  4.43668733e+03],
       [ 1.00000000e+04,  6.66084629e+03],
       [ 1.00000000e+04,  1.00000000e+04]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a40f6af0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a40f6280>
        self       = <MpmathData: modstruve>
        tb         = <traceback object at 0x7fa6a42f8b80>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 2.93606e+121\nMax |rdiff|: 6.43922e+35\nBad results (1 out of 5036) for the following p...235202888 =>        2.9360586115244534e+121 !=         4.5596491027046297e+85  (rdiff          6.439220530770411e+35)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 2.93606e+121
E   Max |rdiff|: 6.43922e+35
E   Bad results (1 out of 5036) for the following points (in output 0):
E               -4216.965034285822               2955.209235202888 =>        2.9360586115244534e+121 !=         4.5596491027046297e+85  (rdiff          6.439220530770411e+35)
        atol       = 1e-300
        data       = array([[-1.00000000e+04,  0.00000000e+00],
       [-1.00000000e+04,  1.00000000e-30],
       [-1.00000000e+04,  5.3948....00000000e+04,  4.43668733e+03],
       [ 1.00000000e+04,  6.66084629e+03],
       [ 1.00000000e+04,  1.00000000e+04]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for modstruve>
        func       = <ufunc 'modstruve'>
        ignore_inf_sign = True
        knownfailure = None
        nan_ok     = True
        nparams    = 2
        param_filter = None
        points     = array([[-1.00000000e+04,  0.00000000e+00],
       [-1.00000000e+04,  1.00000000e-30],
       [-1.00000000e+04,  5.3948....00000000e+04,  4.43668733e+03],
       [ 1.00000000e+04,  6.66084629e+03],
       [ 1.00000000e+04,  1.00000000e+04]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a40f65e0>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a40f65e0>
        rtol       = 5e-10
        vectorized = False
____________________________ TestSystematic.test_riemann_zeta _____________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1860: in test_riemann_zeta
    assert_mpmath_equal(
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a41c2430>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<function zeta at 0x7fa6c35544c0>, <bound method zeta of <mpmath.ctx_mp.MPContext object at 0x7fa6bb4545e0>>, [<scipy.special._mptestutils.Arg object at 0x7fa6a5637c70>])
        d          = <MpmathData: zeta>
        kw         = {'nan_ok': False, 'rtol': 1e-13}
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        dps        = 20
        dps_list   = [20]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a40f6b80>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a40f64c0>
        self       = <MpmathData: zeta>
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
        atol       = 1e-300
        data       = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for zeta>
        func       = <function zeta at 0x7fa6c35544c0>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = False
        nparams    = 1
        param_filter = None
        points     = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a40f6c10>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a40f6c10>
        rtol       = 1e-13
        vectorized = False
scipy/special/_testutils.py:239: in check
    wanted = eval_func_at_params(self.result_func, skip_mask=skip_mask)
        __tracebackhide__ = operator.methodcaller('errisinstance', <class 'AssertionError'>)
        atol       = 1e-300
        data       = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        dtype      = dtype('float64')
        dtypes     = None
        eval_func_at_params = <function FuncData.check.<locals>.eval_func_at_params at 0x7fa6a40f6d30>
        filter     = None
        got        = (array([0.00000000e+00, 8.90512061e+76, 5.40918573e+76, ...,
       1.00000000e+00, 1.00000000e+00, 1.00000000e+00]),)
        idx        = 0
        j          = 0
        param_mask = array([ True,  True,  True, ...,  True,  True,  True])
        params     = [array([-100.        ,  -99.6322642 ,  -99.26588069, ...,   99.26470855,
         99.63167596,  100.        ])]
        rtol       = 1e-13
        self       = <Data for zeta>
        skip_mask  = None
scipy/special/_testutils.py:220: in eval_func_at_params
    got.append(func(*tuple([params[i][j] for i in range(len(params))])))
        func       = <function MpmathData.check.<locals>.<lambda> at 0x7fa6a40f6c10>
        got        = [0.0, 8.905120611274578e+76, 5.409185726103669e+76, 2.1321986264121143e+76, 5.9110460840794665e+75, 7.90241504009859e+74, ...]
        j          = 3251
        params     = [array([-100.        ,  -99.6322642 ,  -99.26588069, ...,   99.26470855,
         99.63167596,  100.        ])]
        self       = <Data for zeta>
        skip_mask  = None
scipy/special/_mptestutils.py:259: in <lambda>
    lambda *a: pytype(self.mpmath_func(*map(mptype, a))),
        a          = (1.0,)
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a40f6b80>
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a40f64c0>
        self       = <MpmathData: zeta>
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/functions/zeta.py:533: in zeta
    return ctx._zeta(s, **kwargs)
        a          = 1
        ctx        = <mpmath.ctx_mp.MPContext object at 0x7fa6bb4545e0>
        d          = 0
        derivative = 0
        kwargs     = {}
        method     = None
        s          = mpf('1.0')
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/ctx_mp_python.py:998: in f
    return ctx.make_mpf(mpf_f(x._mpf_, prec, rounding))
        ctx        = <mpmath.ctx_mp.MPContext object at 0x7fa6bb4545e0>
        kwargs     = {}
        mpc_f      = <function mpc_zeta at 0x7fa6bb5ac430>
        mpf_f      = <function mpf_zeta at 0x7fa6bb5ac3a0>
        name       = 'zeta'
        prec       = 70
        rounding   = 'n'
        x          = mpf('1.0')
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/libmp/gammazeta.py:984: in mpf_zeta
    return mpf_zeta_int(to_int(s), prec, rnd)
        alt        = 0
        bc         = 1
        exp        = 0
        man        = mpz(1)
        prec       = 70
        rnd        = 'n'
        s          = (0, mpz(1), 0, 1)
        sign       = 0
        wp         = 90
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/libmp/gammazeta.py:917: in mpf_zeta_int
    raise ValueError("zeta(1) pole")
E   ValueError: zeta(1) pole
        prec       = 70
        rnd        = 'n'
        s          = 1
        wp         = 90
________________________________ TestSystematic.test_zetac ________________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1869: in test_zetac
    assert_mpmath_equal(sc.zetac,
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa6a40857f0>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<ufunc 'zetac'>, <function TestSystematic.test_zetac.<locals>.<lambda> at 0x7fa6a40f6820>, [<scipy.special._mptestutils.Arg object at 0x7fa6a40858b0>])
        d          = <MpmathData: zetac>
        kw         = {'dps': 45, 'nan_ok': False, 'rtol': 1e-13}
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        dps        = 45
        dps_list   = [45]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a42210d0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a42214c0>
        self       = <MpmathData: zetac>
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
        atol       = 1e-300
        data       = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for zetac>
        func       = <ufunc 'zetac'>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = False
        nparams    = 1
        param_filter = None
        points     = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa66ff9b940>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa66ff9b940>
        rtol       = 1e-13
        vectorized = False
scipy/special/_testutils.py:239: in check
    wanted = eval_func_at_params(self.result_func, skip_mask=skip_mask)
        __tracebackhide__ = operator.methodcaller('errisinstance', <class 'AssertionError'>)
        atol       = 1e-300
        data       = array([[-100.        ],
       [ -99.6322642 ],
       [ -99.26588069],
       ...,
       [  99.26470855],
       [  99.63167596],
       [ 100.        ]])
        dtype      = dtype('float64')
        dtypes     = None
        eval_func_at_params = <function FuncData.check.<locals>.eval_func_at_params at 0x7fa66ff9b670>
        filter     = None
        got        = (array([-1.00000000e+00,  8.90512061e+76,  5.40918573e+76, ...,
        1.31324337e-30,  1.01830299e-30,  7.88860905e-31]),)
        idx        = 0
        j          = 0
        param_mask = array([ True,  True,  True, ...,  True,  True,  True])
        params     = [array([-100.        ,  -99.6322642 ,  -99.26588069, ...,   99.26470855,
         99.63167596,  100.        ])]
        rtol       = 1e-13
        self       = <Data for zetac>
        skip_mask  = None
scipy/special/_testutils.py:220: in eval_func_at_params
    got.append(func(*tuple([params[i][j] for i in range(len(params))])))
        func       = <function MpmathData.check.<locals>.<lambda> at 0x7fa66ff9b940>
        got        = [-1.0, 8.905120611274578e+76, 5.409185726103669e+76, 2.1321986264121143e+76, 5.9110460840794665e+75, 7.90241504009859e+74, ...]
        j          = 3251
        params     = [array([-100.        ,  -99.6322642 ,  -99.26588069, ...,   99.26470855,
         99.63167596,  100.        ])]
        self       = <Data for zetac>
        skip_mask  = None
scipy/special/_mptestutils.py:259: in <lambda>
    lambda *a: pytype(self.mpmath_func(*map(mptype, a))),
        a          = (1.0,)
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa6a42210d0>
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa6a42214c0>
        self       = <MpmathData: zetac>
scipy/special/tests/test_mpmath.py:1870: in <lambda>
    lambda x: mpmath.zeta(x) - 1,
        x          = mpf('1.0')
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/functions/zeta.py:533: in zeta
    return ctx._zeta(s, **kwargs)
        a          = 1
        ctx        = <mpmath.ctx_mp.MPContext object at 0x7fa6bb4545e0>
        d          = 0
        derivative = 0
        kwargs     = {}
        method     = None
        s          = mpf('1.0')
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/ctx_mp_python.py:998: in f
    return ctx.make_mpf(mpf_f(x._mpf_, prec, rounding))
        ctx        = <mpmath.ctx_mp.MPContext object at 0x7fa6bb4545e0>
        kwargs     = {}
        mpc_f      = <function mpc_zeta at 0x7fa6bb5ac430>
        mpf_f      = <function mpf_zeta at 0x7fa6bb5ac3a0>
        name       = 'zeta'
        prec       = 153
        rounding   = 'n'
        x          = mpf('1.0')
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/libmp/gammazeta.py:984: in mpf_zeta
    return mpf_zeta_int(to_int(s), prec, rnd)
        alt        = 0
        bc         = 1
        exp        = 0
        man        = mpz(1)
        prec       = 153
        rnd        = 'n'
        s          = (0, mpz(1), 0, 1)
        sign       = 0
        wp         = 173
/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/mpmath/libmp/gammazeta.py:917: in mpf_zeta_int
    raise ValueError("zeta(1) pole")
E   ValueError: zeta(1) pole
        prec       = 153
        rnd        = 'n'
        s          = 1
        wp         = 173
____________________________ TestSystematic.test_spherical_in _____________________________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/special/tests/test_mpmath.py:1977: in test_spherical_in
    assert_mpmath_equal(lambda n, z: sc.spherical_in(int(n), z),
        mp_spherical_in = <function TestSystematic.test_spherical_in.<locals>.mp_spherical_in at 0x7fa6a40f6e50>
        self       = <scipy.special.tests.test_mpmath.TestSystematic object at 0x7fa66ffdc550>
scipy/special/_mptestutils.py:287: in assert_mpmath_equal
    d.check()
        a          = (<function TestSystematic.test_spherical_in.<locals>.<lambda> at 0x7fa66ff9ba60>, <function exception_to_nan.<locals>....cipy.special._mptestutils.IntArg object at 0x7fa66ffdc7f0>, <scipy.special._mptestutils.Arg object at 0x7fa66ffdc2e0>])
        d          = <MpmathData: wrap>
        kw         = {'atol': 1e-278, 'dps': 200}
scipy/special/_mptestutils.py:274: in check
    raise value
        argarr     = array([[ 0.00000000e+000, -8.98846567e+307],
       [ 0.00000000e+000, -3.84127030e+269],
       [ 0.00000000e+000, -1...000e+002,  1.64158802e+231],
       [ 1.65000000e+002,  3.84127030e+269],
       [ 1.65000000e+002,  8.98846567e+307]])
        dps        = 200
        dps_list   = [200]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa66ff9b9d0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa66ff9b820>
        self       = <MpmathData: wrap>
        tb         = <traceback object at 0x7fa6b32be780>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 1.10313e-273\nMax |rdiff|: 1\nBad results (1 out of 837) for the following points (in o...38786e-27 =>                             0. !=          1.10313444760931e-273  (rdiff                             1.)')
scipy/special/_mptestutils.py:258: in check
    assert_func_equal(self.scipy_func,
        argarr     = array([[ 0.00000000e+000, -8.98846567e+307],
       [ 0.00000000e+000, -3.84127030e+269],
       [ 0.00000000e+000, -1...000e+002,  1.64158802e+231],
       [ 1.65000000e+002,  3.84127030e+269],
       [ 1.65000000e+002,  8.98846567e+307]])
        dps        = 200
        dps_list   = [200]
        j          = 0
        mptype     = <function MpmathData.check.<locals>.mptype at 0x7fa66ff9b9d0>
        old_dps    = 15
        old_prec   = 53
        pytype     = <function MpmathData.check.<locals>.pytype at 0x7fa66ff9b820>
        self       = <MpmathData: wrap>
        tb         = <traceback object at 0x7fa6b32be780>
        tp         = <class 'AssertionError'>
        value      = AssertionError('\nMax |adiff|: 1.10313e-273\nMax |rdiff|: 1\nBad results (1 out of 837) for the following points (in o...38786e-27 =>                             0. !=          1.10313444760931e-273  (rdiff                             1.)')
scipy/special/_testutils.py:83: in assert_func_equal
    fdata.check()
E   AssertionError: 
E   Max |adiff|: 1.10313e-273
E   Max |rdiff|: 1
E   Bad results (1 out of 837) for the following points (in output 0):
E                              10.          5.2250558491838786e-27 =>                             0. !=          1.10313444760931e-273  (rdiff                             1.)
        atol       = 1e-278
        data       = array([[ 0.00000000e+000, -8.98846567e+307],
       [ 0.00000000e+000, -3.84127030e+269],
       [ 0.00000000e+000, -1...000e+002,  1.64158802e+231],
       [ 1.65000000e+002,  3.84127030e+269],
       [ 1.65000000e+002,  8.98846567e+307]])
        distinguish_nan_and_inf = True
        dtype      = None
        fdata      = <Data for <lambda>>
        func       = <function TestSystematic.test_spherical_in.<locals>.<lambda> at 0x7fa66ff9ba60>
        ignore_inf_sign = False
        knownfailure = None
        nan_ok     = True
        nparams    = 2
        param_filter = None
        points     = array([[ 0.00000000e+000, -8.98846567e+307],
       [ 0.00000000e+000, -3.84127030e+269],
       [ 0.00000000e+000, -1...000e+002,  1.64158802e+231],
       [ 1.65000000e+002,  3.84127030e+269],
       [ 1.65000000e+002,  8.98846567e+307]])
        result_columns = None
        result_func = <function MpmathData.check.<locals>.<lambda> at 0x7fa66ff9b310>
        results    = <function MpmathData.check.<locals>.<lambda> at 0x7fa66ff9b310>
        rtol       = 1e-07
        vectorized = False
_________ TestNumericalInverseHermite.test_basic_all_scipy_dists[anglit-shapes1] __________
[gw2] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/stats/tests/test_sampling.py:927: in test_basic_all_scipy_dists
    fni = NumericalInverseHermite(dist)
        dist       = <scipy.stats._distn_infrastructure.rv_frozen object at 0x7fa66fd752e0>
        distname   = 'anglit'
        fail_dists = {'beta', 'gausshyper', 'genhyperbolic', 'geninvgauss', 'invgauss', 'kappa4', ...}
        self       = <scipy.stats.tests.test_sampling.TestNumericalInverseHermite object at 0x7fa66fc76b80>
        shapes     = ()
        slow_dists = {'ksone', 'kstwo', 'levy_stable', 'skewnorm'}
        sup        = <numpy.testing._private.utils.suppress_warnings object at 0x7fa66fd6d340>
scipy/stats/unuran_wrapper.pyx:1920: in scipy.stats.unuran_wrapper.NumericalInverseHermite.__cinit__
    ???
E   RuntimeWarning: [objid: HINV.048] 98 : one or more intervals very short; possibly due to numerical problems with a pole or very flat tail => (serious) round-off error
        DiscreteAliasUrn = <class 'scipy.stats.unuran_wrapper.DiscreteAliasUrn'>
        DiscreteGuideTable = <class 'scipy.stats.unuran_wrapper.DiscreteGuideTable'>
        Method     = <class 'scipy.stats.unuran_wrapper.Method'>
        NumericalInverseHermite = <class 'scipy.stats.unuran_wrapper.NumericalInverseHermite'>
        NumericalInversePolynomial = <class 'scipy.stats.unuran_wrapper.NumericalInversePolynomial'>
        SimpleRatioUniforms = <class 'scipy.stats.unuran_wrapper.SimpleRatioUniforms'>
        TransformedDensityRejection = <class 'scipy.stats.unuran_wrapper.TransformedDensityRejection'>
        UError     = <class 'importlib._bootstrap.UError'>
        UNURANError = <class 'scipy.stats.unuran_wrapper.UNURANError'>
        _URNG      = <class 'scipy.stats.unuran_wrapper._URNG'>
        __all__    = ['UNURANError', 'TransformedDensityRejection', 'DiscreteAliasUrn', 'NumericalInversePolynomial']
        __builtins__ = <builtins>
        __doc__    = None
        __file__   = '/home/rgommers/code/scipy/installdir/lib/python3.9/site-packages/scipy/stats/_unuran/unuran_wrapper.cpython-39-x86_64-linux-gnu.so'
        __loader__ = <_frozen_importlib_external.ExtensionFileLoader object at 0x7fa6bea7c520>
        __name__   = 'scipy.stats._unuran.unuran_wrapper'
        __package__ = 'scipy.stats._unuran'
        __pyx_unpickle_Enum = <built-in function __pyx_unpickle_Enum>
        __pyx_unpickle__URNG = <built-in function __pyx_unpickle__URNG>
        __spec__   = ModuleSpec(name='scipy.stats._unuran.unuran_wrapper', loader=<_frozen_importlib_external.ExtensionFileLoader object at...s/code/scipy/installdir/lib/python3.9/site-packages/scipy/stats/_unuran/unuran_wrapper.cpython-39-x86_64-linux-gnu.so')
        __test__   = {'TransformedDensityRejection.ppf_hat (line 940)': '\n        ppf_hat(u)\n\n        Evaluate the inverse of the CDF of...rm.ppf(0.5)\n        0.0\n        >>> u = np.linspace(0, 1, num=1000)\n        >>> ppf_hat = rng.ppf_hat(u)\n        '}
        _validate_domain = <built-in function _validate_domain>
        argsreduce = <function argsreduce at 0x7fa6bf1f8a60>
        check_random_state = <function check_random_state at 0x7fa6c4ad1b80>
        functools  = <module 'functools' from '/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/functools.py'>
        namedtuple = <function namedtuple at 0x7fa6d62d2670>
        np         = <module 'numpy' from '/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/site-packages/numpy/__init__.py'>
        rv_frozen  = <class 'scipy.stats._distn_infrastructure.rv_frozen'>
        stats      = <module 'scipy.stats' from '/home/rgommers/code/scipy/installdir/lib/python3.9/site-packages/scipy/stats/__init__.py'>
        threading  = <module 'threading' from '/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/threading.py'>
        warnings   = <module 'warnings' from '/home/rgommers/anaconda3/envs/scipy-dev/lib/python3.9/warnings.py'>
______________________ test_moments[vonmises-arg102-True-True-True] _______________________
[gw1] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
scipy/stats/tests/test_continuous_basic.py:292: in test_moments
    check_normalization(distfn, arg, distname)
        arg        = (3.99390425810714,)
        distfn     = <scipy.stats._continuous_distns.vonmises_gen object at 0x7f311291b580>
        distname   = 'vonmises'
        higher_ok  = True
        is_xfailing = True
        k          = array(0.69057153)
        m          = array(-1.23486137e-16)
        normalization_ok = True
        s          = array(-6.31595765e-15)
        sup        = <numpy.testing._private.utils.suppress_warnings object at 0x7f306bd17700>
        v          = array(0.29880424)
scipy/stats/tests/common_tests.py:34: in check_normalization
    npt.assert_allclose(normalization_expect, 1.0, atol=atol, rtol=rtol,
E   AssertionError: 
E   Not equal to tolerance rtol=1e-07, atol=1e-07
E   vonmises
E   Mismatched elements: 1 / 1 (100%)
E   Max absolute difference: 153415.76801168
E   Max relative difference: 153415.76801168
E    x: array(153416.768012)
E    y: array(1.)
        args       = (3.99390425810714,)
        atol       = 1e-07
        distfn     = <scipy.stats._continuous_distns.vonmises_gen object at 0x7f311291b580>
        distname   = 'vonmises'
        norm_moment = 1.0
        normalization_expect = 153416.76801168497
        rtol       = 1e-07
______ installdir/lib/python3.9/site-packages/scipy/spatial/tests/test_hausdorff.py _______
[gw3] linux -- Python 3.9.7 /home/rgommers/anaconda3/envs/scipy-dev/bin/python
worker 'gw3' crashed while running 'installdir/lib/python3.9/site-packages/scipy/spatial/tests/test_hausdorff.py::test_massive_arr_overflow'
================================= short test summary info =================================
FAILED scipy/linalg/tests/test_lapack.py::test_sgesdd_lwork_bug_workaround - pytest.Pyte...
FAILED scipy/sparse/linalg/tests/test_onenormest.py::TestOnenormest::test_onenormest_table_6_t_1
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_airyai_prime - Assertion...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_airybi - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_ci_complex - AssertionEr...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_igam_fac - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_legenp - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_sinpi - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_struvel - AssertionError: 
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_riemann_zeta - ValueErro...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_zetac - ValueError: zeta...
FAILED scipy/special/tests/test_mpmath.py::TestSystematic::test_spherical_in - Assertion...
FAILED scipy/stats/tests/test_sampling.py::TestNumericalInverseHermite::test_basic_all_scipy_dists[anglit-shapes1]
FAILED scipy/stats/tests/test_continuous_basic.py::test_moments[vonmises-arg102-True-True-True]
FAILED scipy/spatial/tests/test_hausdorff.py::test_massive_arr_overflow
=== 15 failed, 47544 passed, 2382 skipped, 237 xfailed, 11 xpassed in 633.30s (0:10:33) ===

Some of these are small test issues; tolerance being slightly exceeded for a couple of special functions. Others indicate problems that need investigating and fixing. A good way of doing this is probably to look at the full set of flags with which each extension is compiled. These tests are not failing in numpy.distutils based CI, so the difference in build flags (e.g. -std=c99 enabled, or optimization level) is likely the explanation.

Issue Analytics

  • State:open
  • Created 2 years ago
  • Comments:14 (14 by maintainers)

github_iconTop GitHub Comments

1reaction
WarrenWeckessercommented, Oct 10, 2022

https://github.com/scipy/scipy/pull/17190 should fix two of the xslow failures in scipy.special.

1reaction
WarrenWeckessercommented, Jul 20, 2022

I created https://github.com/scipy/scipy/issues/16655 to explain the failure indicated by

FAILED scipy/stats/tests/test_continuous_basic.py::test_moments[vonmises-arg103-True-True-True]

The test test_moments calls check_normalization in common_tests.py, and check_normalization computes

normalization_expect = distfn.expect(lambda x: 1, args=args)

which should be 1. As explained in https://github.com/scipy/scipy/issues/16655, vonmises.expect is broken, so this calculation will return nonsense and the test will fail.

Read more comments on GitHub >

github_iconTop Results From Across the Web

Why test suites fail—and how to avoid breakdowns
1. Your test suite has become so large, unwieldy, and unstable that nobody takes the feedback it produces into consideration anymore. It fails...
Read more >
Test Automation - IV - Excelon Development
It would be awful nice to come up with a consistent model for different types of testing tools. Failing that, I'd settle for...
Read more >
How to Fix Flaky Tests - Semaphore CI
Randomly failing tests are the hardest to debug. Here's a framework you can use to fix them and keep your test suite healthy....
Read more >
How to run failed test cases in TestNG & Selenium
Failed tests will inevitably show up in any software testing lifecycle ... How to Run a Failed Test Case using TestNG in Selenium...
Read more >
crackback meth - Il gusto in piazza
A super oily or thick puddle that runs away from heat is not good. ... is used for strengthening cartilage Pipe Test: If...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found