randomize - random double c++

uniform_real_distribution Is Not Uniform

Question

Por favor, ayúdame a entender esto. Después de ejecutar el fragmento:

random_device randomEngine;
mt19937 generatorEngine(randomEngine());
uniform_real_distribution<double> unifRandomValue(
numeric_limits<double>::min(),
numeric_limits<double>::max());
double test[1000];
for (int i{ 0 }; i < 1000; ++i) {
test[i] = unifRandomValue(generatorEngine);
}

¿Por qué cada valor generado en el rango [1.0E306, 1.8E308] ? Esperaba un valor aleatorio distribuido uniformemente desde cerca de 0 hasta el máximo de tipo doble.

¡Gracias por adelantado!

Aquí hay un ejemplo más completo:

const size_t size{ 1000 };
std::random_device randomEngine;
std::mt19937 generatorEngine(randomEngine());
std::uniform_real_distribution<double> unifRandomValue(
std::numeric_limits<double>::min(),
std::numeric_limits<double>::max());
std::array<double, size> test;
for (int i{ 0 }; i < size; ++i) {
test[i] = unifRandomValue(generatorEngine);
}
auto minMaxIt = std::minmax_element(test.begin(), test.end());
// average without overflow
double average{ 0.0 };
for (int i{ 0 }; i < size; ++i) {
average += test[i] / size;
}
std::cout << "min value : " << *minMaxIt.first << std::endl;
std::cout << "max value : " << *minMaxIt.second << std::endl;
std::cout << "average : " << average << endl;
// one possible output (they are all similar)
// min value : 1.73361e+305
// max value : 1.79661e+308
// average : 8.78467e+307

Answer 1

Bueno, esa es la propiedad de un generador uniforme :

El 90% de los valores estarán en el orden de magnitud más alto que especifique.

Pensar más pequeño; considere el rango entero de 0 a 99 inclusive: el 90% de los dibujos tendrá 2 dígitos.