o 灛iP@sddlmZddlmZddlmZdZeefdejfddZeefdejfd d Z eefdejfd d Z eefdejfd dZ eddZ eefdejfddZ eefdejfddZeddZeefdejfddZeefdejfddZdS))jit)core)math n_roundscCs|jtjkrd}d}d} d} nt|jtjkdd}d}d} d } t|D]/} | } | } ||}}t| ||A|A}t| ||A|A}| |}| |}||}||}q&||||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_c2rX/sda-disk/www/egybert/egybert_env/lib/python3.10/site-packages/triton/language/random.py philox_impl s*    r c Cst|}t|}t|}t|}t|}|tj}t|jjdkr>tj}|d?d@tj}|d@tj}ntt|jjdkdtj}tj dd|d}|}|j|dd }|j|dd }|j|dd }|j|dd }t |||||||S) N l@z bitwidth not supported in philox)r)rTbitcast) r to_tensortor constexprrprimitive_bitwidthr r fullr ) seedrrrrr int_dtypeseed_hiseed_lorrrphilox-s&      r/cCst|||\}}}}|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+offsetrretrrrrrandintEs r3cCs|d}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+r1r_0rrrr0Us r0cCst|jtjkst|jtjkr|jtjdd}d}ntt|jtjkp0t|jtjk|jtjdd}d}t |dk| d|}||S)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(rr int32r'r r int64where)xscalerrruint_to_uniform_floatss$*r:cCst|||}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. )r3r:)r+r1rsourcerrrrands r<c Cs@t|||\}}}}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. )r0r:) r+offsetsri1i2i3i4u1u2u3u4rrrrand4xs   rFcCsDtd|}d|}tdt|}|t||t|fS)zBox-Muller transformgHz>g-DT!@g)r maximumrsqrtlogcossin)rBrCthrrrrpair_uniform_to_normals rNc Cs6t|||\}}}}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. )r0r:rN) r+r1rr>r?rrBrCn1rrrrandns  rPc Cs<t|||\}}}}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. )rFrN) r+r1rrBrCrDrErOn2n3n4rrrrandn4xs  rTN) runtime.jitrrr rN_ROUNDS_DEFAULTr(r r/r3r0r:r<rFrNrPrTrrrrs0