U ,d'@sddlZddlmZmZmZmZmZmZmZm Z m Z ddl Z e j de jfde jfgddZe eZejZZe eeddd Ze eedd d Ze eedd d Ze eedddZe eedddZe eedddZe eedddZe eedddZe dejZe dejZe eedddZe eedddZe eedddZ d#dd Z!d$d!d"Z"dS)%N) configcudafloat32float64uint32int64uint64 from_dtypejits0s1T)Zalign)ZforceobjZloopliftcCsxt|}t|}|td}||td?Atd}||td?Atd}||td?A}|||d<|||d<d S) aUse SplitMix64 to generate an xoroshiro128p state from 64-bit seed. This ensures that manually set small seeds don't result in a predictable initial sequence from the random number generator. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: uint64 :param index: offset in states to update :type seed: int64 :param seed: seed value to use when initializing state l|~fq le9z lb&&& r r N)rrr)statesindexseedzr5/tmp/pip-unpacked-wheel-eu7e0c37/numba/cuda/random.pyinit_xoroshiro128p_state-s  rcCs(t|}t|}||>|td|?BS)zLeft rotate x by k bits.@)rr)xkrrrrotlGsrcCsxt|}||d}||d}||}||N}tt|td|A|td>A||d<tt|td||d<|S)zReturn the next random uint64 and advance the RNG in states[index]. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: int64 :param index: offset in states to update :rtype: uint64 r r 7$)rrrr)rrr r resultrrrxoroshiro128p_nextOs   *rcCst|}tdtdf}td}td}tdD]T}tdD]F}||tdt|>@rx|||dN}|||dN}t||qS)z7Convert uint64 to float64 value in the range [0.0, 1.0) r!5)rrrrrrruint64_to_unit_float64sr)cCst|}tt|S)z7Convert uint64 to float32 value in the range [0.0, 1.0))rrr)r(rrruint64_to_unit_float32sr*cCst|}tt||S)zReturn a float32 in range [0.0, 1.0) and advance ``states[index]``. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: int64 :param index: offset in states to update :rtype: float32 )rr*rrrrrrxoroshiro128p_uniform_float32s r,cCst|}tt||S)zReturn a float64 in range [0.0, 1.0) and advance ``states[index]``. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: int64 :param index: offset in states to update :rtype: float64 )rr)rr+rrrxoroshiro128p_uniform_float64s r-r cCsHt|}t||}t||}ttd t|tt|}|S)aReturn a normally distributed float32 and advance ``states[index]``. The return value is drawn from a Gaussian of mean=0 and sigma=1 using the Box-Muller transform. This advances the RNG sequence by two steps. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: int64 :param index: offset in states to update :rtype: float32 @)rr,mathsqrtrlogcosTWO_PI_FLOAT32rru1u2Zz0rrrxoroshiro128p_normal_float32s   (r7cCsHt|}t||}t||}ttd t|tt|}|S)aReturn a normally distributed float32 and advance ``states[index]``. The return value is drawn from a Gaussian of mean=0 and sigma=1 using the Box-Muller transform. This advances the RNG sequence by two steps. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: int64 :param index: offset in states to update :rtype: float64 r.)rr,r/r0rr1r2TWO_PI_FLOAT64r4rrrxoroshiro128p_normal_float64s   (r9cCst|jd}t|}t|}|dkrpt|d|t|D]}t|dq6td|D]}||d||<t||qPdS)Nrr!)shaperrr"r%)rrsubsequence_startn_r#rrrinit_xoroshiro128p_states_cpus    r>cCs.tj|jtd}t||||j||ddS)aInitialize RNG states on the GPU for parallel generators. This initializes the RNG states so that each state in the array corresponds subsequences in the separated by 2**64 steps from each other in the main sequence. Therefore, as long no CUDA thread requests more than 2**64 random numbers, all of the RNG states produced by this function are guaranteed to be independent. The subsequence_start parameter can be used to advance the first RNG state by a multiple of 2**64 steps. :type states: 1D DeviceNDArray, dtype=xoroshiro128p_dtype :param states: array of RNG states :type seed: uint64 :param seed: starting seed for list of generators )r:dtype)streamN)npemptyr:xoroshiro128p_dtyper>Zcopy_to_device)rrr;r@Z states_cpurrrinit_xoroshiro128p_statess rDcCs"tj|t|d}t|||||S)a5Returns a new device array initialized for n random number generators. This initializes the RNG states so that each state in the array corresponds subsequences in the separated by 2**64 steps from each other in the main sequence. Therefore, as long no CUDA thread requests more than 2**64 random numbers, all of the RNG states produced by this function are guaranteed to be independent. The subsequence_start parameter can be used to advance the first RNG state by a multiple of 2**64 steps. :type n: int :param n: number of RNG states to create :type seed: uint64 :param seed: starting seed for list of generators :type subsequence_start: uint64 :param subsequence_start: :type stream: CUDA stream :param stream: stream to run initialization kernel on )r?r@)rZ device_arrayrCrD)r<rr;r@rrrrcreate_xoroshiro128p_states srE)rr)rr)#r/Znumbarrrrrrrr r ZnumpyrAr?rCZxoroshiro128p_typeZENABLE_CUDASIMZ _forceobjZ _loopliftrrrr%r)r*r,r-pir3r8r7r9r>rDrErrrrsB,