src/float/conv/cl_FF_to_float.cc

   1 // cl_FF_to_float().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_ffloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_FF.h"
  13
  14 float cl_float_approx (const cl_FF& obj)
  15 {
  16         union { ffloat eksplicit; float machine_float; } u;
  17         #define val u.eksplicit
  18         val = cl_ffloat_value(obj);
  19         // Der Exponent muß um FF_exp_mid-126 erniedrigt werden.
  20         if (FF_exp_mid>126)
  21           { var uintL exp = (val >> FF_mant_len) & (bit(FF_exp_len)-1); // e
  22             if (exp < FF_exp_mid-126+1)
  23               { // produziere denormalisiertes Float
  24                 val = (val & minus_bit(FF_exp_len+FF_mant_len)) // selbes Vorzeichen
  25                       | (0 << FF_mant_len) // Exponent 0
  26                       | (((val & (bit(FF_mant_len)-1)) | bit(FF_mant_len)) // Mantisse shiften
  27                          >> (FF_exp_mid-126+1 - exp) // shiften
  28                         );
  29               }
  30               else
  31               { val -= (FF_exp_mid - 126) << FF_mant_len; }
  32           }
  33         #undef val
  34         return u.machine_float;
  35 }