Random: Delphi vs. Lazarus?

[tigernut]

Hi Himitsu.

Zwei Dinge verstehe ich nicht.

1. In der Rekursionsformel für den Zufallsgenerator heißt es x[n+1] = (a*x[n] + 1) mod 2^32
   Es wird der Rest der Division durch 2^32 als neue Zufallszahl ausgegeben. Ist das gleichbedeutend mit shr 32? shr ist doch eher äquivalent zu div?
2. Wie liefert Radom eine Kommazahl, wenn wir ausschließlich mit Ganzzahlen rechnen?

[BUG]

Zitat von tigernut:

In der Rekursionsformel für den Zufallsgenerator heißt es x[n+1] = (a*x[n] + 1) mod 2^32
Es wird der Rest der Division durch 2^32 als neue Zufallszahl ausgegeben. Ist das gleichbedeutend mit shr 32? shr ist doch eher äquivalent zu div?

Das mod 32 wird schon durch den Wertebereich von longint erledigt. Das Shiften ist zum Aufbereiten der Ausgabe.

Zitat von tigernut:

Wie liefert Radom eine Kommazahl, wenn wir ausschließlich mit Ganzzahlen rechnen?

Man kann das auf Bitebene regeln. Andererseits kann man einfach 2^32-1 durch die 32bit Zufallszahl teilen (nicht durch 0 teilen).
Satty67 hat sich hier auch ein paar Möglichkeiten überlegt.

[himitsu]

Zitat von tigernut:

Wie liefert Radom eine Kommazahl, wenn wir ausschließlich mit Ganzzahlen rechnen?

Indem dort anders gerechnet wird?

Delphi-Quellcode:

function Random: Extended;
const
  two2neg32: double = ((1.0/$10000) / $10000);  // 2^-32
{$IFDEF PUREPASCAL}
var
  Temp: LongInt;
  F: Extended;
begin
  Temp := RandSeed * $08088405 + 1;
  RandSeed := Temp;
  F  := Temp;
  Result := F * two2neg32;
end;

[gammatester]

Zitat von himitsu:

Delphi-Quellcode:

  //...
  F  := Temp;
  Result := F * two2neg32;
end;

Dieser Code liefert allerdings Zahlen im Intervall [-0.5,0.5). Für das normalerweise von random verwendete Interval [0,1) kann man

Delphi-Quellcode:

  //...
  Result := (Temp + 2147483648.0) / 4294967296.0;
end;

benutzen, und eine ähnliche Formel wenn man mehr als 32 Bit braucht.