
    _Mh                     @    d dl Zd dlmZmZ d dlmZmZ d Zd Z	d Z
dS )    N)assert_allcloseassert_equal)tukeylambda_variancetukeylambda_kurtosisc                     t          d          } t          | t          j        dz  dz  d           t	          d          }t          |dd           t          d          } t          | d	t          j        z
  d           t	          d          }d
t          j        dz  z
  t          j        d	z  dz
  dz  z  dz
  }t          ||d           t          d          } t          | dd           t	          d          }t          |dd           t          d          } t          | dd           t	          d          }t          |dd           dS )z/Compare results with some known exact formulas.r         -q=)atol333333?绽|=g      ?   g?   gUUUUUU?g333333gUUUUUU?N)r   r   nppir   )varkurtdesireds      h/var/www/html/test/jupyter/venv/lib/python3.11/site-packages/scipy/stats/tests/test_tukeylambda_stats.py"test_tukeylambda_stats_known_exactr      sX    q
!
!CCAE2222""DD#E**** s
#
#CCRU////$$DbeAg~"%'A+!11A5GD'.... q
!
!CCu----""DD$U++++ q
!
!CC....""DD$U++++++    c                     t          dd          } t          dd          }g dg dg dg dg d	g d
g dg dg dg dg dg}|D ]@\  }}}t          |          }t          ||fi | t          |          }t          ||fi |  At	          | \  }}}t          |          }t          ||fi | t          |          }t          ||fi |  dS )zACompare results with some values that were computed using mpmath.r   r   )r   rtolr
   )gg;@g2V#H@)gQIg9@gk+\)@)ggxuh~@g@)gMbPg٤{i
@ghYn?)gMbP?gӦNFb:
@g4\>?)g      ?gt@g<n\?)g?gU+2C@gӤB?)gQI?gGI8@gx@As?)r   gڽP?g@᜿)g      $@gds5O?gX\@)g      4@g ?gc@N)dictr   r   r   zip)a10a12datalamvar_expectedkurt_expectedr   r   s           r   test_tukeylambda_stats_mpmathr"   +   sY   
E
"
"
"C
E
"
"
"C 	988;;;999:::999<<<999;;;999;;;;;;D -1 4 4(\="3''\11S111#C((m33s3333 (+Dz$C}
s
#
#CC-----$$DD-//3/////r   c                  2   ddg} t          |           }t          |t          j        t          j        t          j        g                     ddg} t          |           }t          |t          j        t          j        t          j        g                     dS )z;Test values of lambda outside the domains of the functions.g      g      g      пN)r   r   r   arraynaninfr   )r   r   r   s      r   test_tukeylambda_stats_invalidr'   M   sz    ,C
s
#
#Cbh/00111-C$$Drx 01122222r   )numpyr   numpy.testingr   r   scipy.stats._tukeylambda_statsr   r   r   r"   r'    r   r   <module>r,      s        7 7 7 7 7 7 7 7B B B B B B B B ,  ,  ,F0 0 0D3 3 3 3 3r   