o
    k)hM                 	   @   s  d dl Z d dlZd dlZd dl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 d dlmZmZ d dlmZ ddddded	d
 fddddZeejjdk rcdddddddddZndddddddd dZejd!ed"d#d$ Zejd%d& Zd'd( ZG d)d* d*ZG d+d, d,Z G d-d. d.Z!G d/d0 d0Z"G d1d2 d2Z#G d3d4 d4Z$G d5d6 d6Z%ej&j'ed7d8G d9d: d:Z(G d;d< d<Z)d=d> Z*d?d@ Z+dAdB Z,dCdD Z-dEdF Z.dGdH Z/dIdJ Z0dKdL Z1dMdN Z2dS )O    N)	assert_assert_raisesassert_equalassert_warnsassert_no_warningsassert_array_equalassert_array_almost_equalsuppress_warningsIS_WASM)MT19937PCG64)random)g      Y@g333333?)      ?)   r   
   r      g      @)d   r   )      $@)   )binomial	geometrichypergeometric	logseriesmultinomialnegative_binomialpoissonzipf        Z@2fbead005fc63942decb5326d36a1f32fe2c9d32c904ee61e46866b88447c263Z@23ead5dcde35d4cfd4ef2c105e4c3d43304b45dc1b1444b7823b9ee4fa144ebbZ@0d764db64f5c3bad48c8c33551c13b4d07a1e7b470f77629bef6c985cac76fcfZ@7b59bf2f1691626c5815cdcd9a49e1dd68697251d4521575219e4d2a1b8b2c67Z@d754fa5b92943a38ec07630de92362dd2e02c43577fc147417dc5b9db94ccdd3Z@8eb216f7cb2a63cf55605422845caaff002fddc64a7dc8b2d45acd477a49e824Z@70c891d76104013ebd6f6bcf30d403a9074b886ff62e4e6b8eb605bf1a4673b7Z@01f074f97517cd5d21747148ac6ca4074dde7fcb7acbaec0a936606fecacd93f)r   r   r   r   r   r   r   r   Z@8626dd9d052cb608e93d8868de0a7b347258b199493871a1dc56e2a26cacb112Z@8edd53d272e49c4fc8fbbe6c7d08d563d62e482921f3131d0a0e068af30f0db9Z@83496cc4281c77b786c9b7ad88b74d42e01603a55c60577ebab81c3ba8d45657Z@65878a38747c176bc00e930ebafebb69d4e1e16cd3a704e264ea8f5e24f548dbZ@7a984ae6dca26fd25374479e118b22f55db0aedccd5a0f2584ceada33db98605Z@d636d968e6a24ae92ab52fe11c46ac45b0897e98714426764e820a7d77602a61Z@956552176f77e7c9cb20d0118fc9cf690be488d790ed4b4c4747b965e61b0bb4Z@f84ba7feffda41e606e20b28dfc0f1ea9964a74574513d4a4cbc98433a8bfa45module)Zscopeparamsc                 C   s   | j t| j  t| j  fS N)Zparam	INT_FUNCSINT_FUNC_HASHES)Zrequest r#   Clib/python3.10/site-packages/numpy/random/tests/test_randomstate.pyint_func2   s   r%   c                  c   s"    t j } dV  t j|  dS )z:Ensures that the singleton bitgen is restored after a testN)npr   get_bit_generatorset_bit_generator)Zorig_bitgenr#   r#   r$   restore_singleton_bitgen8   s   
r)   c                 C   sn   t | d |d  t| d d |d d  t| d d |d d  t | d |d  t | d |d  d S )Nbit_generatorstatekeypos	has_gaussgauss)r   r   )abr#   r#   r$   assert_mt19937_state_equal@   s
   r2   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S )TestSeedc                 C   s8   t d}t|dd t d}t|dd d S )Nr     i      i  )r   RandomStater   randintselfsr#   r#   r$   test_scalarI   s   

zTestSeed.test_scalarc                 C   sz   t td}t|dd t td}t|dd t dg}t|dd t dg}t|dd d S )Nr   r4   i  r   i  r5   i	  )r   r6   ranger   r7   r&   aranger8   r#   r#   r$   
test_arrayO   s   zTestSeed.test_arrayc                 C   s    t ttjd t ttjd d S )N      r   	TypeErrorr   r6   
ValueErrorr9   r#   r#   r$   test_invalid_scalarY      zTestSeed.test_invalid_scalarc                 C   sX   t ttjdg t ttjdg t ttjdg t ttjg d t ttjg d d S )Nr?   r@   r   )   r   r   )rG   r   rA   rD   r#   r#   r$   test_invalid_array^   s
   zTestSeed.test_invalid_arrayc                 C   sL   t ttjtjg tjd t ttjg dg t ttjg dg dg d S )NdtyperG   r      )      r   )r   rC   r   r6   r&   arrayint64rD   r#   r#   r$   test_invalid_array_shapef   s   
z!TestSeed.test_invalid_array_shapec                 C   sF   t td}tt |d W d    d S 1 sw   Y  d S )Nr     )r   r6   r   r   rB   seed)r9   rsr#   r#   r$   test_cannot_seedn   s   
"zTestSeed.test_cannot_seedc                 C   s   t ttjt d S r    )r   rC   r   r6   r   rD   r#   r#   r$   test_invalid_initializations   s   z$TestSeed.test_invalid_initializationN)
__name__
__module____qualname__r;   r>   rE   rI   rR   rV   rW   r#   r#   r#   r$   r3   H   s    
r3   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 )Nr   intrJ   )r   r   rG   r   )r&   zerosr   r   r   r   )r9   r]   pr#   r#   r$   test_n_zerox   s
   zTestBinomial.test_n_zeroc                 C   s   t ttjdtj d S )NrG   )r   rC   r   r   r&   nanrD   r#   r#   r$   test_p_is_nan   s   zTestBinomial.test_p_is_nanN)rX   rY   rZ   r_   ra   r#   r#   r#   r$   r[   w   s    	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S )TestMultinomialc                 C      t dddg d S )Nr   皙?皙?r   r   rD   r#   r#   r$   
test_basic      zTestMultinomial.test_basicc                 C   s   t dg d d S )Nr   )rd   re           ri   ri   rf   rD   r#   r#   r$   test_zero_probability   rh   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 )Nr@   rO   )r   r   r7   r&   allr9   xr#   r#   r$   test_int_negative_interval   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 )Nr   rG   rG   r   r   r   r   r   r   r   )
r   r   r   r&   uint32shaperP   r   rB   floatr9   r^   r#   r#   r$   	test_size   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 )Nr   皙?rd   皙?r   rC   r   r   rD   r#   r#   r$   test_invalid_prob   s   z!TestMultinomial.test_invalid_probc                 C   s   t ttjdddg d S )Nr@   re   rd   r{   rD   r#   r#   r$   test_invalid_n      zTestMultinomial.test_invalid_nc                 C   st   t d}|t |dd d  }|dd d }td tjd|d}td tjdt |d}t|| d S )Ng      .@rG   rM   i۠iUr   )pvals)r&   r=   sumr   rT   r   ascontiguousarrayr   )r9   r^   r   
non_contigcontigr#   r#   r$   test_p_non_contiguous   s   


z%TestMultinomial.test_p_non_contiguousc                 C   sd   t jg dt jd}||  }d}tjt|d td| W d    d S 1 s+w   Y  d S )N)
Gz?r   &.>r   r   r   r   r   r   r   rJ   z-[\w\s]*pvals array is cast to 64-bit floatingmatchrG   )	r&   rP   float32r   pytestraisesrC   r   r   )r9   rn   r   r   r#   r#   r$   test_multinomial_pvals_float32   s   
"z.TestMultinomial.test_multinomial_pvals_float32c                 C   rc   )Ng      Y@rd   re   rf   rD   r#   r#   r$   test_multinomial_n_float   s   z(TestMultinomial.test_multinomial_n_floatN)rX   rY   rZ   rg   rj   ro   rw   r|   r}   r   r   r   r#   r#   r#   r$   rb      s    
rb   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
dd Zdd Zdd ZdS )TestSetStatec                 C   s$   d| _ t| j | _| j | _d S NiI)rT   r   r6   random_state	get_stater+   rD   r#   r#   r$   setup_method   s   zTestSetState.setup_methodc                 C   s<   | j d}| j | j | j d}tt||k d S )N   )r   tomaxint	set_stater+   r   r&   rl   r9   oldnewr#   r#   r$   rg      s   zTestSetState.test_basicc                 C   s@   | j jdd}| j | j | j jdd}tt||k d S NrM   size)r   standard_normalr   r+   r   r&   rl   r   r#   r#   r$   test_gaussian_reset   s   z TestSetState.test_gaussian_resetc                 C   sR   | j   | j  }| j jdd}| j | | j jdd}tt||k d S r   )r   r   r   r   r   r&   rl   )r9   r+   r   r   r#   r#   r$    test_gaussian_reset_in_media_res   s   

z-TestSetState.test_gaussian_reset_in_media_resc                 C   sz   | j d d }| jjdd}| j| | jjdd}| j| j  | jjdd}tt||k tt||k d S )NrH   r   r   )r+   r   r   r   r   r&   rl   )r9   Z	old_stateZx1Zx2Zx3r#   r#   r$   test_backwards_compatibility   s   z)TestSetState.test_backwards_compatibilityc                 C   s   | j dd d S Nr   )r   r   rD   r#   r#   r$   test_negative_binomial   s   z#TestSetState.test_negative_binomialc                 C   sx   t t }t )}|t}| }tt|dk t	|t
s"J |d dks*J W d    d S 1 s5w   Y  d S )NrG   r*   r   )r   r6   r   r	   recordRuntimeWarningr   r   len
isinstancedict)r9   rU   supwr+   r#   r#   r$   test_get_state_warning   s   
"z#TestSetState.test_get_state_warningc                 C   sl   | j  }d|dd   }tt| j j| tt| j jtj|td | j jdd}|d= tt| j j| d S )N)ZUnknownrG   rJ   FZlegacyr*   )	r   r   r   rC   r   rB   r&   rP   object)r9   r+   	new_stater#   r#   r$   !test_invalid_legacy_state_setting   s   

z.TestSetState.test_invalid_legacy_state_settingc                 C   sj   | j d | j d | j   | j jdd}t|d d tt| j }|jdd}t	|| d S )Nr   r   Fr   r.   rG   )
r   rT   random_sampler   r   r   pickleloadsdumpsr2   )r9   ZpickledZ	rs_unpickZ	unpickledr#   r#   r$   test_pickle   s   
zTestSetState.test_picklec                 C   s<   | j  }| j   | j | | j jdd}t|| d S )NFr   )r   __getstate__r   __setstate__r   r2   )r9   Z
attr_stater+   r#   r#   r$   test_state_setting
  s
   

zTestSetState.test_state_settingc                 C   s   t | jds
J d S )NzRandomState(MT19937))reprr   
startswithrD   r#   r#   r$   	test_repr  r~   zTestSetState.test_reprN)rX   rY   rZ   r   rg   r   r   r   r   r   r   r   r   r   r#   r#   r#   r$   r      s    	

r   c                	   @   s   e Zd ZejZej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ejjedjdk dddd Zdd Zdd ZdS )TestRandintc                 C   s   t t| jdtd d S NrG   rJ   )r   rB   rfuncru   rD   r#   r#   r$   test_unsupported_type     z!TestRandint.test_unsupported_typec                 C   s   | j D ]J}|tju rdnt|j}|tju rdnt|jd }tt| j|d ||d tt| j||d |d tt| j|||d tt| jdd|d qd S )Nr   r   rG   rJ   )	ityper&   booliinfominmaxr   rC   r   )r9   dtlbndubndr#   r#   r$   test_bounds_checking   s   
z TestRandint.test_bounds_checkingc                 C   s   | j D ]T}|tju rdnt|j}|tju rdnt|jd }|d }t| j||d d|d| |}t| j||d d|d| || d }t| j||d d|d| qd S )Nr   r   rG   r4   r   rK   )r   r&   r   r   r   r   r   r   )r9   r   r   r   tgtr#   r#   r$   test_rng_zero_and_extremes)  s   
z&TestRandint.test_rng_zero_and_extremesc                 C   s   | j D ]<}|tju rdnt|j}|tju rdnt|jd }z
| j|||d W q ty? } ztdt	| d }~ww d S )Nr   r   rG   rJ   zMNo error should have been raised, but one was with the following message:

%s)
r   r&   r   r   r   r   r   	ExceptionAssertionErrorstr)r9   r   r   r   er#   r#   r$   test_full_range7  s   
zTestRandint.test_full_rangec                 C   s   t   | jdd  D ] }dD ]}| jd|d|d}t| |k  t| dk qq| jdddtjd}t| dk  t| dk d S )NrG   )rN      r   r   i   r   r   )	r   rT   r   r   r   r   r   r&   r   )r9   r   r   valsr#   r#   r$   test_in_bounds_fuzzE  s   zTestRandint.test_in_bounds_fuzzc              
   C   s   dddddddddd	}| j dd  D ]8}td tjd	kr)| jd
dd|d}n| jd
dd|d }t|	t
j }t|t
|j |k qtd | jd
ddtd	t
j}t| }t|t
tj |k d S )NZ@509aea74d792fb931784c4b0135392c65aec64beee12b0cc167548a2c3d31e71Z@7b07f1a920e46f6d0fe02314155a2330bcfd7635e708da50e536c5ebb631a7d4Z@e577bfed6c935de944424667e3da285012e741892dcb7051a8f1ce68ab05c92fZ@0fbead0b06759df2cfb55e43148822d4a1ff953c7eb19a5b08445a63bb64fa9eZ@001aac3a5acb935a9b186cbe14a1ca064b8bb2dd0b045d48abeacf74d0203404)	r   int16int32rQ   int8uint16rs   uint64uint8rG   rS   littler   r   r4   r   r   )r   r   rT   sys	byteorderr   byteswaphashlibsha256viewr&   r   	hexdigestr   rK   namer   )r9   r   r   valresr#   r#   r$   test_repeatabilityT  s*   



zTestRandint.test_repeatabilitylr   zCannot test with 32-bit C longreasonc              
   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 ]$}td | jdgdgdggg d|d}t||d urX|n|d  q9d S )N)l   Q[ l   dm" l   0 )l   ;i l   vlm l   iXp )i\L_l   ZN iwU')l   <~ l   20l iP?)ifiwLl   @
Vo )l   o\ iwi
1])i{YSl   r=| i)l   0 l   l i&})l   QZ) l   (q l   '}Q )l   WqV1 iMl   't )inMl   b" l   )  )ial    l   S"f )idi98l   nDk )l   Jc i-l    )iMi'*l   l, )N)rO   rM   rM   i90  r@   r   rG   )r5   r   l       r   )r&   rP   r   rT   r   r   )r9   desiredr   rn   r#   r#   r$   .test_repeatability_32bit_boundary_broadcastingt  s8   

z:TestRandint.test_repeatability_32bit_boundary_broadcastingc                 C   s^   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&   rQ   r   r   r   r   r7   r   )r9   r   r   r   r   actualr#   r#   r$   test_int64_uint64_corner_case  s   z)TestRandint.test_int64_uint64_corner_casec                 C   s   | j D ]/}|tju rdnt|j}|tju rdnt|jd }| j|||d}t|jt| qtt	fD ];}|t	u r?dnd}|tu rGdnt|j}|tu rSdnt|jd }| j|||d}t
t|d  tt|| q7d S )Nr   r   rG   rJ   longr   rK   )r   r&   r   r   r   r   r   r   rK   r\   r   hasattrtype)r9   r   r   r   sampleZop_dtyper#   r#   r$   test_respect_dtype_singleton  s   
z(TestRandint.test_respect_dtype_singletonN)rX   rY   rZ   r   r7   r   r&   r   r   r   r   r   r   rs   rQ   r   r   r   r   r   r   r   r   r   markskipifr   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	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"dAdB Z#dCdD Z$dEdF Z%dGdH Z&dIdJ Z'dKdL Z(dMdN Z)dOdP Z*dQdR Z+dSdT Z,dUdV Z-dWdX Z.dYdZ Z/d[d\ Z0d]d^ Z1e2j34d_e56d`dadbe5j7dcgddde Z8dfdg Z9dhdi Z:djdk Z;dldm Z<dndo Z=dpdq Z>drds Z?dtdu Z@dvdw ZAdxdy ZBdzd{ ZCd|d} ZDd~d ZEdd ZFdd ZGdd ZHdd ZIdd ZJdd ZKdd ZLdd ZMdd ZNdd ZOdd ZPdd ZQdd ZRdd ZSdd ZTdd ZUdd ZVdd ZWdd ZXdd ZYdd ZZdS )TestRandomDistc                 C   
   d| _ d S r   rT   rD   r#   r#   r$   r        
zTestRandomDist.setup_methodc                 C   sF   t | j t dd}tddgddgddgg}t||d	d
 d S )NrM   r   _\*?琛?<p?,o?H?Z?   decimal)r   rT   randr&   rP   r   r9   r   r   r#   r#   r$   	test_rand     
zTestRandomDist.test_randc                 C   s*   t | j t  }d}t||dd d S )Nr   r   r   )r   rT   r   r   r   r#   r#   r$   test_rand_singleton  s   z"TestRandomDist.test_rand_singletonc                 C   sl   t | j t dd}tddgddgddgg}t||d	d
 t | j t  }t||d d	d
 d S )NrM   r   !<Oq?4u,?3?LI/DͿ@? @HS\@r   r   r   r   r   rT   Zrandnr&   rP   r   r   r#   r#   r$   
test_randn  s   
zTestRandomDist.test_randnc                 C   F   t | j t jdddd}tddgddgd	d
gg}t|| d S )Nc   rM   r   r      rM   )   )r   rT   r7   r&   rP   r   r   r#   r#   r$   test_randint     
zTestRandomDist.test_randintc                 C   s   t | j t }|t}t jdddd}tt|dk W d    n1 s)w   Y  t	ddgdd	gd
dgg}t
|| t | j t }|t}t jddd}tt|dk W d    n1 siw   Y  t
||d  d S )Nr  r  r  r   rG   r  rM   r  r  r  r     r   )r   rT   r	   r   DeprecationWarningrandom_integersr   r   r&   rP   r   )r9   r   r   r   r   r#   r#   r$   test_random_integers  s$   



z#TestRandomDist.test_random_integersc                 C   s   t | j t | j}|jdd}ttjjdkr-tjddgddgdd	ggtj	d
}ntjddgddgddggtj	d
}t
|| || j | }t
||d  d S )Nr  r   ii4Oi˕+iRKiriFbqirJ   l   K+W+y l   r&&J] l   ;	 l   nBz l   @XHT l   D+l++ r
  )r   rT   r6   r   r&   r   r   r   rP   rQ   r   )r9   rU   r   r   r#   r#   r$   test_tomaxint  s(   


zTestRandomDist.test_tomaxintc                 C   s   t  #}|t}ttdjtdj}tt	|dk W d    n1 s)w   Y  tdj}t
|| t  -}|t}tdj}t|tdj|tdj}tt	|dk W d    n1 slw   Y  t
|| d S )Nr   rG   )r	   r   r  r   r  r&   r   r   r   r   r   rK   r   )r9   r   r   r   r   Ztyperr#   r#   r$   test_random_integers_max_int  s$   



z+TestRandomDist.test_random_integers_max_intc                 C   sp   t  * t dt tttjtdj	 tttjtdj	tdj	 W d    d S 1 s1w   Y  d S )Nerrorr   )
warningscatch_warningssimplefilterr  r   r   r  r&   r   r   rD   r#   r#   r$   test_random_integers_deprecated  s   

"z.TestRandomDist.test_random_integers_deprecatedc                 C   sj   t | j t d}tddgddgddgg}t||dd	 t | j t  }t||d
 dd	 d S )Nr  r   r   r   r   r   r   r   r   r
  )r   rT   r   r&   rP   r   r   r#   r#   r$   test_random_sample,  s   

z!TestRandomDist.test_random_samplec                 C   s4   t | j t dd}tg d}t|| d S )NrN   )r   rM   r   rM   r   rT   choicer&   rP   r   r   r#   r#   r$   test_choice_uniform_replace8  s   z*TestRandomDist.test_choice_uniform_replacec                 C   s<   t | j t jddg dd}tg d}t|| d S )NrN   )皙?r'  皙?r(  r^   )rG   rG   r   r   r$  r   r#   r#   r$   test_choice_nonuniform_replace>  s   z-TestRandomDist.test_choice_nonuniform_replacec                 C   s8   t | j t jdddd}tg d}t|| d S )NrN   rM   Freplace)r   rG   rM   r$  r   r#   r#   r$   test_choice_uniform_noreplaceD     z,TestRandomDist.test_choice_uniform_noreplacec                 C   s>   t | j t jdddg dd}tg d}t|| d S )NrN   rM   F)r(  333333?r   r(  r,  r^   )r   rM   rG   r$  r   r#   r#   r$    test_choice_nonuniform_noreplaceJ  s   z/TestRandomDist.test_choice_nonuniform_noreplacec                 C   s8   t | j t g dd}tg d}t|| d S )N)r0   r1   cdrN   )r2  r3  r2  r3  r$  r   r#   r#   r$   test_choice_nonintegerP  r.  z%TestRandomDist.test_choice_nonintegerc                 C   s2  t j}tt|dd tt|dd tt|ddgddgg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 )Nr@   rM         @rG   r   rN   )rG   r   rM   rN   g      ?r)  )r'  r'  rd   rx   ry   r'  rL   Fr+  rH   )r@   )r@   rG   )rG   r   r   r0  )r   r%  r   rC   )r9   r   r#   r#   r$   test_choice_exceptionsV  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 )Nr(  rz   r   Tr+  Fr0  rG   rJ   r   )r   rM   )r(  r(  r(  r(  r'  rd   r   )rM   r   rN   r   i)r   r   r0   r1   )r   r&   Zisscalarr   r%  rP   emptyr   tuplendimitemr   rt   r=   r7   r   rC   )r9   r^   r0   Zarrr:   r#   r#   r$   test_choice_return_shapei  sL    z'TestRandomDist.test_choice_return_shapec                 C   s,   t g d}g d}tttj||d d S )N)*   rG   r   )NNNr)  )r&   rP   r   rC   r   r%  )r9   r0   r^   r#   r#   r$   test_choice_nan_probabilities  s   z,TestRandomDist.test_choice_nan_probabilitiesc              	   C   s|   t dd }d|dd d< t| j tjdd|d d d d}t| j tjddt |d d d d}t|| d S )Nr   rO   r5  rG   r   rM   r)  )r&   onesr   rT   r%  r   r   )r9   r^   r   r   r#   r#   r$   test_choice_p_non_contiguous  s    z+TestRandomDist.test_choice_p_non_contiguousc                 C   s(   t | j t d}d}t|| d S )Nr   s
   Ui+Wf)r   rT   bytesr   r   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 ] }t | j |g d}t | |}|g d}t|| q#d S )Nc                 S   s
   t g S r    )r&   rP   rn   r#   r#   r$   <lambda>  s   
 z-TestRandomDist.test_shuffle.<locals>.<lambda>c                 S   s   | S r    r#   rB  r#   r#   r$   rC    s    c                 S      t | t jS r    )r&   asarrayastyper   rB  r#   r#   r$   rC        c                 S   rD  r    )r&   rE  rF  r   rB  r#   r#   r$   rC    rG  c                 S   rD  r    )r&   rE  rF  Z	complex64rB  r#   r#   r$   rC    rG  c                 S   s   t | tS r    )r&   rE  rF  r   rB  r#   r#   r$   rC        c                 S   s   dd | D S )Nc                 S      g | ]}||fqS r#   r#   .0ir#   r#   r$   
<listcomp>      ATestRandomDist.test_shuffle.<locals>.<lambda>.<locals>.<listcomp>r#   rB  r#   r#   r$   rC    s    c                 S   s   t dd | D S )Nc                 S   s   g | ]}||gqS r#   r#   rJ  r#   r#   r$   rM    rN  rO  )r&   rE  rB  r#   r#   r$   rC    rN  c                 S   s   t | | gjS r    )r&   ZvstackTrB  r#   r#   r$   rC    rH  c                 S   s*   t dd | D dtfdtfgt jS )Nc                 S   rI  r#   r#   rJ  r#   r#   r$   rM    rN  rO  r0   r1   )r&   rE  r\   r   ZrecarrayrB  r#   r#   r$   rC    s    
c                 S   s(   t dd | D dtdfdt jdfgS )Nc                 S   rI  r#   r#   rJ  r#   r#   r$   rM    rN  rO  r0   rG   r1   )r&   rE  r   r   rB  r#   r#   r$   rC    s
    

rG   r   rM   rN   rO   r      r   	   r   
r   rG   rT  r   r   rN   rO   r   rS  rM   )r   rT   shuffler   )r9   convalistr   r   r#   r#   r$   test_shuffle  s&   
zTestRandomDist.test_shufflec                 C   s   t jt tddd d d}t jt dd d d}| }| }tdD ]2}t| t	t
|j|j  t
|j|j   t| t	t
|j|j  t
|j|j   q+dd }d S )	Nr   )rO   rN   rM   rG   r@   2   c                 S   s   t d}tttj| d S )NrM   )r&   rP   r   rB   r   rV  rm   r#   r#   r$   test_shuffle_invalid_objects  s   
zHTestRandomDist.test_shuffle_masked.<locals>.test_shuffle_invalid_objects)r&   ZmaZmasked_valuesZreshaper<   r=   copyr   rV  r   sorteddatamask)r9   r0   r1   Za_origZb_origrL  r[  r#   r#   r$   test_shuffle_masked  s   "
 
 z"TestRandomDist.test_shuffle_maskedc                 C   s   t | j g d}t |}g d}t|| t | j tg dj}t |}t|t|j t | j d}ttt j| t | j d}ttt j| d}g d}t | j t |}t|| d S )NrR  rU  Zabcdg333333?r   )
rT  r   r   rO   rG   rM   rN   rS  r   r   )	r   rT   Zpermutationr   r&   Z
atleast_2drP  r   
IndexError)r9   rX  r   r   Zarr_2dZ	bad_x_strZbad_x_floatZinteger_valr#   r#   r$   test_permutation  s(   



zTestRandomDist.test_permutationc                 C   J   t | j t jdddd}tddgddgd	d
gg}t||dd d S )Nr(  rz   r  r   gAKč?ghA?g ҇o>gάXf+q?gS1$?g'8zX ?r   r   )r   rT   betar&   rP   r   r   r#   r#   r$   	test_beta  s   zTestRandomDist.test_betac                 C   sl   t | j t jdddd}tddgddgd	d
gg}t|| t | j t dd}d}t|| d S )NgZd;Y@gv/?r  r   %   +   r<  0   .   -   )r   rT   r   r&   rP   r   r   r#   r#   r$   test_binomial  s   

zTestRandomDist.test_binomialc                 C   H   t | j t jddd}tddgddgdd	gg}t||d
d d S )NrZ  r  r   g+Q]uO@g+Q@gY7`ZqP@ge(fG@gvtR@gEIȋR@   r   )r   rT   	chisquarer&   rP   r   r   r#   r#   r$   test_chisquare	     
zTestRandomDist.test_chisquarec                 C   s   t | j tddg}t 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t j| t | j tddg}t |}t||d dd d S )NI<I@EZC@r  r   gi{s?gn,	A?g]^?gpDB?gȧ#y?gr?g1U?gѨUC?g^,%?gB?gډOoX;?gM`!O?r   r   HzG?缉ؗҜr
  )r   rT   r&   rP   	dirichletr   r   rC   )r9   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 )Nrq  rr  rG   rp   r   rq   rr   )
r&   rP   r   r   ru  rs   rt   r   rB   ru   rv   r#   r#   r$   test_dirichlet_size$  s   z"TestRandomDist.test_dirichlet_sizec                 C   s    t ddg}tttj| d S )Nrs  rt  )r&   rP   r   rC   r   ru  )r9   rv  r#   r#   r$   test_dirichlet_bad_alpha0  rF   z'TestRandomDist.test_dirichlet_bad_alphac                 C   sd   t g d}|d d d }t| j tj|dd}t| j tjt |dd}t|| d S )N)rq  g      rr  r   r  r   )r&   rP   r   rT   ru  r   r   )r9   r0   rv  r   r   r#   r#   r$   #test_dirichlet_alpha_non_contiguous5  s   z2TestRandomDist.test_dirichlet_alpha_non_contiguousc                 C   rl  )Ng8EGr?r  r   gNU?gbI0?gVU\@g:3@gU?g"qo?r   r   )r   rT   exponentialr&   rP   r   r   r#   r#   r$   test_exponential?  rp  zTestRandomDist.test_exponentialc                 C   &   t tjddd tttjdd d S Nr   scale       )r   r   r{  r   rC   rD   r#   r#   r$   test_exponential_0G     z!TestRandomDist.test_exponential_0c                 C   rc  )N   M   r  r   g7?g`j?g6"+?gܳ2>?gc>+Y?g_}S?r   r   )r   rT   fr&   rP   r   r   r#   r#   r$   test_fK     
zTestRandomDist.test_fc                 C   rc  )NrO   rM   r  r   gԏGM8@gȌ<@g	":@gv#)@gZQ?@gr[@@   r   )r   rT   gammar&   rP   r   r   r#   r#   r$   
test_gammaS  r  zTestRandomDist.test_gammac                 C   s*   t tjdddd tttjddd d S )Nr   )rt   r  r  )r   r   r  r   rC   rD   r#   r#   r$   test_gamma_0[  s   zTestRandomDist.test_gamma_0c                 C   D   t | j t jddd}tddgddgddgg}t|| d S )	N_c97ݚ?r  r   r   rS     rO   r  )r   rT   r   r&   rP   r   r   r#   r#   r$   test_geometric_     
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 !}|t t ttjtj t ttjtjgd  W d    d S 1 sIw   Y  d S )Nrx   r   ry   )	r   rC   r   r   r	   r   r   r&   r`   r9   r   r#   r#   r$   test_geometric_exceptionsg  s   
"z(TestRandomDist.test_geometric_exceptionsc                 C   rc  )Nr         @r  locr  r   gؕ?gs0 ?g#0g(xg-KD?gůp`@r   r   )r   rT   gumbelr&   rP   r   r   r#   r#   r$   test_gumbelq  r  zTestRandomDist.test_gumbelc                 C   r}  r~  )r   r   r  r   rC   rD   r#   r#   r$   test_gumbel_0y  r  zTestRandomDist.test_gumbel_0c                 C   s   t | j t jddddd}tddgddgddgg}t|| t jdd	d
dd}tg d}t|| t jdd	ddd}tg d}t|| t jd	dd
dd}tg d}t|| t jd	dddd}tg d}t|| d S )Ng333333$@g      @r  r  r   r   rT  rO   r   rM   rN   )rM   rM   rM   rM   r   r  )r  r  r  r  )r   r   r   r   )r   rT   r   r&   rP   r   r   r#   r#   r$   test_hypergeometric}  s&   




z"TestRandomDist.test_hypergeometricc                 C   rc  )Nr  r  r  r  g bO?g¸?gpJm	@g{K*t	@g:^w-gY~o?r   r   )r   rT   laplacer&   rP   r   r   r#   r#   r$   test_laplace  r  zTestRandomDist.test_laplacec                 C   r}  r~  )r   r   r  r   rC   rD   r#   r#   r$   test_laplace_0  r  zTestRandomDist.test_laplace_0c                 C   rc  )Nr  r  r  r  g#K-z?gJ?g92@gر[@g#+^˿g-@r   r   )r   rT   logisticr&   rP   r   r   r#   r#   r$   test_logistic  r  zTestRandomDist.test_logisticc                 C   rc  )Nr  r  r  )meansigmar   gɁ0@go,+4FB@g)ʭ6@g5?gnxNnP@g)U@rm  r   )r   rT   	lognormalr&   rP   r   r   r#   r#   r$   test_lognormal  r  zTestRandomDist.test_lognormalc                 C   s&   t tjddd tttjdd d S )Nr   )r  rG   r  )r   r   r  r   rC   rD   r#   r#   r$   test_lognormal_0  r  zTestRandomDist.test_lognormal_0c                 C   sD   t | j t jddd}tddgddgddgg}t|| d S )Ngƀ@?r  )r^   r   r   r   r  rM   )r   rT   r   r&   rP   r   r   r#   r#   r$   test_logseries  r  zTestRandomDist.test_logseriesc                 C   s   t ddks	J d S )Nr   rG   )r   r   rD   r#   r#   r$   test_logseries_zero  r   z"TestRandomDist.test_logseries_zerovalueri   r@         ?g      @c              	   C   s   t jddp tt t| W d    n1 sw   Y  tt tt |gd  W d    n1 s<w   Y  tt tt |gd d d d  W d    n1 saw   Y  W d    d S W d    d S 1 syw   Y  d S )NignoreZinvalidr   r   )r&   errstater   r   rC   r   r   rP   )r9   r  r#   r#   r$   test_logseries_exceptions  s   ""z(TestRandomDist.test_logseries_exceptionsc                 C   sd   t | j t 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   )rN   rM   rO   rN   r   r   )rO   r   r   r   r   rG   )rM   rN   rM   r   r   rN   )r   rG   rN   rM   r   rN   )rN   rN   r   rO   r   rM   )rN   rM   rN   r   rM   rN   )r   rT   r   r&   rP   r   r   r#   r#   r$   test_multinomial  s   
zTestRandomDist.test_multinomialc           	      C   s  t | j d}ddgddgg}d}t |||}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 ||}tddg}t||dd ddg}ddgddgg}ttt j|| tt j||dd t	t
t j||dd tjddgddggtjd}t }t || |t}t|dksJ W d    n1 sw   Y  td}td}t	t
t j||dd t	t
t jtd| t	t
t j|td t	t
t j|td d S )N)r  r   rG   r   r  g4Vj?gy'@gI1?gSB#@ge<@g 4W(@g6R¿?g9r
v"@gс?g{#@g4xg"f"@r   r   gN}O6?gYMW"@r   r  )Zcheck_validraiser(  rJ   other)r   rG   rG   rM   )r   rT   Zmultivariate_normalr&   rP   r   r   r   r   r   rC   r   r	   r   r   r]   Zeyer7  )	r9   r  Zcovr   r   r   r   r   mur#   r#   r$   test_multivariate_normal  s`   







z'TestRandomDist.test_multivariate_normalc                 C   r  )Nr   g|Pk?r  )nr^   r   iP  iI  i|  ic  i  i  )r   rT   r   r&   rP   r   r   r#   r#   r$   r     r  z%TestRandomDist.test_negative_binomialc                 C   s`   t  #}|t tttjdtj tttjdtjgd  W d    d S 1 s)w   Y  d S )Nr   r   )	r	   r   r   r   rC   r   r   r&   r`   r  r#   r#   r$   !test_negative_binomial_exceptions  s   


"z0TestRandomDist.test_negative_binomial_exceptionsc                 C   s   t | j t jdddd}tddgddgdd	gg}t||d
d t jdddd}tddgddgddgg}t||d
d t | j t jdddd}tddgddgddgg}t||d
d d S )NrO   r  )dfnoncr   gG7@g! ܴ*@g1}z9?@g+驸0@g26E_r#@gI!1@r  r   r   rd   g?gTD?g?T?g+qnl?g~|D?g밬?r   g{01#@gprs'@g~;$@g"dI@g*@g}',@)r   rT   noncentral_chisquarer&   rP   r   r   r#   r#   r$   test_noncentral_chisquare  s(   


z(TestRandomDist.test_noncentral_chisquarec                 C   L   t | j t jddddd}tddgdd	gd
dgg}t||dd d S )NrO   r   rG   r  )dfnumdfdenr  r   g^~?g/+c?g".@g@g7?gE?r  r   )r   rT   noncentral_fr&   rP   r   r   r#   r#   r$   test_noncentral_f0     

z TestRandomDist.test_noncentral_fc                 C   s0   t | j t jddtjd}t|sJ d S )NrO   r   )r  r  r  )r   rT   r  r&   r`   isnan)r9   r   r#   r#   r$   test_noncentral_f_nan9  s   z$TestRandomDist.test_noncentral_f_nanc                 C   rc  )Nr  r  r  r  g&&n@gOG_@g~w@gtw]տgm&5@g$>@r   r   )r   rT   normalr&   rP   r   r   r#   r#   r$   test_normal>  r  zTestRandomDist.test_normalc                 C   r}  r~  )r   r   r  r   rC   rD   r#   r#   r$   test_normal_0F  r  zTestRandomDist.test_normal_0c                 C   sL   t | j t jddd}tddgddgdd	gg}tjj||d
d d S )Nr  r  r0   r   gWI@gԨy@gy]0Ag1~\Aga@gY
7A   )Znulp)r   rT   paretor&   rP   ZtestingZassert_array_almost_equal_nulpr   r#   r#   r$   test_paretoJ  s   
zTestRandomDist.test_paretoc                 C   sD   t | j t jddd}tddgddgddgg}t|| d S )Nr  r  )lamr   r   rG   )r   rT   r   r&   rP   r   r   r#   r#   r$   test_poissonY  r  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 !}|t	 tttjt j
 tttjt j
gd  W d    d S 1 sQw   Y  d S )Nr   r@   r   )r&   r   r   r   rC   r   r   r	   r   r   r`   )r9   ZlambigZlamnegr   r#   r#   r$   test_poisson_exceptionsa  s   
"z&TestRandomDist.test_poisson_exceptionsc                 C   rl  )Nr  r  r  g;'?gܴw*?g5?gFG?g D]?g$`?r   r   )r   rT   powerr&   rP   r   r   r#   r#   r$   
test_powerm  rp  zTestRandomDist.test_powerc                 C   rl  )Nr   r  )r  r   g"u+@gN$MB*@g7v8B4@g55@ge&@g~ӣZ1@r  r   )r   rT   rayleighr&   rP   r   r   r#   r#   r$   test_rayleighu  rp  zTestRandomDist.test_rayleighc                 C   r}  r~  )r   r   r  r   rC   rD   r#   r#   r$   test_rayleigh_0}  r  zTestRandomDist.test_rayleigh_0c                 C   F   t | j t jdd}tddgddgddgg}t||d	d
 d S )Nr  r   gMoDL?g6([9grF=?gM'	- gVgVA:Qy?r   r   )r   rT   Zstandard_cauchyr&   rP   r   r   r#   r#   r$   test_standard_cauchy  r  z#TestRandomDist.test_standard_cauchyc                 C   r  )Nr  r   g܁?gc!z?g&N"@gNI@gF>?g~ME?r   r   )r   rT   Zstandard_exponentialr&   rP   r   r   r#   r#   r$   test_standard_exponential  r  z(TestRandomDist.test_standard_exponentialc                 C   rl  )NrM   r  )rt   r   g)@gľ@gd],q@gPk{@g9B#1@gv @r  r   )r   rT   standard_gammar&   rP   r   r   r#   r#   r$   test_standard_gamma  rp  z"TestRandomDist.test_standard_gammac                 C   r}  )Nr   )rt   r  )r   r   r  r   rC   rD   r#   r#   r$   test_standard_gamma_0  r  z$TestRandomDist.test_standard_gamma_0c                 C   r  )Nr  r   r  r  r  r  r  r	  r   r   )r   rT   r   r&   rP   r   r   r#   r#   r$   test_standard_normal  r  z#TestRandomDist.test_standard_normalc                 C   s0   t | j t  }td}t||dd d S )Nr  r   r   r  r   r#   r#   r$   test_randn_singleton  s   
z#TestRandomDist.test_randn_singletonc                 C   rl  )Nr   r  )r  r   g\VH?ghU%gH$M?gLgVqdzǿg㖀?r   r   )r   rT   
standard_tr&   rP   r   r   r#   r#   r$   test_standard_t  rp  zTestRandomDist.test_standard_tc                 C   r  )Ng{Gz@g(\u$@gףp=
W4@r  )leftmoderightr   g:}\)@ghرTj(@gǨL30@g6ڜA0@g6sh&@g<,@r  r   )r   rT   
triangularr&   rP   r   r   r#   r#   r$   test_triangular  r  zTestRandomDist.test_triangularc                 C   rc  )NGz?gGz%@r  )lowhighr   g>#H@gu@g#@g}d]#@gN@gw)* @r   r   )r   rT   uniformr&   rP   r   r   r#   r#   r$   test_uniform  r  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 )Nru   r   rG   g ؅W4vC)r  r  )
r&   Zfinfor   r   r   r  r   OverflowErrorinf	nextafter)r9   ZfminZfmaxfuncr#   r#   r$   test_uniform_range_bounds  s   z(TestRandomDist.test_uniform_range_boundsc                 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                   @      e Zd Zdd ZdS )zGTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingFloatc                 S      t r    rB   rD   r#   r#   r$   	__float__     zQTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingFloat.__float__N)rX   rY   rZ   r  r#   r#   r#   r$   ThrowingFloat      r  r  c                   @   r  )zITestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingIntegerc                 S   r  r    r  rD   r#   r#   r$   __int__  r  zQTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingInteger.__int__N)rX   rY   rZ   r  r#   r#   r#   r$   ThrowingInteger  r  r  rG   )	r&   ZndarrayrP   r   r   rB   r   r  r   )r9   r  Zthrowing_floatr  Zthrowing_intr#   r#   r$   !test_scalar_exception_propagation  s   
z0TestRandomDist.test_scalar_exception_propagationc                 C   rc  )Nr  p=
ף?r  r  kappar   gLZI@g'N"@gkr?gk=m@gډ?g^=\?r   r   r   rT   vonmisesr&   rP   r   r   r#   r#   r$   test_vonmises  r  zTestRandomDist.test_vonmisesc                 C   s2   t | j t jdddd}tt|  d S )Nri   gg5_PG>@B r  )r   rT   r  r   r&   Zisfiniterl   r9   rr#   r#   r$   test_vonmises_small  s   z"TestRandomDist.test_vonmises_smallc                 C   s<   t | j t jdddd}tg d}t||dd d S )Nri   g    cArM   r  )g A^>?g @R7?g x.r   r   r  r   r#   r#   r$   test_vonmises_large  s   z"TestRandomDist.test_vonmises_largec                 C   s.   t | j t jdtjd}tt| d S )Nri   )r  r  )r   rT   r  r&   r`   r   r  r  r#   r#   r$   test_vonmises_nan  s   z TestRandomDist.test_vonmises_nanc                 C   rc  )Nr  r  r  )r  r  r   gZ@g-g@g8
m?gVP"?g`V?gx^&3?r  r   )r   rT   waldr&   rP   r   r   r#   r#   r$   	test_wald  r  zTestRandomDist.test_waldc                 C   rl  )Nr  r  r  g@6?g]A?g(1ڥR?gU?g2_u?g^}-Q?r   r   )r   rT   weibullr&   rP   r   r   r#   r#   r$   test_weibull  rp  zTestRandomDist.test_weibullc                 C   s:   t | j tt jdddtd ttt jdd d S )Nr   r  r  r  )r0   )r   rT   r   r  r&   r]   r   rC   rD   r#   r#   r$   test_weibull_0  s   zTestRandomDist.test_weibull_0c                 C   r  )	Nr  r  r  B      rG   rM   rm  )r   rT   r   r&   rP   r   r   r#   r#   r$   	test_zipf  r  zTestRandomDist.test_zipfN)[rX   rY   rZ   r   r  r  r  r  r  r  r  r"  r#  r&  r*  r-  r1  r4  r6  r;  r=  r?  rA  rY  r`  rb  re  rk  ro  rw  rx  ry  rz  r|  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r  r   r   Zparametrizer&   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   r  r  r  r#   r#   r#   r$   r     s    .		


3			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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"dAS )BTestBroadcastc                 C   r   )N[r   rD   r#   r#   r$   r   '  r   zTestBroadcast.setup_methodc                 C   s   t | j d S r    )r   rT   rD   r#   r#   r$   set_seed*  s   zTestBroadcast.set_seedc                 C   sl   dg}dg}t j}tg d}|   ||d |}t||dd |   |||d }t||dd d S )Nr   rG   gl?g_x?g=*BN?rM   r  r   )r   r  r&   rP   r	  r   )r9   r  r  r  r   r   r#   r#   r$   r  -  s   zTestBroadcast.test_uniformc                 C      dg}dg}dg}t j}tg d}|   ||d |}t||dd tt||d | |   |||d }t||dd tt|||d  d S )Nr   rG   r@   )g-r@gYl@g-`w?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  	bad_scaler  r   r   r#   r#   r$   r  =     zTestBroadcast.test_normalc                 C      dg}dg}dg}dg}t j}tg d}|   ||d |}t||dd tt||d | tt||d | |   |||d }t||dd tt|||d  tt|||d  d S )	NrG   r   r@   rH   )gnQVf?gVKB?gv`k?rM   r  r   )r   rd  r&   rP   r	  r   r   rC   )r9   r0   r1   bad_aZbad_brd  r   r   r#   r#   r$   re  P      zTestBroadcast.test_betac                 C   V   dg}dg}t j}tg d}|   ||d }t||dd tt||d  d S NrG   r@   )g.UgZ?gi Prq?gEK?rM   r  r   )r   r{  r&   rP   r	  r   r   rC   )r9   r  r  r{  r   r   r#   r#   r$   r|  f     zTestBroadcast.test_exponentialc                 C   r  r  )r   r  r&   rP   r	  r   r   rC   )r9   rt   	bad_shapeZ	std_gammar   r   r#   r#   r$   r  s  r  z!TestBroadcast.test_standard_gammac                 C   r  )	NrG   r   r@   rH   )g.UgZ?gi Prq?gEK?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   rt   r  r  r  r  r   r   r#   r#   r$   r    r  zTestBroadcast.test_gammac                 C   r  )	NrG   r   r@   rH   )gbyʜ?go?gT7@rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  	bad_dfnum	bad_dfdenr  r   r   r#   r#   r$   r    r  zTestBroadcast.test_fc           
   	   C   s  dg}dg}dg}dg}dg}dg}t j}tg d}|   ||d ||}	t|	|dd	 tt|||tjgd sAJ t	t
||d || t	t
||d || t	t
||d || |   |||d |}	t|	|dd	 t	t
|||d | t	t
|||d | t	t
|||d | |   ||||d }	t|	|dd	 t	t
||||d  t	t
||||d  t	t
||||d  d S )
Nr   rM   rN   r   r@   rH   )gr^G"@g/*@gц!@r  r   )r   r  r&   rP   r	  r   rl   r  r`   r   rC   )
r9   r  r  r  r  r  bad_noncZnonc_fr   r   r#   r#   r$   r    s6   $zTestBroadcast.test_noncentral_fc                 C   s:   |    tddg}tjddddd}t||dd d S )Ngaz@gn4M;%?rz   r   r   r  r   )r	  r&   rP   r   r  r   )r9   r   r   r#   r#   r$   test_noncentral_f_small_df  s   z(TestBroadcast.test_noncentral_f_small_dfc                 C   r  )NrG   r@   )g~}N??g$F?g'rS?rM   r  r   )r   rn  r&   rP   r	  r   r   rC   )r9   r  bad_dfrn  r   r   r#   r#   r$   ro    r  zTestBroadcast.test_chisquarec                 C   r  )	NrG   r   r@   rH   )g w "@g2WR@gRY@rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  r  r  Znonc_chir   r   r#   r#   r$   r    r  z'TestBroadcast.test_noncentral_chisquarec                 C   h   dg}dg}t j}tg d}|   ||d }t||dd tt||d  ttt j|d  d S )NrG   r@   )gT@g~Fl@gƚtp?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  tr   r   r#   r#   r$   r       zTestBroadcast.test_standard_tc                 C   r  )Nr   rG   r@   )gW!@gOBg}?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  Z	bad_kappar  r   r   r#   r#   r$   r    r  zTestBroadcast.test_vonmisesc                 C   r  )NrG   r@   )g&8??geFX?gՋː?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r0   r  r  r   r   r#   r#   r$   r    r  zTestBroadcast.test_paretoc                 C   r  r  )r   r  r&   rP   r	  r   r   rC   )r9   r0   r  r  r   r   r#   r#   r$   r  '  r  zTestBroadcast.test_weibullc                 C   r  )NrG   r@   r
  rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r0   r  r  r   r   r#   r#   r$   r  5  r  zTestBroadcast.test_powerc                 C   r  )Nr   rG   r@   )g:BKc?g3k?g~q?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  r  r  r   r   r#   r#   r$   r  C  r  zTestBroadcast.test_laplacec                 C   r  )Nr   rG   r@   )g'Zy?gdSO=?g3L1?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  r  r  r   r   r#   r#   r$   r  V  r  zTestBroadcast.test_gumbelc                 C   s   dg}dg}dg}t j}tg d}|   ||d |}t||dd tt||d | |   |||d }t||dd tt|||d  tt dd	d d S )
Nr   rG   r@   )gy?gJR?g
~&?rM   r  r   r  ri   )	r   r  r&   rP   r	  r   r   rC   r   )r9   r  r  r  r  r   r   r#   r#   r$   r  i  s   zTestBroadcast.test_logisticc                 C   s   dg}dg}dg}t j}tg d}|   ||d |}t||dd tt||d | ttt j|d | |   |||d }t||dd tt|||d  ttt j||d  d S )Nr   rG   r@   )gH"@gܮ @gvv	:@rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  Z	bad_sigmar  r   r   r#   r#   r$   r  }  s   zTestBroadcast.test_lognormalc                 C   r  )NrG   r@   )gݱlo?g/ӂ?gz~R?rM   r  r   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  r  r   r   r#   r#   r$   r    r  zTestBroadcast.test_rayleighc                 C   s  dg}dg}dg}dg}t j}tg d}|   ||d |}t||dd tt||d | tt||d | ttt j|d | ttt j|d | |   |||d }t||dd tt|||d  tt|||d  tt|d	d tt|dd	 d S )
Nr   rG   r   rH   )g@m e?g`I߿?g|2?rM   r  r   ri   )r   r  r&   rP   r	  r   r   rC   )r9   r  r  Zbad_meanr  r  r   r   r#   r#   r$   r     s(   zTestBroadcast.test_waldc                 C   s  dg}dg}dg}dg}dg}|d \}}t j}tg d}	|   ||d ||}
t|
|	dd tt||d || tt||d || tt||d || |   |||d |}
t|
|	dd tt|||d | tt|||d | tt|||d | |   ||||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 )NrG   rM   r   rN   )g,P;bD @g%G @g,~S @r  r   r   ri   g      4@g      9@)r   r  r&   rP   r	  r   r   rC   )r9   r  r  r  Zbad_left_oneZbad_mode_oneZbad_left_twoZbad_mode_twor  r   r   r#   r#   r$   r    sF   
zTestBroadcast.test_triangularc           	      C      dg}dg}dg}dg}dg}t j}tg d}|   ||d |}t|| tt||d | tt||d | tt||d | |   |||d }t|| tt|||d  tt|||d  tt|||d  d S )NrG   r   r@         ?rG   rG   rG   rM   )r   r   r&   rP   r	  r   r   rC   )	r9   r  r^   bad_n	bad_p_one	bad_p_twoZbinomr   r   r#   r#   r$   rk    &   

zTestBroadcast.test_binomialc           	      C   r  )NrG   r   r@   r  )rG   r   rG   rM   )r   r   r&   rP   r	  r   r   rC   )	r9   r  r^   r   r!  r"  Z	neg_binomr   r   r#   r#   r$   r     r#  z$TestBroadcast.test_negative_binomialc                 C   sv   t  j}dg}dg}|d g}t j}tg d}|   ||d }t|| tt	||d  tt	||d  d S )NrG   r@   r   )rG   rG   r   rM   )
r   r6   Z_poisson_lam_maxr   r&   rP   r	  r   r   rC   )r9   Zmax_lamr  Zbad_lam_oneZbad_lam_twor   r   r   r#   r#   r$   r    s   


zTestBroadcast.test_poissonc                 C   s   dg}dg}t j}tg d}|   ||d }t|| tt||d  tjdd tt|tj	 tt|ddtj	g W d    d S 1 sJw   Y  d S )Nr   r   )r   r   rG   rM   r  r  )
r   r   r&   rP   r	  r   r   rC   r  r`   )r9   r0   r  r   r   r   r#   r#   r$   r    s   
"zTestBroadcast.test_zipfc                 C   h   dg}dg}dg}t j}tg d}|   ||d }t|| tt||d  tt||d  d S )Nr   r@   r  rq   rM   )r   r   r&   rP   r	  r   r   rC   )r9   r^   r!  r"  Zgeomr   r   r#   r#   r$   r  ,     
zTestBroadcast.test_geometricc                 C   s  dg}dg}dg}dg}dg}dg}dg}t j}tg d}	|   ||d ||}
t|
|	 tt||d || tt||d || tt||d || tt||d || |   |||d |}
t|
|	 tt|||d | tt|||d | tt|||d | tt|||d | |   ||||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 d S )NrG   r   r@   rH   r   rN   r  rM   r   r      )r   r   r&   rP   r	  r   r   rC   )r9   ZngoodZnbadZnsampleZ	bad_ngoodZbad_nbadZbad_nsample_oneZbad_nsample_twoZ	hypergeomr   r   r#   r#   r$   r  9  sD   


z!TestBroadcast.test_hypergeometricc                 C   r$  )Nr   r   r@   r  rM   )r   r   r&   rP   r	  r   r   rC   )r9   r^   r!  r"  r   r   r   r#   r#   r$   r  a  r%  zTestBroadcast.test_logseriesN)#rX   rY   rZ   r   r	  r  r  re  r|  r  r  r  r  r  ro  r  r  r  r  r  r  r  r  r  r  r  r   r  rk  r   r  r  r  r  r  r#   r#   r#   r$   r  $  sB    #((r  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 )NrN   )r<   seedsrD   r#   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|| q>t	 j
jdkr_tjdkr_t|| d S t|| d S )	Nr   )Threadc                    s&   g | ]\}} t ||fd qS ))targetargs)r   r6   )rK  r:   or)  functionr#   r$   rM  |  s    z-TestThread.check_function.<locals>.<listcomp>c                 S      g | ]}|  qS r#   )startrK  rn   r#   r#   r$   rM  ~  rN  c                 S   r/  r#   )joinr1  r#   r#   r$   rM    rN  rN   win32)Z	threadingr)  r&   r7  r   r(  zipr   r6   ZintprK   itemsizer   platformr   r   )r9   r.  szZout1Zout2r  r:   r,  r#   r-  r$   check_functionu  s   
zTestThread.check_functionc                 C      dd }| j |dd d S )Nc                 S   s   | j dd|d< d S )N'  r   .)r  r+   outr#   r#   r$   
gen_random  rh   z*TestThread.test_normal.<locals>.gen_random)r:  r7  r8  r9   r=  r#   r#   r$   r       zTestThread.test_normalc                 C   r9  )Nc                 S   s   | j tdd|d< d S )Nr   r4   r  .)r{  r&   r>  r;  r#   r#   r$   r=    s   z'TestThread.test_exp.<locals>.gen_randomrB  r>  r?  r@  r#   r#   r$   test_exp  rA  zTestThread.test_expc                 C   r9  )Nc                 S   s   | j ddgd dd|d< d S )Nr   r  r   r:  r   .)r   r;  r#   r#   r$   r=    s   z/TestThread.test_multinomial.<locals>.gen_random)r:  r   r>  r?  r@  r#   r#   r$   r    rA  zTestThread.test_multinomialN)rX   rY   rZ   r   r8  r  rC  r  r#   r#   r#   r$   r'  o  s    r'  c                   @   s,   e Z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 )Nr   rM   rN   rQ  )r&   rP   argOneargTwoargThreetgtShaperD   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 ]}||v r1|tdg}n|| j}t|j| j q"d S r   )r   r{  r  rn  r  r  r  r  r  r   r   r   r   r&   rP   rE  r   rt   rH  )r9   funcs	probfuncsr  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 r3tdg}n| j}|| j|}t|j| j || jd |}t|j| j || j|d }t|j| j q&d S )Nr   r   )r   r  r  rd  r  r  r  r  r  r  r  r  r  r   r   r&   rP   rF  rE  r   rt   rH  )r9   rI  rJ  r  rF  r<  r#   r#   r$   test_two_arg_funcs  s(   z*TestSingleEltArrayInput.test_two_arg_funcsc                 C   s   t jt jt jg}|D ]6}|| j| j| j}t|j| j	 || jd | j| j}t|j| j	 || j| jd | j}t|j| j	 q
d S Nr   )
r   r  r  r   rE  rF  rG  r   rt   rH  )r9   rI  r  r<  r#   r#   r$   test_three_arg_funcs  s   z,TestSingleEltArrayInput.test_three_arg_funcsN)rX   rY   rZ   r   rK  rL  rN  r#   r#   r#   r$   rD    s
    rD  c                 C   sD   t d | \}}}tt |}||ddi}t|jtdk d S )Nr  r   r   r   )r   rT   getattrr   rK   r&   )r%   fnamer+  r   r  r   r#   r#   r$   test_integer_dtype  s
   


rQ  c                 C   sd   t d | \}}}tt |}||ddi}tjdkr| }t|t	j
 }t||k d S )Nr  r   r  r   )r   rT   rO  r   r   r   r   r   r   r&   r   r   r   )r%   rP  r+  r   r  r   r   r#   r#   r$   test_integer_repeat  s   



rR  c                   C   s   t t tjdddgdd W d    n1 sw   Y  t t tjddgddd W d    n1 s9w   Y  t t tjddgddgdd W d    d S 1 s[w   Y  d S )NrG   r/  gffffff?)r   rG   r   r   )r   r   rC   r   r   r#   r#   r#   r$   test_broadcast_size_error  s   "rS  c                  C   s   t jtd} | d |  \}}}|d jjdksJ |d }|| |j	dd}t
|d |d  t|d d	 |d d	  t|d d
 |d d
  t
|d |d  t
|d |d  d S )Nr   rG   r   )r   Fr   r*   r+   r,   r-   r.   r/   )r&   r   r6   r   r   
__reduce__	__class__rX   r   r   r   r   )rU   Zctorr+  Zstate_ar1   Zstate_br#   r#   r$   &test_randomstate_ctor_old_style_pickle  s   

rV  c                 C   sR   t jd}|j}t j| tt jjjjt	|sJ t j
 }||u s'J d S rM  )r&   r   Zdefault_rngr*   r(   r   ZmtrandZ_randZ_bit_generatorr   r'   )r)   Zdef_bgbgZ	second_bgr#   r#   r$   test_hot_swap  s   
rX  c                 C   s   t d}tj| tjjdd}tjd tjjdd}t| t| |d dks.J |d d |d d ks<J |d d |d d ksJJ d S )	Nr   Fr   rG   r*   r   r+   inc)r   r&   r   r(   r   rT   print)r)   rW  r+   r   r#   r#   r$   test_seed_alt_bit_gen  s    r[  c                 C   s^   t j }td}t j| tjtdd t j| W d    d S 1 s(w   Y  d S )Nr   zstate must be for a PCG64r   )	r&   r   r   r   r(   r   r   rC   r   )r)   r+   rW  r#   r#   r$   test_state_error_alt_bit_gen,  s   
"r\  c                 C   s   t jd t jddd}td}|j}t j| t jjdd}|D ]}|| || ks0J q$t jd t jddd}t ||krHJ |j}|d d |d d ksYJ |d d |d d ksgJ d S )	Ń r      @r   Fr   r+   rY  )	r&   r   rT   r7   r   r+   r(   r   rl   )r)   r   rW  r+   Zstate_directZfieldZpcg_valsr   r#   r#   r$   test_swap_worked5  s    r_  c                 C   sL   t jtd t jddd}t jtd}|ddd}t|| d S )Nr]  r   r^  r   )r&   r   r(   r   r7   r6   r   )r)   Zsingleton_valsZrgZnon_singleton_valsr#   r#   r$   %test_swapped_singleton_against_directG  s
   r`  )3r   r   r   r  Znumpyr&   r   Znumpy.testingr   r   r   r   r   r   r   r	   r
   Znumpy.randomr   r   r   r>  r!   r   r   r   r"   Zfixturer%   r)   r2   r3   r[   rb   r   r   r   r  r   r   r'  rD  rQ  rR  rS  rV  rX  r[  r\  r_  r`  r#   r#   r#   r$   <module>   s    ,



/8W %      q    O/E
	