U hP@sddlmZddlmZddlmZdZeefejdddZeefejdd d Z eefejdd d Z eefejdd dZ eddZ eefejdddZ eefejdddZeddZeefejdddZeefejdddZdS))jit)core)math )n_roundscCs|jtjkrd}d}d} d} n$t|jtjkdd}d}d} d } t|D]^} | } | } ||}}t| ||A|A}t| ||A|A}| |}| |}||}||}qL||||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_assertuint64Z static_rangerZumulhi)c0c1c2c3Zk0Zk1rZ PHILOX_KEY_AZ PHILOX_KEY_BZPHILOX_ROUND_AZPHILOX_ROUND_B_ABZ_c0Z_c2rH/var/www/html/venv/lib/python3.8/site-packages/triton/language/random.py philox_impl s*   rc Cst|}t|}t|}t|}t|}|tj}t|jjdkr|tj}|d?d@tj}|d@tj}n4tt|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_tensortor constexprrZprimitive_bitwidthr r fullr) seedr rrrrZ int_dtypeZseed_hiZseed_lorrrphilox-s&      rcCst|||\}}}}|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)roffsetrretrrrrrandintEs 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)rr!rZ_0rrrr Us r cCst|jtjks$t|jtjkr:|jtjdd}d}n>tt|jtjkp`t|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_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. )r#r')rr!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') roffsetsri1i2Zi3i4u1u2u3u4rrrrand4xs r2cCsDtd|}d|}tdt|}|t||t|fS)zBox-Muller transformgHz>g-DT!@g)r maximumrsqrtlogcossin)r.r/thrrrrpair_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:) rr!rr+r,rr.r/n1rrrrandns 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. )r2r:) rr!rr.r/r0r1r;Zn2Zn3Zn4rrrrandn4xs r=N)Z runtime.jitrrr rZN_ROUNDS_DEFAULTrrrr#r r'r)r2r:r<r=rrrrs.