Random: Delphi vs. Lazarus?

[tigernut]

Zitat von gammatester:

Der in Delphi und Pascal ist seit langem der gleiche, ein linearer http://de.wikipedia.org/wiki/Kongruenzgenerator
x[n+1] = a*x[n] + 1 mod 2^32; mit a = $08088405 = 134775813

Edit: der gleiche heißt hier, daß er sich seit TurboPascal 5/6/7 nicht geändert hat, nicht das er der gleiche wie in neueren FPC-Versionen ist.

Hallo,

weisst du wie die Funktion random(pWert:integer):integer arbeitet? Ist das einfach:

random:=zufallszahl mod pWert; Gruß

[himitsu]

Nee

Delphi-Quellcode:

function Random(const ARange: Integer): Integer;
{$IFDEF PUREPASCAL}
var
  Temp: Longint;
begin
  Temp := RandSeed * $08088405 + 1;  // berechne den nächsten Wert
  RandSeed := Temp;  //
  Result := (UInt64(Cardinal(ARange)) * UInt64(Cardinal(Temp))) shr 32;
end;

Es ist im Prinzip eher etwas wie
Result := Random * pWert; // Random liefert ja 0 bis 0.999999 (oder war es 0 bis 1?)

[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.