breaking changes to random

v1
Brett 2024-07-09 14:11:07 -04:00
parent 5e65416684
commit c88f1c3e38
2 changed files with 79 additions and 62 deletions

View File

@ -1,7 +1,7 @@
cmake_minimum_required(VERSION 3.20) cmake_minimum_required(VERSION 3.20)
include(cmake/color.cmake) include(cmake/color.cmake)
set(BLT_VERSION 0.17.29) set(BLT_VERSION 0.18.0)
set(BLT_TEST_VERSION 0.0.1) set(BLT_TEST_VERSION 0.0.1)
set(BLT_TARGET BLT) set(BLT_TARGET BLT)

View File

@ -7,87 +7,104 @@
#ifndef BLT_RANDOM_H #ifndef BLT_RANDOM_H
#define BLT_RANDOM_H #define BLT_RANDOM_H
#include <blt/std/types.h>
#include <random> #include <random>
namespace blt::random { namespace blt::random
{
// https://github.com/avaneev/komihash/tree/main
static inline uint32_t PCG_Hash(uint32_t input) { static inline blt::u32 PCG_Hash(blt::u32 input)
uint32_t state = input * 747796405u + 2891336453u; {
uint32_t word = ((state >> ((state >> 28u) + 4u)) ^ state) * 277803737u; blt::u32 state = input * 747796405u + 2891336453u;
return (word >> 22u) ^ blt::u32 word = ((state >> ((state >> 28u) + 4u)) ^ state) * 277803737u;
word; return (word >> 22u) ^ word;
} }
static inline float randomFloat(uint32_t& seed){ //https://lemire.me/blog/2018/08/15/fast-strongly-universal-64-bit-hashing-everywhere/
static inline blt::u64 murmur64(blt::u64 h)
{
h ^= h >> 33;
h *= 0xff51afd7ed558ccdL;
h ^= h >> 33;
h *= 0xc4ceb9fe1a85ec53L;
h ^= h >> 33;
return h;
}
static inline double pcg_double(blt::u32& seed)
{
seed = PCG_Hash(seed); seed = PCG_Hash(seed);
return (float)seed / (float)std::numeric_limits<uint32_t>::max(); return static_cast<double>(seed) / static_cast<double>(std::numeric_limits<blt::u32>::max());
} }
static inline float pcg_float(blt::u32& seed)
{
return static_cast<float>(pcg_double(seed));
}
/** /**
* @return random float without changing seed * @return random float without changing seed
*/ */
static inline float randomFloat_c(uint32_t seed){ static inline float pcg_float_c(blt::u32 seed)
return randomFloat(seed); {
return pcg_float(seed);
} }
static inline double pcg_double_c(blt::u32 seed)
{
return pcg_double(seed);
}
/** /**
* @param seed seed for random * @param seed seed for random
* @param min inclusive min * @param min inclusive min
* @param max exclusive max * @param max exclusive max
* @return random int between min (inclusive) and max (exclusive) * @return random int between min (inclusive) and max (exclusive)
*/ */
static inline int randomInt(uint32_t& seed, int min = 0, int max = 1){ static inline int pcg_int(blt::u32& seed, int min = 0, int max = 2)
return (int)((randomFloat(seed) * (float)(max - min)) + (float)min); {
return static_cast<int>((pcg_double(seed) * static_cast<double>(max - min)) + static_cast<double>(min));
} }
static inline int randomInt_c(uint32_t seed, int min = 0, int max = 1){ static inline int pcg_int_c(blt::u32 seed, int min = 0, int max = 2)
return randomInt(seed, min, max); {
return pcg_int(seed, min, max);
}
static inline double murmur_double(blt::u64& seed)
{
seed = murmur64(seed);
return static_cast<double>(seed) / static_cast<double>(std::numeric_limits<blt::u64>::max());
}
static inline float murmur_float(blt::u64& seed)
{
return static_cast<float>(murmur_double(seed));
}
static inline float murmur_float_c(blt::u64 seed)
{
return murmur_float(seed);
}
static inline double murmur_double_c(blt::u64 seed)
{
return murmur_double(seed);
} }
/**
* Creates a container class for generating random number distributions
* @tparam T numeric type
* @tparam dist std::uniform_real_distribution or std::uniform_int_distribution
*/
template<typename T, template<typename = T> typename dist = std::uniform_real_distribution>
class random {
private:
std::random_device rd; // obtain a random number from hardware
std::mt19937 gen;
dist<T>* distribution = nullptr;
public:
/**
* Construct the random number generator.
* @param min min value possible to generate. (default: 0)
* @param max max value possible to generate. (default: 1)
* @param seed seed to use in generating random values. (default: 0)
*/
explicit random(T min = (T) 0, T max = (T) 1, long seed = 0): gen(std::mt19937(seed)) {
distribution = new dist(min, max);
}
/**
* Note the min/max are inclusive and defaults to a **uniform** distribution.
* @return random number between the defined min/max or the default of [0,1].
*/
T get() {
return (*distribution)(gen);
}
~random() {
delete distribution;
}
};
template<typename T> template<typename T>
class simplex_noise { static inline T murmur_integral_64(blt::u64& seed, T min = 0, T max = 2)
private: {
return static_cast<T>((murmur_double(seed) * static_cast<double>(max - min)) + static_cast<double>(min));
public: }
simplex_noise() {
template<typename T>
} static inline T murmur_integral_64c(blt::u64 seed, T min = 0, T max = 2)
}; {
return murmur_integral_64(seed, min, max);
}
} }