pnumpy breaks scipy's skewtest

[saimn]

In [1]: import scipy.stats as st
   ...: arr = st.norm.rvs(loc=6000, scale=20, size=(300, 10))
   ...: st.skewtest(arr)
Out[1]: 
SkewtestResult(statistic=array([ 1.21508957, -0.60569814,  0.07676428,  0.10929194, -1.78128941,
        1.3046305 ,  1.03649957, -0.90213437,  1.39273949,  0.87478829]), pvalue=array([0.22433193, 0.54471522, 0.93881107, 0.91297094, 0.07486518,
       0.1920187 , 0.29996914, 0.36698549, 0.16369858, 0.38168911]))

In [2]: import pnumpy; pnumpy.enable()

In [3]: st.skewtest(arr)
Out[3]: 
SkewtestResult(statistic=array([-4.34541497e+03, -4.34723576e+03, -4.34655330e+03, -4.34652077e+03,
       -4.34841135e+03, -4.34532543e+03, -4.34559356e+03, -4.34753220e+03,
        1.39273949e+00,  8.74788294e-01]), pvalue=array([0.        , 0.        , 0.        , 0.        , 0.        ,
       0.        , 0.        , 0.        , 0.16369858, 0.38168911]))

In [4]: pnumpy.__version__
Out[4]: '2.0.23'