a kh@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}tj| |d d }tj| |d d }tj ||d d }tj ||d d }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 F)Zsanitize_overflow) dtypetluint32 static_assertuint64Z static_rangerZumulhimuladd)c0c1c2c3Zk0Zk1rZ PHILOX_KEY_AZ PHILOX_KEY_BZPHILOX_ROUND_AZPHILOX_ROUND_B_ABZ_c0Z_c2rD/var/www/auris/lib/python3.9/site-packages/triton/language/random.py philox_impl s*  rc Cs t|}t|j|tj}t|}t|}t|}t|}t|jjdkrtj }|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_tensorr rZis_inttor constexprprimitive_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!offsetrretrrrrrandintHs r&cCsN|tj}|d}t|jjdkr8|d?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. rr)rr r rrrr")r!r$rZ offset_loZ_0Z offset_hirrrr#Xs  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 int32rr r int64where)xZscalerrruint_to_uniform_float~s$*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 r6cCsDtd|}d|}tdt|}|t||t|fS)zBox-Muller transformgHz>g-DT!@g)r maximumrsqrtlogcossin)r2r3thrrrrpair_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/r0rr2r3n1rrrrandns 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. )r6r>) r!r$rr2r3r4r5r?Zn2Zn3Zn4rrrrandn4xs rAN)Z runtime.jitrrr rZN_ROUNDS_DEFAULTrrr"r&r#r+r-r6r>r@rArrrrs.    %