src/float/conv/cl_DF_to_SF.cc

   1 // cl_DF_to_SF().
   2
   3 // General includes.
   4 #include "cl_sysdep.h"
   5
   6 // Specification.
   7 #include "cl_F.h"
   8
   9
  10 // Implementation.
  11
  12 #include "cl_DF.h"
  13 #include "cl_SF.h"
  14
  15 const cl_SF cl_DF_to_SF (const cl_DF& x)
  16 {
  17         // x entpacken:
  18         var cl_signean sign;
  19         var sintL exp;
  20         #if (cl_word_size==64)
  21         var uint64 mant;
  22         DF_decode(x, { return SF_0; }, sign=,exp=,mant=);
  23         // 52-16=36 Bits wegrunden:
  24         var const int shiftcount = DF_mant_len-SF_mant_len;
  25         if ( ((mant & bit(shiftcount-1)) ==0) // Bit 35 war 0 -> abrunden
  26              || ( ((mant & (bit(shiftcount-1)-1)) ==0) // war 1, Bits 34..0 >0 -> aufrunden
  27                   // round-to-even
  28                   && ((mant & bit(shiftcount)) ==0)
  29            )    )
  30           // abrunden
  31           { mant = mant >> shiftcount; }
  32           else
  33           // aufrunden
  34           { mant = mant >> shiftcount;
  35             mant = mant+1;
  36             if (mant >= bit(SF_mant_len+1))
  37               // Überlauf durchs Runden
  38               { mant = mant>>1; exp = exp+1; } // Mantisse rechts schieben
  39           }
  40         return encode_SF(sign,exp,mant);
  41         #else
  42         var uint32 manthi;
  43         var uint32 mantlo;
  44         DF_decode2(x, { return SF_0; }, sign=,exp=,manthi=,mantlo=);
  45         // 52-16=36 Bits wegrunden:
  46         var const int shiftcount = DF_mant_len-SF_mant_len-32;
  47         if ( ((manthi & bit(shiftcount-1)) ==0) // Bit 35 war 0 -> abrunden
  48              || ( ((manthi & (bit(shiftcount-1)-1)) ==0) // war 1, Bits 34..0 >0 -> aufrunden
  49                   && (mantlo==0)
  50                   // round-to-even
  51                   && ((manthi & bit(shiftcount)) ==0)
  52            )    )
  53           // abrunden
  54           { manthi = manthi >> shiftcount; }
  55           else
  56           // aufrunden
  57           { manthi = manthi >> shiftcount;
  58             manthi = manthi+1;
  59             if (manthi >= bit(SF_mant_len+1))
  60               // Überlauf durchs Runden
  61               { manthi = manthi>>1; exp = exp+1; } // Mantisse rechts schieben
  62           }
  63         return encode_SF(sign,exp,manthi);
  64         #endif
  65 }