Rundungs-Probleme bei Ausgabe als String

[Fritzew]

Zitat von raller09:

Ich hab schon bei mir nach SetRoundMode u. ä. gesucht aber nichts gefunden. nachdem ich den String erstellt habe,m wird auch eine .dll genutzt, die den erstellten String verarbeitet.

Was ist das für eine Dll?

Einfach aus Erfahrung: Nach dem Aufruf einer Fremdle kannst Du nie davon ausgehen das
die FPU noch den gewünschten State hat. Deshalb for dem Aufruf mit
GetFPURoundMode sichern und nachher mit SetFPURoundMode() wieder setzen.

[raller09]

Zitat von Fritzew:

Was ist das für eine Dll?

Einfach aus Erfahrung: Nach dem Aufruf einer Fremdle kannst Du nie davon ausgehen das
die FPU noch den gewünschten State hat. Deshalb for dem Aufruf mit
GetFPURoundMode sichern und nachher mit SetFPURoundMode() wieder setzen.

Das ist ein sehr guter Hinweis!

Ich nutze lpsolve für eine lineare Optimierung (https://sourceforge.net/projects/lpsolve/).
Das ist eine in C geschriebene .dll.

In der Delphi-Zugriffs-Datei steht:

Code:

 *  Important information
 *  Solver library is compiled for a different Control Word, you should change
 *  the Delphi control Word to avoid foating point operation errors.
 *                _control87     Set8087CW
 *  Visual Studio $9001F     ->   639
 *  GCC           $8001F     ->   895

und in der IDE (Main.pas -> in Delphi geschrieben):

Delphi-Quellcode:

initialization
  set8087cw(895);

end.

Ich habe jetzt meine Kommunikation mit der .dll in folgenden Block gepackt:

Delphi-Quellcode:

  //lRoundingMode := GetRoundMode;
  l8087CW       := Get8087CW;
  try
    {
    Hier der gesammte Zugriff auf die .dll
    }
  finally
    ShowMessage( (*
                 'alt: ' + TRttiEnumerationType.GetName( lRoundingMode) +
                 sLineBreak +
                 'neu: ' + TRttiEnumerationType.GetName( GetRoundMode) +
                 sLineBreak +
                *)
                 'alt: ' + IntToStr( l8087CW) +
                 sLineBreak +
                 'neu: ' + IntToStr( Get8087CW)
                 );

    //SetRoundMode( lRoundingMode);

    Set8087CW( l8087CW);
  end;

Das liefert mir:

1.Aufruf
alt: 4991
neu: 610

folgende Aufrufe:
alt: 4991
neu: 4962

Da ich den Wert mittels "Set8087CW" nach der .dll zurückstelle, passt "meine" Rundung/String-Erstellung, ohne das ich irgendwie zusätzlich runden muss...(hat auch tatsächlich nicht korrekt funktioniert)!


Was mich jetzt aber noch stört: Die .dll hat beim 2ten Aufruf einen anderen Wert gesetzt.
Verändert sich dadurch das Verhalten der .dll? (vermutlich, oder?)
Könnt ihr mir hier weiterhelfen, oder muss ich das mit dem Erstellern der .dll klären?


Darf ich im laufenden Programm Set8087CW nutzen? Was wird dadurch gemacht?


Danke,

[Zacherl]

Zitat von raller09:

Darf ich im laufenden Programm Set8087CW nutzen? Was wird dadurch gemacht?

Darfst du schon. Das X87 Control Word besteht aus mehreren Bits, die unterschiedliche Bedeutungen haben:
7jKfh.png

Nähere Infos dazu entnimmst du am besten der "Intel Architecture Software Developer's Manual; Vol I: Basic Architecture".

Was RC angeht, sind hier die relevanten Werte:

Code:

Round to nearest even       = $00B
Round down toward infinity  = $01B
Round up toward infinity    = $10B
Round toward zero (trunc)   = $11B

Zitat von raller09:

Was mich jetzt aber noch stört: Die .dll hat beim 2ten Aufruf einen anderen Wert gesetzt.
Verändert sich dadurch das Verhalten der .dll? (vermutlich, oder?)

Kann gut sein, dass nur die Exception Masks geändert wurden. Ich schaue mal nach .. Edit: Ja genau z.b. zwischen 4991 und 4962 ändern sich nur die Exception Masks.

[hoika]

Hallo,

Zitat:

Erst wenn ich die Rundungs-Funktion mit "lAnzahlNachkommaStellen + 1" aufrufe, scheint es so, dass ich ein reproduzierbares Ergebnis bekomme...

Das ist genau der richtige Ansatz!

[Luckie]

Zitat von hoika:

Hallo,

Zitat:

Erst wenn ich die Rundungs-Funktion mit "lAnzahlNachkommaStellen + 1" aufrufe, scheint es so, dass ich ein reproduzierbares Ergebnis bekomme...

Das ist genau der richtige Ansatz!

Ich wollte gerade sagen. Das kann es nicht sein, dass man für eine solche Funktion der Maßen tief in den Prozessorinstruktionen eingreifen muss. Wenn das der Fall ist, dann macht man irgendwas verkehrt.

[raller09]

Zitat von Zacherl:

Darfst du schon. Das X87 Control Word besteht aus mehreren Bits, die unterschiedliche Bedeutungen haben.
Nähere Infos dazu entnimmst du am besten der "Intel Architecture Software Developer's Manual; Vol I: Basic Architecture".

Was RC angeht, sind hier die relevanten Werte:

Code:

Round to nearest even       = $00B
Round down toward infinity  = $01B
Round up toward infinity    = $10B
Round toward zero (trunc)   = $11B

Habe gesten nachmittag auch das Dokument noch gefunden. Ist dort recht gut erklärt.

Zitat von Zacherl:

Kann gut sein, dass nur die Exception Masks geändert wurden. Ich schaue mal nach .. Edit: Ja genau z.b. zwischen 4991 und 4962 ändern sich nur die Exception Masks.

Hier mal die Bit-Schreibweise meiner Werte:

Code:

                                      keine Berücksichtigung
Dec     Hex        C RC   Exception   von mir
4991    137F    0001 0011 0111 1111   Delphi vor erster Start
 610     262           10 0110 0010   nach 1. Aufruf
 610     262           10 0110 0010   folgende Aufrufe

                   I   PC             Schutz mit try ...
Dec     Hex        C RC   Exception   finally
4991    137F    0001 0011 0111 1111   Delphi vor erster Start
 610     262           10 0110 0010   nach 1. Aufruf
4962    1362    0001 0011 0110 0010   folgende Aufrufe

                   I   PC             Schutz mit try...
Dec     Hex        C RC   Exception   finally + Setzen "37F"
4991    137F    0001 0011 0111 1111   Delphi vor erster Start
 895     37F           11 0111 1111   manuell gesetzt
 610     262           10 0110 0010   nach 1. Aufruf
 895     37F           11 0111 1111   manuell gesetzt
 866     362           11 0110 0010   folgende Aufrufe

                   I   PC             Schutz mit try...
Dec     Hex        C RC   Exception   finally + Setzen "262"
4991    137F    0001 0011 0111 1111   Delphi vor erster Start
 610     262           10 0110 0010   manuell gesetzt
 610     262           10 0110 0010   nach 1. Aufruf
 610     262           10 0110 0010   manuell gesetzt
 610     262           10 0110 0010   folgende Aufrufe

Fazit:

• die Dll setzt sich immer die gleiche Exception Mask
• Der Rundungs-Modus (RC) wird nicht geändert
• Infinity wird zwar geändert, spielt aber entsprechend der Dokumentation aus dem Link keine Rolle:

  Zitat:

  8.1.6 Infinity Control Flag
  The infinity control flag (bit 12 of the x87 FPU control word) is provided for compatibility with the Intel 287 Math
  Coprocessor; it is not meaningful for later version x87 FPU coprocessors or IA-32 processors. See Section 4.8.3.3,
  “Signed Infinities,” for information on how the x87 FPUs handle infinity values.

  Ist der Wert also jetzt egal?
• es gibt ggf. eine Abweichung bei Precision Control (PC), mögliche Werte
  Code:

  Precision                           PC Field
  Single Precision (24 bits)          00B
  Reserved                            01B
  Double Precision (53 bits)          10B
  Double Extended Precision (64 bits) 11B

  -> erster Aufruf läuft immer auf "Double Precision", folgende auf "Double Extended Precision", falls der vorherige Wert nicht übereinstimmt. Ich habe mir nicht die genauen Gründe angesehen. Das werde ich mit den Erstellern der dll abklären. Wenn ich dort einfach auch "262" setze, dann wird der Wert nicht geändert.
• Durch die .dll wird die Eigenschaft des "x87 FPU Control Word" für das Programm verstellt, ich muss/sollte also mir vorher den Wert merken und ich nachher zurücksetzen.

Zitat von Luckie:

Zitat von hoika:

Zitat:

Erst wenn ich die Rundungs-Funktion mit "lAnzahlNachkommaStellen + 1" aufrufe, scheint es so, dass ich ein reproduzierbares Ergebnis bekomme...

Das ist genau der richtige Ansatz!

Ich wollte gerade sagen. Das kann es nicht sein, dass man für eine solche Funktion der Maßen tief in den Prozessorinstruktionen eingreifen muss. Wenn das der Fall ist, dann macht man irgendwas verkehrt.

Das wäre mir auch lieber. Aber es hat sich leider bei mir gezeigt, dass auch mit der Rundung andere Ergebnisse (an anderer Stelle) von den "Format" und "ToStr"- Funktionen ausgegeben werden (also wird das Round... schon irgendwie anders arbeiten).
Daher ist das "try ... finally" zum "Merken und Zurücksetzen" schon notwendig, denke ich.


Vielen Dank für eure Hilfe!

[Fritzew]

Also wenn es sogar im Header der Dll so vermerkt ist, solltest Du meiner Meinung nach bei Aufrufen in die Dll so verfahren:

// Nur so hingetipt Pseudo Code

Delphi-Quellcode:

Procedure callToDll ;
begin
  Set8087CW(DllCW);  // DllCw = CtrlWord das die Dll erwartet
  CallToDll;
  Set8087CW(MYCW); // Das was Du erwartest
end;

Einfach um sicher zu gehen das die Dll auch sauber arbeitet und davon ausgeht das immer "Ihr" Wert gesetzt ist

[Luckie]

Ach das steht so in der Doku der DLL? Was macht die um Gottes Willen für komische Sachen.