çgPddlmZddlmZddlmZdZeefdejfdZeefdejfdZ 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} | } | } ||}}t j| ||z |z }t j| ||z |z }| |z}| |z}||z}||z}U||||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 )dtypetluint32 static_assertuint64 static_rangerumulhi)c0c1c2c3k0k1r PHILOX_KEY_A PHILOX_KEY_BPHILOX_ROUND_APHILOX_ROUND_B_AB_c0_c2s R/var/www/html/ai-engine/env/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 % * W W ,  ,  r2r>c tj|}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_tensortor constexprr primitive_bitwidthr r fullr!) seedrrrrr int_dtypeseed_hiseed_los r philoxr2-ss <  D b  B b  B b  B b  B 7729  D |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)r.offsetrretrs r randintr7Es"T6844LCAq Jr"c4|dz}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&)r2)r.r5r_0s r r4r4Us% !B $BH 5 55r"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*r rint64where)xscales r uint_to_uniform_floatr@ss |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. )r7r@)r.r5rsources r randrCs#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. )r4r@) r.offsetsri1i2i3i4u1u2u3u4s r rand4xrNs_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)rJrKthrs r pair_uniform_to_normalrWs` 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. )r4r@rW) r.r5rrFrGrrJrKn1s r randnrZsOT6844LBAq 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. )rNrW) r.r5rrJrKrLrMrYn2n3n4s r randn4xr_sQD&(33NBB #B + +FB #B + +FB r2r>r"N) runtime.jitrrr rN_ROUNDS_DEFAULTr+r!r2r7r4r@rCrNrWrZr_r"r rdsAQ",@:JCC2<CCCC.3C  BL    5E 6 6bl 6 6 6 6:&0@ ) ) ) ) ) )3C  BL    *...1A  ",    3C  BL      r"