| OLD | NEW |
| 1 /* | 1 /* |
| 2 * Copyright (c) 2015 The WebRTC project authors. All Rights Reserved. | 2 * Copyright (c) 2015 The WebRTC project authors. All Rights Reserved. |
| 3 * | 3 * |
| 4 * Use of this source code is governed by a BSD-style license | 4 * Use of this source code is governed by a BSD-style license |
| 5 * that can be found in the LICENSE file in the root of the source | 5 * that can be found in the LICENSE file in the root of the source |
| 6 * tree. An additional intellectual property rights grant can be found | 6 * tree. An additional intellectual property rights grant can be found |
| 7 * in the file PATENTS. All contributing project authors may | 7 * in the file PATENTS. All contributing project authors may |
| 8 * be found in the AUTHORS file in the root of the source tree. | 8 * be found in the AUTHORS file in the root of the source tree. |
| 9 */ | 9 */ |
| 10 #include "webrtc/rtc_base/random.h" | 10 #include "webrtc/rtc_base/random.h" |
| 11 | 11 |
| 12 #include <math.h> | 12 #include <math.h> |
| 13 | 13 |
| 14 #include "webrtc/rtc_base/checks.h" | 14 #include "webrtc/rtc_base/checks.h" |
| 15 #include "webrtc/rtc_base/safe_conversions.h" |
| 15 | 16 |
| 16 namespace webrtc { | 17 namespace webrtc { |
| 17 | 18 |
| 18 Random::Random(uint64_t seed) { | 19 Random::Random(uint64_t seed) { |
| 19 RTC_DCHECK(seed != 0x0ull); | 20 RTC_DCHECK(seed != 0x0ull); |
| 20 state_ = seed; | 21 state_ = seed; |
| 21 } | 22 } |
| 22 | 23 |
| 23 uint32_t Random::Rand(uint32_t t) { | 24 uint32_t Random::Rand(uint32_t t) { |
| 24 // Casting the output to 32 bits will give an almost uniform number. | 25 // Casting the output to 32 bits will give an almost uniform number. |
| 25 // Pr[x=0] = (2^32-1) / (2^64-1) | 26 // Pr[x=0] = (2^32-1) / (2^64-1) |
| 26 // Pr[x=k] = 2^32 / (2^64-1) for k!=0 | 27 // Pr[x=k] = 2^32 / (2^64-1) for k!=0 |
| 27 // Uniform would be Pr[x=k] = 2^32 / 2^64 for all 32-bit integers k. | 28 // Uniform would be Pr[x=k] = 2^32 / 2^64 for all 32-bit integers k. |
| 28 uint32_t x = NextOutput(); | 29 uint32_t x = NextOutput(); |
| 29 // If x / 2^32 is uniform on [0,1), then x / 2^32 * (t+1) is uniform on | 30 // If x / 2^32 is uniform on [0,1), then x / 2^32 * (t+1) is uniform on |
| 30 // the interval [0,t+1), so the integer part is uniform on [0,t]. | 31 // the interval [0,t+1), so the integer part is uniform on [0,t]. |
| 31 uint64_t result = x * (static_cast<uint64_t>(t) + 1); | 32 uint64_t result = x * (static_cast<uint64_t>(t) + 1); |
| 32 result >>= 32; | 33 result >>= 32; |
| 33 return result; | 34 return result; |
| 34 } | 35 } |
| 35 | 36 |
| 36 uint32_t Random::Rand(uint32_t low, uint32_t high) { | 37 uint32_t Random::Rand(uint32_t low, uint32_t high) { |
| 37 RTC_DCHECK(low <= high); | 38 RTC_DCHECK(low <= high); |
| 38 return Rand(high - low) + low; | 39 return Rand(high - low) + low; |
| 39 } | 40 } |
| 40 | 41 |
| 41 int32_t Random::Rand(int32_t low, int32_t high) { | 42 int32_t Random::Rand(int32_t low, int32_t high) { |
| 42 RTC_DCHECK(low <= high); | 43 RTC_DCHECK(low <= high); |
| 43 // We rely on subtraction (and addition) to be the same for signed and | 44 const int64_t low_i64{low}; |
| 44 // unsigned numbers in two-complement representation. Thus, although | 45 return rtc::dchecked_cast<int32_t>( |
| 45 // high - low might be negative as an int, it is the correct difference | 46 Rand(rtc::dchecked_cast<uint32_t>(high - low_i64)) + low_i64); |
| 46 // when interpreted as an unsigned. | |
| 47 return Rand(high - low) + low; | |
| 48 } | 47 } |
| 49 | 48 |
| 50 template <> | 49 template <> |
| 51 float Random::Rand<float>() { | 50 float Random::Rand<float>() { |
| 52 double result = NextOutput() - 1; | 51 double result = NextOutput() - 1; |
| 53 result = result / 0xFFFFFFFFFFFFFFFEull; | 52 result = result / 0xFFFFFFFFFFFFFFFEull; |
| 54 return static_cast<float>(result); | 53 return static_cast<float>(result); |
| 55 } | 54 } |
| 56 | 55 |
| 57 template <> | 56 template <> |
| (...skipping 19 matching lines...) Expand all Loading... |
| 77 double u2 = static_cast<double>(NextOutput()) / 0xFFFFFFFFFFFFFFFFull; | 76 double u2 = static_cast<double>(NextOutput()) / 0xFFFFFFFFFFFFFFFFull; |
| 78 return mean + standard_deviation * sqrt(-2 * log(u1)) * cos(2 * kPi * u2); | 77 return mean + standard_deviation * sqrt(-2 * log(u1)) * cos(2 * kPi * u2); |
| 79 } | 78 } |
| 80 | 79 |
| 81 double Random::Exponential(double lambda) { | 80 double Random::Exponential(double lambda) { |
| 82 double uniform = Rand<double>(); | 81 double uniform = Rand<double>(); |
| 83 return -log(uniform) / lambda; | 82 return -log(uniform) / lambda; |
| 84 } | 83 } |
| 85 | 84 |
| 86 } // namespace webrtc | 85 } // namespace webrtc |
| OLD | NEW |