mersenne_twister Class

Artykuł
02/27/2013

Generuje losowe przez algorytm Mersenne twister.Zachowane zgodności TR1.Użyj mersenne_twister_engine Class zamiast.

template<class UIntType,
   int W, int N, int M, int R,
   UIntType A, int U, int S,
   UIntType B, int T, UIntType C, int L>
   class mersenne_twister {
public:
   typedef mersenne_twister<UIntType, W, N, M, R,
      P, U, S, B, T, C, L> _MyT;
   typedef UIntType result_type;
   static const int word_size = W;
   static const int state_size = N;
   static const int shift_size = M;
   static const int mask_bits = R;
   static const int UIntType parameter_a = A;
   static const int output_u = U;
   static const int output_s = S;
   static const UIntType output_b = B;
   static const int output_t = T;
   static const UIntType output_c = C;
   static const int output_l = L;
   static const UIntType default_seed = 5489U;
   explicit mersenne_twister(unsigned long x0 = default_seed);
   mersenne_twister(const mersenne_twister& right);
   mersenne_twister(mersenne_twister& right);
   template<class Gen>
      mersenne_twister(Gen& gen);
   void seed(unsigned long x0 = default_seed);
   template<class Gen>
      void seed(Gen& gen);
   result_type min() const;
   result_type max() const;
   result_type operator()();
   };

Parametry

UIntType
Typ wyniku liczba całkowita bez znaku.
W
Parametr w silnika.
N
Parametr n silnika.
R
Parametr silnika R.
A
Parametr silnika A.
U
Parametr u silnika.
S
Parametr silnika S.
B
B parametr silnika.
T
Parametr silnika T.
C
Parametr c silnika.
L
Parametr l silnika.

Uwagi

Decribes klasa szablonu proste silnika.Przechowuje ona dużą wartość integralną z W * (N - 1) + R bitów.Wyodrębnia on W – bitów w czasie, z tym dużą wartość, i gdy wykorzystane wszystkie bity on ruch dużą wartość przesunięcia i mieszanie bitów, aby nowy zestaw bitów, aby wyodrębnić z.Państwo silnika jest ostatnim NW-bitowe wartości używane, jeśli operator() została wywołana co najmniej N razy inaczej MW-bitowe wartości, które zostały użyte i ostatnią N - M wartości materiału siewnego.

Argument szablonu UIntType musi być wystarczająco duży, aby pomieścić do wartości 2W - 1.Wartości innych argumentów szablon musi spełniać następujące wymagania:

0 < M <= N

0 <= R, U, S, T, L <= W

0 <= A, B, C <= 2W

W * (N - 1) + Rmusi być Mersenne prime

Generator – ruch dużą wartość, przechowującą wykonując następujący kod:

    for (int i = 0; i < N; ++i)
        {
        temp = (x[i] & LMASK) << (W - 1) | (x[i + 1] & HMASK) >> 1;
        if (temp & 1)
            y[i] = (temp >> 1) ^ A ^ x[(i + R) % N];
        else
            y[i] = (temp >> 1) ^ x[(i + R) % N];
        }
        for (int i = 0; i < N; ++i)
            x[i] = y[i];

gdy LMASK jest niepodpisany W-bitowe wartości z jego niską R bity są ustawione na 1, a pozostałe bity, jego wartość 0, i HMASK jest dopełnienie LMASK.

Generator posiada bieżący indeks idx równe 0.Rozpakowuje bitów wykonując następujący kod:

    temp = x[idx++];
    temp = temp ^ (temp >> U);
    temp = temp ^ ((temp << S) & B);
    temp = temp ^ ((temp << T) & C);
    temp = temp ^ (temp >> L);

Gdy idx osiągnie N generator – ruch wartość przechowywana i ustawia idx do 0.