splitmix64_next_range Function

private function splitmix64_next_range(gen, lo, hi, err) result(r)

Uniform integer in [lo, hi] from a splitmix64_t generator.

Arguments

Type IntentOptional Attributes Name
type(splitmix64_t), intent(inout) :: gen

Generator, advanced by the draw.
integer(kind=default_int), intent(in) :: lo

Inclusive lower bound.
integer(kind=default_int), intent(in) :: hi

Inclusive upper bound; must not be below lo.
type(error_t), intent(inout), optional :: err

Set to ERROR_VALIDATION when hi < lo.

Return Value integer(kind=default_int)


Calls

proc~~splitmix64_next_range~~CallsGraph proc~splitmix64_next_range splitmix64_next_range proc~error_set error_t%error_set proc~splitmix64_next_range->proc~error_set proc~splitmix64_uniform_below splitmix64_uniform_below proc~splitmix64_next_range->proc~splitmix64_uniform_below proc~u64_add u64_add proc~splitmix64_next_range->proc~u64_add proc~rejection_mask~2 rejection_mask proc~splitmix64_uniform_below->proc~rejection_mask~2 proc~splitmix64_next_u64 splitmix64_t%splitmix64_next_u64 proc~splitmix64_uniform_below->proc~splitmix64_next_u64 proc~u64_less u64_less proc~splitmix64_uniform_below->proc~u64_less proc~u64_shr u64_shr proc~rejection_mask~2->proc~u64_shr proc~splitmix64_next splitmix64_t%splitmix64_next proc~splitmix64_next_u64->proc~splitmix64_next proc~splitmix64_next->proc~u64_add proc~mix64 mix64 proc~splitmix64_next->proc~mix64 proc~mix64->proc~u64_shr proc~u64_mul u64_mul proc~mix64->proc~u64_mul

Called by

proc~~splitmix64_next_range~~CalledByGraph proc~splitmix64_next_range splitmix64_next_range interface~next_range next_range interface~next_range->proc~splitmix64_next_range

Variables

Type Visibility Attributes Name Initial
integer(kind=int64), private :: span_minus_1

Source Code

   function splitmix64_next_range(gen, lo, hi, err) result(r)
      !! Uniform integer in `[lo, hi]` from a splitmix64_t generator.
      type(splitmix64_t), intent(inout) :: gen
         !! Generator, advanced by the draw.
      integer(default_int), intent(in) :: lo
         !! Inclusive lower bound.
      integer(default_int), intent(in) :: hi
         !! Inclusive upper bound; must not be below `lo`.
      type(error_t), intent(inout), optional :: err
         !! Set to `ERROR_VALIDATION` when `hi < lo`.
      integer(default_int) :: r

      integer(int64) :: span_minus_1

      r = lo
      if (hi < lo) then
         if (present(err)) then
            call err%set(ERROR_VALIDATION, "pic_random_dist: next_range requires hi >= lo")
         end if
         return
      end if
      if (hi == lo) return

      ! hi - lo in unsigned 64-bit, which is exact even when the span is wider
      ! than the signed range. Negation is spelled as `not(x) + 1` because
      ! negating the most negative value directly is signed overflow.
      span_minus_1 = u64_add(int(hi, int64), u64_add(not(int(lo, int64)), 1_int64))
      r = int(u64_add(int(lo, int64), splitmix64_uniform_below(gen, span_minus_1)), default_int)
   end function splitmix64_next_range