10
Ho eseguito il test con l'implementazione twister Mersenne della libreria standard gcc C++. Supera sia il generatore a congruenza lineare che il C rand, che è probabilmente un LCG. Anche il A boost documentation sembra dare un risultato simile, ma favorisce ancora di più il twister di Mersenne. Qualcuno può spiegarlo?Perché Mersenne è più veloce del generatore a congruenza lineare?
#include <cstdlib>
#include <iostream>
#include <chrono>
#include <random>
class Timer
{
private:
std::chrono::high_resolution_clock::time_point start_time;
std::chrono::high_resolution_clock::time_point stop_time;
public:
void start()
{
start_time = std::chrono::high_resolution_clock::now();
}
void stop()
{
stop_time = std::chrono::high_resolution_clock::now();
}
double measure()
{
using namespace std::chrono;
return duration_cast<microseconds>
(stop_time - start_time).count()/1000000.0;
}
};
template<typename T>
class Random
{
private:
T generator;
public:
Random()
: generator
(std::chrono::high_resolution_clock::now().time_since_epoch().count())
{
}
int generate_integer(int begin, int end)
{
return std::uniform_int_distribution<int>(begin, end - 1)(generator);
}
};
int main()
{
constexpr int n = 300000000;
Random<std::minstd_rand> mr;
Random<std::mt19937> mt;
Timer t;
for (int j = 0; j < 3; ++j)
{
t.start();
for (int i = 0; i < n; ++i)
{
static_cast<volatile void>(mr.generate_integer(0, 10));
}
t.stop();
std::cout << "minstd " << t.measure() << std::endl;
t.start();
for (int i = 0; i < n; ++i)
{
static_cast<volatile void>(mt.generate_integer(0, 10));
}
t.stop();
std::cout << "mersenne " << t.measure() << std::endl;
t.start();
for (int i = 0; i < n; ++i)
{
static_cast<volatile void>(std::rand() % 10);
}
t.stop();
std::cout << "rand " << t.measure() << std::endl;
}
}
risultato
minstd 4.70876
mersenne 1.55853
rand 4.11873
minstd 4.53199
mersenne 1.55928
rand 4.15159
minstd 4.5374
mersenne 1.55667
rand 4.13715
Ti aspettavi un risultato diverso? – emlai
@zenith Certo, un LCG è solo alcune operazioni aritmetiche, mentre mt19937 è di diverse pagine di codice. – xiver77
Hai attivato le ottimizzazioni? Per buoni risultati è necessario attivare ottimizzazioni complete e forzare un effetto collaterale (accumulando un checksum?) Dell'attività temporizzata come stampare un risultato di checksum * dopo * che il timer si è fermato per impedire al compilatore di ottimizzare l'attività temporizzata. – Galik