XhddlmZddlmZddlmZdZeefdejfdZeefdejfdZ eefdejfd Z eefdejfd Z ed Z eefdejfd Z eefdejfd ZedZeefdejfdZeefdejfdZdS))jit)core)math n_roundsc |jtjkr d}d}d} d} n0tj|jtjkdd}d}d} d } tj|D]} | } | } ||}}t j| ||z |z }t j| ||z |z }tj| |d }tj| |d }tj ||d }tj ||d }||||fS) zV Run `n_rounds` rounds of Philox for state (c0, c1, c2, c3) and key (k0, k1). lyn<l.vlS$lW =z"dtype not supported in philox_impll|~fq l;' :=[ llB;: lW$* R F)sanitize_overflow) dtypetluint32 static_assertuint64 static_rangerumulhimuladd)c0c1c2c3k0k1r PHILOX_KEY_A PHILOX_KEY_BPHILOX_ROUND_APHILOX_ROUND_B_AB_c0_c2s h/var/www/tools.fuzzalab.pt/emblema-extractor/venv/lib/python3.11/site-packages/triton/language/random.py philox_implr$ s&  x29%/ %/ '1'1 RY.0TUUU%7 %7 '9'9 _X & & ? ?  rS [C 2 % * [C 2 % * VAse 4 4 4 VAse 4 4 4 VB  > > > VB  > > > r2r>c tj|}tj|j|tj}tj|}tj|}tj|}tj|}tj|jjdkrTtj }|dz dztj }|dztj }nZtjtj|jjdkdtj}tj dd|}|}||d }||d }||d }||d }t|||||||S) N l@z bitwidth not supported in philox)r)r Tbitcast) r to_tensorrr is_inttor constexprprimitive_bitwidthr fullr$) seedrrrrr int_dtypeseed_hiseed_los r#philoxr6-s <  DTZ&&(())) 7729  D b  B b  B b  B b  B |BH/00B66I BJ*,00;;*$((33 bh&ABBbHJlmmmI '%)444 y$ ' 'B y$ ' 'B y$ ' 'B y$ ' 'B r2r2w B BBr%c2t|||\}}}}|S)aq Given a :code:`seed` scalar and an :code:`offset` block, returns a single block of random :code:`int32`. If you need multiple streams of random numbers, using `randint4x` is likely to be faster than calling `randint` 4 times. :param seed: The seed for generating random numbers. :param offset: The offsets to generate random numbers for. ) randint4x)r2offsetrretrs r#randintr;Hs"T6844LCAq Jr%c|tj}|dz}tj|jjdkr#|dz tj}n|}t ||||||S)aU Given a :code:`seed` scalar and an :code:`offset` block, returns four blocks of random :code:`int32`. This is the maximally efficient entry point to Triton's Philox pseudo-random number generator. :param seed: The seed for generating random numbers. :param offsets: The offsets to generate random numbers for. r)r')r.r r r/r r0r6)r2r9r offset_lo_0 offset_his r#r8r8Xsr ")$$I QB |FL344r99r\%%bi00  $ 9b"h ? ??r%c2tj|jtjks'tj|jtjkr$|tjd}d}ntjtj|jtjkp&tj|jtjk|tjd}d}tj |dk| dz |}||zS)zo Numerically stable function to convert a random uint into a random float uniformly sampled in [0, 1). Tr*g=g;r)r) r r/r r int32r.rrint64where)xscales r#uint_to_uniform_floatrF~s |AGry())'R\!'RX:M-N-N' DD4D ( (  ag&:;;`r|AGWYW_L_?`?`aaa DD4D ( (& QQ""A u9r%cBt|||}t|S)z Given a :code:`seed` scalar and an :code:`offset` block, returns a block of random :code:`float32` in :math:`U(0, 1)`. :param seed: The seed for generating random numbers. :param offsets: The offsets to generate random numbers for. )r;rF)r2r9rsources r#randrIs#T68 , ,F  ( ((r%ct|||\}}}}t|}t|}t|} t|} ||| | fS)a Given a :code:`seed` scalar and an :code:`offsets` block, returns 4 blocks of random :code:`float32` in :math:`U(0, 1)`. :param seed: The seed for generating random numbers. :param offsets: The offsets to generate random numbers for. )r8rF) r2offsetsri1i2i3i4u1u2u3u4s r#rand4xrTs_tWh77NBB r " "B r " "B r " "B r " "B r2r>r%ctjd|}d|z}tjdtj|z}|tj|z|tj|zfS)zBox-Muller transformgHz>g-DT!@g)r maximumrsqrtlogcossin)rPrQthrs r#pair_uniform_to_normalr]s` FB  B R B $"%&&A tx|| Q"- --r%ct|||\}}}}t|}t|}t||\}}|S)a Given a :code:`seed` scalar and an :code:`offset` block, returns a block of random :code:`float32` in :math:`\mathcal{N}(0, 1)`. :param seed: The seed for generating random numbers. :param offsets: The offsets to generate random numbers for. )r8rFr]) r2r9rrLrMrrPrQn1s r#randnr`sOT6844LBAq r " "B r " "B "2r * *EB Ir%ct|||\}}}}t||\}}t||\} } ||| | fS)a  Given a :code:`seed` scalar and an :code:`offset` block, returns 4 blocks of random :code:`float32` in :math:`\mathcal{N}(0, 1)`. :param seed: The seed for generating random numbers. :param offsets: The offsets to generate random numbers for. )rTr]) r2r9rrPrQrRrSr_n2n3n4s r#randn4xresQD&(33NBB #B + +FB #B + +FB r2r>r%N) runtime.jitrrr rN_ROUNDS_DEFAULTr/r$r6r;r8rFrIrTr]r`rer%r#rjsAQ",@:JCC2<CCCC43C  BL    5E@@bl@@@@J&0@ ) ) ) ) ) )3C  BL    *...1A  ",    3C  BL      r%