src/float/sfloat/elem/cl_SF_fround.cc

   1 // fround().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cln/sfloat.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_SF.h"
  13
  14 namespace cln {
  15
  16 const cl_SF fround (const cl_SF& x)
  17 {
  18 // Methode:
  19 // x = 0.0 oder e<0 -> Ergebnis 0.0
  20 // 0<=e<=16 -> letzte (17-e) Bits der Mantisse wegrunden,
  21 //             Exponent und Vorzeichen beibehalten.
  22 // e>16 -> Ergebnis x
  23       var uintL uexp = SF_uexp(x); // e + SF_exp_mid
  24       if (uexp < SF_exp_mid) // x = 0.0 oder e<0 ?
  25         { return SF_0; }
  26         else
  27         { if (uexp > SF_exp_mid+SF_mant_len) // e > 16 ?
  28             { return x; }
  29             else
  30             if (uexp > SF_exp_mid+1) // e>1 ?
  31               { var cl_uint bitmask = // Bitmaske: Bit 16-e gesetzt, alle anderen gelöscht
  32                   bit(SF_mant_len+SF_mant_shift + SF_exp_mid-uexp);
  33                 var cl_uint mask = // Bitmaske: Bits 15-e..0 gesetzt, alle anderen gelöscht
  34                   bitmask - bit(SF_mant_shift);
  35                 if ( ((x.word & bitmask) ==0) // Bit 16-e =0 -> abrunden
  36                      || ( ((x.word & mask) ==0) // Bit 16-e =1 und Bits 15-e..0 >0 -> aufrunden
  37                           // round-to-even, je nach Bit 17-e :
  38                           && ((x.word & (bitmask<<1)) ==0)
  39                    )    )
  40                   // abrunden
  41                   { mask |= bitmask; // Bitmaske: Bits 16-e..0 gesetzt, alle anderen gelöscht
  42                     return cl_SF_from_word(x.word & ~mask);
  43                   }
  44                   else
  45                   // aufrunden
  46                   { return cl_SF_from_word(
  47                       (x.word | mask) // alle diese Bits 15-e..0 setzen (Bit 16-e schon gesetzt)
  48                       + bit(SF_mant_shift) // letzte Stelle erhöhen, dabei evtl. Exponenten incrementieren
  49                       );
  50                   }
  51               }
  52             elif (uexp == SF_exp_mid+1) // e=1 ?
  53               // Wie bei 1 < e <= 16, nur daß Bit 17-e stets gesetzt ist.
  54               { if ((x.word & bit(SF_mant_len+SF_mant_shift-1)) ==0) // Bit 16-e =0 -> abrunden
  55                   // abrunden
  56                   { return cl_SF_from_word(x.word & ~(bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))); }
  57                   else
  58                   // aufrunden
  59                   { return cl_SF_from_word(
  60                       (x.word | (bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))) // alle diese Bits 16-e..0 setzen
  61                       + bit(SF_mant_shift) // letzte Stelle erhöhen, dabei evtl. Exponenten incrementieren
  62                       );
  63                   }
  64               }
  65             else // e=0 ?
  66               // Wie bei 1 < e <= 16, nur daß Bit 16-e stets gesetzt
  67               // und Bit 17-e stets gelöscht ist.
  68               { if ((x.word & (bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))) ==0)
  69                   // abrunden von +-0.5 zu 0.0
  70                   { return SF_0; }
  71                   else
  72                   // aufrunden
  73                   { return cl_SF_from_word(
  74                       (x.word | (bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))) // alle Bits 15-e..0 setzen
  75                       + bit(SF_mant_shift) // letzte Stelle erhöhen, dabei Exponenten incrementieren
  76                       );
  77               }   }
  78         }
  79 }
  80
  81 }  // namespace cln