a
    ;e                    @   sH  d dl Z d dlZd dlZd dlZd dlmZ d dlmZm	Z	m
Z
mZmZmZmZmZmZmZ d dlmZmZmZmZ ee Zd ddddd	d
dddddddddddg dddddddddgZejdddgddd ZG dd dZG d d! d!ZG d"d# d#ZG d$d% d%ZG d&d' d'ZG d(d) d)Z G d*d+ d+Z!G d,d- d-Z"ej#j$ed.d/G d0d1 d1Z%G d2d3 d3Z&ej#'d4ed5d6 Z(d7d8 Z)d9d: Z*d;d< Z+ej#'d=d>d>ggej#'dddgd?d@ Z,ej#'dAdBdCgdDdE Z-ej#'dAdBdCgej#'dFdGdHgej#'dIej.ejgdJdK Z/dLdM Z0dS )N    N)LinAlgError)
assert_assert_raisesassert_equalassert_allcloseassert_warnsassert_no_warningsassert_array_equalassert_array_almost_equalsuppress_warningsIS_WASM)	GeneratorMT19937SeedSequenceRandomState
   Z@bb1636883c2707b51c5b7fc26c6927af4430f2e0785a8c7bc886337f919f9edf	   )
key_sha256posZ@ff682ac12bb140f2d72fba8d3506cf4e46817a0db27aae1683867629031d8d55iV  )seedstepsinitialjumpedi$<i8  Z@16b791a1e04886ccbbb4d448d6ff791267dc458ae599475d08d5cced29d11614i7  Z@a0110a2cf23b56be0feaed8f787a7fc84bef0cb5623003d75b26bdfa1c18002ci  )i|2iEm:iS43i5  i  Z@d306cf01314d51bd37892d874308200951a35265ede54d200f1e065004c3e9eai  Z@0e00ab449f01a5195a83b4aee0dfbc2ce8d46466a640b92e33977d2e42f777f8i  moduleTF)Zscopeparamsc                 C   s   | j S N)param)request r   ]/var/www/sunrise/env/lib/python3.9/site-packages/numpy/random/tests/test_generator_mt19937.pyendpoint'   s    r    c                   @   s<   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd ZdS )TestSeedc                 C   s<   t td}t|dd t td}t|dd d S )Nr           D  )r   r   r   integersselfsr   r   r   test_scalar-   s    zTestSeed.test_scalarc                 C   s   t ttd}t|dd t ttd}t|dd t tdg}t|dd t tdg}t|dd d S )Nr   r"   i  r   r#   r$   r%   )r   r   ranger   r&   nparanger'   r   r   r   
test_array3   s    zTestSeed.test_arrayc                 C   s    t td}t|dd d S )Nr      iz)r   r   r   
random_rawr'   r   r   r   test_seedsequence=   s    zTestSeed.test_seedsequencec                 C   s   t ttd t ttd d S )N      r   	TypeErrorr   
ValueErrorr(   r   r   r   test_invalid_scalarA   s    zTestSeed.test_invalid_scalarc                 C   s0   t ttdg t ttdg t ttg d d S )Nr2   r3   )r/           r4   r7   r   r   r   test_invalid_arrayF   s    zTestSeed.test_invalid_arrayc                 C   s   t ttt d S r   )r   r6   r   r   r7   r   r   r   test_noninstantized_bitgenL   s    z#TestSeed.test_noninstantized_bitgenN)	__name__
__module____qualname__r*   r.   r1   r8   r;   r<   r   r   r   r   r!   ,   s   
r!   c                   @   s   e Zd Zdd Zdd ZdS )TestBinomialc                 C   sB   t jddd}dD ]*}ttd|dk tt||| qd S )N   intdtype)r         ?r/   r   )r,   zerosr   randombinomialr	   )r(   rF   pr   r   r   test_n_zeroQ   s    zTestBinomial.test_n_zeroc                 C   s   t ttjdtj d S )Nr/   )r   r6   rG   rH   r,   nanr7   r   r   r   test_p_is_nanZ   s    zTestBinomial.test_p_is_nanN)r=   r>   r?   rJ   rL   r   r   r   r   r@   P   s   	r@   c                   @   sL   e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dS )TestMultinomialc                 C   s   t dddg d S )Nd   皙?皙?rG   multinomialr7   r   r   r   
test_basic`   s    zTestMultinomial.test_basicc                 C   s   t dg d d S )NrN   )rO   rP           rT   rT   rQ   r7   r   r   r   test_zero_probabilityc   s    z%TestMultinomial.test_zero_probabilityc                 C   sZ   t dtdd  kodk n   tddd}t td|k t t|dk  d S )Nr3      )r   rG   r&   r,   allr(   xr   r   r   test_int_negative_intervalf   s    $z*TestMultinomial.test_int_negative_intervalc                 C   s   ddg}t td|tdjd t td|tdjd t td|tdjd t td|ddgjd t td|djd t td|tdjd tttjd|t	d d S )NrE   r/   r/   rA   rA   rA   rA   rA   rA   rA   )
r   rG   rR   r,   uint32shapearrayr   r5   floatr(   rI   r   r   r   	test_sizel   s    zTestMultinomial.test_sizec                 C   s,   t ttjdddg t ttjdddg d S )NrN   皙?rO   皙?r   r6   rG   rR   r7   r   r   r   test_invalid_probz   s    z!TestMultinomial.test_invalid_probc                 C   s2   t ttjdddg t ttjdgd ddg d S )Nr3   rP   rO   r   rh   r7   r   r   r   test_invalid_n~   s    zTestMultinomial.test_invalid_nc                 C   sx   t d}|t |dd d  }|dd d }ttd}|jd|d}ttd}|jdt |d}t|| d S )Ng      .@r/      ۠iUrN   )pvals)r,   r-   sumr   r   rR   ascontiguousarrayr	   )r(   rI   rm   rG   
non_contigcontigr   r   r   test_p_non_contiguous   s    
z%TestMultinomial.test_p_non_contiguousc                 C   sn   t jg dt jd}||  }ttd}d}tjt|d |	d| W d    n1 s`0    Y  d S )N)
Gz?rs   &.>rt   rt   rt   rt   rt   rt   rt   rC   rl   z-[\w\s]*pvals array is cast to 64-bit floatingmatchr/   )
r,   ra   float32rn   r   r   pytestraisesr6   rR   )r(   rZ   rm   rG   rv   r   r   r   test_multinomial_pvals_float32   s    
z.TestMultinomial.test_multinomial_pvals_float32N)r=   r>   r?   rS   rU   r[   rd   ri   rj   rr   rz   r   r   r   r   rM   _   s   
rM   c                   @   s   e Zd Zdd Zdd Zejdddgdd	 Zejd
g dejdddgejdg ddd Z	dd Z
dd Zdd ZdS )TestMultivariateHypergeometricc                 C   s
   d| _ d S )Ni_ r   r7   r   r   r   setup_method   s    z+TestMultivariateHypergeometric.setup_methodc                 C   s   t ttjdd t ttjg dd t ttjg dd t ttjg dd t ttjg d t ttjdd	gd
dd ttj}|j}||jj	 }t ttj|d d	gd
dd d S )Nr      )rA   rk   r~   r3   )r3   rA   rk   rA   r/   iɚ;e   rW   	marginalsrN   count)
r   r6   rG   multivariate_hypergeometricr,   iinfoint64maxrD   itemsize)r(   Z
int64_infoZ	max_int64Zmax_int64_indexr   r   r   test_argument_validation   s0    z7TestMultivariateHypergeometric.test_argument_validationmethodr   r   c                 C   s
  t t| j}|jg dd|d}t|g d |jg d|d}t|g  |jg dd|d}t|tjdtjd |jg dd|d}t|g d |jg d	d
|d}t|g d g d}|j|t||d}t|| |jg ddd
|d}t|g dgd
  d S )Nr   r   r   r   r   r/   sizer   )r/   r   rC   r/   rA   rk   )r   r   r   rk   )rk   r   r   )r/   r/   r   r/   r/   rk   r~   rW      )	r   r   r   r   r	   r,   emptyr   rn   )r(   r   rG   rZ   colorsr   r   r   test_edge_cases   s(    

z.TestMultivariateHypergeometric.test_edge_casesnsample)      -   7   r   )rW   rA   rk   iI c                 C   s   t t| j}tg d}|j||||d}t|trF|f|j }n
||j }t	|j| t
|dk  t
||k  t|jddtj||td t|tr|dkrt|jdd|| |  dd	d
 d S )N)r   rW      r   r   r   r3   axis)Z
fill_valuerD   順 MbP?{Gzt?)rtolZatol)r   r   r   r,   ra   r   
isinstancerB   r`   r   r   rX   r	   rn   fullr   mean)r(   r   r   r   rG   r   sampleexpected_shaper   r   r   test_typical_cases   s&    


z1TestMultivariateHypergeometric.test_typical_casesc                 C   sZ   t t| j}|jg ddddd}tg dg dg dg dg dg}t|| d S )Nr   rW   r   r   )rA   r/   rA   r/   r/   rk   rA   r   rk   r   r   r   r   r,   ra   r	   r(   rG   r   expectedr   r   r   test_repeatability1   s    
z2TestMultivariateHypergeometric.test_repeatability1c                 C   sZ   t t| j}|jg ddddd}tg dg dg dg dg d	g}t|| d S )
Nr      2   r   rW   r   r   )r         )      r   )r         )r      r   r   r   r   r   r   test_repeatability2  s    
z2TestMultivariateHypergeometric.test_repeatability2c                 C   sZ   t t| j}|jg ddddd}tg dg dg dg dg d	g}t|| d S )
Nr   r   rW   r   r   )rA   rk   r   )rW   rk   r~   )rA   rW   rW   )r/   rW      r   r   r   r   r   test_repeatability3  s    
z2TestMultivariateHypergeometric.test_repeatability3N)r=   r>   r?   r}   r   rx   markparametrizer   r   r   r   r   r   r   r   r   r{      s   !
"r{   c                   @   s,   e Zd Zdd Zdd Zdd Zdd Zd	S )
TestSetStatec                 C   sT   d| _ tt| j | _| jj| _| jj| _| jd | jd d | jd d f| _d S )NIbit_generatorstatekeyr   )r   r   r   rgr   r   Zlegacy_stater7   r   r   r   r}     s    

zTestSetState.setup_methodc                 C   s<   | j jdd}| j| j_| j jdd}tt||k d S Nrk   r   )r   standard_normalr   r   r   r,   rX   )r(   oldnewr   r   r   test_gaussian_reset&  s    
z TestSetState.test_gaussian_resetc                 C   sL   | j   | jj}| j jdd}|| j_| j jdd}tt||k d S r   )r   r   r   r   r   r,   rX   )r(   r   r   r   r   r   r    test_gaussian_reset_in_media_res-  s    
z-TestSetState.test_gaussian_reset_in_media_resc                 C   s   | j dd d S NrE   )r   negative_binomialr7   r   r   r   test_negative_binomial8  s    z#TestSetState.test_negative_binomialN)r=   r>   r?   r}   r   r   r   r   r   r   r   r     s   	r   c                   @   s\  e Zd ZejZeejej	ej
ejejejejejg	Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Zdd Zdd Zdd Zdd Zejddeg dfdeg dfdeg dfgdd Z d d! Z!d"d# Z"d$d% Z#d&d' Z$d(d) Z%d*d+ Z&d,d- Z'ejj(ejd.d/d0ejd1fd/d2ej	d3fd4d5ej
d6fd7d8ejd9fgd:d; Z)d<S )=TestIntegersc                 C   s   t t| jd|td d S Nr/   r    rD   )r   r5   rfuncrb   )r(   r    r   r   r   test_unsupported_typeE  s    z"TestIntegers.test_unsupported_typec              	   C   s@  | j D ]2}|tu rdn
t|j}|tu r0dnt|jd }|rL|d n|}tt| j|d |||d tt| j||d ||d tt| j||||d tt| jdd||d tt| j|d g|||d tt| j|g|d g||d tt| j|g|g||d tt| jddg||d tt| j|d g|g||d qd S Nr   rA   r/   r   	itypeboolr,   r   minr   r   r6   r   r(   r    dtlbndubndr   r   r   test_bounds_checkingH  s>    z!TestIntegers.test_bounds_checkingc              	   C   s   | j D ]}|tu rdn
t|j}|tu r.dnt|j|  }tt| j|d gd |gd ||d tt| j|gd |d gd ||d tt| j||gd ||d tt| jdgd d||d qd S r   r   r   r   r   r   test_bounds_checking_arraya  s    
z'TestIntegers.test_bounds_checking_arrayc              	   C   s0  | j D ]"}|tu rdn
t|j}|tu r0dnt|jd }|rL|d n|}| }|d }t| j||| d||d| t| j|g|| d||d| |}t| j||| d||d| t| j||| gd||d| || d }t| j||| d||d| t| j|g|| gd||d| qd S )Nr   rA   r/   r"   r   r    rD   r   r   r,   r   r   r   r   r   )r(   r    r   r   r   Zis_opentgtr   r   r   test_rng_zero_and_extremeso  sL    z'TestIntegers.test_rng_zero_and_extremesc                 C   s  d}| j D ]}|tu rdn
t|j}|tu r4dnt|jd }|rP|d n|}|d }t| j|g|d g||d| t| j|g| |d g| |d| t| j|g| |d g| ||d| |}t| j|g|d g||d| t| j|g| |d g| |d| t| j|g| |d g| ||d| || d }t| j|g|d g||d| t| j|g| |d g| |d| t| j|g| |d g| ||d| q
d S )Nr"   r   rA   r/   r   rD   rC   r   )r(   r    r   r   r   r   r   r   r   r    test_rng_zero_and_extremes_array  sj    z-TestIntegers.test_rng_zero_and_extremes_arrayc                 C   s   | j D ]}|tu rdn
t|j}|tu r.dnt|jd }|rJ|d n|}z| j||||d W q ty } ztdt	| W Y d }~qd }~0 0 qd S Nr   rA   r/   r   zMNo error should have been raised, but one was with the following message:

%s
r   r   r,   r   r   r   r   	ExceptionAssertionErrorstrr(   r    r   r   r   er   r   r   test_full_range  s    
zTestIntegers.test_full_rangec                 C   s   | j D ]}|tu rdn
t|j}|tu r.dnt|jd }|rJ|d n|}z| j|gd |g||d W q ty } ztdt	| W Y d }~qd }~0 0 qd S r   r   r   r   r   r   test_full_range_array  s    
z"TestIntegers.test_full_range_arrayc              	   C   s   t t }| jdd  D ]F}dD ]<}| jd|| d||d}t| |k  t| dk q q| jdd| d|td}t| dk  t| dk d S )Nr/   )r~   r      rA   i   r   r   )r   r   r   r   r   r   r   r   )r(   r    rG   r   r   valsr   r   r   test_in_bounds_fuzz  s    
z TestIntegers.test_in_bounds_fuzzc           
      C   s   | j D ]}|tu rdn
t|j}|tu r.dnt|jd }|rJ|d n|}d}ttd}|j|||||d}ttd}|j|g|g|||d}ttd}|j|g| |g| |||d}	t	|| t	||	 qd S )Nr   rA   r/   r"     r   )
r   r   r,   r   r   r   r   r   r&   r	   )
r(   r    r   r   r   r   rG   ZscalarZscalar_arrayra   r   r   r   test_scalar_array_equiv  s*    

z$TestIntegers.test_scalar_array_equivc              
   C   s   dddddddddd	}| j dd  D ]v}ttd}tjd	krZ|jd
d| d||d}n|jd
d| d||d }t|	 }t
|t|j |k q&ttd}|jd
d| d|tdtj}t|	 }t
|ttj |k d S )NZ@053594a9b82d656f967c54869bc6970aa0358cf94ad469c81478459c6a90eee3Z@54de9072b6ee9ff7f20b58329556a46a447a8a29d67db51201bf88baa6e4e5d4Z@d3a0d5efb04542b25ac712e50d21f39ac30f312a5052e9bbb1ad3baa791ac84bZ@14e224389ac4580bfbdccb5697d6190b496f91227cf67df60989de3d546389b1Z@0e203226ff3fbbd1580f15da4621e5f7164d0d8d6b51696dd42d004ece2cbec1)	r   int16int32r   int8uint16r_   uint64uint8r/   r   littler   r   r"   r   rA   )r   r   r   sys	byteorderr&   byteswaphashlibsha256	hexdigestr   r,   rD   namer   viewr   )r(   r    r   r   rG   valresr   r   r   test_repeatability  s:    


zTestIntegers.test_repeatabilityc                 C   s   | j D ]}|ttjfv rdn
t|j}|ttjfv r:dnt|jd }|rV|d n|}ttd}|j	||d||d}ttd}|j	|gd |||d}t
|| ttd}|j	|gd |gd ||d}t
|| qd S )Nr   rA   r/   r   r"   r   r   )r   r   r,   bool_r   r   r   r   r   r&   r	   )r(   r    r   r   r   rG   r   Zval_bcr   r   r   test_repeatability_broadcasting  s$    
"
z,TestIntegers.test_repeatability_broadcastingzbound, expectedr$   )iui,YQiXgi$#Pl   bA i{El   /Yi r:   )ui,YQXgi%#Pl   bA i|El   0Yi        )r  r  l   bA l   1Yi i{,miuHl   HM c                 C   sJ   d t |fD ]8}ttd}|j||d}t||d ur:|n|d  qd S )Nr   r   r   )lenr   r   r&   r   )r(   boundr   r   rG   rZ   r   r   r   !test_repeatability_32bit_boundary%  s    
z.TestIntegers.test_repeatability_32bit_boundaryc              
   C   s   t g dg dg dgg dg dg dgg dg dg d	gg d
g dg dgg dg dg dgg}dD ]J}ttd}|jdgdgdggg d|d}t||d ur|n|d  qrd S )N)i
`l   Y/ i0b)i)H{TiP-iTq)l   C iZ'l   !g )l   k40 l   n2c l   DBJ )i|(iy+l   _&A )ip@GiG^l   J|xK )l   
f l   Dmv l   E"g )i>{l   wT2L iJQ)iQ//itLol   k )i\siJi)i=l   }D l   E9l )l   
+x i#Jl   b )i0i7Ol   rn )ixiiAi.t)l   q/q i;-/i)N)rW   rk   rk   i90  r3   r   r/   )r$   r:   r  r   )r,   ra   r   r   r&   r	   )r(   desiredr   rG   rZ   r   r   r   .test_repeatability_32bit_boundary_broadcasting4  s8    
z;TestIntegers.test_repeatability_32bit_boundary_broadcastingc              
   C   s4  t jddd fdddft jdi}|D ]}|| D ]}|\}}|| }t |gd g}t |gd }tttj||||d	 tttj||||d	 tttj||||d	 tttj||||d	 t j|gd gtd
}	t j|gd td
}
tttj|	|||d	 tttj||
||d	 tttj|	|
||d	 q2q$d S )Nr   rA   A   )r3               r   r   r   r   ))r   l            )l         r  r  r  )        r  r   r   rC   )	r,   r   r   ra   r   r6   rG   r&   object)r(   r    ZconfigsrD   configlowhighZlow_aZhigh_aZlow_oZhigh_or   r   r   &test_int64_uint64_broadcast_exceptionsK  s@    
z3TestIntegers.test_int64_uint64_broadcast_exceptionsc                 C   sd   t j}t t jj}t t t jj}t t t jjd | }tj||||d}t|| d S r   )r,   r   r   r   r   rG   r&   r   )r(   r    r   r   r   r   actualr   r   r   test_int64_uint64_corner_caseg  s    z*TestIntegers.test_int64_uint64_corner_casec                 C   s  | j D ]x}|tu rdn
t|j}|tu r.dnt|jd }|rJ|d n|}|tu r\tjn|}| j||||d}t|j	| qtt
fD ]v}|tu rdn
t|j}|tu rdnt|jd }|r|d n|}| j||||d}t|drJ tt|| qd S )Nr   rA   r/   r   rD   )r   r   r,   r   r   r   r  r   r   rD   rB   hasattrtyper(   r    r   r   r   r   r   r   r   test_respect_dtype_singleton~  s    
z)TestIntegers.test_respect_dtype_singletonc                 C   s   | j D ]}|tu rdn
t|j}|tu r.dnt|jd }|rJ|d n|}|tu r\tjn|}| j|g|g||d}t|j	| | j|gd |gd ||d}t|j	| qd S r   )
r   r   r,   r   r   r   r  r   r   rD   r  r   r   r   test_respect_dtype_array  s    
z%TestIntegers.test_respect_dtype_arrayc                 C   s   | j D ]}| jddd||d}|jdks,J |j|ks:J | jddd||djdksXJ ttjddddjd ttjddddjd ttjddddjd qd S )Nr   rk   r   r~   r   r   r   r   )r   r   r`   rD   r   rG   r&   )r(   r    r   r   r   r   r   test_zero_size  s    
zTestIntegers.test_zero_sizec                 C   sR   t jdkrdnd}tt" tjddd|d W d    n1 sD0    Y  d S )Nbigz<i4z>i4r      r   r   )r   r   rx   ry   r6   rG   r&   )r(   Zother_byteord_dtr   r   r   test_error_byteorder  s    z!TestIntegers.test_error_byteorderzsample_size,high,dtype,chi2maxi@KL rW   g     @_@r   g     b@i i	  g     ȩ@i  g     d@c           
      C   sN   t j|||d}tj|dd\}}|| }|| d |  }	|	|k sJJ d S )Nr   T)Zreturn_countsrA   )rG   r&   r,   uniquern   )
r(   Zsample_sizer  rD   Zchi2maxsamplesvaluescountsr   Zchi2r   r   r   $test_integers_small_dtype_chisquared  s
    
z1TestIntegers.test_integers_small_dtype_chisquaredN)*r=   r>   r?   rG   r&   r   r   r,   r   r   r   r   r   r_   r   r   r   r   r   r   r   r   r   r   r   r   r   r  rx   r   r   ra   r  r
  r  r  r  r  r  r"  slowr(  r   r   r   r   r   >  sN   #	
r   c                   @   sd  e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
ejdejejfejejfgdd Zdd Zdd Zdd Zdd Zdd Zdd Zd d! Zd"d# Zd$d% Zd&d' Zd(d) Zd*d+ Zd,d- Z d.d/ Z!d0d1 Z"d2d3 Z#d4d5 Z$d6d7 Z%d8d9 Z&d:d; Z'd<d= Z(d>d? Z)d@dA Z*dBdC Z+dDdE Z,ejdFe-e.gejdGdHe/g dIg dJgfdKe/g dLg dMgfdNe/g dOg dPgfgdQdR Z0dSdT Z1dUdV Z2ejdWdXdYgdZd[ Z3d\d] Z4d^d_ Z5d`da Z6dbdc Z7ddde Z8dfdg Z9dhdi Z:djdk Z;dldm Z<dndo Z=ejj>dpdq Z?ejdrg dsg dtg dug dvg dwg dxdygg dzgd{d| Z@d}d~ ZAdd ZBdd ZCdd ZDdd ZEdd ZFdd ZGdd ZHdd ZIdd ZJdd ZKdd ZLdd ZMdd ZNdd ZOdd ZPdd ZQejdeRdyddejSdgdd ZTdd ZUejjVeWddejdg ddd ZXejddKgdggfdgdNggfgdd ZYejdg ddd ZZdd Z[dd Z\dd Z]dd Z^dd Z_dd Z`dd ZaddÄ Zbddń ZcddǄ ZdddɄ Zedd˄ Zfdd̈́ Zgddτ Zhddф Ziddӄ ZjddՄ Zkddׄ Zlddل Zmddۄ Zndd݄ Zodd߄ Zpdd Zqdd Zrdd Zsdd Ztdd Zudd Zvdd Zwdd Zxdd Zydd Zzdd Z{dd Z|dd Z}dd Z~dd Zejddd gdd Zejddej dejdgejdg ddd	 Zd
d Zdd Zdd Zdd ZdHS (  TestRandomDistc                 C   s
   d| _ d S )Nr   r|   r7   r   r   r   r}     s    zTestRandomDist.setup_methodc                 C   sH   t t| j}|jdddd}tddgddgd	d
gg}t|| d S )Nc   rk   rA   r   )   %   ir   r   r   r&   r,   ra   r	   r(   rG   r  r	  r   r   r   test_integers  s    zTestRandomDist.test_integersc                 C   sR   t t| j}|jdddtjd}tjddgddgd	d
ggtjd}t|| d S )Nr   r,  r-  r   r      F   D   r   r0  rC   )r   r   r   r&   r,   r_   ra   r	   r4  r   r   r   test_integers_masked  s    "z#TestRandomDist.test_integers_maskedc                 C   sJ   t t| j}|jddddd}tddgdd	gd
dgg}t|| d S )Nr+  r,  r-  T)r   r    r.  r/  r0  &   r2  ir3  r4  r   r   r   test_integers_closed  s    z#TestRandomDist.test_integers_closedc                 C   s:   t jtdjtdjdd}tdj}t|| d S )NlTr    )rG   r&   r,   r   r   r   )r(   r  r	  r   r   r   test_integers_max_int  s
    z$TestRandomDist.test_integers_max_intc                 C   sn   t t| j}|d}tddgddgddgg}t||dd	 t t| j}| }t||d
 dd	 d S )Nr-  gG0Ը?gæ?g2?gu	A?gH
>H?g;dE?r   decimalr  r   r   r   rG   r,   ra   r
   r4  r   r   r   test_random  s    

zTestRandomDist.test_randomc                 C   sF   t t| j}|d}tddgddgddgg}t||dd	 d S )
Nr-  CԸ?gP?gtX?g81
A?grO>H?g88dE?r   r?  rA  r4  r   r   r   test_random_float  s    

z TestRandomDist.test_random_floatc                 C   s2   t t| j}|jtjd}d}t||dd d S )NrC   rC  r   r?  )r   r   r   rG   r,   rw   r
   r4  r   r   r   test_random_float_scalar  s    z'TestRandomDist.test_random_float_scalarzdtype, uint_view_typec                 C   sN   t t| j}|jd|d}t||d@ }d|  k rDdk sJn J d S )Nr   rC   r/   i$^  i,e  )r   r   r   rG   r,   Zcount_nonzeror   )r(   rD   Zuint_view_typerG   r   Znum_ones_in_lsbr   r   r   test_random_distribution_of_lsb  s    z.TestRandomDist.test_random_distribution_of_lsbc                 C   s   t ttjdd d S )Nr   rC   )r   r5   rG   r7   r   r   r   test_random_unsupported_type  s    z+TestRandomDist.test_random_unsupported_typec                 C   s<   t t| j}|dd}tjg dtjd}t|| d S )Nr~   )r   r   rA   rA   rC   r   r   r   choicer,   ra   r   r	   r4  r   r   r   test_choice_uniform_replace  s    z*TestRandomDist.test_choice_uniform_replacec                 C   sD   t t| j}|jddg dd}tjg dtjd}t|| d S )Nr~   )皙?rK  皙?rL  rI   )r   r/   r   r/   rC   rH  r4  r   r   r   test_choice_nonuniform_replace  s    z-TestRandomDist.test_choice_nonuniform_replacec                 C   sl   t t| j}|jdddd}tjg dtjd}t|| |jddddd}tjdtjd}t|| d S )Nr~   rk   Freplacer   rC   )rP  shuffle)	r   r   r   rI  r,   ra   r   r	   r-   r4  r   r   r   test_choice_uniform_noreplace"  s    
z,TestRandomDist.test_choice_uniform_noreplacec                 C   sF   t t| j}|jdddg dd}tjg dtjd}t|| d S )Nr~   rk   F)rL  333333?rE   rL  rP  rI   )r   rA   rk   rC   rH  r4  r   r   r    test_choice_nonuniform_noreplace+  s    z/TestRandomDist.test_choice_nonuniform_noreplacec                 C   s:   t t| j}|g dd}tg d}t|| d S )N)abcdr~   )rV  rV  rX  rX  r   r   r   rI  r,   ra   r	   r4  r   r   r   test_choice_noninteger1  s    z%TestRandomDist.test_choice_nonintegerc                 C   s\   t t| j}|ddgddgddgddggd}tddgddgddgg}t|| d S )	Nr   r/   rA   rk   r~   rW   r   r   rZ  r4  r   r   r   )test_choice_multidimensional_default_axis7  s    $z8TestRandomDist.test_choice_multidimensional_default_axisc                 C   s^   t t| j}|jddgddgddgddggddd	}tdgdgdgdgg}t|| d S )
Nr   r/   rA   rk   r~   rW   r   r   r   rZ  r4  r   r   r   (test_choice_multidimensional_custom_axis=  s    (z7TestRandomDist.test_choice_multidimensional_custom_axisc                 C   s  t j}tt|dd tt|dd tt|g d tt|g ddddgddggd tt|ddgdg d	d tt|ddgdd
dgd tt|ddgdddgd tt|g dddd tt|g dddd tt|g dddd tt|g dddd tt|g dddg dd d S )Nr3   rk         @)r/   rA   rk   r~         ?rM  r/   rA   )rK  rK  rO   re   rf   rK  r   r~   FrO  r9   )r3   )r3   r/   )r/   r   r   rT  )rG   rI  r   r6   )r(   r   r   r   r   test_choice_exceptionsC  s"    z%TestRandomDist.test_choice_exceptionsc              	   C   s   ddg}t ttjddd t ttjddd t ttjdd|d t ttjdd|d t ttjddgdd t tjd gddd u  tddg}tjdtd	}||d
< t tj|dd|u  t }t ttjd|dd  t ttjd|dd  t ttjd|d|d  t ttjd|d|d  t ttjddg|dd  t tjd g|ddj	d
k tddg}tjdtd	}||d
< t tj||dd
 |u  d}g d}ttjd|ddj| ttjd|ddj| ttjd|d|dj| ttjd|d|dj| ttjtd|ddj| ttjd
d
ddjd ttjd
dd
djd ttjddd
djd ttjd
d
djd ttjg ddjd ttjddgddjd tttjg d d S )NrL  rg   rA   TrO  FrT  r/   rC   r   r   )rL  rL  rL  rL  rK  rO   r   r  r   r  r  r   rV  rW  )r   r,   ZisscalarrG   rI  ra   r   r  tuplendimitemr   r`   r-   r&   r   r6   )r(   rI   rV  arrr)   r   r   r   test_choice_return_shapeU  sL     z'TestRandomDist.test_choice_return_shapec                 C   s,   t g d}g d}tttj||d d S )N)*   r/   rA   )NNNrM  )r,   ra   r   r6   rG   rI  )r(   rV  rI   r   r   r   test_choice_nan_probabilities  s    z,TestRandomDist.test_choice_nan_probabilitiesc              	   C   s   t dd }d|dd d< tt| j}|jdd|d d d d}tt| j}|jddt |d d d d}t|| d S )Nr   rW   r^  r/   rA   rk   rM  )r,   onesr   r   r   rI  ro   r	   )r(   rI   rG   rp   rq   r   r   r   test_choice_p_non_contiguous  s     z+TestRandomDist.test_choice_p_non_contiguousc                 C   s   t dd }tdd}|jt jks*J tjdddd}|jt jksJJ tjdd|d}|jt jksjJ tjdd|dd}|jt jksJ d S )Nr~   g      @rA   FrO  rM  )rI   rP  )r,   rh  rG   rI  rD   r   )r(   rI   r  r   r   r   test_choice_return_type  s    z&TestRandomDist.test_choice_return_typec                 C   sZ   d}t t| j}|jdddd}tjdkr4| }t|	t
j }t||k d S )NZ@4266599d12bfcfb815213303432341c06b4349f5455890446578877bb322e222'  r#  FrO  r   )r   r   r   rI  r   r   r   r   r   r   r,   r   r   r   )r(   Zchoice_hashrG   r  r   r   r   r   test_choice_large_sample  s    
z'TestRandomDist.test_choice_large_samplec                 C   s*   t t| j}|d}d}t|| d S )Nr   s
   	8%)r   r   r   bytesr   r4  r   r   r   
test_bytes  s    
zTestRandomDist.test_bytesc                 C   s   dd dd dd dd dd dd dd d	d d
d dd dd fD ]B}t t| j}|g d}|| |}|g d}t|| qFd S )Nc                 S   s
   t g S r   )r,   ra   rZ   r   r   r   <lambda>      z-TestRandomDist.test_shuffle.<locals>.<lambda>c                 S   s   | S r   r   ro  r   r   r   rp    rq  c                 S   s   t | t jS r   )r,   asarrayastyper   ro  r   r   r   rp    rq  c                 S   s   t | t jS r   )r,   rr  rs  rw   ro  r   r   r   rp    rq  c                 S   s   t | t jS r   )r,   rr  rs  Z	complex64ro  r   r   r   rp    rq  c                 S   s   t | tS r   )r,   rr  rs  r  ro  r   r   r   rp    rq  c                 S   s   dd | D S )Nc                 S   s   g | ]}||fqS r   r   .0ir   r   r   
<listcomp>  rq  ATestRandomDist.test_shuffle.<locals>.<lambda>.<locals>.<listcomp>r   ro  r   r   r   rp    rq  c                 S   s   t dd | D S )Nc                 S   s   g | ]}||gqS r   r   rt  r   r   r   rw    rq  rx  )r,   rr  ro  r   r   r   rp    rq  c                 S   s   t | | gjS r   )r,   vstackTro  r   r   r   rp    rq  c                 S   s*   t dd | D dtfdtfgt jS )Nc                 S   s   g | ]}||fqS r   r   rt  r   r   r   rw    rq  rx  rV  rW  )r,   rr  rB   r   Zrecarrayro  r   r   r   rp    s   c                 S   s(   t dd | D dtdfdt jdfgS )Nc                 S   s   g | ]}||fqS r   r   rt  r   r   r   rw    rq  rx  rV  r/   rW  )r,   rr  r  r   ro  r   r   r   rp    s   

r/   rA   rk   r~   rW   r   r   r   r   r   
r~   r/   r   r   r   rW   rk   r   rA   r   )r   r   r   rQ  r	   )r(   convrG   alistr  r	  r   r   r   test_shuffle  s$    
zTestRandomDist.test_shufflec                 C   s   t t| j}tdd}|j|dd tg dg dg dg dg}t|| t t| j}tdd}|j|d	d t|| d S )
Nr   r~   r~   r/   r   r   rk   r/   rA   r~   r   rW   r   r      r   r   r   r   r   r   r3   )	r   r   r   r,   r-   reshaperQ  ra   r	   r4  r   r   r   test_shuffle_custom_axis  s    

z'TestRandomDist.test_shuffle_custom_axisc                 C   sT   t t| j}tg d}dD ],}tg d}|j||d t|| q"d S )N)r   r   )r   r/   r   )r   r   r   r,   ra   r  rQ  r	   )r(   rG   r	  r   r  r   r   r   test_shuffle_custom_axis_empty  s    z-TestRandomDist.test_shuffle_custom_axis_emptyc                 C   s^   t ddd}| }tt| j}|j|dd tt| j}||j t	|| d S )Nr   rA   r   r/   r   )
r,   r-   r  copyr   r   r   rQ  rz  r	   )r(   y1y2rG   r   r   r   test_shuffle_axis_nonsquare  s    z*TestRandomDist.test_shuffle_axis_nonsquarec                 C   s   t jt tddd d d}t jt dd d d}| }| }tdD ]d}t| t	t
|j|j  t
|j|j   t| t	t
|j|j  t
|j|j   qVd S )Nr   )rW   r~   rk   r/   r3   r   )r,   maZmasked_valuesr  r+   r-   r  rG   rQ  r   sorteddatamask)r(   rV  rW  Za_origZb_origrv  r   r   r   test_shuffle_masked  s    "
 
 z"TestRandomDist.test_shuffle_maskedc                 C   s   t t| j}td}ttj|j|d tdd}ttj|j|d tt	|j|t
ddd  g dg dg}tt|j|d td}tt	|j| td	}ttj|j|d d S )
Nr   r/   r   rk   rk   rk   rA   r   )r~   rW   r   r-  )r   r   r   r,   r-   r   	AxisErrorrQ  r  r5   sliceNotImplementedErrorra   rh  r(   rG   rd  r   r   r   test_shuffle_exceptions  s    


z&TestRandomDist.test_shuffle_exceptionsc                 C   s\   t t| j}td}d|j_tjt	dd |
| W d    n1 sN0    Y  d S )NrW   F	read-onlyru   )r   r   r   r,   rF   flags	writeablerx   ry   r6   rQ  )r(   rG   rV  r   r   r   test_shuffle_not_writeable  s
    
z)TestRandomDist.test_shuffle_not_writeablec           	      C   s   t t| j}g d}||}g d}t|| t t| j}tg dj}||}t|t|j d}ttj	|j| d}ttj	|j| t t| j}d}g d}||}t|| d S )Nr|  r}  Zabcdg333333?r   )
rk   r   r   r   r   r~   rA   rW   r/   r   )
r   r   r   permutationr	   r,   Z
atleast_2drz  r   r  )	r(   rG   r  r  r	  Zarr_2dZ	bad_x_strZbad_x_floatZinteger_valr   r   r   test_permutation
  s$    



zTestRandomDist.test_permutationc                 C   s   t dd}t g dg dg dg dg}tt| j}|j|dd}t|| tt| j}|j|d	d}t|| d S )
Nr   r  r  r  r  r  r/   r   r3   )	r,   r-   r  ra   r   r   r   r  r	   )r(   rV  r	  rG   r  r   r   r   test_permutation_custom_axis#  s    

z+TestRandomDist.test_permutation_custom_axisc                 C   sh   t t| j}td}ttj|j|d tdd}ttj|j|d tt	|j|t
ddd  d S )Nr   r/   r   r  rk   rA   )r   r   r   r,   r-   r   r  r  r  r5   r  r  r   r   r   test_permutation_exceptions0  s    
z*TestRandomDist.test_permutation_exceptionsrD   zaxis, expectedN)rk   r   r   r   r   r  )r   r~   rA   rW   r/   r   r   )r   r/   rA   r   r   r  )r   r   r   rk   r~   rW   r/   )rW   rk   r~   r   rA   r/   )r  r   r   r   r   r   c                 C   s   t t| j}tddd|}|j|||d t|| t t| j}tddd|}|j||d}|j	|ksJ t|| d S )Nr   rA   r   r   outr   )
r   r   r   r,   r-   r  rs  permutedr	   rD   )r(   rD   r   r   rG   rZ   yr   r   r   test_permuted8  s    	
zTestRandomDist.test_permutedc                 C   s   t t| j}tddd}| }|d d d d df }|j|d|d}tg dg dg}t	|| ||d d d d df< t	|| d S )	N   rA   r  rk   r/   r  )r   r   rk   r   )r   r   r  r   )
r   r   r   r,   r-   r  r  r  ra   r	   )r(   rG   Zx0x1rZ   r  r   r   r   r   test_permuted_with_stridesL  s    

z)TestRandomDist.test_permuted_with_stridesc                 C   s   t g }t|g  d S r   )rG   r  r	   )r(   r  r   r   r   test_permuted_emptyY  s    
z"TestRandomDist.test_permuted_emptyoutshaper   rW   c                 C   s^   t g d}t j||jd}tjtdd tj||d W d    n1 sP0    Y  d S )Nr   rC   z
same shaperu   r  )	r,   ra   rF   rD   rx   ry   r6   rG   r  )r(   r  rV  r  r   r   r   "test_permuted_out_with_wrong_shape]  s    z1TestRandomDist.test_permuted_out_with_wrong_shapec                 C   s\   t jdt jd}t d}tjtdd  tj|d|d W d    n1 sN0    Y  d S )N)rk   rW   rC   zCannot castru   r/   r  )	r,   rF   r   rh  rx   ry   r5   rG   r  )r(   r  rZ   r   r   r   !test_permuted_out_with_wrong_typed  s    
z0TestRandomDist.test_permuted_out_with_wrong_typec                 C   sT   t d}d|j_tjtdd  tj|d|d W d    n1 sF0    Y  d S )N)rA   rW   Fr  ru   r/   r  )	r,   rF   r  r  rx   ry   r6   rG   r  rY   r   r   r   test_permuted_not_writeablej  s    
z*TestRandomDist.test_permuted_not_writeablec                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )NrL  rg   r-  r   gW<$=g^
=g/:Ћ?gG!Gc>g'`2?gU?r   r?  )r   r   r   betar,   ra   r
   r4  r   r   r   	test_betap  s    zTestRandomDist.test_betac                 C   sp   t t| j}|jdddd}tddgddgdd	gg}t|| t t| j}|dd}d}t|| d S )
NgZd;Y@gv/?r-  r   rf  r0  0   ,   r   )r   r   r   rH   r,   ra   r	   r4  r   r   r   test_binomialy  s    

zTestRandomDist.test_binomialc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nr   r-  r   gmE~@@g
1ςC@grӳL@gm%xL@g<kQK@gtK@r   r?  )r   r   r   	chisquarer,   ra   r
   r4  r   r   r   test_chisquare  s    
zTestRandomDist.test_chisquarec                 C   s   t t| j}tddg}|j|dd}tddgddggd	d
gddggddgddggg}t||dd tddg}tt|j| t t| j}tddg}||}t||d dd d S )NI<I@EZC@r-  r   g28\h?gɎG/?gAp?g|);?gLX?gfOh?g됷W?g(ސP?gz$Q?gmWt ?g搝K?g4hz?r   r?  HzG?缉ؗҜr  )	r   r   r   r,   ra   	dirichletr
   r   r6   )r(   rG   alphar  r	  Z	bad_alphar   r   r   test_dirichlet  s(    

zTestRandomDist.test_dirichletc                 C   s   t ddg}tt|t djd tt|t djd tt|t djd tt|ddgjd tt|djd tt|t djd tttj|t	d d S )Nr  r  r/   r\   rA   r]   r^   )
r,   ra   r   rG   r  r_   r`   r   r5   rb   rc   r   r   r   test_dirichlet_size  s    z"TestRandomDist.test_dirichlet_sizec                 C   s   t ddg}tttj| tttjddgg tttjdgdgg tttjdgdggdgdggg tttjt ddgddgg d S )Nr  r  rW   r/   )r,   ra   r   r6   rG   r  )r(   r  r   r   r   test_dirichlet_bad_alpha  s    "z'TestRandomDist.test_dirichlet_bad_alphac                 C   sh   t g d}|d d d }tt| j}|j|dd}tt| j}|jt |dd}t|| d S )N)r  g      r  rA   r-  r   )r,   ra   r   r   r   r  ro   r
   )r(   rV  r  rG   rp   rq   r   r   r   #test_dirichlet_alpha_non_contiguous  s    z2TestRandomDist.test_dirichlet_alpha_non_contiguousc                 C   sx   d}|t ddg }tt| j}|j|dd}t ddgddggddgddggddgddggg}t||dd d S )	Nrt         ?r   r-  r   rT   r   r?  )r,   ra   r   r   r   r  r
   )r(   Zepsr  rG   r  r   r   r   r   test_dirichlet_small_alpha  s     z)TestRandomDist.test_dirichlet_small_alphac                 C   sT   t g d}||  }tt| j}|j|dd}|jdd}t||dd d S )N)g{Gz?{Gz?Q?i -1r   r   r   r   )r   )	r,   ra   rn   r   r   r   r  r   r   )r(   r  Z
exact_meanrG   r   Zsample_meanr   r   r   %test_dirichlet_moderately_small_alpha  s    z4TestRandomDist.test_dirichlet_moderately_small_alphar  )rW   r   r   r   )rE   r   r   r   )r/   rW   r   r         ?r   r   r   )g{Gz?r  r   r   )gh㈵>r   r   r   )gMb`?gQ?r   r   r  r   r   r   rT   r   c                 C   s*   t |}t|}t||dk d d S )Nr   rT   )r,   ra   rG   r  r   )r(   r  r  r   r   r   &test_dirichlet_multiple_zeros_in_alpha  s    

z5TestRandomDist.test_dirichlet_multiple_zeros_in_alphac                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Ng8EGr?r-  r   gvM?g?g#c?gQ:N?gAP?gE_9@r   r?  )r   r   r   exponentialr,   ra   r
   r4  r   r   r   test_exponential  s    
zTestRandomDist.test_exponentialc                 C   s&   t tjddd tttjdd d S Nr   scale       )r   rG   r  r   r6   r7   r   r   r   test_exponential_0  s    z!TestRandomDist.test_exponential_0c                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )Nr   M   r-  r   g}(Ҍ?gi?g=r?g]hDA?gjӨ?g8k*u?r   r?  )r   r   r   fr,   ra   r
   r4  r   r   r   test_f  s    
zTestRandomDist.test_fc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )NrW   rk   r-  r   g,שn'@gt@g<e2@g1a3@g3r
y-2@gs1-2@r   r?  )r   r   r   gammar,   ra   r
   r4  r   r   r   
test_gamma  s    
zTestRandomDist.test_gammac                 C   s*   t tjdddd tttjddd d S )Nr   )r`   r  r  )r   rG   r  r   r6   r7   r   r   r   test_gamma_0
  s    zTestRandomDist.test_gamma_0c                 C   sF   t t| j}|jddd}tddgddgddgg}t|| d S )N_c97ݚ?r-  r   r/   r  r   r   )r   r   r   	geometricr,   ra   r	   r4  r   r   r   test_geometric  s    
zTestRandomDist.test_geometricc                 C   s   t ttjd t ttjdgd  t ttjd t ttjdgd  tjdd6 t ttjtj t ttjtjgd  W d    n1 s0    Y  d S )Nre   r   rf   ignoreinvalid)r   r6   rG   r  r,   errstaterK   r7   r   r   r   test_geometric_exceptions  s    z(TestRandomDist.test_geometric_exceptionsc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )Nr         @r-  locr  r   gA7G@g1Fiҿgei@g_Cg0cPg`(_տr   r?  )r   r   r   gumbelr,   ra   r
   r4  r   r   r   test_gumbel  s    
zTestRandomDist.test_gumbelc                 C   s&   t tjddd tttjdd d S r  )r   rG   r  r   r6   r7   r   r   r   test_gumbel_0'  s    zTestRandomDist.test_gumbel_0c                 C   s   t t| j}|jddddd}tddgddgddgg}t|| |jdd	d
dd}tg d}t|| |jdd	ddd}tg d}t|| |jd	dd
dd}tg d}t|| |jd	dddd}tg d}t|| d S )Ng333333$@g      @r   r-  r   r   r   rW   r   rk   r~   )rk   rk   rk   rk   r   r   )r   r   r   r   )r   r   r   r   )r   r   r   hypergeometricr,   ra   r	   r4  r   r   r   test_hypergeometric+  s&    




z"TestRandomDist.test_hypergeometricc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )Nr  r  r-  r  g6@	ga@bDA"?g*[r{g>K?g(?g<?r   r?  )r   r   r   laplacer,   ra   r
   r4  r   r   r   test_laplaceE  s    
zTestRandomDist.test_laplacec                 C   s&   t tjddd tttjdd d S r  )r   rG   r  r   r6   r7   r   r   r   test_laplace_0M  s    zTestRandomDist.test_laplace_0c                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )Nr  r  r-  r  gbZgHʊg$>?g>&gQh@g%?n?gĉ8|?r   r?  )r   r   r   logisticr,   ra   r
   r4  r   r   r   test_logisticQ  s    
zTestRandomDist.test_logisticc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )Nr  r  r-  )r   sigmar   gܒx?gk*+@g<Ҿ?g
@g5a@g%jz(@r   r?  )r   r   r   	lognormalr,   ra   r
   r4  r   r   r   test_lognormalY  s    
zTestRandomDist.test_lognormalc                 C   s&   t tjddd tttjdd d S )Nr   )r  r/   r  )r   rG   r  r   r6   r7   r   r   r   test_lognormal_0a  s    zTestRandomDist.test_lognormal_0c                 C   sF   t t| j}|jddd}tddgddgdd	gg}t|| d S )
Ngƀ@?r-  )rI   r   r   r   rk      rW   r/   )r   r   r   	logseriesr,   ra   r	   r4  r   r   r   test_logseriese  s    
zTestRandomDist.test_logseriesc                 C   s$   t t| j}|ddks J d S )Nr   r/   )r   r   r   r  r(   rG   r   r   r   test_logseries_zerom  s    z"TestRandomDist.test_logseries_zerovaluer3   r  g      @c              	   C   s   t t| j}tjdd tt || W d    n1 sF0    Y  tt& |t	|gd  W d    n1 s0    Y  tt0 |t	|gd d d d  W d    n1 s0    Y  W d    n1 s0    Y  d S )Nr  r  r   rA   )
r   r   r   r,   r  rx   ry   r6   r  ra   )r(   r  rG   r   r   r   test_logseries_exceptionsq  s    (4z(TestRandomDist.test_logseries_exceptionsc                 C   sf   t t| j}|jddgd dd}tg dg dgg dg d	gg d
g dgg}t|| d S )Nr   UUUUUU?r   r-  r   )r/   rW   r/   r   r~   rk   )r~   rA   r   rA   r~   rA   )rW   rk   rA   r   rk   r/   )r~   r~   r   rA   rk   r   )r   rk   r/   rW   rk   rA   )rW   rW   rk   r/   rA   r~   )r   r   r   rR   r,   ra   r	   r4  r   r   r   test_multinomial~  s    
zTestRandomDist.test_multinomialzfp errors don't work in wasmreasonr   )svdeighcholeskyc                 C   s  t t| j}d}ddgddgg}d}|j||||d}tddgdd	ggd
dgddggddgddggg}t||dd |j|||d}tddg}t||dd ddg}ddgddgg}tt|j||dd ddgddgg}t	t
|j|| t	t
|j||dd tt|j||dd t|j||dd tt|j||dd tt|j||ddd ddgddgg}|dv r|j||d|d}t|d |d  d!d ntt|j||dd tjdd"gd"dggtjd#}t <}	|j|||d |	t
}
t|
dksJ W d    n1 s0    Y  td}td}tt|j||d$d tt|jtd%| tt|j|td tt|j|td& d S )'N)r  r   r/   r   r-  r   gtg$&@g&بO=g!$@gDo?g\&@g,R{?g
f!@g?ժ?gT^o#@g`R?gP'@r   r?  gl?g$("@rA   raise)check_validr  r  r  )r  r   )r  r  r   ).r   ).r/   r   rL  rC   other)rA   r/   r/   rk   )r   r   r   multivariate_normalr,   ra   r
   r   r6   r   RuntimeWarningr   r   rw   r   recordr  rF   Zeyer   )r(   r   rG   r   covr   r  r	  r%  supwmur   r   r   test_multivariate_normal  s    




2




z'TestRandomDist.test_multivariate_normalz	mean, covy      ?      ?y                c                 C   sL   t t| j}tjtdd ||| W d    n1 s>0    Y  d S )Nzmust not be complexru   )r   r   r   rx   ry   r5   r  )r(   r   r  rG   r   r   r   )test_multivariate_normal_disallow_complex  s    z8TestRandomDist.test_multivariate_normal_disallow_complexc           	      C   s   t t| j}d}tddg}tddgddgg}|j|||f|d}|| }|j| |d  }tt|	ddk sJ tt|| dk sJ d S )Nr"   r/   rA   r   r9   rL  rO   )
r   r   r   r,   ra   r  rz  rX   absr   )	r(   r   rG   Zn_sr   r  r)   Zs_centerZcov_empr   r   r   $test_multivariate_normal_basic_stats  s    z3TestRandomDist.test_multivariate_normal_basic_statsc                 C   sH   t t| j}|jdddd}tddgddgd	d
gg}t|| d S )NrN   g|Pk?r-  )nrI   r   i  i  i  i  iX  i  )r   r   r   r   r,   ra   r	   r4  r   r   r   r     s    
z%TestRandomDist.test_negative_binomialc                 C   sZ   t jdd: tttjdt j tttjdt jgd  W d    n1 sL0    Y  d S )Nr  r  rN   r   )r,   r  r   r6   rG   r   rK   r7   r   r   r   !test_negative_binomial_exceptions  s
    

z0TestRandomDist.test_negative_binomial_exceptionsc                 C   s8   t t tdd}W d    n1 s*0    Y  d S Nr/   r   )r   r6   rG   r   rY   r   r   r   #test_negative_binomial_p0_exception  s    
z2TestRandomDist.test_negative_binomial_p0_exceptionc                 C   sT   t jdd4 tttjdd tttjdgdg W d    n1 sF0    Y  d S )Nr  r  r  rL  )r,   r  r   r6   rG   r   r7   r   r   r   .test_negative_binomial_invalid_p_n_combination  s    z=TestRandomDist.test_negative_binomial_invalid_p_n_combinationc                 C   s   t t| j}|jdddd}tddgddgdd	gg}t||d
d |jdddd}tddgddgddgg}t||d
d t t| j}|jdddd}tddgddgddgg}t||d
d d S )NrW   r-  )dfnoncr   gs؀3J?gugو/@gNJm+@gA}-4@gq=&@gjhn@r   r?  rE   rO   gZօDN?gE!X&?gt%T?gQ)V?g5Zlf?gΎ{?r   gmK?gZ?g%yKY@g@W@gĈ{>@g&BB=@)r   r   r   noncentral_chisquarer,   ra   r
   r4  r   r   r   test_noncentral_chisquare  s(    


z(TestRandomDist.test_noncentral_chisquarec                 C   sN   t t| j}|jddddd}tddgdd	gd
dgg}t||dd d S )NrW   rA   r/   r-  )dfnumdfdenr  r   gHP
?gn?g%f?gch?gzƎ?g]?:?r   r?  )r   r   r   noncentral_fr,   ra   r
   r4  r   r   r   test_noncentral_f  s    

z TestRandomDist.test_noncentral_fc                 C   s2   t t| j}|jddtjd}t|s.J d S )NrW   rA   )r  r  r  )r   r   r   r  r,   rK   isnanr(   rG   r  r   r   r   test_noncentral_f_nan  s    z$TestRandomDist.test_noncentral_f_nanc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )Nr  r  r-  r  g?`xg\P@g/7u g0aj?gi%?gog@r   r?  )r   r   r   normalr,   ra   r
   r4  r   r   r   test_normal$  s    
zTestRandomDist.test_normalc                 C   s&   t tjddd tttjdd d S r  )r   rG   r  r   r6   r7   r   r   r   test_normal_0,  s    zTestRandomDist.test_normal_0c                 C   sN   t t| j}|jddd}tddgddgdd	gg}tjj||d
d d S )Nr  r-  rV  r   gá?gƫh@gi`>?g0"&Agc@g=fAr   )Znulp)r   r   r   paretor,   ra   ZtestingZassert_array_almost_equal_nulpr4  r   r   r   test_pareto0  s    
	zTestRandomDist.test_paretoc                 C   sF   t t| j}|jddd}tddgddgddgg}t|| d S )Nr  r-  )lamr   r   )r   r   r   poissonr,   ra   r	   r4  r   r   r   test_poisson>  s    
zTestRandomDist.test_poissonc                 C   s   t dj}d}tttj| tttj|gd  tttj| tttj|gd  t jdd6 tttjt j tttjt jgd  W d    n1 s0    Y  d S )Nr   r3   r   r  r  )	r,   r   r   r   r6   rG   r   r  rK   )r(   ZlambigZlamnegr   r   r   test_poisson_exceptionsF  s    z&TestRandomDist.test_poisson_exceptionsc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nr  r-  r  gZ >g?gÄ(=g҇?gŚE?gߢH?r   r?  )r   r   r   powerr,   ra   r
   r4  r   r   r   
test_powerQ  s    
zTestRandomDist.test_powerc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nr   r-  )r  r   gry@gtQP0@gAI_@g@u81@g	G0@g"wc5@r   r?  )r   r   r   rayleighr,   ra   r
   r4  r   r   r   test_rayleighY  s    
zTestRandomDist.test_rayleighc                 C   s&   t tjddd tttjdd d S r  )r   rG   r%  r   r6   r7   r   r   r   test_rayleigh_0a  s    zTestRandomDist.test_rayleigh_0c                 C   sH   t t| j}|jdd}tddgddgddgg}t||d	d
 d S )Nr-  r   gH1gW|3
g\?gLeg5P	gFA?r   r?  )r   r   r   Zstandard_cauchyr,   ra   r
   r4  r   r   r   test_standard_cauchye  s    
z#TestRandomDist.test_standard_cauchyc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nr-  invr   ga\5?gr7k?g8PC*?g5$[?gW3T?g*?r   r?  )r   r   r   standard_exponentialr,   ra   r
   r4  r   r   r   test_standard_exponentialm  s    
z(TestRandomDist.test_standard_exponentialc                 C   s   t ttjtjd d S NrC   )r   r5   rG   r*  r,   r   r7   r   r   r   #test_standard_expoential_type_erroru  s    z2TestRandomDist.test_standard_expoential_type_errorc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nrk   r-  r`   r   gΈ̏&?gTϿ?g1n21@gDb<z@gBv@gc^@r   r?  r   r   r   standard_gammar,   ra   r
   r4  r   r   r   test_standard_gammax  s    
z"TestRandomDist.test_standard_gammac                 C   s4   t t| j}|jdtjd}d}t||dd d S )Nrk   rC   g   d@r   r?  )r   r   r   r0  r,   rw   r
   r4  r   r   r   !test_standard_gammma_scalar_float  s    z0TestRandomDist.test_standard_gammma_scalar_floatc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nrk   r-  r.  gf&?gMO?g1@g):@g@g;i@rW   r?  r/  r4  r   r   r   test_standard_gamma_float  s    
z(TestRandomDist.test_standard_gamma_floatc                 C   s   t jdt jd}tt| j}|jd|t jd t jddgddgd	d
ggt jd}t||dd tt| j}|jd|dt jd t||dd d S )Nr-  rC         $@r  rD   gM»L$@ga {@gN]"@gPS")@g_+@gmB@rW   r?  r  r   rD   )	r,   rF   rw   r   r   r   r0  ra   r
   )r(   r  rG   r	  r   r   r   test_standard_gammma_float_out  s    
z-TestRandomDist.test_standard_gammma_float_outc                 C   s   t ttjddd d S )Nr  r   rC   )r   r5   rG   r0  r7   r   r   r    test_standard_gamma_unknown_type  s    
z/TestRandomDist.test_standard_gamma_unknown_typec                 C   s6   t d}tttjdd|d tttjdd|d d S )Nr   r4  r   )r   r  )r   r/   )r,   rF   r   r6   rG   r0  )r(   r  r   r   r   test_out_size_mismatch  s    
z%TestRandomDist.test_out_size_mismatchc                 C   s&   t tjddd tttjdd d S )Nr   )r`   r  )r   rG   r0  r   r6   r7   r   r   r   test_standard_gamma_0  s    z$TestRandomDist.test_standard_gamma_0c                 C   sH   t t| j}|jdd}tddgddgddgg}t||d	d
 d S )Nr-  r   g
bYg? ?gGH{Lgz\?gc`+?gT}?r   r?  )r   r   r   r   r,   ra   r
   r4  r   r   r   test_standard_normal  s    
z#TestRandomDist.test_standard_normalc                 C   s   t ttjtjd d S r,  )r   r5   rG   r   r,   r   r7   r   r   r   %test_standard_normal_unsupported_type  s    z4TestRandomDist.test_standard_normal_unsupported_typec                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )Nr   r-  )r  r   gV1g	|ڙ(?g-?g NhڿgMc?g!uLM r   r?  )r   r   r   
standard_tr,   ra   r
   r4  r   r   r   test_standard_t  s    
zTestRandomDist.test_standard_tc                 C   sN   t t| j}|jddddd}tddgdd	gd
dgg}t||dd d S )Ng{Gz@g(\u$@gףp=
W4@r-  )leftmoderightr   gEpw@g4DC+@gŋ@g/,@gkuR*@gHNPcr+@r   r?  )r   r   r   
triangularr,   ra   r
   r4  r   r   r   test_triangular  s    

zTestRandomDist.test_triangularc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )NGz?gGz%@r-  )r  r  r   gm2@gJaD@g: @g
Q @g@T{e@gb@r   r?  r   r   r   uniformr,   ra   r
   r4  r   r   r   test_uniform  s    
zTestRandomDist.test_uniformc                 C   s   t dj}t dj}tj}tt|t j d tt|dt j tt||| tt|t j gdg tt|dgt jg tjt 	|d|d d d S )Nrb   r   r/   g ؅W4vC)r  r  )
r,   Zfinfor   r   rG   rF  r   OverflowErrorinf	nextafter)r(   ZfminZfmaxfuncr   r   r   test_uniform_range_bounds  s    z(TestRandomDist.test_uniform_range_boundsc                 C   sn   t j}|dd}t|d |dtjgdtjg}t|dtjg |dgdggddg}t|dtd  d S )Nr  rT   g
ףp= @r^   )rG   rF  r   r,   pirF   )r(   rK  resultr   r   r   test_uniform_zero_range  s    

z&TestRandomDist.test_uniform_zero_rangec                 C   sH   t j}tt|dd tt|ddgddg tt|ddgddggd d S )NrA   r/   r   rk   )rG   rF  r   r6   )r(   rK  r   r   r   test_uniform_neg_range  s    z%TestRandomDist.test_uniform_neg_rangec                 C   sj   G dd dt j}t d|}tttj|| G dd dt j}t d|}tttj|dd d S )Nc                   @   s   e Zd Zdd ZdS )zGTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingFloatc                 S   s   t d S r   r5   r7   r   r   r   	__float__  s    zQTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingFloat.__float__N)r=   r>   r?   rR  r   r   r   r   ThrowingFloat  s   rS  r  c                   @   s   e Zd Zdd ZdS )zITestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingIntegerc                 S   s   t d S r   rQ  r7   r   r   r   __int__  s    zQTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingInteger.__int__N)r=   r>   r?   rT  r   r   r   r   ThrowingInteger  s   rU  r/   )	r,   Zndarrayra   r   r   r5   rG   rF  r  )r(   rS  Zthrowing_floatrU  Zthrowing_intr   r   r   !test_scalar_exception_propagation  s    
z0TestRandomDist.test_scalar_exception_propagationc                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )NrD  p=
ף?r-  r  kappar   gLA?gxyw@gzVR?gpE ?g8&Aпgv @r   r?  )r   r   r   vonmisesr,   ra   r
   r4  r   r   r   test_vonmises  s    
zTestRandomDist.test_vonmisesc                 C   s4   t t| j}|jdddd}tt|  d S )NrT   gg5_PG>i@B rX  )r   r   r   rZ  r   r,   isfiniterX   r(   rG   rr   r   r   test_vonmises_small  s    z"TestRandomDist.test_vonmises_smallc                 C   s0   t t| j}|jdtjd}tt| d S )NrT   )r  rY  )r   r   r   rZ  r,   rK   r   r  r]  r   r   r   test_vonmises_nan  s    z TestRandomDist.test_vonmises_nanrY  g     @  4&kCc                 C   sr   t t| j}t|j}|jj}|jd|dd}||j_|jd|dd}|dk r\t|| nt	||ksnJ d S )Nr   r   r   g    .A)
r   r   r   r   r   r   rZ  r   r,   rX   )r(   rY  rG   rsr   Zrandom_state_valsZgen_valsr   r   r   test_vonmises_large_kappa  s    
z(TestRandomDist.test_vonmises_large_kappar  g      gg	@)rt   gư>r/   g     @@ra  c                 C   sF   t t| j}|||d}tt|tj ko>t|tjk d S )Nr   )r   r   r   rZ  r   r,   rX   rM  )r(   r  rY  rG   r^  r   r   r   test_vonmises_large_kappa_range$  s    z.TestRandomDist.test_vonmises_large_kappa_rangec                 C   sL   t t| j}|jdddd}tddgddgd	d
gg}t||dd d S )NrD  rW  r-  )r   r  r   g*س2?gt	@g;S@g\ F@gxF @g6b?r   r?  )r   r   r   waldr,   ra   r
   r4  r   r   r   	test_wald+  s    
zTestRandomDist.test_waldc                 C   sJ   t t| j}|jddd}tddgddgdd	gg}t||d
d d S )NrD  r-  r  gnJZ?g*/F?gs(RY?g]'%?g,qD?g8 ?r   r?  )r   r   r   weibullr,   ra   r
   r4  r   r   r   test_weibull3  s    
zTestRandomDist.test_weibullc                 C   s<   t t| j}t|jdddtd tt|jdd d S )Nr   r   r  r  )rV  )	r   r   r   r   rg  r,   rF   r   r6   r  r   r   r   test_weibull_0;  s    zTestRandomDist.test_weibull_0c                 C   sF   t t| j}|jddd}tddgddgddgg}t|| d S )	NrD  r-  r  r/   r   ic  ib  rA   )r   r   r   zipfr,   ra   r	   r4  r   r   r   	test_zipf@  s    
zTestRandomDist.test_zipf)r=   r>   r?   r}   r5  r9  r;  r>  rB  rD  rE  rx   r   r   r,   rw   r_   float64r   rF  rG  rJ  rN  rR  rU  r[  r\  r]  r`  re  rg  ri  rj  rl  rn  r  r  r  r  r  r  r  r  r  r  rB   r  ra   r  r  r  r  r  r  r  r  r  r  r  r  r  r  r)  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  rJ  rK   r  r  skipifr   r  r  r  r   r
  r  r  r  r  r  r  r  r  r!  r"  r$  r&  r'  r(  r+  r-  r1  r2  r3  r7  r8  r9  r:  r;  r<  r>  rC  rG  rL  rO  rP  rV  r[  r_  r`  rc  rM  rd  rf  rh  ri  rk  r   r   r   r   r*    s6  


		.			
	


	
I$

			$r*  c                	   @   sT  e Zd Zdd Zdd Zdd Zdd Zd	d
 Zdd Zdd Z	dd Z
dd Zdd Zdd Zdd Zdd Zdd Zdd Zdd  Zd!d" Zd#d$ Zd%d& Zd'd( Zd)d* Zd+d, Zd-d. Zd/d0 Zd1d2 Zd3d4 Zd5d6 Zd7d8 Zd9d: Zd;d< Z d=d> Z!d?d@ Z"e#j$%dAdBe&'dBdBge&'dBggdBggggdCdD Z(dEdF Z)dGdH Z*dIS )JTestBroadcastc                 C   s
   d| _ d S )Ni[r|   r7   r   r   r   r}   L  s    zTestBroadcast.setup_methodc                 C   s   t t| j}dg}dg}|j}tg d}t t| j}||d |}t||dd t t| j}|||d }t||dd d S )Nr   r/   )g4
7^?gF]
"?g .?rk   r   r?  rE  )r(   rG   r  r  rF  r	  r  r   r   r   rG  O  s    zTestBroadcast.test_uniformc                 C   s   dg}dg}dg}t t| j}tg d}t t| j}||d |}t||dd tt|j|d | t t| j}|j}|||d }t||dd tt|||d  d S )Nr   r/   r3   )gw-Dؿgn0_x?gʊ;<.?rk   r   r?  )	r   r   r   r,   ra   r  r
   r   r6   )r(   r  r  	bad_scalerG   r	  r  r  r   r   r   r  ^  s    zTestBroadcast.test_normalc           	      C   s   dg}dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt||d | tt| j}|||d }t||dd d S )	Nr/   rA   r3   r9   )gf8?g%eo?g$q??rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )	r(   rV  rW  bad_aZbad_br	  rG   r  r  r   r   r   r  p  s    zTestBroadcast.test_betac                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S Nr/   r3   )g`ʄ?gE]?g?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   r(   r  ro  r	  rG   r  r   r   r   r    s    zTestBroadcast.test_exponentialc                 C   s\   dg}dg}t g d}tt| j}|j}||d }t||dd tt||d  d S rq  )	r,   ra   r   r   r   r0  r
   r   r6   )r(   r`   	bad_shaper	  rG   Z	std_gammar  r   r   r   r1    s    z!TestBroadcast.test_standard_gammac           	      C   s   dg}dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt||d | tt| j}|j}|||d }t||dd tt|||d  tt|||d  d S )	Nr/   rA   r3   r9   )g3ʄ?gJE]?g?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )	r(   r`   r  rs  ro  r	  rG   r  r  r   r   r   r    s"    zTestBroadcast.test_gammac           	      C   s   dg}dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt||d | tt| j}|j}|||d }t||dd tt|||d  tt|||d  d S )	Nr/   rA   r3   r9   )gM3B?gybi@g0?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )	r(   r  r  	bad_dfnum	bad_dfdenr	  rG   r  r  r   r   r   r    s"    zTestBroadcast.test_fc              	   C   s  dg}dg}dg}dg}dg}dg}t g d}tt| j}|j}	|	|d ||}
t|
|dd	 t t |	||t j	gd sJ t
t|	|d || t
t|	|d || t
t|	|d || tt| j}|j}	|	||d |}
t|
|dd	 t
t|	||d | t
t|	||d | t
t|	||d | tt| j}|j}	|	|||d }
t|
|dd	 t
t|	|||d  t
t|	|||d  t
t|	|||d  d S )
NrA   rk   r~   r   r3   r9   )gŜ,n1 @g!6)@g9?r   r?  )r,   ra   r   r   r   r  r
   rX   r  rK   r   r6   )r(   r  r  r  rt  ru  bad_noncr	  rG   Znonc_fr  r   r   r   r    s:    $zTestBroadcast.test_noncentral_fc                 C   s@   t t| j}tddg}|jddddd}t||dd d S )NgZ\#?gDqw?rg   rA   r   r   r?  )r   r   r   r,   ra   r  r
   )r(   rG   r	  r  r   r   r   test_noncentral_f_small_df  s    z(TestBroadcast.test_noncentral_f_small_dfc                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S )Nr/   r3   )gC{?g7Q#?g%h1T@rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   r(   r  bad_dfr	  rG   r  r   r   r   r    s    zTestBroadcast.test_chisquarec           	      C   s   dg}dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt||d | tt| j}|j}|||d }t||dd tt|||d  tt|||d  d S )	Nr/   rA   r3   r9   )g-S?g[@gr3/.?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )	r(   r  r  ry  rv  r	  rG   Znonc_chir  r   r   r   r    s"    z'TestBroadcast.test_noncentral_chisquarec                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S )Nr/   r3   )g>Qg{g,e?rk   r   r?  )	r,   ra   r   r   r   r=  r
   r   r6   rx  r   r   r   r>  	  s    zTestBroadcast.test_standard_tc                 C   s   dg}dg}dg}t g d}tt| j}||d |}t||dd tt|j|d | tt| j}|||d }t||dd tt|j||d  d S )NrA   r/   r3   )gD()@g,)ʸ@gPork   r   r?  )	r,   ra   r   r   r   rZ  r
   r   r6   )r(   r  rY  Z	bad_kappar	  rG   r  r   r   r   r[    s    zTestBroadcast.test_vonmisesc                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S )Nr/   r3   )g/h?gcE?g/gV?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   r(   rV  rp  r	  rG   r  r   r   r   r  #  s    zTestBroadcast.test_paretoc                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S rq  )	r,   ra   r   r   r   rg  r
   r   r6   rz  r   r   r   rh  -  s    zTestBroadcast.test_weibullc                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S )Nr/   r3   )gT?gʥ?gKܮc?rk   r   r?  )	r,   ra   r   r   r   r#  r
   r   r6   rz  r   r   r   r$  7  s    zTestBroadcast.test_powerc                 C   s   dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt| j}|j}|||d }t||dd tt|||d  d S )Nr   r/   r3   )gh&BgsgR#?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )r(   r  r  ro  r	  rG   r  r  r   r   r   r  A  s    zTestBroadcast.test_laplacec                 C   s   dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt| j}|j}|||d }t||dd tt|||d  d S )Nr   r/   r3   )gs4?g$Cw%T?gտrk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )r(   r  r  ro  r	  rG   r  r  r   r   r   r  S  s    zTestBroadcast.test_gumbelc                 C   s   dg}dg}dg}t g d}tt| j}||d |}t||dd tt|j|d | tt| j}|||d }t||dd tt|j||d  t	|dd	d d S )
Nr   r/   r3   )g$DgXPg,?rk   r   r?  r  rT   )
r,   ra   r   r   r   r  r
   r   r6   r   )r(   r  r  ro  r	  rG   r  r   r   r   r  e  s    zTestBroadcast.test_logisticc                 C   s   dg}dg}dg}t g d}tt| j}|j}||d |}t||dd tt||d | tt| j}|||d }tt|j||d  d S )Nr   r/   r3   )g c?gu@gVQ?rk   r   r?  )	r,   ra   r   r   r   r  r
   r   r6   )r(   r   r  Z	bad_sigmar	  rG   r  r  r   r   r   r  v  s    zTestBroadcast.test_lognormalc                 C   sZ   dg}dg}t g d}tt| j}||d }t||dd tt|j|d  d S )Nr/   r3   )gs(?g4mK?gVk+?rk   r   r?  )	r,   ra   r   r   r   r%  r
   r   r6   rr  r   r   r   r&    s    zTestBroadcast.test_rayleighc                 C   s   dg}dg}dg}dg}t g d}tt| j}||d |}t||dd tt|j|d | tt|j|d | tt| j}|||d }t||dd tt|j||d  tt|j||d  d S )	NrE   r/   r   r9   )gk7Z?g*}z9?gc-	?rk   r   r?  )	r,   ra   r   r   r   re  r
   r   r6   )r(   r   r  Zbad_meanro  r	  rG   r  r   r   r   rf    s    zTestBroadcast.test_waldc                 C   s  dg}dg}dg}dg}dg}|d \}}t g d}tt| j}	|	j}
|
|d ||}t||dd tt|
|d || tt|
|d || tt|
|d || tt| j}	|	j}
|
||d |}t||dd tt|
||d | tt|
||d | tt|
||d | tt| j}	|	j}
|
|||d }t||dd tt|
|||d  tt|
|||d  tt|
|||d  tt|
dd	d
 tt|
ddd
 tt|
ddd d S )Nr/   rk   rA   r~   )gR8T>?gK?g;>h@r   r?  r4  rT   g      4@g      9@)	r,   ra   r   r   r   rB  r
   r   r6   )r(   r?  rA  r@  Zbad_left_oneZbad_mode_oneZbad_left_twoZbad_mode_twor	  rG   rB  r  r   r   r   rC    sJ    
zTestBroadcast.test_triangularc           
      C   s   dg}dg}dg}dg}dg}t g d}tt| j}|j}||d |}	t|	| tt||d | tt||d | tt||d | tt| j}|||d }	t|	| tt|||d  tt|||d  tt|||d  d S )Nr/   rE   r3   r  r   r   r/   rk   )	r,   ra   r   r   r   rH   r	   r   r6   )
r(   r	  rI   bad_n	bad_p_one	bad_p_twor	  rG   Zbinomr  r   r   r   r    s&    

zTestBroadcast.test_binomialc           
      C   s   dg}dg}dg}dg}dg}t jg dt jd}tt| j}|j}||d |}	t|	| tt	||d | tt	||d | tt	||d | tt| j}|j}|||d }	t|	| tt	|||d  tt	|||d  tt	|||d  d S )Nr/   rE   r3   r  )r   rA   r/   rC   rk   )
r,   ra   r   r   r   r   r   r	   r   r6   )
r(   r	  rI   r|  r}  r~  r	  rG   Z	neg_binomr  r   r   r   r     s(    

z$TestBroadcast.test_negative_binomialc           	      C   sx   dg}dg}t g d}tt| j}|j}|d g}|j}||d }t|| tt	||d  tt	||d  d S )Nr/   r3   )r   r   rk   rA   rk   )
r,   ra   r   r   r   Z_poisson_lam_maxr   r	   r   r6   )	r(   r  Zbad_lam_oner	  rG   Zmax_lamZbad_lam_twor   r  r   r   r   r!    s    

zTestBroadcast.test_poissonc                 C   s   dg}dg}t g d}tt| j}|j}||d }t|| tt||d  t j	dd2 tt|t j
 tt|ddt j
g W d    n1 s0    Y  d S )NrA   r   )r/   r   r/   rk   r  r  )r,   ra   r   r   r   rj  r	   r   r6   r  rK   )r(   rV  rp  r	  rG   rj  r  r   r   r   rk  	  s    
zTestBroadcast.test_zipfc                 C   sn   dg}dg}dg}t g d}tt| j}|j}||d }t|| tt||d  tt||d  d S )NrE   r3   r  r   rk   )	r,   ra   r   r   r   r  r	   r   r6   )r(   rI   r}  r~  r	  rG   r  r  r   r   r   r  	  s    
zTestBroadcast.test_geometricc                 C   s2  dg}dg}dg}dg}dg}dg}dg}t g d}tt| j}	|	|d ||}
t|
| tt|	j|d || tt|	j|d || tt|	j|d || tt|	j|d || tt| j}	|	||d |}
t|
| tt|	j||d | tt|	j||d | tt|	j||d | tt|	j||d | tt| j}	|	j}||||d }
t|
| tt||||d  tt||||d  tt||||d  tt||||d  tt|ddd	 tt|ddd	 tt|ddd tt|ddd
 tt|ddd	 tt|ddd tt|dddgd d S )Nr/   rA   r3   r9   r~   r{  rk   r   r   r   i   @i  l        r   i    r"   )	r,   ra   r   r   r   r  r	   r   r6   )r(   ZngoodZnbadr   Z	bad_ngoodZbad_nbadZbad_nsample_oneZbad_nsample_twor	  rG   r  Z	hypergeomr   r   r   r  (	  sJ    


z!TestBroadcast.test_hypergeometricc                 C   sn   dg}dg}dg}t g d}tt| j}|j}||d }t|| tt||d  tt||d  d S )NrE   rA   r3   )r/   r/   r/   rk   )	r,   ra   r   r   r   r  r	   r   r6   )r(   rI   r}  r~  r	  rG   r  r  r   r   r   r  U	  s    
zTestBroadcast.test_logseriesc                 C   sn  t t| j}|jddgdgd dd}tjg dg dgg d	g d
gg dg dggtjd}t|| t t| j}|ddgdgd }tjg dg dgtjd}t|| t t| j}|ddgdgd gd }tjg dg dgtjd}t|| t t| j}|dgdggdgd gd }tjg dg dgg dg d
ggtjd}t|| d S )NrW   r   r  r   r-  r   )r   r   rA   r/   rA   r   )rA   rk   r   r~   rA   rk   )r/   r   r/   r   rA   r/   )r   rA   rA   r/   r~   r~   )r   rA   r   r/   rA   r   )rk   rA   rk   rk   r~   rW   rC   rA   )r   r   rA   r/   r/   r/   )r~   rA   rk   rk   rW   rk   )r   r   r   rR   r,   ra   r   r	   r4  r   r   r   r  b	  sV    





 
zTestBroadcast.test_multinomialr	  r   c                 C   sJ  t t| j}tdgd }|||}t|tr:t n|j	}|d }|j	|ksVJ t
||g}|||}t||j	d d d }|j	|ksJ t
|g|gg}|||}t||j	d d }|j	|d ksJ |j||d| d}|j	d| d ks
J tt  |j||dd}W d    n1 s<0    Y  d S )Nr_  r~   )r~   r3   r-  r   )r/   r/   r/   r/   r/   r/   )r   r   r   r,   ra   rR   r   rB   ra  r`   ry  Zbroadcast_shapesrx   ry   r6   )r(   r	  rG   rm   r  Zn_shaper   r   r   r   test_multinomial_pval_broadcast	  s$    z-TestBroadcast.test_multinomial_pval_broadcastc                 C   sF   t t| j}dgd dgd g}tt|jd| tt|jdd d S )Nr  r   r_  r/   rE   )r   r   r   r   r6   rR   )r(   rG   rm   r   r   r   test_invalid_pvals_broadcast	  s    z*TestBroadcast.test_invalid_pvals_broadcastc                 C   s   t t| j}|tdddgd }|jdks6J |dtd}|jdksVJ |tddtd	}|jd	ks~J d S )
N)r   r   r   i8r  r   )r   r   r   r   r   )r   r   r   )rk   r   r   )rk   r   r   r~   )r   r   r   rR   r,   r   r`   r  r   r   r   test_empty_outputs	  s    z TestBroadcast.test_empty_outputsN)+r=   r>   r?   r}   rG  r  r  r  r1  r  r  r  rw  r  r  r>  r[  r  rh  r$  r  r  r  r  r&  rf  rC  r  r   r!  rk  r  r  r  r  rx   r   r   r,   ra   r  r  r  r   r   r   r   rn  I  sN   
#




(-

rn  zcan't start threadr  c                   @   s4   e Zd Zdd Zdd Zdd Zdd Zd	d
 ZdS )
TestThreadc                 C   s   t d| _d S )Nr~   )r+   seedsr7   r   r   r   r}   	  s    zTestThread.setup_methodc                    s   ddl m  tt| jf| }tt| jf| } fddt| j|D }dd |D  dd |D  t| j|D ]\}}tt|| q|t	 j
jdkrtjdkrt|| n
t|| d S )	Nr   )Threadc                    s(   g | ] \}} t t||fd qS ))targetargs)r   r   )ru  r)   or  functionr   r   rw  	  s   z-TestThread.check_function.<locals>.<listcomp>c                 S   s   g | ]}|  qS r   )startru  rZ   r   r   r   rw  	  rq  c                 S   s   g | ]}|  qS r   )joinr  r   r   r   rw  	  rq  r~   win32)	threadingr  r,   r   r  r  zipr   r   ZintprD   r   r   platformr
   r	   )r(   r  szZout1Zout2tr)   r  r   r  r   check_function	  s    
zTestThread.check_functionc                 C   s   dd }| j |dd d S )Nc                 S   s   | j dd|d< d S )Nrk  r   .)r  r   r  r   r   r   
gen_random	  s    z*TestThread.test_normal.<locals>.gen_random)rk  r  r  r(   r  r   r   r   r  	  s    zTestThread.test_normalc                 C   s   dd }| j |dd d S )Nc                 S   s   | j tdd|d< d S )NrN   r"   r  .)r  r,   rh  r  r   r   r   r  	  s    z'TestThread.test_exp.<locals>.gen_randomr  r  r  r  r   r   r   test_exp	  s    zTestThread.test_expc                 C   s   dd }| j |dd d S )Nc                 S   s   | j ddgd dd|d< d S )Nr   r  r   rk  r   .)rR   r  r   r   r   r  	  s    z/TestThread.test_multinomial.<locals>.gen_random)rk  r   r  r  r  r   r   r   r  	  s    zTestThread.test_multinomialN)r=   r>   r?   r}   r  r  r  r  r   r   r   r   r  	  s
   r  c                   @   s4   e Zd Zdd Zdd Zdd Zdd Zd	d
 ZdS )TestSingleEltArrayInputc                 C   s4   t dg| _t dg| _t dg| _d| _d S )NrA   rk   r~   r{  )r,   ra   argOneargTwoargThreetgtShaper7   r   r   r   r}   	  s    z$TestSingleEltArrayInput.setup_methodc                 C   s   t jt jt jt jt jt jt jt jt j	t j
t jt jf}t jt jf}|D ]6}||v rb|tdg}n
|| j}t|j| j qDd S r   )rG   r  r0  r  r=  r  rg  r#  r%  r   rj  r  r  r,   ra   r  r   r`   r  )r(   funcs	probfuncsrK  r  r   r   r   test_one_arg_funcs	  s    
z*TestSingleEltArrayInput.test_one_arg_funcsc                 C   s   t jt jt jt jt jt jt jt jt j	t j
t jt jt jt jf}t jt jf}|D ]v}||v rftdg}n| j}|| j|}t|j| j || jd |}t|j| j || j|d }t|j| j qLd S )NrE   r   )rG   rF  r  r  r  r  r  rZ  r  r  r  r  re  rH   r   r,   ra   r  r  r   r`   r  )r(   r  r  rK  r  r  r   r   r   test_two_arg_funcs	  s&    z*TestSingleEltArrayInput.test_two_arg_funcsc              	   C   s   t jt jt jt jt jt jt jt jt j	g	}t
j}t dg}t dg}|D ]f}|||||d}t|j| j ||d |||d}t|j| j |||d ||d}t|j| j qJd S )Nr/   r   r   )r,   r  r   r   r   r   r   r_   r   r   rG   r&   ra   r   r`   r  )r(   r    r   rK  r  r  r   r  r   r   r   r5  
  s    z%TestSingleEltArrayInput.test_integersc                 C   s   t jt jt jg}|D ]l}|| j| j| j}t|j| j	 || jd | j| j}t|j| j	 || j| jd | j}t|j| j	 qd S )Nr   )
rG   r  rB  r  r  r  r  r   r`   r  )r(   r  rK  r  r   r   r   test_three_arg_funcs$
  s    z,TestSingleEltArrayInput.test_three_arg_funcsN)r=   r>   r?   r}   r  r  r5  r  r   r   r   r   r  	  s
   r  r  c                 C   s   | d }| d }t |}|| |jd d }tjdkrB| }t|}|jd d | d d ksjJ | | d d ksJ |	 }|jd d }tjdkr| }t|}|jd d | d	 d ksJ | | d	 d ksJ d S )
Nr   r   r   r   r   r   r   r   r   )
r   r0   r   r   r   r   r   r   r   r   )r  r   r   Zmt19937r   r   r   r   r   r   test_jumped3
  s"    




r  c                  C   sn  t d} t d}d}tj| |ddjdks2J tt  tj| ||d W d    n1 sb0    Y  tt  tj| |dd W d    n1 s0    Y  tt  tj| |dd W d    n1 s0    Y  t d}tt tj||d W d    n1 s0    Y  tt tj|dd W d    n1 sT0    Y  tt tj|dd W d    n1 s0    Y  t 	|}tt tj||d	 W d    n1 s0    Y  tt$ tj
d
ddgdd W d    n1 s0    Y  tt$ tj
d
dgddd W d    n1 sR0    Y  tt( tj
d
dgddgdd W d    n1 s0    Y  tt( tjddgddgdd W d    n1 s0    Y  tjddd}tjddd}tjddd}t|||jdks(J tjtdd" tj|||dd W d    n1 s`0    Y  d S )Nrk   )r~   rk   )r   r~   rA   )rW   r~   rk   r   )r/   rk   )r~   r/   r/   )rk   r  r/   rS  gffffff?)rA   r/   rA   rW   zOutput size \(6, 5, 1, 1\) isru   )r   rW   r/   r/   )r,   rh  rG   r  r`   rx   ry   r6   r0  r   rH   rR   r  r  )r  r  r   r`   r  rV  rW  rX  r   r   r   test_broadcast_size_errorO
  sD    

...
...
.4488r  c                  C   sb   t d} t d}tj| |dd tt  tj| |dd W d    n1 sT0    Y  d S )Nrk   r   rA   )r,   rh  rG   r  rx   ry   r6   )r  r  r   r   r   test_broadcast_size_scalarz
  s
    

r  c                  C   s8   g g dg} t td}t|j|  | dg g gks4J d S r  )r   r   r   rQ  )seqgenr   r   r   test_ragged_shuffle
  s    
r  r  r9   c                 C   s   t td}|rdnd}tjt|d |j| |d W d    n1 sJ0    Y  |r\dnd}tjt|d  |jd| |d W d    n1 s0    Y  tjt|d" |jdg| |d W d    n1 s0    Y  d S )	Nr   zhigh < 0z	high <= 0ru   r=  z
low > highzlow >= highr3   )r   r   rx   ry   r6   r&   )r  r    r  msgr   r   r   !test_single_arg_integer_exception
  s    ,.r  rD   Zf4Zf8c                 C   s   t jdd| d}g d}tjtdd  tj||| d W d    n1 sL0    Y  tjtdd$ tj|||j| d W d    n1 s0    Y  d S )	Nr   F)orderrD   r   zSupplied output arrayru   r5  r6  )r,   r   rx   ry   r6   rG   r0  r`   )rD   r  r`   r   r   r   test_c_contig_req_out
  s    .r  r  r  Cdistc                 C   sH   t jd||d}| ||d}||u s(J | |||jd}||u sDJ d S )Nr   )rD   r  r5  )r  rD   r   )r,   r   r`   )r  r  rD   r  Zvariatesr   r   r   test_contig_req_out
  s
    r  c                  C   sp   t jt jd} | d |  \}}}|d d dks@J ||d d  }||j_|jj}||kslJ d S )Nr   r/   )	PCG64DXSM)r,   rG   r   r  r   
__reduce__r   r   )r   Zctorr  Zstate_arW  Zstate_br   r   r   $test_generator_ctor_old_style_pickle
  s    
r  )1r   r   rx   numpyr,   Znumpy.linalgr   Znumpy.testingr   r   r   r   r   r   r	   r
   r   r   Znumpy.randomr   r   r   r   rG   ZJUMP_TEST_DATAZfixturer    r!   r@   rM   r{   r   r   r*  rn  r   rm  r  r  r   r  r  r  r  r  r  r   r  r  r   r   r   r   <module>   s   0

$7 "                i/U
+

	