10
Testowałem przy użyciu implementacji Twister Mersenne w bibliotece standardowej gcc C++. Przewyższa zarówno liniowy generator kongruencji, jak i C rand, co jest najprawdopodobniej LCG. A boost documentation również wydaje się dawać podobny wynik, ale faworyzując twister Mersenne jeszcze bardziej. Czy ktoś może to wyjaśnić?Dlaczego tornado Mersenne jest szybsze niż liniowy generator kongruencji?
#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;
}
}
wynik
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
Czy oczekiwałeś innego wyniku? – emlai
@zenith Pewnie, LCG to tylko kilka operacji arytmetycznych, a mt19937 ma kilka stron kodu. – xiver77
Czy włączono optymalizacje? Aby uzyskać dobre wyniki, należy włączyć pełne optymalizacje i wymusić efekt uboczny (sumowanie sumy kontrolnej?) Działania czasowego, np. Wydrukowanie wyniku sumy kontrolnej * po * zatrzymaniu stopera, aby uniemożliwić kompilatorowi optymalizację działania czasowego. – Galik