mersenne_twister Class

Článek
03/04/2013

Vygeneruje náhodné posloupnosti algoritmus twister Mersenne.Zachována kompatibilita TR1.Jako náhradu použijte mersenne_twister_engine Class.

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
Výsledek typu celé číslo bez znaménka.
W
Parametr w motoru.
N
Parametr n motoru.
R
Parametr r motoru.
A
Parametr a motoru.
U
Parametr u motoru.
S
Parametr s motorem.
B
Parametr b motoru.
T
Parametr t motoru.
C
Parametr c motoru.
L
Parametr l motoru.

Poznámky

Popsáno třída Šablona jednoduchého motoru.Obsahuje velký integrální hodnota s W * (N - 1) + R bitů.Extrahuje W bitů v době od této velké hodnoty a použije všechny bity je twists velké hodnoty posunu a tak, aby nová sada bitů extrahovat z promíchání bitů.Stav motoru je poslední NW-bit hodnoty použít, pokud operator() byla volána alespoň N jinak časy MW-bit hodnoty, které byly použity a poslední N - M hodnoty osiva.

Argument šablony UIntType musí být dostatečně velké pro uložení hodnot až 2W - 1.Hodnoty jiné šablony argumenty musí splňovat následující požadavky:

0 < M <= N

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

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

W * (N - 1) + Rmusí být Mersenne prime

Generátor twists velká hodnota obsahující provedením následující kód:

    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];

kde LMASK je nepodepsaný W-bit hodnotu jeho nízké R bity hodnotu 1 a zbytek bitů nastaven na hodnotu 0, a HMASK je doplňkem z LMASK.

Generátor uchovává aktuální index idx 0 inicializován.Extrahuje bitů provedením následující kód:

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

Při idx dosáhne N twists zdroje tepla uložené a nastaví idx zpět na 0.