Single runden

[Scurra]

Hallo zusammen,

ich stehe gerade vor einem Problem: Ich möchte Gleitkommazahlen (vom Typ Single, um genauer zu sein), exakt (d. h. kaufmännisch) runden. Das Runden soll für verschiedene Genauigkeiten (Nachkommastellen) möglich sein. Dazu habe ich mir folgende Funktion geschrieben (ich arbeite hier mit dem Datentyp Extended, weil ich bei der Verwendung von Single noch mehr Probleme habe, obwohl ich in meinem Programm immer Werte vom Typ Single an die Funktion übergebe):

Delphi-Quellcode:

function SetPrecision(aValue: Extended; aPrecision: Integer): Extended;
var
  multi : Single;
begin
  if aPrecision < 0 then Result := aValue
  else
  begin
    multi := IntPower(10, aPrecision);
    if aValue >= 0 then aValue := Trunc(aValue*multi + 0.5)
    else aValue := Trunc(aValue*multi - 0.5);
    Result := aValue/multi;
  end;
end;

Das funktioniert auch fast immer. Allerdings gibt es noch Probleme, wie mir mit dUnit gerade aufgefallen ist. Wenn an der entscheidenen Stelle fürs Runden eine 5 steht, dann kommt es vor, dass falsch gerundet wird, da die Zahlen intern wohl anders gespeichert werden (die Zahl 2.5 ist nicht nur 2.5, sondern so etwas wie 2.4999845). Dann wird natürlich 2.5 nicht auf 3, sondern auf 2 gerundet.

Um dieses Problem zu lösen, hatte ich die Idee, die Funktion SetPrecision zunächst noch einmal rekursiv aufzurufen, etwa mit:

Delphi-Quellcode:

function SetPrecision(aValue: Extended; aPrecision: Integer; firstTime: Boolean = True): Extended;
var
  multi : Single;
begin
  if firstTime then aValue := SetPrecision(aValue, aPrecision+1, false);

  if aPrecision < 0 then Result := aValue
  else
  begin
    multi := IntPower(10, aPrecision);
    if aValue >= 0 then //aValue := Trunc(aValue*multi + 0.5)
    begin
      aValue := aValue*multi;
      aValue := aValue +0.5;
      if Frac(aValue) <> 0 then aValue := Trunc(aValue);
    end
    else aValue := Trunc(aValue*multi - 0.5);
    Result := aValue/multi;
  end;
end;

Das funktioniert leider immer noch nicht. Ich habe den Debugger verwendet, um nachvollziehen zu können, was schief läuft (daher z. B. auch die Zeile if Frac(aValue) <> 0 then aValue := Trunc(aValue); Genau in dieser Zeile tritt das Problem auf. Nehmen wir als Beispiel den Aufruf Aufruf von SetPrecision(x, 4) mit der Variable x = 44.99995 als Single. Laut Debugger ist im zweiten Durchlauf direkt vor Ausführen der oben genannten Zeile aValue = 450000. Intern wird da aber wohl wieder irgendein Rest mitgezogen, so dass die if-Abfrage ein true liefert und nach Trunc(aValue) habe ich 449999. Am Ende steht also der Wert 44.9999 dran, obwohl es 45 sein sollte.

Hat jemand eine Idee, wie ich das Problem in den Griff bekommen kann?

P.S.: Ja, ich habe schon im Internet nach Lösungen gesucht. Wie ich gesehen habe, gibt es schon einige Beiträge dazu, aber mein Problem konnte ich leider nicht lösen.

[Sherlock]

An welcher Stelle versagen denn da die Bordmittel? Es gibt ja einige Rundungsmethoden in Delphi. Exaktes Runden ist übrigens ein Klasse Oxymoron.

Sherlock

[gammatester]

Du vermischt zwei verschiedene Aspekte beim Runden: 2.5 ist exakt als single darstellbar! Aber es wird round-to-even benutzt und deshalb nicht auf die ungerade 3 aufgerundet sondern auf die gerade 2 abgerundet. 3.5 ist auch exakt darstellbar, wird aber auf 4 gerundet.

[himitsu]

RoundTo? :stupid:


Den Datentypen Extended solltest du besser nicht verwenden.
Reicht Double denn wirklich nicht aus?


http://docwiki.embarcadero.com/RADSt...eitkommawerten
http://docwiki.embarcadero.com/RADSt...e_%28Delphi%29
http://docwiki.embarcadero.com/Libra...ystem.Extended
http://docwiki.embarcadero.com/RADSt...e_Datenformate -> Auf der 64-Bit-Intel-Plattform und ARM-Platform ist der Typ Extended ein Alias für den Typ Double ...
http://docwiki.embarcadero.com/RADSt...ommaarithmetik
http://docwiki.embarcadero.com/RADSt...d_%28Delphi%29

[BadenPower]

Zitat:

Zitat von himitsu (Beitrag 1293945)

RoundTo? :stupid:

SimpleRoundTo()

Dann braucht er nicht noch SetRoundMode, wenn er kaufmännisch runden will.

[Scurra]

Zitat:

Exaktes Runden ist übrigens ein Klasse Oxymoron.

Ich hätte wohl besser "richtiges Runden" schreiben sollen ;)

Zitat:

Zitat von himitsu (Beitrag 1293945)

Reicht Double denn wirklich nicht aus?

http://docwiki.embarcadero.com/RADSt...eitkommawerten
http://docwiki.embarcadero.com/RADSt...e_%28Delphi%29
http://docwiki.embarcadero.com/Libra...ystem.Extended
http://docwiki.embarcadero.com/RADSt...e_Datenformate -> Auf der 64-Bit-Intel-Plattform und ARM-Platform ist der Typ Extended ein Alias für den Typ Double ...
http://docwiki.embarcadero.com/RADSt...ommaarithmetik
http://docwiki.embarcadero.com/RADSt...d_%28Delphi%29

Ok, dann ist es wohl besser, wenn ich Double verwende. Diese Genauigkeit genügt mir eigentlich auch schon. Das Runden sollte am besten für Single und Double gleiche Ergebnisse liefern. Gemäß dem obersten Link ist SimpleRoundTo genau die Funktion, die ich haben möchte.

Allerdings gibt es auch hier ein für mich unerwünschtes Verhalten (siehe die jeweiligen Ergebnisse als Kommentar):

Delphi-Quellcode:

procedure TForm1.Button1Click(Sender: TObject);
var
  x1: Single;
  x2: Double;
  y1: Single;
  y2: Double;
begin
  x1 := 44.99995;
  x2 := 44.99995;
  ShowMessage(FloatToStr(SimpleRoundTo(x1, -4))); // liefert 45
  ShowMessage(FloatToStr(SimpleRoundTo(x2, -4))); // liefert 44.9999

  y1 := 4.99995;
  y2 := 4.99995;
  ShowMessage(FloatToStr(SimpleRoundTo(y1, -4))); // liefert 4.999899864
  ShowMessage(FloatToStr(SimpleRoundTo(y2, -4))); // liefert 5
end;

Oder liegt das nur daran, dass ich es von zu Hause aus gerade nicht mit Delphi testen kann, sondern nur mit Lazarus?

Zitat:

Du vermischt zwei verschiedene Aspekte beim Runden: 2.5 ist exakt als single darstellbar! Aber es wird round-to-even benutzt und deshalb nicht auf die ungerade 3 aufgerundet sondern auf die gerade 2 abgerundet. 3.5 ist auch exakt darstellbar, wird aber auf 4 gerundet.

Ich sehe zwar in meinem Code nicht so etwas wie round-to-even, aber sollte dann 49.99995 nicht erst recht auf 45 gerundet werden, wenn man auf 4 Nachkommastellen rundet (9.5 wird zur geraden Zahl, also zu 10 gerunden)?

Danke schon einmal für eure Hilfe. Bisher waren mir nur die Rundungsfunktionen wie Trunc, Ceil und Round bekannt. Beim googlen bin ich fast nur auf Beiträge gestoßen, wo das kaufmännische Runden selbst in Funktionen implementiert wurde. Das wundert mich, wenn es doch schon bereitgestellte Funktionen dafür gibt.

[Sir Rufo]

Ich stell mal die entscheidende Frage:

Der Wert der da herauskommt hat welche Bedeutung? Etwa ein Währungsbetrag?

Dann bist du mit deiner Rechnerei eh auf dem Holzweg, denn der Wert sollte dann final vom Typ

Delphi-Quellcode:

Currency

sein.

[Scurra]

Zitat:

Der Wert der da herauskommt hat welche Bedeutung?

Ich habe ein Dokument, in dem sich Formen (Rechtecke, Bilder o. ä.) befinden. Ich möchte beispielsweise die Positionen der Formen oder die Seitenränder (z. B. in cm) o. ä. speichern und später wieder abrufen können. Da die Werte, die ich auslesen kann, aber immer x-Stellen lang sind und ich nicht so eine hohe Genauigkeit brauche, runde ich die Zahlen lieber.

Natürlich sollte immer nach der gleichen Vorschrift gerundet werden, aber solche Grenzfälle mit .5 bestehen meistens (meistens wird wohl auch bei .5 nach meiner Erwartung gerundet) meine Testcases... Dass da etwas nicht richtig funktioniert, ist mir aber erst durch Testcases (mit dUnit) aufgefallen.

[BadenPower]

Wenn Dir ein Double reicht, dann kannst Du folgendes machen, um die richtigen gerundeten Werte zu erhalten:

Delphi-Quellcode:

var
  x2: Double;
  y2: Double;
begin
  SetPrecisionMode(pmDouble);
  SetRoundMode(rmTruncate);

  x2 := 44.99995;
  y2 := 4.99995;

  ShowMessage(FloatToStr(SimpleRoundTo(x2, -4)));
  ShowMessage(FloatToStr(SimpleRoundTo(y2, -4)));
end;

Deine Gleitkomma-Variablen müssen allerdings vom Typ Double sein.

[Scurra]

Funktioniert das wirklich?

Delphi-Quellcode:

ShowMessage(FloatToStr(SimpleRoundTo(x2, -4)));

Das liefert bei mir 44,9999 (anstatt 45).

Was mir sonst gerade einfällt: Es ist zwar umständlich, aber falls das mit dem Runden einfach nicht klappt, muss man eben "von Hand" runden: Man frägt ab, ob die fürs Runden entscheidende Ziffer (z. B. 5. Ziffer hinter dem Punkt/Komma) eine 1, 2, 3 oder 4 ist, dann soll abgerundet werden, bei 5, 6, 7, 8 oder 9 soll aufgerundet werden und bei 0 soll die Zahl nur abgeschnitten werden. Eine einfachere und elegantere Lösung würde ich bevorzugen :D