U
    d?                 	   @   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 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dk rd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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/d0 d0Z!G d1d2 d2Z"G d3d4 d4Z#G d5d6 d6Z$G d7d8 d8Z%d9d: Z&d;d< Z'd=d> Z(dS )?    N)assert_assert_raisesassert_equalassert_warnsassert_no_warningsassert_array_equalassert_array_almost_equalsuppress_warnings)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)param	INT_FUNCSINT_FUNC_HASHES)request r$   G/tmp/pip-unpacked-wheel-c1z6a0q8/numpy/random/tests/test_randomstate.pyint_func2   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_equal8   s
    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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_scalarA   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   r1   i  r   i  r2   i	  )r   r3   ranger   r4   nparanger5   r$   r$   r%   
test_arrayG   s    zTestSeed.test_arrayc                 C   s    t ttjd t ttjd d S )N      r   	TypeErrorr   r3   
ValueErrorr6   r$   r$   r%   test_invalid_scalarQ   s    zTestSeed.test_invalid_scalarc                 C   s\   t ttjdg t ttjdg t ttjdg t ttjdddg t ttjdddg d S )Nr=   r>   r      r   r?   rB   r$   r$   r%   test_invalid_arrayV   s
    zTestSeed.test_invalid_arrayc                 C   sR   t ttjtjg tjd t ttjdddgg t ttjdddgdddgg d S )NdtyperD   r            r   )r   rA   r   r3   r:   arrayint64rB   r$   r$   r%   test_invalid_array_shape^   s    z!TestSeed.test_invalid_array_shapec              	   C   s0   t td}tt |d W 5 Q R X d S )Nr     )r   r3   r   r   r@   seed)r6   rsr$   r$   r%   test_cannot_seedf   s    
zTestSeed.test_cannot_seedc                 C   s   t ttjt d S r   )r   rA   r   r3   r
   rB   r$   r$   r%   test_invalid_initializationk   s    z$TestSeed.test_invalid_initializationN)
__name__
__module____qualname__r8   r<   rC   rF   rN   rR   rS   r$   r$   r$   r%   r0   @   s   
r0   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   intrG   )r   r   rD   r   )r:   zerosr   r   r   r   )r6   rY   pr$   r$   r%   test_n_zerop   s    zTestBinomial.test_n_zeroc                 C   s   t ttjdtj d S )NrD   )r   rA   r   r   r:   nanrB   r$   r$   r%   test_p_is_nany   s    zTestBinomial.test_p_is_nanN)rT   rU   rV   r[   r]   r$   r$   r$   r%   rW   o   s   	rW   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 )Nr   皙?皙?r   r   rB   r$   r$   r%   
test_basic   s    zTestMultinomial.test_basicc                 C   s   t ddddddg d S )Nr   r_   r`           ra   rB   r$   r$   r%   test_zero_probability   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 )Nr>   rK   )r   r   r4   r:   allr6   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   rD   rD   r   r   r   r   r   r   r   )
r   r   r   r:   uint32shaperL   r   r@   floatr6   rZ   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   皙?r_   皙?r   rA   r   r   rB   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>   r`   r_   ru   rB   r$   r$   r%   test_invalid_n   s    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      .@rD   rI   i۠iUr   )pvals)r:   r;   sumr   rP   r   ascontiguousarrayr   )r6   rZ   rx   
non_contigcontigr$   r$   r%   test_p_non_contiguous   s    


z%TestMultinomial.test_p_non_contiguousc                 C   s^   t jddddddddddg
t jd}||  }d}tjt|d td| W 5 Q R X d S )NgGz?g&.>rG   z-[\w\s]*pvals array is cast to 64-bit floating)matchrD   )	r:   rL   float32ry   pytestraisesrA   r   r   )r6   rh   rx   r~   r$   r$   r%   test_multinomial_pvals_float32   s       z.TestMultinomial.test_multinomial_pvals_float32N)rT   rU   rV   rb   rd   ri   rq   rv   rw   r}   r   r$   r$   r$   r%   r^   ~   s   
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
dd Zdd Zdd ZdS )TestSetStatec                 C   s$   d| _ t| j | _| j | _d S NiI)rP   r   r3   random_state	get_stater(   rB   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:   rf   r6   oldnewr$   r$   r%   rb      s    zTestSetState.test_basicc                 C   s@   | j jdd}| j | j | j jdd}tt||k d S NrI   size)r   standard_normalr   r(   r   r:   rf   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:   rf   )r6   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 )NrE   r   r   )r(   r   r   r   r   r:   rf   )r6   	old_statex1Zx2Zx3r$   r$   r%   test_backwards_compatibility   s    z)TestSetState.test_backwards_compatibilityc                 C   s   | j dd d S Nr   )r   r   rB   r$   r$   r%   test_negative_binomial   s    z#TestSetState.test_negative_binomialc              	   C   sb   t t }t F}|t}| }tt|dk t	|t
sDt|d dksTtW 5 Q R X d S )NrD   r'   r   )r   r3   r   r	   recordRuntimeWarningr   r   len
isinstancedictAssertionError)r6   rQ   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)UnknownrD   rG   Flegacyr'   )	r   r   r   rA   r   r@   r:   rL   object)r6   r(   Z	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+   rD   )
r   rP   random_sampler   r   r   pickleloadsdumpsr/   )r6   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   r/   )r6   Z
attr_stater(   r$   r$   r%   test_state_setting   s
    

zTestSetState.test_state_settingc                 C   s   t | jdstd S )NzRandomState(MT19937))reprr   
startswithr   rB   r$   r$   r%   	test_repr  s    zTestSetState.test_reprN)rT   rU   rV   r   rb   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 NrD   rG   )r   r@   rfuncro   rB   r$   r$   r%   test_unsupported_type  s    z!TestRandint.test_unsupported_typec                 C   s   | j D ]}|tjkrdn
t|j}|tjkr2d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   rD   rG   )	ityper:   bool_iinfominmaxr   rA   r   )r6   dtlbndubndr$   r$   r%   test_bounds_checking  s    
z TestRandint.test_bounds_checkingc                 C   s   | j D ]}|tjkrdn
t|j}|tjkr2d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   rD   r1   r   rH   )r   r:   r   r   r   r   r   r   )r6   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krdn
t|j}|tjkr2dnt|jd }z| j|||d W q tk
r } ztdt	| W 5 d }~X Y qX qd S )Nr   r   rD   rG   zMNo error should have been raised, but one was with the following message:

%s)
r   r:   r   r   r   r   r   	Exceptionr   str)r6   r   r   r   er$   r$   r%   test_full_range,  s    
zTestRandint.test_full_rangec                 C   s   t   | jdd  D ]@}dD ]6}| 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 )NrD   )rJ      r   r   i   r   r   )	r   rP   r   r   r   r   r   r:   r   )r6   r   r   valsr$   r$   r%   test_in_bounds_fuzz:  s    zTestRandint.test_in_bounds_fuzzc              
   C   s   dddddddddd	}| j dd  D ]p}td tjd	krR| 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)	boolint16int32rM   int8uint16rm   uint64uint8rD   rO   littler   r   r1   r   r   )r   r   rP   sys	byteorderr   byteswaphashlibsha256viewr:   r   	hexdigestr   rH   namer   )r6   r   r   valresr$   r$   r%   test_repeatabilityI  s*    



zTestRandint.test_repeatabilitylr   zCannot test with 32-bit C long)reasonc                 C   s   t dddgdddgddd	ggd
ddgdddgdddggdddgdddgdddggdddgdd d!gd"d#d$ggd%d&d'gd(d)d*gd+d,d-ggg}d.D ]J}td/ | jd0gd1gd2ggd3d4d5g|d6}t||d k	r|n|d1  qd S )7Nl   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)rK   rI   rI   i90  r>   r   rD   r2   r   l       r   )r:   rL   r   rP   r   r   )r6   desiredr   rh   r$   r$   r%   .test_repeatability_32bit_boundary_broadcastingi  s6    
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:   rM   r   r   r   r   r4   r   )r6   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krdn
t|j}|tjkr2dnt|jd }| j|||d}t|jt| qt	t
tjjfD ]f}|t	krdn
t|j}|t	krdnt|jd }| j|||d}tt|d  tt|| qtd S )Nr   r   rD   rG   rH   )r   r:   r   r   r   r   r   r   rH   r   rX   compatlongr   hasattrtype)r6   r   r   r   sampler$   r$   r%   test_respect_dtype_singleton  s    
z(TestRandint.test_respect_dtype_singletonN)rT   rU   rV   r   r4   r   r:   r   r   r   r   r   r   rm   rM   r   r   r   r   r   r   r   r   r   markZ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   s
   d| _ d S r   rP   rB   r$   r$   r%   r     s    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 )NrI   r   _\*?琛?<p?,o?H?Z?   decimal)r   rP   randr:   rL   r   r6   r   r   r$   r$   r%   	test_rand  s    
zTestRandomDist.test_randc                 C   s*   t | j t  }d}t||dd d S )Nr   r   r   )r   rP   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 )NrI   r   !<Oq?4u,?3?LI/DͿ@? @HS\@r   r   r   r   r   rP   Zrandnr:   rL   r   r   r$   r$   r%   
test_randn  s    
zTestRandomDist.test_randnc                 C   sF   t | j t jdddd}tddgddgd	d
gg}t|| d S )Nc   rI   r   r      rI   )   )r   rP   r4   r:   rL   r   r   r$   r$   r%   test_randint  s    
zTestRandomDist.test_randintc              	   C   s   t | j t 0}|t}t jdddd}tt|dk W 5 Q R X 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 5 Q R X t
||d  d S )Nr  r  r  r   rD   r  rI   r  r  r  r     r   )r   rP   r	   r   DeprecationWarningrandom_integersr   r   r:   rL   r   )r6   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dkrXt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rG   l   K+W+y l   r&&J] l   ;	 l   nBz l   @XHT l   D+l++ r   )r   rP   r3   r   r:   r   rX   r   rL   rM   r   )r6   rQ   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 5 Q R X tdj}t
|| t  P}|t}tdj}t|tdj|tdj}tt	|dk W 5 Q R X t
|| d S )Nr   rD   )r	   r   r  r   r  r:   r   r   r   r   r   rH   r   )r6   r   r   r   r   Ztyperr$   r$   r%   test_random_integers_max_int  s     



z+TestRandomDist.test_random_integers_max_intc              	   C   sZ   t  H t dt tttjtdj	 tttjtdj	tdj	 W 5 Q R X d S )Nerrorr   )
warningscatch_warningssimplefilterr  r   r   r  r:   r   r   rB   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   rP   r   r:   rL   r   r   r$   r$   r%   test_random_sample  s    

z!TestRandomDist.test_random_samplec                 C   s8   t | j t dd}tddddg}t|| d S )NrJ   r   rI   r   rP   choicer:   rL   r   r   r$   r$   r%   test_choice_uniform_replace*  s    z*TestRandomDist.test_choice_uniform_replacec                 C   sD   t | j t jddddddgd}tddddg}t|| d S )NrJ   皙?皙?rZ   rD   r   r  r   r$   r$   r%   test_choice_nonuniform_replace0  s    z-TestRandomDist.test_choice_nonuniform_replacec                 C   s:   t | j t jdddd}tdddg}t|| d S )NrJ   rI   Freplacer   rD   r  r   r$   r$   r%   test_choice_uniform_noreplace6  s    z,TestRandomDist.test_choice_uniform_noreplacec                 C   sD   t | j t jdddddddgd}tddd	g}t|| d S )
NrJ   rI   Fr  333333?r   r  rZ   r   rD   r  r   r$   r$   r%    test_choice_nonuniform_noreplace<  s    z/TestRandomDist.test_choice_nonuniform_noreplacec                 C   s@   t | j t ddddgd}tddddg}t|| d S )Nr-   r.   cdrJ   r  r   r$   r$   r%   test_choice_nonintegerB  s    z%TestRandomDist.test_choice_nonintegerc              	   C   sD  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|ddddgdddgddggd tt|ddgdd	d	d
gd tt|ddgdddgd tt|ddgdd	d	gd tt|dddgddd tt|dddgddd tt|dddgddd tt|dddgddd tt|dddgdddddgd d S )Nr>   rI         @rD   r   rJ   g      ?r  r  r_   rr   rs   Fr  rE   )r>   )r>   rD   r   r!  )r   r  r   rA   )r6   r   r$   r$   r%   test_choice_exceptionsH  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 k tddg}tjdtd	}||d
< t tj|dd|k 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
 |k d}ddddddg}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  rt   r   Tr  Fr!  rD   rG   r   )r   rI   r  r_   r   )rI   r   rJ   r   i)r   r   r-   r.   )r   r:   Zisscalarr   r  rL   emptyr   tuplendimitemr   rn   r;   r4   r   rA   )r6   rZ   r-   Zarrr7   r$   r$   r%   test_choice_return_shape[  sL     z'TestRandomDist.test_choice_return_shapec                 C   s0   t dddg}d d d g}tttj||d d S )N*   rD   r   r  )r:   rL   r   rA   r   r  )r6   r-   rZ   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   rK   r&  rD   r   rI   r  )r:   onesr   rP   r  rz   r   )r6   rZ   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   rP   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 |ddddddddddg
}t | |}|ddddddddddg
}t|| qFd S )Nc                 S   s
   t g S r   )r:   rL   rh   r$   r$   r%   <lambda>      z-TestRandomDist.test_shuffle.<locals>.<lambda>c                 S   s   | S r   r$   r3  r$   r$   r%   r4    r5  c                 S   s   t | t jS r   )r:   asarrayastyper   r3  r$   r$   r%   r4    r5  c                 S   s   t | t jS r   )r:   r6  r7  r   r3  r$   r$   r%   r4    r5  c                 S   s   t | t jS r   )r:   r6  r7  Z	complex64r3  r$   r$   r%   r4    r5  c                 S   s   t | tS r   )r:   r6  r7  r   r3  r$   r$   r%   r4    r5  c                 S   s   dd | D S )Nc                 S   s   g | ]}||fqS r$   r$   .0ir$   r$   r%   
<listcomp>  s     ATestRandomDist.test_shuffle.<locals>.<lambda>.<locals>.<listcomp>r$   r3  r$   r$   r%   r4    r5  c                 S   s   t dd | D S )Nc                 S   s   g | ]}||gqS r$   r$   r8  r$   r$   r%   r;    s     r<  )r:   r6  r3  r$   r$   r%   r4    r5  c                 S   s   t | | gjS r   )r:   ZvstackTr3  r$   r$   r%   r4    r5  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$   r8  r$   r$   r%   r;    s     r<  r-   r.   )r:   r6  rX   r   Zrecarrayr3  r$   r$   r%   r4    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$   r8  r$   r$   r%   r;    s     r<  r-   rD   r.   )r:   r6  r   r   r3  r$   r$   r%   r4    s   
rD   r   rI   rJ   rK   r      r   	   r   )r   rP   shuffler   )r6   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 ]d}t| t	t
|j|j  t
|j|j   t| t	t
|j|j  t
|j|j   qVdd }d S )	Nr   )rK   rJ   rI   rD   r>   2   c                 S   s   t d}tttj| d S )NrI   )r:   rL   r   r@   r   rA  rg   r$   r$   r%   test_shuffle_invalid_objects  s    
zHTestRandomDist.test_shuffle_masked.<locals>.test_shuffle_invalid_objects)r:   maZmasked_valuesZreshaper9   r;   copyr   rA  r   sorteddatamask)r6   r-   r.   Za_origZb_origr:  rF  r$   r$   r%   test_shuffle_masked  s     "
 
 z"TestRandomDist.test_shuffle_maskedc                 C   s  t | j ddddddddd	d
g
}t |}d
dd	dddddddg
}t|| t | j tddddddddd	d
g
j}t |}t|t|j t | j d}ttt j| t | j d}ttt j| d}d	d
ddddddddg
}t | j t |}t|| d S )NrD   r   rI   rJ   rK   r   r?  r   r@  r   Zabcdg333333?r   )	r   rP   Zpermutationr   r:   Z
atleast_2dr=  r   
IndexError)r6   rC  r   r   Zarr_2dZ	bad_x_strZbad_x_floatZinteger_valr$   r$   r%   test_permutation  s(    

 

zTestRandomDist.test_permutationc                 C   sJ   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  rt   r  r   gAKč?ghA?g ҇o>gάXf+q?gS1$?g'8zX ?r   r   )r   rP   betar:   rL   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   rP   r   r:   rL   r   r   r$   r$   r%   test_binomial  s    

zTestRandomDist.test_binomialc                 C   sH   t | j t jddd}tddgddgdd	gg}t||d
d d S )NrE  r  r   g+Q]uO@g+Q@gY7`ZqP@ge(fG@gvtR@gEIȋR@   r   )r   rP   	chisquarer:   rL   r   r   r$   r$   r%   test_chisquare  s    
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   rP   r:   rL   	dirichletr   r   rA   )r6   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 )NrZ  r[  rD   rj   r   rk   rl   )
r:   rL   r   r   r^  rm   rn   r   r@   ro   rp   r$   r$   r%   test_dirichlet_size  s    z"TestRandomDist.test_dirichlet_sizec                 C   s    t ddg}tttj| d S )Nr\  r]  )r:   rL   r   rA   r   r^  )r6   r_  r$   r$   r%   test_dirichlet_bad_alpha"  s    z'TestRandomDist.test_dirichlet_bad_alphac                 C   sf   t dddg}|d d d }t| j tj|dd}t| j tjt |dd}t|| d S )NrZ  g      r[  r   r  r   )r:   rL   r   rP   r^  rz   r   )r6   r-   r_  r{   r|   r$   r$   r%   #test_dirichlet_alpha_non_contiguous'  s    z2TestRandomDist.test_dirichlet_alpha_non_contiguousc                 C   sH   t | j t jddd}tddgddgdd	gg}t||d
d d S )Ng8EGr?r  r   gNU?gbI0?gVU\@g:3@gU?g"qo?r   r   )r   rP   exponentialr:   rL   r   r   r$   r$   r%   test_exponential1  s    
zTestRandomDist.test_exponentialc                 C   s&   t tjddd tttjdd d S Nr   scale       )r   r   rd  r   rA   rB   r$   r$   r%   test_exponential_09  s    z!TestRandomDist.test_exponential_0c                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )N   M   r  r   g7?g`j?g6"+?gܳ2>?gc>+Y?g_}S?r   r   )r   rP   fr:   rL   r   r   r$   r$   r%   test_f=  s    
zTestRandomDist.test_fc                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )NrK   rI   r  r   gԏGM8@gȌ<@g	":@gv#)@gZQ?@gr[@@   r   )r   rP   gammar:   rL   r   r   r$   r$   r%   
test_gammaE  s    
zTestRandomDist.test_gammac                 C   s*   t tjdddd tttjddd d S )Nr   )rn   rh  ri  )r   r   rp  r   rA   rB   r$   r$   r%   test_gamma_0M  s    zTestRandomDist.test_gamma_0c                 C   sD   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   r?     rK   rk  )r   rP   r   r:   rL   r   r   r$   r$   r%   test_geometricQ  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 6}|t t ttjtj t ttjtjgd  W 5 Q R X d S )Nrr   r   rs   )	r   rA   r   r   r	   r   r   r:   r\   r6   r   r$   r$   r%   test_geometric_exceptionsY  s    
z(TestRandomDist.test_geometric_exceptionsc                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )Nrs         @r  locrh  r   gؕ?gs0 ?g#0g(xg-KD?gůp`@r   r   )r   rP   gumbelr:   rL   r   r   r$   r$   r%   test_gumbelc  s    
zTestRandomDist.test_gumbelc                 C   s&   t tjddd tttjdd d S rf  )r   r   r{  r   rA   rB   r$   r$   r%   test_gumbel_0k  s    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d
d
d
d
g}t|| t jdd	ddd}tddddg}t|| t jd	dd
dd}td	d	d	d	g}t|| t jd	dddd}td	d	d	d	g}t|| d S )Ng333333$@g      @ro  r  r   r   r@  rK   r   rI   rJ   r   rk  )r   rP   r   r:   rL   r   r   r$   r$   r%   test_hypergeometrico  s&    




z"TestRandomDist.test_hypergeometricc                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )Nrs  rx  r  ry  g bO?g¸?gpJm	@g{K*t	@g:^w-gY~o?r   r   )r   rP   laplacer:   rL   r   r   r$   r$   r%   test_laplace  s    
zTestRandomDist.test_laplacec                 C   s&   t tjddd tttjdd d S rf  )r   r   r  r   rA   rB   r$   r$   r%   test_laplace_0  s    zTestRandomDist.test_laplace_0c                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )Nrs  rx  r  ry  g#K-z?gJ?g92@gر[@g#+^˿g-@r   r   )r   rP   logisticr:   rL   r   r   r$   r$   r%   test_logistic  s    
zTestRandomDist.test_logisticc                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )Nrs  rx  r  )meansigmar   gɁ0@go,+4FB@g)ʭ6@g5?gnxNnP@g)U@rW  r   )r   rP   	lognormalr:   rL   r   r   r$   r$   r%   test_lognormal  s    
zTestRandomDist.test_lognormalc                 C   s&   t tjddd tttjdd d S )Nr   )r  rD   ri  )r   r   r  r   rA   rB   r$   r$   r%   test_lognormal_0  s    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  )rZ   r   r   r   rt  rI   )r   rP   r   r:   rL   r   r   r$   r$   r%   test_logseries  s    
zTestRandomDist.test_logseriesc                 C   s   t ddkstd S )Nr   rD   )r   r   r   rB   r$   r$   r%   test_logseries_zero  s    z"TestRandomDist.test_logseries_zerovaluerc   r>         ?g      @c              
   C   s   t jdd tt t| W 5 Q R X tt tt |gd  W 5 Q R X tt& tt |gd d d d  W 5 Q R X W 5 Q R X d S )Nignoreinvalidr   r   )r:   errstater   r   rA   r   r   rL   )r6   r  r$   r$   r%   test_logseries_exceptions  s     z(TestRandomDist.test_logseries_exceptionsc                 C   s   t | j t jddgd dd}tddddd	d	gdd	d
d	d	dggddddddgd	dddddggddd	dd	dgdddd	ddggg}t|| d S )Nr   UUUUUU?r   r  r   rJ   rI   rK   r   r   rD   r   )r   rP   r   r:   rL   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s6tW 5 Q R X 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)rs  r   rD   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  rG   other)r   rD   rD   rI   )r   rP   Zmultivariate_normalr:   rL   r   r   r   r   r   rA   r   r	   r   r   r   rY   Zeyer(  )	r6   r  Zcovr   r   r   r   r   mur$   r$   r%   test_multivariate_normal  sd    




   z'TestRandomDist.test_multivariate_normalc                 C   sF   t | j t jdddd}tddgddgd	d
gg}t|| d S )Nr   g|Pk?r  )nrZ   r   iP  iI  i|  ic  i  i  )r   rP   r   r:   rL   r   r   r$   r$   r%   r     s    
z%TestRandomDist.test_negative_binomialc              	   C   sJ   t  :}|t tttjdtj tttjdtjgd  W 5 Q R X d S )Nr   r   )	r	   r   r   r   rA   r   r   r:   r\   rv  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 )NrK   r  )dfnoncr   gG7@g! ܴ*@g1}z9?@g+驸0@g26E_r#@gI!1@ro  r   r   r_   g?gTD?g?T?g+qnl?g~|D?g밬?r   g{01#@gprs'@g~;$@g"dI@g*@g}',@)r   rP   noncentral_chisquarer:   rL   r   r   r$   r$   r%   test_noncentral_chisquare  s(    


z(TestRandomDist.test_noncentral_chisquarec                 C   sL   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 )NrK   r   rD   r  )dfnumdfdenr  r   g^~?g/+c?g".@g@g7?gE?ro  r   )r   rP   noncentral_fr:   rL   r   r   r$   r$   r%   test_noncentral_f"  s    

z TestRandomDist.test_noncentral_fc                 C   s0   t | j t jddtjd}t|s,td S )NrK   r   )r  r  r  )r   rP   r  r:   r\   isnanr   )r6   r   r$   r$   r%   test_noncentral_f_nan+  s    z$TestRandomDist.test_noncentral_f_nanc                 C   sJ   t | j t jdddd}tddgddgd	d
gg}t||dd d S )Nrs  rx  r  ry  g&&n@gOG_@g~w@gtw]տgm&5@g$>@r   r   )r   rP   normalr:   rL   r   r   r$   r$   r%   test_normal0  s    
zTestRandomDist.test_normalc                 C   s&   t tjddd tttjdd d S rf  )r   r   r  r   rA   rB   r$   r$   r%   test_normal_08  s    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 )Nrs  r  r-   r   gWI@gԨy@gy]0Ag1~\Aga@gY
7A   )Znulp)r   rP   paretor:   rL   ZtestingZassert_array_almost_equal_nulpr   r$   r$   r%   test_pareto<  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 )Nrs  r  )lamr   r   rD   )r   rP   r   r:   rL   r   r   r$   r$   r%   test_poissonK  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 6}|t	 tttjt j
 tttjt j
gd  W 5 Q R X d S )Nr   r>   r   )r:   r   r   r   rA   r   r   r	   r   r   r\   )r6   ZlambigZlamnegr   r$   r$   r%   test_poisson_exceptionsS  s    
z&TestRandomDist.test_poisson_exceptionsc                 C   sH   t | j t jddd}tddgddgdd	gg}t||d
d d S )Nrs  r  r  g;'?gܴw*?g5?gFG?g D]?g$`?r   r   )r   rP   powerr:   rL   r   r   r$   r$   r%   
test_power_  s    
zTestRandomDist.test_powerc                 C   sH   t | j t jddd}tddgddgdd	gg}t||d
d d S )Nr   r  )rh  r   g"u+@gN$MB*@g7v8B4@g55@ge&@g~ӣZ1@ro  r   )r   rP   rayleighr:   rL   r   r   r$   r$   r%   test_rayleighg  s    
zTestRandomDist.test_rayleighc                 C   s&   t tjddd tttjdd d S rf  )r   r   r  r   rA   rB   r$   r$   r%   test_rayleigh_0o  s    zTestRandomDist.test_rayleigh_0c                 C   sF   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   rP   Zstandard_cauchyr:   rL   r   r   r$   r$   r%   test_standard_cauchys  s    
z#TestRandomDist.test_standard_cauchyc                 C   sF   t | j t jdd}tddgddgddgg}t||d	d
 d S )Nr  r   g܁?gc!z?g&N"@gNI@gF>?g~ME?r   r   )r   rP   Zstandard_exponentialr:   rL   r   r   r$   r$   r%   test_standard_exponential{  s    
z(TestRandomDist.test_standard_exponentialc                 C   sH   t | j t jddd}tddgddgdd	gg}t||d
d d S )NrI   r  )rn   r   g)@gľ@gd],q@gPk{@g9B#1@gv @ro  r   )r   rP   standard_gammar:   rL   r   r   r$   r$   r%   test_standard_gamma  s    
z"TestRandomDist.test_standard_gammac                 C   s&   t tjddd tttjdd d S )Nr   )rn   ri  )r   r   r  r   rA   rB   r$   r$   r%   test_standard_gamma_0  s    z$TestRandomDist.test_standard_gamma_0c                 C   sF   t | j t jdd}tddgddgddgg}t||d	d
 d S )Nr  r   r   r   r   r   r   r   r   r   )r   rP   r   r:   rL   r   r   r$   r$   r%   test_standard_normal  s    
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   sH   t | j t 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\VH?ghU%gH$M?gLgVqdzǿg㖀?r   r   )r   rP   
standard_tr:   rL   r   r   r$   r$   r%   test_standard_t  s    
zTestRandomDist.test_standard_tc                 C   sL   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 )Ng{Gz@g(\u$@gףp=
W4@r  )leftmoderightr   g:}\)@ghرTj(@gǨL30@g6ڜA0@g6sh&@g<,@ro  r   )r   rP   
triangularr:   rL   r   r   r$   r$   r%   test_triangular  s    

zTestRandomDist.test_triangularc                 C   sJ   t | j t 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  )lowhighr   g>#H@gu@g#@g}d]#@gN@gw)* @r   r   )r   rP   uniformr:   rL   r   r   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 )Nro   r   rD   g ؅W4vC)r  r  )
r:   Zfinfor   r   r   r  r   OverflowErrorinf	nextafter)r6   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                   @   s   e Zd Zdd ZdS )zGTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingFloatc                 S   s   t d S r   r@   rB   r$   r$   r%   	__float__  s    zQTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingFloat.__float__N)rT   rU   rV   r  r$   r$   r$   r%   ThrowingFloat  s   r  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   r  rB   r$   r$   r%   __int__  s    zQTestRandomDist.test_scalar_exception_propagation.<locals>.ThrowingInteger.__int__N)rT   rU   rV   r  r$   r$   r$   r%   ThrowingInteger  s   r  rD   )	r:   ZndarrayrL   r   r   r@   r   r  r   )r6   r  Zthrowing_floatr  Zthrowing_intr$   r$   r%   !test_scalar_exception_propagation  s    
z0TestRandomDist.test_scalar_exception_propagationc                 C   sJ   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  p=
ף?r  r  kappar   gLZI@g'N"@gkr?gk=m@gډ?g^=\?r   r   r   rP   vonmisesr:   rL   r   r   r$   r$   r%   test_vonmises  s    
zTestRandomDist.test_vonmisesc                 C   s2   t | j t jdddd}tt|  d S )Nrc   gg5_PG>@B r  )r   rP   r  r   r:   isfiniterf   r6   rr$   r$   r%   test_vonmises_small  s    z"TestRandomDist.test_vonmises_smallc                 C   s>   t | j t jdddd}tdddg}t||dd	 d S )
Nrc   g    cArI   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 )Nrc   )r  r  )r   rP   r  r:   r\   r   r  r  r$   r$   r%   test_vonmises_nan  s    z TestRandomDist.test_vonmises_nanc                 C   sJ   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  r  r  )r  rh  r   gZ@g-g@g8
m?gVP"?g`V?gx^&3?ro  r   )r   rP   waldr:   rL   r   r   r$   r$   r%   	test_wald  s    
zTestRandomDist.test_waldc                 C   sH   t | j t jddd}tddgddgdd	gg}t||d
d d S )Nr  r  r  g@6?g]A?g(1ڥR?gU?g2_u?g^}-Q?r   r   )r   rP   weibullr:   rL   r   r   r$   r$   r%   test_weibull   s    
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   rk  r  ri  )r-   )r   rP   r   r  r:   rY   r   rA   rB   r$   r$   r%   test_weibull_0  s    zTestRandomDist.test_weibull_0c                 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  r  B      rD   rI   rW  )r   rP   r   r:   rL   r   r   r$   r$   r%   	test_zipf  s    
zTestRandomDist.test_zipfN)[rT   rU   rV   r   r   r   r   r	  r  r  r  r  r  r  r  r  r"  r%  r'  r,  r.  r0  r2  rD  rL  rN  rP  rV  rY  r`  ra  rb  rc  re  rj  rn  rq  rr  ru  rw  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   s
   d| _ d S )N[r   rB   r$   r$   r%   r     s    zTestBroadcast.setup_methodc                 C   s   t | j d S r   )r   rP   rB   r$   r$   r%   set_seed  s    zTestBroadcast.set_seedc                 C   sn   dg}dg}t j}tdddg}|   ||d |}t||dd |   |||d }t||dd d S )	Nr   rD   l?_x?=*BN?rI   ro  r   )r   r  r:   rL   r  r   )r6   r  r  r  r   r   r$   r$   r%   r    s    zTestBroadcast.test_uniformc                 C   s   dg}dg}dg}t j}tdddg}|   ||d |}t||dd	 tt||d | |   |||d }t||dd	 tt|||d  d S )
Nr   rD   r>   g-r@gYl@g-`w?rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   rz  rh  	bad_scaler  r   r   r$   r$   r%   r  /  s     zTestBroadcast.test_normalc                 C   s   dg}dg}dg}dg}t j}tdddg}|   ||d |}t||d	d
 tt||d | tt||d | |   |||d }t||d	d
 tt|||d  tt|||d  d S )NrD   r   r>   rE   gnQVf?gVKB?gv`k?rI   ro  r   )r   rO  r:   rL   r  r   r   rA   )r6   r-   r.   bad_aZbad_brO  r   r   r$   r$   r%   rP  B  s&    zTestBroadcast.test_betac                 C   sX   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  d S 	NrD   r>   g.UgZ?gi Prq?gEK?rI   ro  r   )r   rd  r:   rL   r  r   r   rA   )r6   rh  r  rd  r   r   r$   r$   r%   re  X  s    zTestBroadcast.test_exponentialc                 C   sX   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  d S r  )r   r  r:   rL   r  r   r   rA   )r6   rn   	bad_shapeZ	std_gammar   r   r$   r$   r%   r  e  s    z!TestBroadcast.test_standard_gammac                 C   s   dg}dg}dg}dg}t j}tdddg}|   ||d |}t||d	d
 tt||d | tt||d | |   |||d }t||d	d
 tt|||d  tt|||d  d S )NrD   r   r>   rE   g.UgZ?gi Prq?gEK?rI   ro  r   )r   rp  r:   rL   r  r   r   rA   )r6   rn   rh  r  r  rp  r   r   r$   r$   r%   rq  r  s&    zTestBroadcast.test_gammac                 C   s   dg}dg}dg}dg}t j}tdddg}|   ||d |}t||d	d
 tt||d | tt||d | |   |||d }t||d	d
 tt|||d  tt|||d  d S )NrD   r   r>   rE   gbyʜ?go?gT7@rI   ro  r   )r   rm  r:   rL   r  r   r   rA   )r6   r  r  	bad_dfnum	bad_dfdenrm  r   r   r$   r$   r%   rn    s&    zTestBroadcast.test_fc           
   	   C   s  dg}dg}dg}dg}dg}dg}t j}tddd	g}|   ||d ||}	t|	|d
d tt|||tjgd st	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   rI   rJ   r   r>   rE   gr^G"@g/*@gц!@ro  r   )r   r  r:   rL   r  r   rf   r  r\   r   r   rA   )
r6   r  r  r  r  r  bad_noncZnonc_fr   r   r$   r$   r%   r    s<    $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;%?rt   r   r   ro  r   )r  r:   rL   r   r  r   )r6   r   r   r$   r$   r%   test_noncentral_f_small_df  s    z(TestBroadcast.test_noncentral_f_small_dfc                 C   sX   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  d S )	NrD   r>   g~}N??g$F?g'rS?rI   ro  r   )r   rX  r:   rL   r  r   r   rA   )r6   r  bad_dfrX  r   r   r$   r$   r%   rY    s    zTestBroadcast.test_chisquarec                 C   s   dg}dg}dg}dg}t j}tdddg}|   ||d |}t||d	d
 tt||d | tt||d | |   |||d }t||d	d
 tt|||d  tt|||d  d S )NrD   r   r>   rE   g w "@g2WR@gRY@rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   r  r  r  r  Znonc_chir   r   r$   r$   r%   r    s&    z'TestBroadcast.test_noncentral_chisquarec                 C   sj   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  ttt j|d  d S )	NrD   r>   gT@g~Fl@gƚtp?rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   r  r  tr   r   r$   r$   r%   r    s    zTestBroadcast.test_standard_tc                 C   s   dg}dg}dg}t j}tdddg}|   ||d |}t||dd	 tt||d | |   |||d }t||dd	 tt|||d  d S )
Nr   rD   r>   gW!@gOBg}?rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   r  r  Z	bad_kappar  r   r   r$   r$   r%   r    s     zTestBroadcast.test_vonmisesc                 C   sj   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  ttt j|d  d S )	NrD   r>   g&8??geFX?gՋː?rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   r-   r  r  r   r   r$   r$   r%   r    s    zTestBroadcast.test_paretoc                 C   sj   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  ttt j|d  d S r  )r   r  r:   rL   r  r   r   rA   )r6   r-   r  r  r   r   r$   r$   r%   r    s    zTestBroadcast.test_weibullc                 C   sj   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  ttt j|d  d S )	NrD   r>   r  r  r  rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   r-   r  r  r   r   r$   r$   r%   r  '  s    zTestBroadcast.test_powerc                 C   s   dg}dg}dg}t j}tdddg}|   ||d |}t||dd	 tt||d | |   |||d }t||dd	 tt|||d  d S )
Nr   rD   r>   g:BKc?g3k?g~q?rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   rz  rh  r  r  r   r   r$   r$   r%   r  5  s     zTestBroadcast.test_laplacec                 C   s   dg}dg}dg}t j}tdddg}|   ||d |}t||dd	 tt||d | |   |||d }t||dd	 tt|||d  d S )
Nr   rD   r>   g'Zy?gdSO=?g3L1?rI   ro  r   )r   r{  r:   rL   r  r   r   rA   )r6   rz  rh  r  r{  r   r   r$   r$   r%   r|  H  s     zTestBroadcast.test_gumbelc                 C   s   dg}dg}dg}t j}tdddg}|   ||d |}t||dd	 tt||d | |   |||d }t||dd	 tt|||d  tt d
dd
 d S )Nr   rD   r>   gy?gJR?g
~&?rI   ro  r   r  rc   )	r   r  r:   rL   r  r   r   rA   r   )r6   rz  rh  r  r  r   r   r$   r$   r%   r  [  s"    zTestBroadcast.test_logisticc                 C   s   dg}dg}dg}t j}tdddg}|   ||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   rD   r>   gH"@gܮ @gvv	:@rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   r  r  Z	bad_sigmar  r   r   r$   r$   r%   r  o  s$    zTestBroadcast.test_lognormalc                 C   sX   dg}dg}t j}tdddg}|   ||d }t||dd tt||d  d S )	NrD   r>   gݱlo?g/ӂ?gz~R?rI   ro  r   )r   r  r:   rL   r  r   r   rA   )r6   rh  r  r  r   r   r$   r$   r%   r    s    zTestBroadcast.test_rayleighc                 C   s  dg}dg}dg}dg}t j}tdddg}|   ||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   rD   r   rE   g@m e?g`I߿?g|2?rI   ro  r   rc   )r   r  r:   rL   r  r   r   rA   )r6   r  rh  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dddg}	|   ||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 )NrD   rI   r   rJ   g,P;bD @g%G @g,~S @ro  r   r   rc   g      4@g      9@)r   r  r:   rL   r  r   r   rA   )r6   r  r  r  Zbad_left_oneZbad_mode_oneZbad_left_twoZbad_mode_twor  r   r   r$   r$   r%   r    sL    
zTestBroadcast.test_triangularc           	      C   s   dg}dg}dg}dg}dg}t j}tdddg}|   ||d |}t|| tt||d | tt||d | tt||d | |   |||d }t|| tt|||d  tt|||d  tt|||d  d S )NrD   r   r>         ?rI   )r   r   r:   rL   r  r   r   rA   )	r6   r  rZ   bad_n	bad_p_one	bad_p_twoZbinomr   r   r$   r$   r%   rV    s&    

zTestBroadcast.test_binomialc           	      C   s   dg}dg}dg}dg}dg}t j}tdddg}|   ||d |}t|| tt||d | tt||d | tt||d | |   |||d }t|| tt|||d  tt|||d  tt|||d  d S )NrD   r   r>   r  r   rI   )r   r   r:   rL   r  r   r   rA   )	r6   r  rZ   r  r  r  Z	neg_binomr   r   r$   r$   r%   r     s&    

z$TestBroadcast.test_negative_binomialc                 C   sx   t  j}dg}dg}|d g}t j}tdddg}|   ||d }t|| tt	||d  tt	||d  d S )NrD   r>   r   r   rI   )
r   r3   Z_poisson_lam_maxr   r:   rL   r  r   r   rA   )r6   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dddg}|   ||d }t|| tt||d  tjdd( tt|tj	 tt|ddtj	g W 5 Q R X d S )Nr   r   rD   rI   r  r  )
r   r   r:   rL   r  r   r   rA   r  r\   )r6   r-   r  r   r   r   r$   r$   r%   r    s    
zTestBroadcast.test_zipfc                 C   sj   dg}dg}dg}t j}tdddg}|   ||d }t|| tt||d  tt||d  d S )Nr   r>   r  r   rI   )r   r   r:   rL   r  r   r   rA   )r6   rZ   r  r  Zgeomr   r   r$   r$   r%   ru    s    
zTestBroadcast.test_geometricc                 C   s  dg}dg}dg}dg}dg}dg}dg}t j}tdddg}	|   ||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 )NrD   r   r>   rE   r   rJ   rI   r   r      )r   r   r:   rL   r  r   r   rA   )r6   ZngoodZnbadZnsampleZ	bad_ngoodZbad_nbadZbad_nsample_oneZbad_nsample_twoZ	hypergeomr   r   r$   r$   r%   r~  +  sD    


z!TestBroadcast.test_hypergeometricc                 C   sj   dg}dg}dg}t j}tdddg}|   ||d }t|| tt||d  tt||d  d S )Nr   r   r>   rD   rI   )r   r   r:   rL   r  r   r   rA   )r6   rZ   r  r  r   r   r   r$   r$   r%   r  S  s    
zTestBroadcast.test_logseriesN)#rT   rU   rV   r   r  r  r  rP  re  r  rq  rn  r  r  rY  r  r  r  r  r  r  r  r|  r  r  r  r  r  rV  r   r  r  ru  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 )
TestThreadc                 C   s   t d| _d S )NrJ   )r9   seedsrB   r$   r$   r%   r   c  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rtjdkrt|| n
t|| d S )	Nr   )Threadc                    s&   g | ]\}} t ||fd qS ))targetargs)r   r3   )r9  r7   or  functionr$   r%   r;  m  s   z-TestThread.check_function.<locals>.<listcomp>c                 S   s   g | ]}|  qS r$   )startr9  rh   r$   r$   r%   r;  o  s     c                 S   s   g | ]}|  qS r$   )joinr  r$   r$   r%   r;  p  s     rJ   win32)	threadingr  r:   r(  r   r  zipr   r3   ZintprH   itemsizer   platformr   r   )r6   r  szZout1Zout2r  r7   r   r$   r  r%   check_functionf  s    
zTestThread.check_functionc                 C   s   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}  s    z*TestThread.test_normal.<locals>.gen_random)r  r  r  r6   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r   r1   rg  .)rd  r:   r/  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   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    s    zTestThread.test_multinomialN)rT   rU   rV   r   r  r  r  r  r$   r$   r$   r%   r  a  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   rI   rJ   r>  )r:   rL   argOneargTwoargThreetgtShaperB   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}||krb|tdg}n
|| j}t|j| j qDd S r   )r   rd  r  rX  r  r  r  r  r  r   r   r   r   r:   rL   r  r   rn   r  )r6   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}||krftdg}n| j}|| j|}t|j| j || jd |}t|j| j || j|d }t|j| j qLd S )Nr   r   )r   r  r  rO  rp  rm  r  r  r  r{  r  r  r  r   r   r:   rL   r  r  r   rn   r  )r6   r  r  r  r  r  r$   r$   r%   test_two_arg_funcs  s2          z*TestSingleEltArrayInput.test_two_arg_funcsc                 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   )
r   r  r  r   r  r  r  r   rn   r  )r6   r  r  r  r$   r$   r%   test_three_arg_funcs  s    z,TestSingleEltArrayInput.test_three_arg_funcsN)rT   rU   rV   r   r  r  r  r$   r$   r$   r%   r    s   r  c                 C   sD   t d | \}}}tt |}||ddi}t|jtdk d S )Nr  r   r   r   )r   rP   getattrr   rH   r:   )r&   fnamer  r   rm  r   r$   r$   r%   test_integer_dtype  s
    


r"  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   rP   r   r   r   r   r   r   r   r:   r   r   r   )r&   r!  r  r   rm  r   r   r$   r$   r%   test_integer_repeat  s    



r#  c                	   C   s   t t tjdddgdd W 5 Q R X t t tjddgddd W 5 Q R X t t tjddgddgdd W 5 Q R X d S )NrD   r   gffffff?)r   rD   r   r   )r   r   rA   r   r   r$   r$   r$   r%   test_broadcast_size_error  s    r$  ))r   r   r   r  Znumpyr:   r   Znumpy.testingr   r   r   r   r   r   r   r	   Znumpy.randomr
   r   r   r/  r!   r   rX   r   r"   Zfixturer&   r/   r0   rW   r^   r   r   r   r  r  r  r"  r#  r$  r$   r$   r$   r%   <module>   s|   (


/5W "      q    O/E