a
    ³ÐOc-@  ã                   @   sÀ  d dl Z d dlmZ d dlZd dlZd dlmZmZm	Z	m
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 zd dlZdZW n ey¢   dZY n0 zd dlZdZW n eyÊ   dZY n0 ejjdkrÜdZe j  e j !e"¡¡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/G d!d"„ d"ƒZ0G d#d$„ d$e0ƒZ1G d%d&„ d&e0ƒZ2G d'd(„ d(e0ƒZ3G d)d*„ d*e0ƒZ4G d+d,„ d,e0ƒZ5G d-d.„ d.ƒZ6dS )/é    N)Újoin)Úassert_equalÚassert_allcloseÚassert_array_equalÚassert_raises)	Ú	GeneratorÚMT19937ÚPCG64Ú	PCG64DXSMÚPhiloxÚRandomStateÚSeedSequenceÚSFC64Údefault_rng)Ú	interfaceFTé   c                 C   sh   | D ]^}t | | tƒr*t| | || ƒ qt | | tjƒrNt| | || ƒ q| | || ksJ ‚qd S ©N)Ú
isinstanceÚdictÚassert_state_equalÚnpZndarrayr   )ÚactualÚtargetÚkey© r   úU/var/www/html/gsl.tools/lib/python3.9/site-packages/numpy/random/tests/test_direct.pyr   '   s    r   c                 C   s   | t  d¡? d  t j¡S )Né   g      p>)r   Úuint32ZastypeÚfloat32)Úur   r   r   Úuint32_to_float321   s    r    c                 C   s\   t  | ¡} t j| t  d¡? t jd�}t  d¡}t j| |@ t jd�}t  ||g¡ ¡ }t|ƒS )Né    ©Údtypeì   ÿÿ )r   Úuint64Úarrayr   Zcolumn_stackZravelr    )ÚxÚupperÚlowerZjoinedr   r   r   Úuniform32_from_uint645   s    

r*   c                 C   s0   t  | ¡t  d¡? } t  | t  d¡@ ¡} t| ƒS )Né   r$   )r   r%   r   r    ©r'   r   r   r   Úuniform32_from_uint53>   s    r-   c                 C   s   t | ƒS r   )r    r,   r   r   r   Úuniform32_from_uint32D   s    r.   c                 C   s8   |dkrt | ƒS |dkr t| ƒS |dkr0t| ƒS t‚d S )Né@   é5   r!   )r*   r-   r.   ÚNotImplementedError©r'   Úbitsr   r   r   Úuniform32_from_uintH   s    r4   c                 C   s$   |dv rt | ƒS |dkr t| ƒS d S )N)r/   é?   r0   r!   )Úuniform_from_uint64Úuniform_from_uint32r2   r   r   r   Úuniform_from_uintS   s    r8   c                 C   s   | t  d¡? d S )Né   g       <)r   r%   r,   r   r   r   r6   Z   s    r6   c                 C   s`   t  t| ƒd ¡}tdt| ƒdƒD ]8}| | d? }| |d  d? }|d | d ||d < q"|S )Né   r   é   r   é   g      �Ag      @C)r   ÚemptyÚlenÚrange)r'   ÚoutÚiÚaÚbr   r   r   r7   ^   s    r7   c                 C   s   |   tj¡d S )Nç      ð?)Úviewr   Údoubler,   r   r   r   Úuniform_from_dsfmtg   s    rG   c           
      C   sÜ   |dv rt | ƒ}n|dkr$t| ƒ}nt| ƒ}g }d}d }}t|ƒ|k rÐd}|dks\|dkršd||  d }d||d   d }|| ||  }|d7 }qLt d	t |¡ | ¡}	| |	| ¡ | |	| ¡ q<|d |… S )
N)r/   r5   r!   r   ç        r:   rD   g       @r   g       À)r6   r7   rG   r>   r   ÚsqrtÚlogÚappend)
r'   Únr3   ZdoublesÚgaussÚlocÚx1Zx2Zr2Úfr   r   r   Úgauss_from_uintk   s&    


rQ   c                  C   sž   ddl m} m}m} ttdƒddd�}| d¡ tf i |j¤Ž}t|j|jƒ t|j	|j	ƒ t
t| ƒ t
t|ƒ |ƒ }t
t|jdƒ t| d¡ƒdksšJ ‚d S )Nr   )ÚISeedSequenceÚISpawnableSeedSequenceÚSeedlessSeedSequenceé
   )r   r:   r<   )Z	spawn_keyZ	pool_size)Znumpy.random.bit_generatorrR   rS   rT   r   r?   ÚspawnÚstater   Zn_children_spawnedr   Ú	TypeErrorr1   Zgenerate_stater>   )rR   rS   rT   Ús1Ús2Údummyr   r   r   Útest_seedsequence„   s    


r\   c                   @   sÒ   e Zd ZejZi  ZZedd„ ƒZ	e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ejjedd �d!d"„ ƒZejjed#d �d$d%„ ƒZd&d'„ Zd(S ))ÚBasec                 C   s*   t | _d| _tj| _t| _g | _g | _	d S )Nr/   )
r	   Úbit_generatorr3   r   r%   r#   rX   Úseed_error_typeÚinvalid_init_typesÚinvalid_init_values©Úclsr   r   r   Úsetup_class›   s    zBase.setup_classc                 C   s–   t |ƒ�z}| ¡ }| d¡}dd„ |dd … D ƒ}g }|D ]"}| t| d¡d  ¡ dƒ¡ q:|tj|| jd�dœW  d   ƒ S 1 sˆ0    Y  d S )	Nú,c                 S   s   g | ]}t | ¡ d ƒ‘qS )r   )ÚintÚstrip)Ú.0Úsr   r   r   Ú
<listcomp>©   ó    z"Base._read_csv.<locals>.<listcomp>r   éÿÿÿÿr   r"   )ÚseedÚdata)	ÚopenÚreadlineÚsplitrK   rf   rg   r   r&   r#   )rc   ÚfilenameÚcsvrm   rn   Úliner   r   r   Ú	_read_csv¤   s    

 zBase._read_csvc                 C   s„   | j | jd Ž }| d¡}t|| jd ƒ | j | jd Ž }| ¡ }t|| jd d ƒ | j | jd Ž }| d¡}t|| jd ƒ d S )Nrm   éè  rn   r   )r^   Údata1Ú
random_rawr   Údata2©Úselfr^   Zuintsr   r   r   Útest_raw¯   s    

zBase.test_rawc                 C   sF   | j | jd Ž }|jdd�}|d u s(J ‚|jddd�}|d u sBJ ‚d S )Nrm   F)Úoutputrv   )r^   rw   rx   rz   r   r   r   Útest_random_raw¼   s
    zBase.test_random_rawc                 C   sx   d}t | j| jd Ž ƒ}| |¡}t|t| jd || jƒƒ t | j| jd Ž ƒ}| d¡}t|t| jd || jƒƒ d S )Né   rm   rn   )r   r^   rw   Ústandard_normalr   rQ   r3   ry   )r{   rL   ÚrsrM   r   r   r   Útest_gauss_invÃ   s    
ÿ
ÿzBase.test_gauss_invc                 C   sœ   t | j| jd Ž ƒ}t| jd | jƒ}| t|ƒ¡}t||ƒ t|j	t
jƒ t | j| jd Ž ƒ}t| jd | jƒ}| t|ƒ¡}t||ƒ t|j	t
jƒ d S )Nrm   rn   )r   r^   rw   r8   r3   Úrandomr>   r   r   r#   r   Úfloat64ry   ©r{   r�   ÚvalsZuniformsr   r   r   Útest_uniform_doubleÏ   s    

zBase.test_uniform_doublec                 C   s¨   t | j| jd Ž ƒ}t| jd | jƒ}|jt|ƒtjd�}t	||ƒ t
|jtjƒ t | j| jd Ž ƒ}t| jd | jƒ}|jt|ƒtjd�}t	||ƒ t
|jtjƒ d S )Nrm   rn   r"   )r   r^   rw   r4   r3   rƒ   r>   r   r   r   r   r#   ry   r…   r   r   r   Útest_uniform_floatÜ   s    

zBase.test_uniform_floatc                 C   sL   t | j| jd Ž ƒ}dt|ƒv s$J ‚t|ƒd› ¡  dd¡t|ƒv sHJ ‚d S ©Nrm   r   z#xÚXr'   )r   r^   rw   ÚreprÚidr(   Úreplace©r{   r�   r   r   r   Ú	test_repré   s    zBase.test_reprc                 C   sd   t | j| jd Ž ƒ}dt|ƒv s$J ‚t| jjƒt|ƒv s<J ‚t|ƒd› ¡  dd¡t|ƒvs`J ‚d S r‰   )r   r^   rw   ÚstrÚ__name__rŒ   r(   r�   rŽ   r   r   r   Útest_strî   s    zBase.test_strc           	      C   s–   dd l }| j| jd Ž }|j}| |¡}| |¡}|j}tt|ƒ d¡t|ƒ d¡ƒ ||usbJ ‚t	||ƒ t
dƒ}| | |¡¡}t|j|jƒ d S )Nr   rm   rv   éd   )Úpickler^   rw   rW   ÚdumpsÚloadsr   r   r€   r   r   r   )	r{   r”   r^   rW   Z
bitgen_pklZreloadedZreloaded_stateÚssZaar   r   r   Útest_pickleô   s    

ÿ
zBase.test_picklec                 C   sF   | j | jd Ž }t t¡� dh|_W d   ƒ n1 s80    Y  d S )Nrm   Ú1)r^   rw   ÚpytestÚraisesrX   rW   ©r{   r^   r   r   r   Útest_invalid_state_type  s    zBase.test_invalid_state_typec                 C   sR   | j | jd Ž }|j}d|d< t t¡� ||_W d   ƒ n1 sD0    Y  d S )Nrm   ZotherBitGeneratorr^   )r^   rw   rW   rš   r›   Ú
ValueError)r{   r^   rW   r   r   r   Útest_invalid_state_value
  s
    zBase.test_invalid_state_valuec              	   C   sH   | j }| jD ]6}t t¡� ||Ž  W d   ƒ q1 s80    Y  qd S r   )r^   r`   rš   r›   rX   ©r{   r^   Ústr   r   r   Útest_invalid_init_type  s    
zBase.test_invalid_init_typec              	   C   sL   | j }| jD ]:}t ttf¡� ||Ž  W d   ƒ q1 s<0    Y  qd S r   )r^   ra   rš   r›   rž   ÚOverflowErrorr    r   r   r   Útest_invalid_init_values  s    
zBase.test_invalid_init_valuesc                 C   s`   | j | jd Ž }| d¡ | dd¡ t t¡� | dd¡ W d   ƒ n1 sR0    Y  d S )Nrm   r   rF   Úint32)r^   rw   Z
_benchmarkrš   r›   rž   rœ   r   r   r   Útest_benchmark  s
    
zBase.test_benchmarkzcffi not available)Úreasonc                 C   s:   | j | jd Ž }|j}t|tƒs$J ‚|j}||u s6J ‚d S ©Nrm   )r^   rw   Úcffir   r   )r{   r^   Zcffi_interfaceZother_cffi_interfacer   r   r   Ú	test_cffi$  s
    zBase.test_cffizctypes not availablec                 C   s:   | j | jd Ž }|j}t|tƒs$J ‚|j}||u s6J ‚d S r¨   )r^   rw   Úctypesr   r   )r{   r^   Zctypes_interfaceZother_ctypes_interfacer   r   r   Útest_ctypes,  s
    zBase.test_ctypesc                 C   s,   | j | jd Ž }|j}| ¡ }t||ƒ d S r¨   )r^   rw   rW   Ú__getstate__r   )r{   r^   rW   Z	alt_stater   r   r   Útest_getstate4  s    zBase.test_getstateN) r‘   Ú
__module__Ú__qualname__r   r%   r#   ry   rw   Úclassmethodrd   ru   r|   r~   r‚   r‡   rˆ   r�   r’   r˜   r�   rŸ   r¢   r¤   r¦   rš   ÚmarkZskipifÚMISSING_CFFIrª   ÚMISSING_CTYPESr¬   r®   r   r   r   r   r]   —   s0   




r]   c                   @   s    e Zd Zedd„ ƒZdd„ ZdS )Ú
TestPhiloxc                 C   sb   t | _d| _tj| _|  ttdƒ¡| _	|  ttdƒ¡| _
t| _g | _ddd d dd d fg| _d S )	Nr/   z./data/philox-testset-1.csvz./data/philox-testset-2.csv)r   Nr   ©rl   r:   i  r   )r   r^   r3   r   r%   r#   ru   r   Úpwdrw   ry   rX   r_   r`   ra   rb   r   r   r   rd   <  s    ÿÿzTestPhilox.setup_classc                 C   sF   | j | jd Ž }|j}| j |d d |d d d�}t|j|jƒ d S )Nrm   rW   Úcounterr   )r¸   r   )r^   rw   rW   r   )r{   r^   rW   Zkeyedr   r   r   Útest_set_keyI  s    
ÿzTestPhilox.test_set_keyN)r‘   r¯   r°   r±   rd   r¹   r   r   r   r   rµ   ;  s   
rµ   c                   @   s(   e Zd Zedd„ ƒZdd„ Zdd„ ZdS )Ú	TestPCG64c                 C   s^   t | _d| _tj| _|  ttdƒ¡| _	|  ttdƒ¡| _
ttf| _dd gfdg| _dg| _d S )Nr/   z./data/pcg64-testset-1.csvz./data/pcg64-testset-2.csv©gš™™™™™	@©r   Nr¶   )r	   r^   r3   r   r%   r#   ru   r   r·   rw   ry   rž   rX   r_   r`   ra   rb   r   r   r   rd   R  s    
zTestPCG64.setup_classc                 C   s¢   t | j| jd Ž ƒ}|jj}d}|j |¡ | d¡}||j_|j dd | ¡ | d¡}||j_|j ddd  | ¡ | d¡}||ks’J ‚||ksžJ ‚d S ©Nrm   l÷ÿÿÿ        PAOi_n¼ž rU   r:   é€   ©r   r^   rw   rW   ÚadvanceÚintegers©r{   r�   rW   ÚstepZval_negZval_posZval_bigr   r   r   Útest_advance_symmetry]  s    


zTestPCG64.test_advance_symmetryc                 C   sl   t |  d¡ƒ}|j}|jd }d}|d |ks2J ‚| tdd„ dD ƒƒ¡ |jd }d}|d |kshJ ‚d S )Nì   ï%”@
 rW   ì	   Ír%4=KÆv±uÅ1¹/Ø c                 s   s   | ]}d | V  qdS ©r:   Nr   ©rh   rA   r   r   r   Ú	<genexpr>r  rk   z/TestPCG64.test_advange_large.<locals>.<genexpr>©é`   r/   r!   r+   r   é   r:   r   l	   ÿ@H}.|z³KLf�be ©r   r^   rW   rÀ   Úsum©r{   r�   ZpcgrW   Zinitial_stateZadvanced_stater   r   r   Útest_advange_largel  s    

zTestPCG64.test_advange_largeN©r‘   r¯   r°   r±   rd   rÄ   rÐ   r   r   r   r   rº   Q  s   

rº   c                   @   s(   e Zd Zedd„ ƒZdd„ Zdd„ ZdS )ÚTestPCG64DXSMc                 C   s^   t | _d| _tj| _|  ttdƒ¡| _	|  ttdƒ¡| _
ttf| _dd gfdg| _dg| _d S )Nr/   z./data/pcg64dxsm-testset-1.csvz./data/pcg64dxsm-testset-2.csvr»   r¼   r¶   )r
   r^   r3   r   r%   r#   ru   r   r·   rw   ry   rž   rX   r_   r`   ra   rb   r   r   r   rd   y  s    
zTestPCG64DXSM.setup_classc                 C   s¢   t | j| jd Ž ƒ}|jj}d}|j |¡ | d¡}||j_|j dd | ¡ | d¡}||j_|j ddd  | ¡ | d¡}||ks’J ‚||ksžJ ‚d S r½   r¿   rÂ   r   r   r   rÄ   „  s    


z#TestPCG64DXSM.test_advance_symmetryc                 C   sl   t |  d¡ƒ}|j}|j}d}|d d |ks2J ‚| tdd„ dD ƒƒ¡ |jd }d}|d |kshJ ‚d S )NrÅ   rÆ   rW   c                 s   s   | ]}d | V  qdS rÇ   r   rÈ   r   r   r   rÉ   ™  rk   z3TestPCG64DXSM.test_advange_large.<locals>.<genexpr>rÊ   l	   oF‚	œ[UmO1X·9}Ð rÍ   rÏ   r   r   r   rÐ   “  s    
z TestPCG64DXSM.test_advange_largeNrÑ   r   r   r   r   rÒ   x  s   

rÒ   c                   @   s(   e Zd Zedd„ ƒZdd„ Zdd„ ZdS )ÚTestMT19937c                 C   sP   t | _d| _tj| _|  ttdƒ¡| _	|  ttdƒ¡| _
t| _g | _dg| _d S )Nr!   z./data/mt19937-testset-1.csvz./data/mt19937-testset-2.csvr¶   )r   r^   r3   r   r   r#   ru   r   r·   rw   ry   rž   r_   r`   ra   rb   r   r   r   rd      s    zTestMT19937.setup_classc                 C   s”   t t| jt tjg¡ƒ t t| jt tj g¡ƒ t t| jt tjtj g¡ƒ t t| jt dtjg¡ƒ t t| jtjgƒ t t| jdtjgƒ d S )Nr   )r   rX   r^   r   r&   Úpi)r{   r   r   r   Útest_seed_float_array«  s    z!TestMT19937.test_seed_float_arrayc                 C   sˆ   t | j| jd Ž ƒ}|j}|j}| d¡}|d |d d |d d f}||_| d¡}t||ƒ |d }||_| d¡}t||ƒ d S )Nrm   i   r^   rW   r   Úpos)r   rH   )r   r^   rw   rW   rÁ   r   )r{   r�   r^   rW   ZdesiredÚtupr   r   r   r   Útest_state_tuple³  s    

ÿ


zTestMT19937.test_state_tupleN)r‘   r¯   r°   r±   rd   rÕ   rØ   r   r   r   r   rÓ   Ÿ  s   

rÓ   c                   @   s   e Zd Zedd„ ƒZdS )Ú	TestSFC64c                 C   s^   t | _d| _tj| _|  ttdƒ¡| _	|  ttdƒ¡| _
ttf| _dd gfdg| _dg| _d S )Nr/   z./data/sfc64-testset-1.csvz./data/sfc64-testset-2.csvr»   r¼   r¶   )r   r^   r3   r   r%   r#   ru   r   r·   rw   ry   rž   rX   r_   r`   ra   rb   r   r   r   rd   Ä  s    ÿÿ
zTestSFC64.setup_classN)r‘   r¯   r°   r±   rd   r   r   r   r   rÙ   Ã  s   rÙ   c                   @   s   e Zd Zdd„ Zdd„ ZdS )ÚTestDefaultRNGc                 C   s4   dddddgffD ]}t |Ž }t|jtƒsJ ‚qd S )Nr   r   )éÒ  rÛ   i.  )r   r   r^   r	   )r{   ÚargsÚrgr   r   r   Ú	test_seedÓ  s    zTestDefaultRNG.test_seedc                 C   sB   t ƒ }t|ƒ}|j|u sJ ‚t|ƒ}||u s0J ‚|j|u s>J ‚d S r   )r   r   r^   )r{   ÚbgrÝ   Zrg2r   r   r   Útest_passthroughØ  s    zTestDefaultRNG.test_passthroughN)r‘   r¯   r°   rÞ   rà   r   r   r   r   rÚ   Ò  s   rÚ   )7ÚosÚos.pathr   ÚsysÚnumpyr   Znumpy.testingr   r   r   r   rš   Znumpy.randomr   r   r	   r
   r   r   r   r   r   Znumpy.random._commonr   r©   r³   ÚImportErrorr«   r´   ÚflagsÚoptimizeÚpathÚdirnameÚabspathÚ__file__r·   r   r    r*   r-   r.   r4   r8   r6   r7   rG   rQ   r\   r]   rµ   rº   rÒ   rÓ   rÙ   rÚ   r   r   r   r   Ú<module>   sP   ,


		 %''$