B \$@s:ddlmZmZddlZddlmZmZmZmZm Z m Z m Z m Z ddl Zejdej fdej fgddZe eZe dd Ze d d Ze d d Ze de dfZe ddZe ddZe ddZe ddZe ddZedejZedejZe ddZe ddZe dd Z d%d!d"Z!d&d#d$Z"dS)')print_functionabsolute_importN)cudafloat32float64uint32int64uint64 from_dtypejits0s1T)ZaligncCsxt|}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)rr r)statesindexseedzr0lib/python3.7/site-packages/numba/cuda/random.pyinit_xoroshiro128p_state%s  rcCs(t|}t|}||>|td|?BS)zLeft rotate x by k bits.@)r r)xkrrrrotl?srcCsxt|}||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$)rr rr)rrr r resultrrrxoroshiro128p_nextGs   *r lzKW`u l"U .XC cCst|}td}td}xdtdD]X}xRtdD]F}t|tdt|>@rl|||dN}|||dN}t||q0Wq"W|||d<|||d<dS)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 rrr r N)rr rangeXOROSHIRO128P_JUMPrr )rrr r ibrrrxoroshiro128p_jump`s  r'cCs,t|}|td?tdtdtd>S)z7Convert uint64 to float64 value in the range [0.0, 1.0) r"5)r rr)rrrruint64_to_unit_float64ysr*cCst|}tt|S)z7Convert uint64 to float32 value in the range [0.0, 1.0))r rr*)rrrruint64_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 )rrrrrxoroshiro128p_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 g@)rr,mathsqrtrlogcosTWO_PI_FLOAT32)rru1u2z0rrrxoroshiro128p_normal_float32s   (r6cCsHt|}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 g@)rr,r.r/rr0r1TWO_PI_FLOAT64)rrr3r4r5rrrxoroshiro128p_normal_float64s   (r8cCs||jd}t|}t|}|dkrxt|d|xt|D]}t|dq8Wx,td|D]}||d||<t||qVWdS)Nrr")shaper rr#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 intializes 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 )r9dtype)streamN)npemptyr9xoroshiro128p_dtyper=Zcopy_to_device)rrr:r?Z states_cpurrrinit_xoroshiro128p_statess rCcCs"tj|t|d}t|||||S)a4Returns a new device array initialized for n random number generators. This intializes 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_arrayrBrC)r;rr:r?rrrrcreate_xoroshiro128p_statessrD)rr)rr)#Z __future__rrr.Znumbarrrrrr r r Znumpyr@r>rBZxoroshiro128p_typerrr r$r'r*r+r,r-Zpir2r7r6r8r=rCrDrrrrs*(