Einzelnen Beitrag anzeigen

Benutzerbild von sx2008
sx2008

Registriert seit: 15. Feb 2008
Ort: Baden-Württemberg
2.332 Beiträge
 
Delphi 2007 Professional
 
#3

AW: WELL 1024a Zufallszahlengenerator

  Alt 29. Mai 2013, 14:07
Wie kommen die nur auf solche Periodenlängen?
Das hängt von der Grösse der inneren Zustands ab.
Bei Well1024a ist das STATE : array[0..32-1] of Cardinal ; 32 * sizeof(Cardinal) * 8 = 1024 Bit.
Allerdings werden nicht alle Zustände durchlaufen (insbesondere Zustand 0) daher kommt die -1.

Beim Mersenne-Twister sind es 624 * sizeof(Integer) * 8 = 19968 Bit wobei wohl ein Bit nicht genutzt werden kann und dann halt "nur" 2^219937-1 Zustände bleiben.

Der Vorteil des Well gegenüber Mersenne-Twister ist, das der Well-Generator schon nach ~ 100 Schritten "eingelaufen" ist, während Mersenne-Twister dazu ~ 70000 Schritte braucht.

"Einlaufen" heisst, man muss entsprechend viele Zufallszahlen abrufen und wegwerfen, bevor man mit guten gleichverteilten Zahlen rechnen kann.

Wie sieht die Performance beider PRNG aus?
Hab's nicht gemessen, aber der Well1024 dürfte etwas schneller sein weil er keine for-Schleifen enthält.
fork me on Github
  Mit Zitat antworten Zitat