o 8ήc'@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)align)forceobjloopliftcCsxt|}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)statesindexseedzr8/tmp/pip-target-vg8gfxp4/lib/python/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  *r"cCst|}tdtdf}td}td}tdD]*}tdD]#}||tdt|>@r<|||dN}|||dN}t||qq|||d<|||d<d S) zAdvance the RNG in ``states[index]`` by 2**64 steps. :type states: 1D array, dtype=xoroshiro128p_dtype :param states: array of RNG states :type index: int64 :param index: offset in states to update lzKW`u l"U .XC rrr r N)rrrangerr")rrjumpr r ibrrrxoroshiro128p_jumpes     r)cCs,t|}|td?tdtdtd>S)z7Convert uint64 to float64 value in the range [0.0, 1.0) r$5)rrrrrrruint64_to_unit_float64s$r-cCst|}tt|S)z7Convert uint64 to float32 value in the range [0.0, 1.0))rrr-r,rrruint64_to_unit_float32s r.cCt|}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.r"rrrrrxoroshiro128p_uniform_float32 r1cCr/)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-r"r0rrrxoroshiro128p_uniform_float64r2r3r#cCHt|}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 @)rr1mathsqrtrlogcosTWO_PI_FLOAT32rru1u2z0rrrxoroshiro128p_normal_float32   (r?cCr4)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 r5)rr1r6r7rr8r9TWO_PI_FLOAT64r;rrrxoroshiro128p_normal_float64r@rBcCsx|jd}t|}t|}|dkr8t|d|t|D]}t|dqtd|D]}||d||<t||q(dSdS)Nrr$)shaperrr%r))rrsubsequence_startn_r'rrrinit_xoroshiro128p_states_cpus     rGcCs.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 )rCdtype)streamN)npemptyrCxoroshiro128p_dtyperGcopy_to_device)rrrDrI states_cpurrrinit_xoroshiro128p_statess rOcCs"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 )rHrI)r device_arrayrLrO)rErrDrIrrrrcreate_xoroshiro128p_states srQ)rr)#r6numbarrrrrrrr r numpyrJrHrLxoroshiro128p_typeENABLE_CUDASIM _forceobj _loopliftrrr"r)r-r.r1r3pir:rAr?rBrGrOrQrrrrsD,