src/float/dfloat/elem/cl_DF_minus.cc

   1 // binary operator -
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/dfloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_DF.h"
  13
  14 namespace cln {
  15
  16
  17 const cl_DF operator- (const cl_DF& x1, const cl_DF& x2)
  18 {
  19 // Methode:
  20 // (- x1 x2) = (+ x1 (- x2))
  21 #ifdef FAST_DOUBLE
  22       double_to_DF(DF_to_double(x1) - DF_to_double(x2), return ,
  23                    TRUE, TRUE, // Overflow und subnormale Zahl abfangen
  24                    FALSE, // kein Underflow mit Ergebnis +/- 0.0 möglich
  25                           // (nach Definition der subnormalen Zahlen)
  26                    FALSE, FALSE // keine Singularität, kein NaN als Ergebnis möglich
  27                   );
  28 #else
  29 #if (cl_word_size==64)
  30       var dfloat x2_ = TheDfloat(x2)->dfloat_value;
  31       if (DF_uexp(x2_) == 0)
  32         { return x1; }
  33         else
  34         { return x1 + allocate_dfloat(x2_ ^ bit(63)); }
  35 #else
  36       var uint32 x2_semhi = TheDfloat(x2)->dfloat_value.semhi;
  37       var uint32 x2_mlo = TheDfloat(x2)->dfloat_value.mlo;
  38       if (DF_uexp(x2_semhi) == 0)
  39         { return x1; }
  40         else
  41         { return x1 + allocate_dfloat(x2_semhi ^ bit(31), x2_mlo); }
  42 #endif
  43 #endif
  44 }
  45
  46 }  // namespace cln