o cZh@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}tj| |d d }tj| |d d }tj ||d d }tj ||d d }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 F)Zsanitize_overflow) dtypetluint32 static_assertuint64Z static_rangerZumulhimuladd)c0c1c2c3Zk0Zk1rZ PHILOX_KEY_AZ PHILOX_KEY_BZPHILOX_ROUND_AZPHILOX_ROUND_B_ABZ_c0Z_c2rE/var/www/auris/lib/python3.10/site-packages/triton/language/random.py philox_impl s*   rc Cs t|}t|j|tj}t|}t|}t|}t|}t|jjdkrFtj }|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)rTZbitcast) r Z to_tensorr rZis_inttor constexprZprimitive_bitwidthr fullr) seedrrrrrZ int_dtypeZseed_hiZseed_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 offsetrretrrrrrandintFs r%cCs|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 r#rZ_0rrrr"Vs r"cCst|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). Trg=g;rr) r rrr Zint32rr r Zint64where)xscalerrruint_to_uniform_floatts$*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. )r%r))r r#rsourcerrrrands 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. )r"r)) r offsetsri1i2Zi3i4u1u2u3u4rrrrand4xs   r4cCsDtd|}d|}tdt|}|t||t|fS)zBox-Muller transformgHz>g-DT!@g)r maximumrsqrtlogcossin)r0r1thrrrrpair_uniform_to_normals r<c 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. )r"r)r<) r r#rr-r.rr0r1n1rrrrandns  r>c 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. )r4r<) r r#rr0r1r2r3r=Zn2Zn3Zn4rrrrandn4xs  r?N)Z runtime.jitrrr rZN_ROUNDS_DEFAULTrrr!r%r"r)r+r4r<r>r?rrrrs0