Finalize CLN 1.3.7 release.

[cln.git] / src / float / dfloat / elem / cl_DF_plus.cc

diff --git a/src/float/dfloat/elem/cl_DF_plus.cc b/src/float/dfloat/elem/cl_DF_plus.cc
index bae86a1623ad49ab82ae10936fecd15be65ab6c9..63077e87157ad16d43a6a1f92938346f864459b1 100644 (file)
--- a/src/float/dfloat/elem/cl_DF_plus.cc
+++ b/src/float/dfloat/elem/cl_DF_plus.cc
@@ -1,7 +1,7 @@
 // binary operator +
 
 // General includes.
-#include "cl_sysdep.h"
+#include "base/cl_sysdep.h"
 
 // Specification.
 #include "cln/dfloat.h"
@@ -9,14 +9,12 @@
 
 // Implementation.
 
-#include "cl_DF.h"
-#include "cl_F.h"
-#include "cl_ieee.h"
-#include "cl_xmacros.h"
+#include "float/dfloat/cl_DF.h"
+#include "float/cl_F.h"
+#include "base/cl_xmacros.h"
 
 namespace cln {
 
-NEED_IEEE_FLOATS()
 
 const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
 {
@@ -27,11 +25,11 @@ const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
 // Also e1 >= e2.
 // Falls e1 - e2 >= 52 + 3, Ergebnis x1.
 // Schiebe beide Mantissen um 3 Bits nach links (Vorbereitung der Rundung:
-//   Bei e1-e2=0,1 ist keine Rundung nötig, bei e1-e2>1 ist der Exponent des
+//   Bei e1-e2=0,1 ist keine Rundung nötig, bei e1-e2>1 ist der Exponent des
 //   Ergebnisses =e1-1, =e1 oder =e1+1. Brauche daher 1 Schutzbit und zwei
-//   Rundungsbits: 00 exakt, 01 1.Hälfte, 10 exakte Mitte, 11 2.Hälfte.)
+//   Rundungsbits: 00 exakt, 01 1.Hälfte, 10 exakte Mitte, 11 2.Hälfte.)
 // Schiebe die Mantisse von x2 um e0-e1 Bits nach rechts. (Dabei die Rundung
-// ausführen: Bit 0 ist das logische Oder der Bits 0,-1,-2,...)
+// ausführen: Bit 0 ist das logische Oder der Bits 0,-1,-2,...)
 // Falls x1,x2 selbes Vorzeichen haben: Addiere dieses zur Mantisse von x1.
 // Falls x1,x2 verschiedenes Vorzeichen haben: Subtrahiere dieses von der
 //   Mantisse von x1. <0 -> (Es war e1=e2) Vertausche die Vorzeichen, negiere.
@@ -41,9 +39,9 @@ const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
 #ifdef FAST_DOUBLE
       double_to_DF(DF_to_double(x1) + DF_to_double(x2), return ,
                    TRUE, TRUE, // Overflow und subnormale Zahl abfangen
-                   FALSE, // kein Underflow mit Ergebnis +/- 0.0 möglich
+                   FALSE, // kein Underflow mit Ergebnis +/- 0.0 möglich
                           // (nach Definition der subnormalen Zahlen)
-                   FALSE, FALSE // keine Singularität, kein NaN als Ergebnis möglich
+                   FALSE, FALSE // keine Singularität, kein NaN als Ergebnis möglich
                   );
 #else
 #if (cl_word_size==64)
@@ -87,10 +85,10 @@ const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
         { mant1 = mant1 + mant2; }
       // mant1 = Ergebnis-Mantisse >0, sign1 = Ergebnis-Vorzeichen,
       // exp1 = Ergebnis-Exponent.
-      // Außerdem: Bei expdiff=0,1 sind die zwei letzten Bits von mant1 Null,
+      // Außerdem: Bei expdiff=0,1 sind die zwei letzten Bits von mant1 Null,
       // bei expdiff>=2 ist mant1 >= 2^(DF_mant_len+2).
       // Stets ist mant1 < 2^(DF_mant_len+5). (Daher werden die 2 Rundungsbits
-      // nachher um höchstens eine Position nach links geschoben werden.)
+      // nachher um höchstens eine Position nach links geschoben werden.)
       // [Knuth, S.201, leicht modifiziert:
       //   N1. m>=1 -> goto N4.
       //   N2. [Hier m<1] m>=1/2 -> goto N5.
@@ -129,8 +127,8 @@ const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
          // aufrunden
          { mant1 = mant1+1;
            if (mant1 >= bit(DF_mant_len+1))
-             // Bei Überlauf während der Rundung nochmals rechts schieben
-             // (Runden ist hier überflüssig):
+             // Bei Überlauf während der Rundung nochmals rechts schieben
+             // (Runden ist hier überflüssig):
              { mant1 = mant1>>1; exp1 = exp1+1; } // Mantisse rechts schieben
          }
       }// Runden fertig
@@ -205,10 +203,10 @@ const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
         }
       // mant1 = Ergebnis-Mantisse >0, sign1 = Ergebnis-Vorzeichen,
       // exp1 = Ergebnis-Exponent.
-      // Außerdem: Bei expdiff=0,1 sind die zwei letzten Bits von mant1 Null,
+      // Außerdem: Bei expdiff=0,1 sind die zwei letzten Bits von mant1 Null,
       // bei expdiff>=2 ist mant1 >= 2^(DF_mant_len+2).
       // Stets ist mant1 < 2^(DF_mant_len+5). (Daher werden die 2 Rundungsbits
-      // nachher um höchstens eine Position nach links geschoben werden.)
+      // nachher um höchstens eine Position nach links geschoben werden.)
       // [Knuth, S.201, leicht modifiziert:
       //   N1. m>=1 -> goto N4.
       //   N2. [Hier m<1] m>=1/2 -> goto N5.
@@ -252,8 +250,8 @@ const cl_DF operator+ (const cl_DF& x1, const cl_DF& x2)
            if (mantlo1==0)
              { manthi1 = manthi1+1;
                if (manthi1 >= bit(DF_mant_len-32+1))
-                 // Bei Überlauf während der Rundung nochmals rechts schieben
-                 // (Runden ist hier überflüssig):
+                 // Bei Überlauf während der Rundung nochmals rechts schieben
+                 // (Runden ist hier überflüssig):
                  { manthi1 = manthi1>>1; exp1 = exp1+1; } // Mantisse rechts schieben
          }   }
       }// Runden fertig