U
|
P ±dù ã @ sD d dl Zd dlZd dlmZ d dlmZ e Zejj Z
|
ejj Z dd Z eje ddZdZd d
|
Zej e ¡Zd d Ze Ze Zejefkstejejks¢teed dZeedde d eed dZ ee dde d ej j Z ej j Z!ejdddd Z"ee"dde!j# ejdddd Z$e$ddde!j# dS )é N)ÚPCG64)Útimeitc C sâ t | ¡}t| d d D ]Â}dt| d }dt| d }|| || }|dks^|dkrdt| d }dt| d }|| || }qNt dt |¡ | ¡}|| |d| < d| d | k r|| |d| d < q|S )Né é g @g ð?g g À)ÚnpÚemptyÚrangeÚnext_dÚsqrtÚlog)ÚnÚstateÚoutÚiÚx1Zx2Úr2Úf© r ú]d:\z\workplace\vscode\pyvenv\venv\Lib\site-packages\numpy/random/_examples/numba/extending.pyÚnormals s
|
r T)Znopythoni' c C s
|
t ttS )N)Únormalsjr Ú
|
state_addrr r r r Ú numbacall s r c C s t jtdS )N)Úsize)ÚrgÚnormalr r r r r Ú numpycall% s r iè )Únumberz.2fz
|
secs for z% PCG64 (Numba/PCG64) gaussian randomsz% PCG64 (NumPy/PCG64) gaussian randomsc C sr || }}||d? O }||d? O }||d? O }||d? O }||d? O }t ||@ }||krjt ||@ }qT| | S )Nr r é é é )Únext_u32)ÚlbÚubr ÚmaskÚdeltaÚvalr r r Ú bounded_uint8 s r' iC iC$ c C s2 t j|t jd}t|D ]}t| ||||< qd S )N)Zdtype)r r Zuint32r r' )r" r# r r r r r r r Ú bounded_uintsK s r( i )%Únumpyr ZnumbaÚnbZ numpy.randomr r Zbit_genZcffiZ next_doubler Z state_addressr r Zjitr r r ÚrandomÚ Generatorr r Zr1r ÚshapeÚAssertionErrorÚt1ÚprintÚt2ÚctypesZ next_uint32r! r Z ctypes_stater' Úvaluer( r r r r Ú<module> s8
|
|
|
|
|
