§
    ¢éUgÂZ  ã                   ó  — d dl Z d dlZd dlmZmZmZmZmZmZm	Z	m
Z
mZmZmZmZmZmZ d dlZd dlmZmZmZmZ d dlmZ d dlmZ d dlmZ e j                             dg d¢¦  «        d	„ ¦   «         Z e j                             dd
dg¦  «        d„ ¦   «         Z!d„ Z"d„ Z#d„ Z$d„ Z%d„ Z&d„ Z'd„ Z(d„ Z)e j                             d ej*        e+e,fgdz  Ž ¦  «        d„ ¦   «         Z-d„ Z.d„ Z/d„ Z0d„ Z1d„ Z2d„ Z3d„ Z4e j                             dg d ¢¦  «        d!„ ¦   «         Z5d"„ Z6d#„ Z7d$„ Z8d%„ Z9d&„ Z: G d'„ d(¦  «        Z; G d)„ d*¦  «        Z<e j                             d+g d,¢g d-¢g¦  «        d.„ ¦   «         Z=d/„ Z> G d0„ d1¦  «        Z? G d2„ d3¦  «        Z@dS )4é    N)Ú	betabinomÚ
betanbinomÚ	hypergeomÚ
nhypergeomÚ	bernoulliÚ	boltzmannÚskellamÚzipfÚzipfianÚbinomÚnbinomÚnchypergeom_fisherÚnchypergeom_walleniusÚrandint)Úassert_almost_equalÚassert_equalÚassert_allcloseÚsuppress_warnings)r   )Úroot_scalar)Úquadzk, M, n, N, expected, rtol))é   é
   é   é   gÏó<Ïó<ï?çVçž¯Ò<)ék   é'  é¸  é×   gŒâÙÿÿÿï?r   )r   r   r   r   gŸåèS©;ç»½×Ùß|Ë=c                 óX   — t          j        | |||¦  «        }t          |||¬¦  «         d S ©N©Úrtol)r   Úcdfr   ©ÚkÚMÚnÚNÚexpectedr$   Úps          úe/var/www/surfInsights/venv3-11/lib/python3.11/site-packages/scipy/stats/tests/test_discrete_distns.pyÚtest_hypergeom_cdfr.      s4   € õ 	Œ�a˜˜A˜qÑ!Ô!€AÝ�A�x dÐ+Ñ+Ô+Ð+Ð+Ð+ó    )é   r   r   r   gäûòÿÿÿï?r   )é}   r   r   r   g«™Çö!˜:<r    c                 óX   — t          j        | |||¦  «        }t          |||¬¦  «         d S r"   )r   Úsfr   r&   s          r-   Útest_hypergeom_sfr4   "   s4   € õ 	Œ�Q˜˜1˜aÑ Ô €AÝ�A�x dÐ+Ñ+Ô+Ð+Ð+Ð+r/   c                  óä  — d} d}d}d}t          j        | |||¦  «        }t          j        || z
  |||z
  |¦  «        }t          j        || z
  ||||z
  ¦  «        }t          j        | |||¦  «        }t          ||d¬¦  «         t          ||d¬¦  «         t          ||d¬¦  «         d} d}d}d}t          j        | |||¦  «        }t          j        | ||z  ¦  «        }	t          ||	d¬¦  «         d S )Nr   é2   r   é   )Údecimalé   é   )r   Úlogpmfr   r   )
r'   r*   ÚKr)   Úlogpmf1Úlogpmf2Úlogpmf3Úlogpmf4Úhypergeom_logpmfÚbernoulli_logpmfs
             r-   Útest_hypergeom_logpmfrC   ,   s   € ð 	
€AØ
€AØ
€AØ	€AÝÔ˜q ! Q¨Ñ*Ô*€GÝÔ˜q 1™u a¨¨Q©°Ñ2Ô2€GÝÔ˜q 1™u a¨¨A°©EÑ2Ô2€GÝÔ˜q ! Q¨Ñ*Ô*€GÝ˜ °"Ð5Ñ5Ô5Ð5Ý˜ °"Ð5Ñ5Ô5Ð5Ý˜ °"Ð5Ñ5Ô5Ð5ð 	
€AØ
€AØ	€AØ	€AÝ Ô'¨¨1¨a°Ñ3Ô3ÐÝ Ô'¨¨1¨Q©3Ñ/Ô/ÐÝÐ(Ð*:ÀBÐGÑGÔGÐGÐGÐGr/   c                  óÒ   — d\  } }}d}t          j        || ||¦  «        }t          j        || |||z   dz
  ¦  «        | |z
  |dz
  z
  z  | ||z   dz
  z
  z  }t          ||d¬¦  «         d S )N)é-   é   é   é   r9   ç»½×Ùß|Û=r#   )r   Úpmfr   r   )r(   r)   Úrr'   ÚNHGÚHGs         r-   Útest_nhypergeom_pmfrN   F   s�   € à�G€A€qˆ!Ø	€AÝ
Œ.˜˜A˜q !Ñ
$Ô
$€CÝ	Œ�q˜!˜Q  !¡ A¡Ñ	&Ô	&¨!¨a©%°1°Q±3©-Ñ	8¸AÀÀ1ÁÀQÁ¹KÑ	H€BÝ�B˜ %Ð(Ñ(Ô(Ð(Ð(Ð(r/   c                  óì   — d} d}d}t          j        |dz   ¦  «        }t          j        || ||¦  «        }t          j        || ||¦  «        }t          |g d¢d¬¦  «         t          |g d¢d¬¦  «         d S )	NrG   r   r   r9   )ç’$I’$I²?gÛ¶mÛ¶mË?ç·mÛ¶mÛÖ?rQ   ç‚vIhÂ%<=r#   )rP   ç’$I’$IÒ?g%I’$I’ä?ç      ð?)ÚnpÚaranger   rJ   r%   r   )r(   r)   rK   ÚsupportrJ   r%   s         r-   Útest_nhypergeom_pmfcdfrX   O   sŠ   € à	€AØ	€AØ	€AÝŒi˜˜!™‰nŒn€GÝ
Œ.˜ ! Q¨Ñ
*Ô
*€CÝ
Œ.˜ ! Q¨Ñ
*Ô
*€CÝ�CÐ1Ð1Ð1¸Ð>Ñ>Ô>Ð>Ý�CÐ0Ð0Ð0°uÐ=Ñ=Ô=Ð=Ð=Ð=r/   c                  ó|   — d} d}d}t          j        g d¢g d¢g| ||¦  «        }t          |g d¢g d¢gd¬	¦  «         d S )
Nr   r   r   )r   r9   é   r   )r9   rZ   r   r   )r9   r   r   r9   )r   r   r9   r   rR   r#   )r   rJ   r   )r(   r)   rK   rJ   s       r-   Útest_nhypergeom_r0r[   [   s[   € à
€AØ	€AØ	€AÝ
Œ.˜,˜,˜,¨¨¨Ð5°q¸!¸QÑ
?Ô
?€CÝ�C˜,˜,˜,¨¨¨Ð5¸EÐBÑBÔBÐBÐBÐBr/   c                  ó\   — t          j        dg d¢dgdggd¬¦  «        } | j        dk    sJ ‚d S )Né   )r:   rG   é	   r7   rF   ©r   r9   rZ   r   ©Úsize)r   ÚrvsÚshape)Úxs    r-   Útest_nhypergeom_rvs_shapere   d   sA   € õ 	Œ�r˜9˜9˜9¨ t¨b¨T l¸ÐFÑFÔF€AØŒ7�lÒ"Ð"Ð"Ð"Ð"Ð"r/   c                  ó@  — t           j                             d¦  «         t          j        dddd¬¦  «        } t           j                             d¦  «         t           j                             d¬¦  «        }t          j        |ddd¦  «        }t          | |¦  «         d S )Nr   r]   r:   é   éd   r`   )rU   ÚrandomÚseedr   rb   ÚuniformÚppfr   )rd   r,   Úys      r-   Útest_nhypergeom_accuracyrn   l   s‡   € õ „I‡N‚N�1ÑÔÐÝŒ�r˜1˜b sÐ+Ñ+Ô+€AÝ„I‡N‚N�1ÑÔÐÝ
Œ	×Ò˜sÐÑ#Ô#€AÝŒ�q˜"˜a Ñ$Ô$€AÝ��AÑÔÐÐÐr/   c                  óh  — t          j        dd¦  «        } d}t          j        | d|¦  «        }| dk    }t	          ||¦  «         t          j        d¦  «        }d}t          j        | ||¦  «        }g d¢}t          ||d	¬
¦  «         t          j        | ||¦  «        }g d¢}t          ||d	¬
¦  «         d S )Néýÿÿÿr   r9   g°rh‘í|¿?r   rZ   r   )r   r   r   ç’$I’$Iâ?rS   g’$I’$IÂ?r   r   rR   r#   )r   r   r   rq   gÛ¶mÛ¶më?r9   r9   r9   )rU   rV   r   rJ   r   Úlogr   r%   )r'   r*   r,   r+   ÚlamÚcs         r-   Útest_boltzmann_upper_boundru   w   sÂ   € Ý
Œ	�"�aÑÔ€Aà	€AÝŒ�a˜ Ñ"Ô"€AØ�AŠv€HÝ��HÑÔÐå
Œ&�‰)Œ)€CØ	€AÝŒ�a˜˜aÑ Ô €AØ-Ð-Ð-€HÝ�A�x eÐ,Ñ,Ô,Ð,åŒ�a˜˜aÑ Ô €AØ+Ð+Ð+€HÝ�A�x eÐ,Ñ,Ô,Ð,Ð,Ð,r/   c                  óÜ   — d} t          j        | dz   ¦  «        }t          | dd¦  «                             |¦  «        }t          j        d| dz   z  | dz   ¦  «        }t          ||¦  «         d S )Né   r9   )rU   rV   r   rJ   Úrepeatr   )r)   r'   r,   r+   s       r-   Útest_betabinom_a_and_b_unityry   Š   sn   € ð 	€AÝ
Œ	�!�a‘%ÑÔ€AÝ�!�Q˜ÑÔ×Ò˜qÑ!Ô!€AÝŒy˜˜a !™e™ a¨!¡eÑ,Ô,€HÝ˜˜8Ñ$Ô$Ð$Ð$Ð$r/   Údtypesr   c                 ó    — | \  }}} |d¦  «         |d¦  «         |d¦  «        }}}t          t          j        |||d¬¦  «        d¦  «         d S )Nr   rZ   r   r'   ©Úmomentsgˆa†aæ¿)r   r   Ústats)rz   Ún_typeÚa_typeÚb_typer)   ÚaÚbs          r-   Ú-test_betabinom_stats_a_and_b_integers_gh18026r„   ”   s_   € ð $Ñ€FˆF�FØˆf�R‰jŒj˜&˜& ™)œ) V V¨A¡Y¤Yˆ!€q€AÝ•I”O A q¨!°SÐ9Ñ9Ô9Ð;NÑOÔOÐOÐOÐOr/   c                  óî   — d} d}t          j        d¦  «        }t          d| |¦  «                             |¦  «        }t	          | | |z   z  ¦  «                             |¦  «        }t          ||¦  «         d S )Ngffffff@g)\�Âõ(ä?rZ   r9   )rU   rV   r   rJ   r   r   )r‚   rƒ   r'   r,   r+   s        r-   Útest_betabinom_bernoullir†   �   sq   € à€AØ€AÝ
Œ	�!‰Œ€AÝ�!�Q˜ÑÔ×Ò˜qÑ!Ô!€AÝ˜˜a !™e™Ñ%Ô%×)Ò)¨!Ñ,Ô,€HÝ˜˜8Ñ$Ô$Ð$Ð$Ð$r/   c                  ó\   — d\  } }}t          t          j        | ||¬¦  «        d¦  «         d S )N)gÍÌÌÌÌÌì?r   r9   ©Ú
confidencer)   r,   ©r   r   )r   r   Úinterval©Úalphar)   r,   s      r-   Útest_issue_10317rŽ   §   s4   € Ø�K€Eˆ1ˆaÝ•”¨E°Q¸!Ð<Ñ<Ô<¸fÑEÔEÐEÐEÐEr/   c                  ó\   — d\  } }}t          t          j        | ||¬¦  «        d¦  «         d S )N)gffffffî?r   r   rˆ   rŠ   )r   r   r‹   rŒ   s      r-   Útest_issue_11134r�   ¬   s4   € Ø�K€Eˆ1ˆaÝ•”¨5°A¸Ð;Ñ;Ô;¸VÑDÔDÐDÐDÐDr/   c                  óV  — t           j                             d¦  «         t          t	          j        t           j                             d¦  «        dd¦  «        d¦  «         t          t	          j        ddd¦  «        d¦  «         t          t	          j        ddd¦  «        d¦  «         d S )Nr   r   ç      à?éÿÿÿÿr9   )rU   ri   rj   r   r   rl   Úrand© r/   r-   Útest_issue_7406r–   ±   s†   € Ý„I‡N‚N�1ÑÔÐÝ•”�2œ9Ÿ>š>¨"Ñ-Ô-¨q°#Ñ6Ô6¸Ñ:Ô:Ð:õ •”˜1˜a Ñ%Ô% rÑ*Ô*Ð*Ý•”˜1˜a Ñ%Ô% qÑ)Ô)Ð)Ð)Ð)r/   c                  ód  — d} t           j                             dd¬¦  «        }d}t          j        ||| ¦  «        }t          |d¦  «         t          j        ddd¦  «        }t          j        ||| ¦  «        }t          |d¦  «         d}t          j        ||| ¦  «        }t          ||¦  «         d S )	Nr   rh   r   r`   r“   ç{®Gáz„?g®Gáz®ï?r9   )rU   ri   r   r   rl   r   Úlinspace)r,   r)   rd   rl   s       r-   Útest_issue_5122rš   º   s­   € Ø	€AÝ
Œ	×Ò˜# BÐÑ'Ô'€Aà	€AÝ
Œ)�A�q˜!Ñ
Ô
€CÝ��bÑÔÐå
Œ�D˜$ Ñ#Ô#€AÝ
Œ)�A�q˜!Ñ
Ô
€CÝ��aÑÔÐà	€AÝ
Œ)�A�q˜!Ñ
Ô
€CÝ��aÑÔÐÐÐr/   c            	      óŽ   — t          t          dt          j        dd¦  «        ¦  «                             d¦  «        d¦  «         d S )Néè  rp   iœÿÿÿr˜   r   )r   r   rU   Úlogspacerl   r•   r/   r-   Útest_issue_1603rž   Ë   s<   € Ý•�t�Rœ[¨¨TÑ2Ô2Ñ3Ô3×7Ò7¸Ñ=Ô=¸qÑAÔAÐAÐAÐAr/   c                  óˆ   — d} t          j        ddd¦  «        }t          t          j        |d|z  | ¦  «        dd¬¦  «         d S )Nr’   r   é   r7   rZ   r˜   )Úatol)rU   r�   r   r   r%   )r,   rd   s     r-   Útest_issue_5503r¢   Ï   sG   € Ø€AÝ
Œ�A�r˜2ÑÔ€AÝ•E”I˜a  1¡ aÑ(Ô(¨#°DÐ9Ñ9Ô9Ð9Ð9Ð9r/   zx, n, p, cdf_desired)	)i,  rœ   ç333333Ó?g’24û½à?)r   r   r£   gUž…7i(à?)i0u  i † r£   gÎ­Ì‘Çà?)ià“ i@B r£   g8@¨�
à?)iÀÆ- é€–˜ r£   gfø‡(Gà?)i€ÃÉé áõr£   gÓwàtg à?)i�êÉr¥   r£   g¥Vç}Öˆï?)ipœÉr¥   r£   g¢°�ÿzÌ�?)i0 Ér¥   r£   g¡éÉnåC:c                 óN   — t          t          j        | ||¦  «        |¦  «         d S ©N©r   r   r%   )rd   r)   r,   Úcdf_desireds       r-   Útest_issue_5503pt2rª   Õ   s(   € õ •E”I˜a  AÑ&Ô&¨Ñ4Ô4Ð4Ð4Ð4r/   c                  óN   — t          t          j        ddd¦  «        d¦  «         d S )NrZ   l    J)£gê-�™—q=gB.¸Ã+ní?r¨   r•   r/   r-   Útest_issue_5503pt3r¬   ä   s'   € å•E”I˜a ¨Ñ1Ô1Ð3IÑJÔJÐJÐJÐJr/   c                  óN   — t          t          j        ddd¦  «        d¦  «         d S )Néú   r6   gAAà?g_[CË®i³8)r   r   r3   r•   r/   r-   Útest_issue_6682r¯   é   s)   € õ •F”I˜c 2 wÑ/Ô/Ð1FÑGÔGÐGÐGÐGr/   c                  óv   — t          j        g d¢dd¦  «        } dt          j         dg}t	          | |¦  «         d S )N)r   r“   r9   r   r’   g¬¥ÂÐÑGÞ¿g1DÊ&+´À)r   ÚlogcdfrU   Úinfr   )ÚresultÚ	references     r-   Útest_issue_19747rµ   ð   s?   € åŒ]˜:˜:˜: q¨#Ñ.Ô.€FØ�rœv˜g {Ð3€IÝ�F˜IÑ&Ô&Ð&Ð&Ð&r/   c                  óZ   — d} d}d}t          t          j        || |¦  «        d¦  «         d S )Nrœ   r˜   iä  ç        )r   r   rJ   )r)   r,   r'   s      r-   Ú%test_boost_divide_by_zero_issue_15101r¸   ÷   s5   € Ø€AØ€AØ€AÝ•E”I˜a  AÑ&Ô&¨Ñ,Ô,Ð,Ð,Ð,r/   c                  ób   — g d¢} t          j        d| | ¦  «        }g d¢}t          ||¦  «         d S )N)	r9   r   rh   rœ   éˆ  iº  iì  i‚  ip  r   )	gâ–>l¦ïä?g‹ëõz½oá?gVr�Ržsà?g°�ˆŠ$à?gVWêCWà?g4y�Bà?g4ö€.à?g.ÚLròà?gDáøÅêà?)r	   r%   r   )Úmur%   Úcdf_expecteds      r-   Útest_skellam_gh11474r½   þ   sL   € à	9Ð	9Ð	9€BÝ
Œ+�a˜˜RÑ
 Ô
 €CðPð Pð P€Lõ �C˜Ñ&Ô&Ð&Ð&Ð&r/   c                   ó$  — e Zd Zd„ Zd„ Zd„ Zej                             d¦  «          ej	         ej
        ddd¦  «        ej                             dd	d¦  «        f¦  «        j        Zej                             d
e¦  «        d„ ¦   «         Zd„ ZdS )ÚTestZipfianc                 óô  — d}d}t          j        dd¦  «        }t          t          j        |||¦  «        t          j        ||¦  «        ¦  «         t          t          j        |||¦  «        t          j        ||¦  «        ¦  «         t          t          j        |||¦  «        t          j        ||¦  «        ¦  «         t          t          j        ||d¬¦  «        t          j        |d¬¦  «        ¦  «         d S )Ng      @r¤   r9   é   Úmsvkr|   )	rU   rV   r   r   rJ   r
   r%   r3   r~   )Úselfr‚   r*   r'   s       r-   Útest_zipfian_asymptoticz#TestZipfian.test_zipfian_asymptotic  sÖ   € àˆØˆÝŒI�a˜ÑÔˆÝ�œ A q¨!Ñ,Ô,­d¬h°q¸!©n¬nÑ=Ô=Ð=Ý�œ A q¨!Ñ,Ô,­d¬h°q¸!©n¬nÑ=Ô=Ð=Ý�œ
 1 a¨Ñ+Ô+­T¬W°Q¸©]¬]Ñ;Ô;Ð;Ý�œ a¨°FÐ;Ñ;Ô;Ýœ
 1¨fÐ5Ñ5Ô5ñ	7ô 	7ð 	7ð 	7ð 	7r/   c                 ó  — d\  }}d}t          j        d|dz   ¦  «        }t          t          j        |||¦  «        t          j        |||¦  «        d¬¦  «         t          t          j        |||¦  «        t          j        |||¦  «        d¬¦  «         t          t          j        |||¦  «        t          j        |||¦  «        d¬¦  «         t          t          j        ||d¬¦  «        t          j        ||d¬¦  «        d¬¦  «         d S )N)gGœ¡úÿÿï?gÜ1¯  ð?é   r9   g�íµ ÷Æ >r#   rÂ   r|   )rU   rV   r   r   rJ   r%   r3   r~   )rÃ   Úalt1Úagt1r*   r'   s        r-   Útest_zipfian_continuityz#TestZipfian.test_zipfian_continuity  s%  € ð ,‰
ˆˆdØˆÝŒI�a˜˜Q™ÑÔˆÝ�œ A t¨QÑ/Ô/µ´¸QÀÀaÑ1HÔ1HØ!ð	#ñ 	#ô 	#ð 	#å�œ A t¨QÑ/Ô/µ´¸QÀÀaÑ1HÔ1HØ!ð	#ñ 	#ô 	#ð 	#å�œ
 1 d¨AÑ.Ô.µ´
¸1¸dÀAÑ0FÔ0FØ!ð	#ñ 	#ô 	#ð 	#å�œ d¨A°vÐ>Ñ>Ô>Ýœ d¨A°vÐ>Ñ>Ô>ÀTð	Kñ 	Kô 	Kð 	Kð 	Kð 	Kr/   c                 óþ  — t           j                             d¦  «         t           j                             ddd¬¦  «        }t           j                             d¦  «        dz  dz   }t           j                             ddd¬¦  «        }g d¢}g d¢}t          t          j        |||¦  «        dd …         |dd …         d	¬
¦  «         t          t          j        |||¦  «        dd …         |dd …         d¬
¦  «         d S )Nr   r9   rw   r   r`   rh   )
g÷Ä¹Ô¨Š€?gØ¡–Ÿkïþ>ç	]x�[ï?gÿÇ™ûÊ>g	¤W¿³4?gÂ€
-ä?g˜$fhnåý=g9´’±[.Ä>g±`ý€Ö”Ø>gn²màì ?)
gÉaæñj¯ì?gZ›k=èþï?rË   gõ…�óþÿï?g90×þï?gÉ¸˜UŽÿï?rT   gJ_9ïÿï?g2h2-ûÿï?gB0ýÅþï?ç�íµ ÷Æ°>r#   g-Cëâ6
?)	rU   ri   rj   r   r”   r   r   rJ   r%   )rÃ   r'   r‚   r)   rJ   r%   s         r-   Útest_zipfian_RzTestZipfian.test_zipfian_R(  s  € õ 	Œ	�Š�qÑÔÐÝŒI×Ò˜a ¨"ÐÑ-Ô-ˆÝŒI�NŠN˜2ÑÔ˜rÑ! AÑ%ˆÝŒI×Ò˜a ¨2ÐÑ.Ô.ˆð+ð +ð +ˆðFð Fð Fˆõ 	�œ A q¨!Ñ,Ô,¨Q¨R¨RÔ0°#°a°b°b´'ÀÐEÑEÔEÐEÝ�œ A q¨!Ñ,Ô,¨Q¨R¨RÔ0°#°a°b°b´'ÀÐEÑEÔEÐEÐEÐEr/   r   éþÿÿÿr9   r   rZ   é(   za, nc                 ó¢  ‡— t           j        d„ ¦   «         Št           j        ˆfd„¦   «         }t          j        |dz   ¦  «        } ||||¦  «        }t          j        |¦  «        }t          j        ||¬¦  «        }t          j        ||z
  dz  |¬¦  «        }|dz  }	t          j        ||z
  |	z  dz  |¬¦  «        }
t          j        ||z
  |	z  dz  |¬¦  «        dz
  }t          t          j        |||¦  «        |¦  «         t          t          j        |||¦  «        |¦  «         t          t          j	        ||d	¬
¦  «        |||
|g¦  «         d S )Nc                 ób   — dt          j        d| dz   ¦  «        |z  z                       ¦   «         S )z$Naive implementation of harmonic sumr9   )rU   rV   Úsum)r)   Úss     r-   ÚHnsz+TestZipfian.test_zipfian_naive.<locals>.HnsH  s0   € ð •b”i  1 Q¡3Ñ'Ô'¨Ñ*Ñ*×/Ò/Ñ1Ô1Ð1r/   c                 óJ   •— | dk     s| |k    rdS d| |z  z   ‰||¦  «        z  S )z#Naive implementation of zipfian pmfr9   r·   r•   )r'   r‚   r)   rÔ   s      €r-   Úpzipz,TestZipfian.test_zipfian_naive.<locals>.pzipM  s8   ø€ ð �1Šuˆu˜˜Aš˜Ø�rà˜1˜a™4‘x # # a¨¡)¤)Ñ+Ð+r/   r9   )ÚweightsrZ   r’   r   r   Úmvskr|   )
rU   Ú	vectorizerV   ÚcumsumÚaverager   r   rJ   r%   r~   )rÃ   r‚   r)   rÖ   r'   rJ   r%   ÚmeanÚvarÚstdÚskewÚkurtosisrÔ   s               @r-   Útest_zipfian_naivezTestZipfian.test_zipfian_naiveD  so  ø€ õ 
Œð	2ð 	2ñ 
Œð	2õ 
Œð	,ð 	,ð 	,ð 	,ñ 
Œð	,õ ŒI�a˜‘c‰NŒNˆØˆd�1�a˜‰mŒmˆÝŒi˜‰nŒnˆÝŒz˜! SÐ)Ñ)Ô)ˆÝŒj˜!˜d™( Q™°Ð4Ñ4Ô4ˆØ�3‰hˆÝŒz˜A˜d™F C™<¨!Ñ+°SÐ9Ñ9Ô9ˆÝ”:  $¡¨™|¨aÑ/¸Ð=Ñ=Ô=ÀÑAˆÝ�œ A q¨!Ñ,Ô,¨cÑ2Ô2Ð2Ý�œ A q¨!Ñ,Ô,¨cÑ2Ô2Ð2Ý�œ a¨°FÐ;Ñ;Ô;Ø˜s D¨(Ð3ñ	5ô 	5ð 	5ð 	5ð 	5r/   c                 ó  — t          j        dd¦  «        }|                     t           j        ¦  «        }t	          dd¦  «        }|                     |¦  «        }|                     |¦  «        }t          ||¦  «         d S )Nr   rœ   éo   r]   )rU   rV   ÚastypeÚint32r   rJ   r   ©rÃ   r'   Úk_int32ÚdistrJ   Úpmf_k_int32s         r-   Útest_pmf_integer_kzTestZipfian.test_pmf_integer_kb  sk   € ÝŒI�a˜ÑÔˆØ—(’(�2œ8Ñ$Ô$ˆÝ�s˜BÑÔˆØ�hŠh�q‰kŒkˆØ—h’h˜wÑ'Ô'ˆÝ�S˜+Ñ&Ô&Ð&Ð&Ð&r/   N)Ú__name__Ú
__module__Ú__qualname__rÄ   rÉ   rÍ   rU   ri   rj   Úvstackr�   r   ÚTÚnaive_testsÚpytestÚmarkÚparametrizerá   rê   r•   r/   r-   r¿   r¿     sá   € € € € € ð	7ð 	7ð 	7ðKð Kð KðFð Fð Fð0 „I‡N‚N�1ÑÔÐØ�"”)˜[˜Rœ[¨¨Q°Ñ3Ô3ØœY×.Ò.¨q°"°bÑ9Ô9ð;ñ <ô <Ü<=ð ð „[×Ò˜V [Ñ1Ô1ð5ð 5ñ 2Ô1ð5ð:'ð 'ð 'ð 'ð 'r/   r¿   c                   óJ  — e Zd Zej                             d¦  «         dZdZej                             dee¬¦  «        Z	ej                             dee¬¦  «        Z
e	e
z   Z ej        deej        ¬¦  «        Z ej        dee
z
  ¦  «        Z ej        ee	¦  «        Z ej        eeej        ¬¦  «        Z ej        j        ej        Ž dz  Zej                             ddd	g¦  «        d
„ ¦   «         Zd„ Zd„ Zd„ Zej                             ddd	g¦  «        d„ ¦   «         ZdS )ÚTestNCHrZ   ©rZ   r   r   rh   r9   r`   r   Ú	dist_namer   r   c           	      óì   — t           t          dœ}||         }| j        | j        | j        | j        f\  }}}}t          |                     ||||d¬¦  «        t          j        ||||¦  «        ¦  «         d S )N©r   r   r9   )Úodds)	r   r   rd   r*   Úm1r)   r   rJ   r   )rÃ   r÷   Údistsrè   rd   r*   rû   r)   s           r-   Útest_nch_hypergeomzTestNCH.test_nch_hypergeomy  s‡   € õ
 (:Ý*?ðAð Aˆà�YÔˆØ”f˜dœf d¤g¨t¬vÐ5‰ˆˆ1ˆb�!Ý˜Ÿš  A r¨1°1˜Ñ5Ô5Ý!œ a¨¨B°Ñ2Ô2ñ	4ô 	4ð 	4ð 	4ð 	4r/   c           	      óŠ  — | j         | j        | j        | j        | j        f\  }}}}}t
          j        d„ ¦   «         } ||||||¦  «        \  }}}	t          t          j	        |||||¦  «        |¦  «         t          t          j
        ||||d¬¦  «        |¦  «         t          t          j
        ||||d¬¦  «        |	¦  «         d S )Nc                 ó   ‡‡‡‡‡‡‡— |‰z
  Št          j        d‰‰z
  ¦  «        Št          j        ‰‰¦  «        Šˆˆˆˆfd„Šˆˆˆfd„} |d¦  «        } |d¦  «        } |d¦  «        } ‰| ¦  «        |z  }	||z  }
||z  ||z  dz  z
  }|	|
|fS )Nr   c                 ó`   •— t          ‰| ¦  «        }t          ‰‰| z
  ¦  «        }||z  ‰| z  z  S r§   )Úspecial_binom)rd   Út1Út2rû   Úm2r)   Úws      €€€€r-   ÚfzFTestNCH.test_nchypergeom_fisher_naive.<locals>.pmf_mean_var.<locals>.f�  s7   ø€ Ý" 2 qÑ)Ô)�Ý" 2 q¨1¡uÑ-Ô-�Ø˜B‘w  A¡‘~Ð%r/   c                 ó`   •‡ — t          ˆˆ fd„t          ‰‰dz   ¦  «        D ¦   «         ¦  «        S )Nc              3   ó:   •K  — | ]} ‰|¦  «        |‰z  z  V — Œd S r§   r•   )Ú.0rm   r  r'   s     €€r-   ú	<genexpr>zYTestNCH.test_nchypergeom_fisher_naive.<locals>.pmf_mean_var.<locals>.P.<locals>.<genexpr>–  s3   øè è € Ð@Ð@¨˜1˜1˜Q™4œ4  1¡™9Ð@Ð@Ð@Ð@Ð@Ð@r/   r9   )rÒ   Úrange)r'   r  ÚxlÚxus   `€€€r-   ÚPzFTestNCH.test_nchypergeom_fisher_naive.<locals>.pmf_mean_var.<locals>.P•  s8   øø€ ÝÐ@Ð@Ð@Ð@Ð@­e°B¸¸Q¹Ñ.?Ô.?Ð@Ñ@Ô@Ñ@Ô@Ð@r/   r9   rZ   ©rU   ÚmaximumÚminimum)rd   r*   rû   r)   r  r  ÚP0ÚP1ÚP2rJ   rÜ   rÝ   r  r  r  r  s     ```       @@@@r-   Úpmf_mean_varz;TestNCH.test_nchypergeom_fisher_naive.<locals>.pmf_mean_var‰  sñ   øøøøøøø€ ð �R‘ˆBÝ”˜A˜q ™tÑ$Ô$ˆBÝ”˜A˜rÑ"Ô"ˆBð&ð &ð &ð &ð &ð &ð &ð &ð
Að Að Að Að Að Að Að ��1‘”ˆBØ��1‘”ˆBØ��1‘”ˆBØ�!�A‘$”$˜‘)ˆCØ˜‘7ˆDØ�r‘'˜R "™W q™LÑ(ˆCØ˜˜c�>Ð!r/   Úmr|   Úv)rd   r*   rû   r)   rú   rU   rÙ   r   r   rJ   r~   )
rÃ   rd   r*   rû   r)   rú   r  rJ   rÜ   rÝ   s
             r-   Útest_nchypergeom_fisher_naivez%TestNCH.test_nchypergeom_fisher_naive…  së   € à œF D¤F¨D¬G°T´V¸T¼YÐFÑˆˆ1ˆb�!�Tå	Œð	"ð 	"ñ 
Œð	"ð, &˜ a¨¨B°°4Ñ8Ô8‰ˆˆT�3ÝÕ*Ô.¨q°!°R¸¸DÑAÔAÀ3ÑGÔGÐGÝÕ*Ô0°°B¸¸4ÈÐMÑMÔMØñ	ô 	ð 	åÕ*Ô0°°B¸¸4ÈÐMÑMÔMØñ	ô 	ð 	ð 	ð 	r/   c           
      óì  ‡‡— t           j                             d¦  «         d}d}t           j                             d||¬¦  «        }t           j                             d||¬¦  «        }||z   }t          j        d||j        ¬¦  «        }t          j        d||z
  ¦  «        }t          j        ||¦  «        }t          j        |||j        ¬¦  «        }	t          j        j        |	j        Ž dz  }
d„ Št           j	        ˆfd„¦   «         Št          ¦   «         5 }|                     t          d	¬
¦  «         t          t          j        ||||
¦  «         ‰||||
¦  «        d¬¦  «         d d d ¦  «         n# 1 swxY w Y   t           j	        ˆfd„¦   «         }t          ¦   «         5 }|                     t          d	¬
¦  «         t          t          j        ||||
d¬¦  «         |||||
¦  «        d¬¦  «         d d d ¦  «         n# 1 swxY w Y   t           j	        ˆfd„¦   «         } ||	||||
¦  «        }t          j        |	||||
¦  «        }d\  }}t          j        ||z
  ¦  «        ||t          j        |¦  «        z  z   k     }|                     ¦   «         t          j        |¦  «        dz  k    sJ ‚t+          ||          ||          ||          |
|          ¦  «        D ]”\  }}}}
||z
  } ‰||||
¦  «        \  }}t          j        ||dz   ¦  «        }	 ||	||||
¦  «                             ¦   «         dk     sJ ‚t          t          j        |	||||
¦  «                             ¦   «         d¦  «         Œ•d S )NrZ   rö   rh   r9   r`   r   c                 ón   — | |z
  }t          j        d||z
  ¦  «        }t          j        ||¦  «        }||fS )Nr   r  )r*   rû   r)   r  r  r  r  s          r-   rW   z9TestNCH.test_nchypergeom_wallenius_naive.<locals>.support¶  s:   € Ø�R‘ˆBÝ”˜A˜q ™tÑ$Ô$ˆBÝ”˜A˜rÑ"Ô"ˆBØ�r�6ˆMr/   c                 óx   •‡‡‡‡— | ‰z
  Š ‰| ‰‰‰¦  «        \  }}ˆˆˆˆfd„}t          |||f¬¦  «        j        S )Nc                 ó2   •— | ‰z  d‰| z
  ‰z  z
  ‰z  z   dz
  S ©Nr9   r•   )Úurû   r  r)   r  s    €€€€r-   ÚfunzCTestNCH.test_nchypergeom_wallenius_naive.<locals>.mean.<locals>.funÁ  s(   ø€ Ø˜‘t˜q A a¡C¨¡8™|¨aÑ/Ñ/°!Ñ3Ð3r/   )Úbracket)r   Úroot)	r*   rû   r)   r  r  r  r  r  rW   s	    ```   @€r-   rÜ   z6TestNCH.test_nchypergeom_wallenius_naive.<locals>.mean¼  sm   øøøøø€ à�R‘ˆBØ�W˜Q  A qÑ)Ô)‰FˆB�ð4ð 4ð 4ð 4ð 4ð 4ð 4ð 4õ ˜s¨R°¨HÐ5Ñ5Ô5Ô:Ð:r/   z!invalid value encountered in mean)Úmessageg{®Gáz”?r#   c                 óŠ   •— | |z
  } ‰| |||¦  «        }|||z
  z  }||z
  ||z   |z
  z  }| |z  |z  | dz
  ||z  ||z  z   z  z  S r  r•   )	r*   rû   r)   r  r  r  r‚   rƒ   rÜ   s	           €r-   Úvariancez:TestNCH.test_nchypergeom_wallenius_naive.<locals>.varianceÌ  sl   ø€ à�R‘ˆBØ��Q˜˜A˜qÑ!Ô!ˆAØ�R˜!‘V‘ˆAØ�1‘�q˜2‘v ‘zÑ"ˆAØ�Q‘3�q‘5˜Q˜q™S R¨¡T¨B¨q©D¡[Ñ1Ñ2Ð2r/   r  r|   gš™™™™™©?c                 ót   •‡ ‡‡‡‡‡	— |‰z
  Š	 ‰
|‰‰‰¦  «        \  }}ˆˆ	ˆˆˆ fd„Šˆˆˆ	ˆfd„} |‰ ¦  «        S )Nc                 ój   •— ‰‰‰z
  z  ‰‰‰z
  z
  z   }d| ‰|z  z  z
  ‰z  d| d|z  z  z
  ‰‰z
  z  z  }|S r  r•   )ÚtÚDÚresrû   r  r)   r  rd   s      €€€€€r-   Ú	integrandzHTestNCH.test_nchypergeom_wallenius_naive.<locals>.pmf.<locals>.integrandâ  sQ   ø€ Ø�r˜A‘v‘J "¨¨!©¡*Ñ-�Ø˜˜Q˜q™S™‘z A‘o¨¨1¨q°©s©8©°q¸±sÑ(;Ñ;�Ø�
r/   c                 óŽ   •— t          ‰| ¦  «        }t          ‰‰| z
  ¦  «        }t          ‰dddd¬¦  «        }||z  |d         z  S )Nr   r9   g¼‰Ø—²Òœ<)ÚepsrelÚepsabs)r  r   )rd   r  r  Úthe_integralr*  rû   r  r)   s       €€€€r-   r  z@TestNCH.test_nchypergeom_wallenius_naive.<locals>.pmf.<locals>.fç  sZ   ø€ Ý" 2 qÑ)Ô)�Ý" 2 q¨1¡uÑ-Ô-�Ý# I¨q°!Ø+0¸ð @ñ  @ô  @�à˜B‘w ¨a¤Ñ0Ð0r/   r•   )rd   r*   rû   r)   r  r  r  r  r*  r  rW   s   ` ```   @@€r-   rJ   z5TestNCH.test_nchypergeom_wallenius_naive.<locals>.pmfÝ  s�   øøøøøøø€ à�R‘ˆBØ�W˜Q  A qÑ)Ô)‰FˆB�ðð ð ð ð ð ð ð ð ð
1ð 1ð 1ð 1ð 1ð 1ð 1ð 1ð �1�Q‘4”4ˆKr/   )rÌ   rÌ   r’   )rU   ri   rj   r   rb   rc   r  r  r”   rÙ   r   ÚfilterÚRuntimeWarningr   r   rÜ   r~   rJ   ÚabsrÒ   ÚprodÚziprV   )rÃ   rc   Úmax_mrû   r  r*   r)   r  r  rd   r  Úsupr$  rJ   Úpmf0Úpmf1r¡   r$   ÚirÜ   rW   s                      @@r-   Ú test_nchypergeom_wallenius_naivez(TestNCH.test_nchypergeom_wallenius_naive§  sR  øø€ õ 	Œ	�Š�qÑÔÐØˆØˆÝŒY×Ò˜q %¨eÐÑ4Ô4ˆÝŒY×Ò˜q %¨eÐÑ4Ô4ˆØ�‰GˆÝŒK˜˜1 1¤7Ð+Ñ+Ô+ˆÝŒZ˜˜1˜R™4Ñ Ô ˆÝŒZ˜˜2ÑÔˆÝŒK˜˜B R¤XÐ.Ñ.Ô.ˆÝŒIŒN˜AœGÐ$ QÑ&ˆð	ð 	ð 	õ 
Œð	;ð 	;ð 	;ð 	;ñ 
Œð	;õ Ñ Ô ð 	: CØ�JŠJ•~ØBð ñ Dô Dð DåÕ1Ô6°q¸"¸aÀÑCÔCØ ˜D  B¨¨1Ñ-Ô-°Dð:ñ :ô :ð :ð	:ð 	:ð 	:ñ 	:ô 	:ð 	:ð 	:ð 	:ð 	:ð 	:ð 	:øøøð 	:ð 	:ð 	:ð 	:õ 
Œð	3ð 	3ð 	3ð 	3ñ 
Œð	3õ Ñ Ô ð 	 CØ�JŠJ•~ØBð ñ Dô Dð DåÝ%Ô+¨A¨r°1°aÀÐEÑEÔEØ�˜˜B  1Ñ%Ô%Øðñ ô ð ð	ð 	ð 	ñ 	ô 	ð 	ð 	ð 	ð 	ð 	ð 	øøøð 	ð 	ð 	ð 	õ 
Œð	ð 	ð 	ð 	ñ 
Œð	ð$ ˆs�1�a˜˜Q Ñ"Ô"ˆÝ$Ô(¨¨A¨r°1°aÑ8Ô8ˆà‰
ˆˆdÝŒF�4˜$‘;ÑÔ $¨­b¬f°T©l¬lÑ):Ñ":Ò:ˆØ�uŠu‰wŒw�œ ™œ¨!Ñ+Ò+Ð+Ð+Ð+õ ˜q ! œu b¨!¨¤f¨a°°¬e°Q¸°r´UÑ;Ô;ð 		Pð 		P‰KˆAˆr�1�aà�R‘ˆBØ�W˜Q  A qÑ)Ô)‰FˆB�Ý”	˜"˜b 1™fÑ%Ô%ˆAð �3�q˜!˜R  AÑ&Ô&×*Ò*Ñ,Ô,¨rÒ1Ð1Ð1Ð1ÝÕ1Ô5°a¸¸BÀÀ1ÑEÔE×IÒIÑKÔKÈQÑOÔOÐOÐOð		Pð 		Ps&   ÄAE3Å3E7Å:E7Æ#AHÈHÈHc           	      ó`  — d}d}d}d}t          j        d¦  «        }t          j        g d¢¦  «        }d}d}t          t	          j        |||||¦  «        |d	d	¬
¦  «         t          t	          j        ||||¦  «        |d	¬¦  «         t          t	          j        ||||¦  «        |d¬¦  «         d S )Nr6   rÆ   rw   g      @rÁ   )g9”T€ÆÕå;g»��rþü�<gÝD,P=g³4À VŒ=g'±	,îñð=gC�åõ‹ŒG>gyzžtj…”>gÁÙë«t•×>gÃ+k…†?g ‰Me´GD?gqNZVÉ’o?gœcø2nž‘?g~Z€ÄN¬?g¼EŸ}?YÀ?gõÛÀ¢ßöÊ?gËÎ‡aQÏ?gè€]?…É?gr”…)‚‚º?gÌáÙp{¡?g3Í9	·†y?gµ Ñmç>?gJ®ò•´�-@gñihÞ@rR   )r$   r¡   r#   g•dyáý¥=)rU   rV   Úarrayr   r   rJ   rÜ   rÝ   )	rÃ   r(   r)   r*   rú   r5  rJ   rÜ   rÝ   s	            r-   Útest_wallenius_against_mpmathz%TestNCH.test_wallenius_against_mpmath  sñ   € ð ˆØˆØˆØˆåŒi˜‰mŒmˆÝŒhð /ð /ð /ñ 0ô 0ˆð* "ˆØ ˆõ 	Õ-Ô1°#°q¸!¸QÀÑEÔEÀsØ"¨ð	0ñ 	0ô 	0ð 	0åÕ-Ô2°1°a¸¸DÑAÔAØ 5ð	*ñ 	*ô 	*ð 	*åÕ-Ô1°!°Q¸¸4Ñ@Ô@Ø %ð	)ñ 	)ô 	)ð 	)ð 	)ð 	)r/   c                 óŽ   — t           t          dœ}||         }|                     dddgdggg d¢d¬¦  «        }|j        dk    sJ ‚d S )	Nrù   r6   rÆ   r   rw   )r’   rT   g       @r_   r`   )r   r   rb   rc   )rÃ   r÷   rü   rè   rd   s        r-   Útest_rvs_shapezTestNCH.test_rvs_shape-  sb   € õ (:Ý*?ðAð Aˆà�YÔˆØ�HŠH�R˜˜r˜d R D˜\¨?¨?¨?ÀˆHÑNÔNˆØŒw˜,Ò&Ð&Ð&Ð&Ð&Ð&r/   N)rë   rì   rí   rU   ri   rj   rc   r4  r   rû   r  r*   rb   r)   r  r  r  r  rd   r”   rú   rñ   rò   ró   rý   r  r9  r<  r>  r•   r/   r-   rõ   rõ   k  sœ  € € € € € Ø„I‡N‚N�1ÑÔÐØ€EØ€EØ	Œ×	Ò	˜1˜e¨%Ð	Ñ	0Ô	0€BØ	Œ×	Ò	˜1˜e¨%Ð	Ñ	0Ô	0€BØ
ˆR‰€AØˆŒ�A�q˜qœwÐ'Ñ'Ô'€AØ	ˆŒ�A�q˜‘tÑ	Ô	€BØ	ˆŒ�A�rÑ	Ô	€BØˆŒ�B˜ ¤Ð*Ñ*Ô*€AØˆ2Œ9Œ>˜1œ7Ð# AÑ%€Dð „[×Ò˜[Ø2Ð4KÐLñNô Nð4ð 4ñNô Nð4ð ð  ð  ðDZPð ZPð ZPðx()ð ()ð ()ðT „[×Ò˜[Ø2Ð4KÐLñNô Nð'ð 'ñNô Nð'ð 'ð 'r/   rõ   zmu, q, expected)r   éx   g@	›#á¸)iÜ  r   g"q³V§UÀc                 óf   — d}|||| z   z  }}t          t          j        |||¦  «        |¦  «         d S )Nrw   )r   r   r±   )r»   Úqr+   ra   r)   r,   s         r-   Útest_nbinom_11465rB  :  s@   € ð
 €DØ��t˜B‘w‘€q€Aõ •F”M ! Q¨Ñ*Ô*¨HÑ5Ô5Ð5Ð5Ð5r/   c                  ó  — t          j        ddd¦  «        } d}t          |¦  «                             | ¦  «        }| dz  dk    }t	          |d         |¦  «         t	          |d         d|z
  ¦  «         t          ||          d¦  «         d S )Nr   r9   rg   gš™™™™™é?r“   )rU   r™   r   rJ   r   r   )rd   r,   rJ   r8  s       r-   Útest_gh_17146rD  G  sŠ   € õ 	Œ�A�q˜"ÑÔ€AØ€AÝ
�A‰,Œ,×
Ò
˜1Ñ
Ô
€CØ	
ˆQ‰�!Š€AÝ�C˜”G˜QÑÔÐÝ�C˜”F˜A˜a™CÑ Ô Ð Ý��a�R”˜!ÑÔÐÐÐr/   c                   óÂ   — e Zd Zej                             dg d¢g d¢g d¢g¦  «        d„ ¦   «         Zej                             dg d¢g d¢g d	¢g¦  «        d
„ ¦   «         ZdS )ÚTestBetaNBinomzx, n, a, b, ref)r   g    ÐSAr   rw   g:[Y³³»))rh   r6   r   rw   g ôÍødg?)r   rœ   r   rw   gL&¸[šô>c                 óT   — t          t          j        ||||¦  «        |d¬¦  «         d S )NrI   r#   )r   r   rJ   )rÃ   rd   r)   r‚   rƒ   Úrefs         r-   Útest_betanbinom_pmfz"TestBetaNBinom.test_betanbinom_pmfT  s/   € õ  	�
œ q¨!¨Q°Ñ2Ô2°C¸eÐDÑDÔDÐDÐDÐDr/   zn, a, b, ref)r   rº   r6   g—Í:ÑèoÀ?)r   r^   r^   gÍ>["”°@)rh   rœ   r   gÑ Å>ÿ[ù?c                 óV   — t          t          j        |||d¬¦  «        |d¬¦  «         d S )Nr'   r|   g [n‡ë<r#   )r   r   r~   )rÃ   r)   r‚   rƒ   rH  s        r-   Útest_betanbinom_kurtosisz'TestBetaNBinom.test_betanbinom_kurtosisg  s@   € õ* 	�
Ô(¨¨A¨q¸#Ð>Ñ>Ô>Ø %ð	)ñ 	)ô 	)ð 	)ð 	)ð 	)r/   N)rë   rì   rí   rñ   rò   ró   rI  rK  r•   r/   r-   rF  rF  S  sÄ   € € € € € Ø„[×ÒÐ.ØEÐEÐEØBÐBÐBØHÐHÐHðJñKô KðEð Eñ	Kô KðEð „[×Ò˜^ØCÐCÐCØ:Ð:Ð:Ø?Ð?Ð?ðAñBô Bð)ð )ñ	Bô Bð)ð )ð )r/   rF  c                   ó   — e Zd Zd„ ZdS )ÚTestZipfc                 ó   — t          j        dd¦  «        }|                     t           j        ¦  «        }t	          d¦  «        }|                     |¦  «        }|                     |¦  «        }t          ||¦  «         d S )Nr   rœ   r^   )rU   rV   rä   rå   r
   rJ   r   ræ   s         r-   Útest_gh20692zTestZipf.test_gh20692�  sg   € åŒI�a˜ÑÔˆØ—(’(�2œ8Ñ$Ô$ˆÝ�A‰wŒwˆØ�hŠh�q‰kŒkˆØ—h’h˜wÑ'Ô'ˆÝ�S˜+Ñ&Ô&Ð&Ð&Ð&r/   N)rë   rì   rí   rO  r•   r/   r-   rM  rM  €  s#   € € € € € ð'ð 'ð 'ð 'ð 'r/   rM  )Arñ   Ú	itertoolsÚscipy.statsr   r   r   r   r   r   r	   r
   r   r   r   r   r   r   ÚnumpyrU   Únumpy.testingr   r   r   r   Úscipy.specialr  Úscipy.optimizer   Úscipy.integrater   rò   ró   r.   r4   rC   rN   rX   r[   re   rn   ru   ry   ÚproductÚintÚfloatr„   r†   rŽ   r�   r–   rš   rž   r¢   rª   r¬   r¯   rµ   r¸   r½   r¿   rõ   rB  rD  rF  rM  r•   r/   r-   ú<module>rZ     s¡  ðØ €€€Ø Ð Ð Ð ð"ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð "ð
 Ð Ð Ð ðð ð ð ð ð ð ð ð ð ð ð ð 1Ð 0Ð 0Ð 0Ð 0Ð 0Ø &Ð &Ð &Ð &Ð &Ð &Ø  Ð  Ð  Ð  Ð  Ð  ð
 „×ÒÐ5ð:ð :ð :ñ;ô ;ð,ð ,ñ;ô ;ð,ð „×ÒÐ5ð6ð:ð;ñ<ô <ð
,ð ,ñ<ô <ð
,ð
Hð Hð Hð4)ð )ð )ð	>ð 	>ð 	>ðCð Cð Cð#ð #ð #ðð ð ð-ð -ð -ð&%ð %ð %ð „×Ò˜Ð#4 9Ô#4¸¸U°|°nÀQÑ6FÐ#GÑHÔHðPð Pñ IÔHðPð%ð %ð %ðFð Fð Fð
Eð Eð Eð
*ð *ð *ðð ð ð"Bð Bð Bð:ð :ð :ð „×ÒÐ/ð 
2ð 
2ð 
2ñ 
ô 
ð5ð 5ñ
ô 
ð5ðKð Kð Kð
Hð Hð Hð'ð 'ð 'ð-ð -ð -ð'ð 'ð 'ð['ð ['ð ['ð ['ð ['ñ ['ô ['ð ['ð|L'ð L'ð L'ð L'ð L'ñ L'ô L'ð L'ð^ „×ÒÐ*Ø;Ð;Ð;Ø7Ð7Ð7ð9ñ:ô :ð6ð 6ñ:ô :ð6ð	ð 	ð 	ð*)ð *)ð *)ð *)ð *)ñ *)ô *)ð *)ðZ'ð 'ð 'ð 'ð 'ñ 'ô 'ð 'ð 'ð 'r/   