FMod "streikt"

[Rollo62]

Ich hätte dazu auch noch ne Frage:
Wie bekommt man denn eine 150-stellige Zahl überhaupt hin ?
Extended sollte doch sofort auf 18-20 Stellen kürzen, oder bleibt 1.2E150 etwa so bestehen ?

Delphi-Quellcode:

    LExt := 1.2E-150;
    LExt := LExt + 0.1E-150;  // Das geht tatsächlich
    LExt := Sin(LExt);        // Das aber wohl nicht

    LExt := 1.2E-150;
    LExt := LExt + 0.1E-149;  // Das geht tatsächlich

    LExt := 1.2E150;
    LExt := LExt + 0.1E150;   // Das geht tatsächlich
    LExt := Sin(LExt);        // Das wohl auch nicht

Ich benutze solche Extremen normalweise nicht, bin aber auch etwas schockiert das es nur "teilweise" funktioniert.
Anscheinend kann man mit Mantisse/Exponent rechnen, aber mit keine Funktion.
Da hoffe ich mal das ich sowas nie brauchen werde.

Rollo

[himitsu]

Zitat:

Zitat von Rollo62 (Beitrag 1375402)

Ich hätte dazu auch noch ne Frage:
Wie bekommt man denn eine 150-stellige Zahl überhaupt hin ?
Extended sollte doch sofort auf 18-20 Stellen kürzen, oder bleibt 1.2E150 etwa so bestehen ?

Siehe Post #2.

Das Stichwort heißt "signifikante Dezimalstellen"

Eine Fließkommazahl besteht aus dem Vorzeichen (+ oder -), der Mantisse und dem Exponent.
Die Mantisse ist der Teil, in welchem die 18-20 "signifikanten" Dezimalstellen stecken und mit dem Exponent kann man ihn hin und her schieben, also bei Extended um fast 5000 Dezimalstellen vor oder hinter das Komma. (bei Double etwa 300 Dezimalstellen)

3,4xxxxxxxxxxxxxxxxx * 10^-4932 bis 1,1xxxxxxxxxxxxxxxxxx * 10^4932
3,4xxxxxxxxxxxxxxxxx / 10^4932 bis 1,1xxxxxxxxxxxxxxxxxx * 10^4932
das im Positiven und nochmal im degativen Bereich, zuzüglich der 0 und den Sonterwerten NAN, Infinity und NegInfinity.

https://www.h-schmidt.net/FloatConverter/IEEE754de.html
https://de.wikipedia.org/wiki/Mantisse

[Sherlock]

Zur Frage wie man 150 Stellen überhaupt hinbekommt: Man muss sich mit eigenen Typen behelfen. Hier mal ein Beispiel:
http://www.delphiforfun.org/Programs...g_integers.htm

Sherlock

[gammatester]

Zitat:

Zitat von Rollo62 (Beitrag 1375402)

Ich hätte dazu auch noch ne Frage:
Wie bekommt man denn eine 150-stellige Zahl überhaupt hin ?
Extended sollte doch sofort auf 18-20 Stellen kürzen, oder bleibt 1.2E150 etwa so bestehen ?

Delphi-Quellcode:

    LExt := 1.2E150;
    LExt := LExt + 0.1E150;   // Das geht tatsächlich
    LExt := Sin(LExt);        // Das wohl auch nicht

Ich benutze solche Extremen normalweise nicht, bin aber auch etwas schockiert das es nur "teilweise" funktioniert.

Ich glaube Du solltest Dir erst einmal klar werden, was Du unter "Stellen" verstehst. Bei Ganzzahl-Werten sind die "Stellen" tatsächlich so etwas wie die Anzahl der einzelnen (Dezimal-)Ziffern.

Bei Fließkomma-Zahlen (double) hat man sowas wie m*10^e, wobei im wesentlichen m eine Ganzzahl zwischen 0 und ca 10^16 ist. Hier hat ein double also ca 16 (signifikante "Stellen").

Bei Fließkomma-Rechnungen wird immer auf diese 16 Stellen gerundet und das hat manchmal komische Konsequenzen, zB wird 2E16 + 1 = 20000000000000000 + 1 = 20000000000000001 zu 20000000000000000 gerundet, d.h. Du kannst beliebig oft 1 zu 2E16 addieren ohne das sich der Double-Wert ändert.

Daß sin(1e150) nicht funktioniert, ist ein Delphi-Bug (EMBA wird sagen 'Feature').

[Rollo62]

Ja ich bin von signifikanten Stellen ausgegangen. Das sin (e150) nicht funktioniert finde ich schon erschreckend. Ich hatte mit dem erstbesten online js Rechner getestet und der konnte etwas sinnvolles zurückgeben.

Das heisst für mich erstmal Achtung bei fliesskomma Rechnungen.
Mit Rundungsfehlern hätte ich gerechnet aber solche logikfehler eher nicht.

[Johann Steiner]

Hallo Leute,
die Ursache für die Entstehung des Problems war das März-Preisrätsel 2017 in Bild der Wissenschaft.
Dabei hat mich nicht die "Lösung des Problems" interessiert, sondern wie das Endergebnis ausschaut und ob bei jeder diversen Losziehung das gleiche Ergebnis herauskommt.
Die Berechnung war kein Problem, Double genügt, aber das gesamte Ergebnis als String darzustellen war mir unmöglich.
Als FMode bei e^21 ausstieg war ich schon sehr erstaunt (nicht einmal Single wird überschritten!), zumal in der Hilfe doch die Funktion mit Extended-Werten umgehen sollte können.
So wollte ich das machen:
function ExtToString(source:Extended): string;
const
CharDigits = '0123456789';
var Rest:integer;
begin
Result := '';
repeat
Rest:= Round(FMod(Source, 10))+1;
Result := CharDigits[Rest] + Result;
ASource := Int(Source / 10);
until (Source = 0);
end;
Nach euren Hinweisen habe ich die Unit UBigIntsV5 von Gary Darby verwendet und es klappt hervorragend.

Nochmals herzlichen Dank an alle!