U hd(@sddlZddlmZmZmZmZmZmZmZm Z m Z ddl Z e j de jfde jfgddZe eZejZZej Ze eeeddd Ze eeedd d Ze eeedd d Ze eeedddZe eeedddZe eeedddZe eeedddZe eeedddZe dejZe dejZe eeedddZe eeedddZ e eeedddZ!d#dd Z"d$d!d"Z#dS)%N) configcudafloat32float64uint32int64uint64 from_dtypejits0s1T)align)ZforceobjZloopliftZnopythoncCsxt|}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)statesindexseedzrJ/home/sam/Atlas/atlas_env/lib/python3.8/site-packages/numba/cuda/random.pyinit_xoroshiro128p_state.s  rcCs(t|}t|}||>|td|?BS)zLeft rotate x by k bits.@)rr)xkrrrrotlHsrcCsxt|}||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_nextPs   *r cCst|}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+r rrrrrxoroshiro128p_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*r r,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   (r8cCsHt|}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-r0r1rr2r3TWO_PI_FLOAT64r5rrrxoroshiro128p_normal_float64s   (r:cCst|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<rAZ states_cpurrrinit_xoroshiro128p_statess rEcCs"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@rA)rZ device_arrayrDrE)r=rr<rArrrrcreate_xoroshiro128p_states srF)rr)rr)$r0Znumbarrrrrrrr r numpyrBr@rDZxoroshiro128p_typeZENABLE_CUDASIMZ _forceobjZ _loopliftZ _nopythonrrr r&r*r+r-r.pir4r9r8r:r?rErFrrrrsD,