src/float/conv/cl_DF_to_LF.cc

   1 // cl_DF_to_LF().
   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_LF.h"
  14 #include "cl_LF_impl.h"
  15 #include "cl_DS.h"
  16
  17 const cl_LF cl_DF_to_LF (const cl_DF& x, uintC len)
  18 {
  19         // x entpacken:
  20         var cl_signean sign;
  21         var sintL exp;
  22         #if (cl_word_size==64)
  23         var uint64 mant;
  24         DF_decode(x, { return encode_LF0(len); }, sign=,exp=(sintL),mant=);
  25         #else
  26         var uint32 manthi;
  27         var uint32 mantlo;
  28         DF_decode2(x, { return encode_LF0(len); }, sign=,exp=(sintL),manthi=,mantlo=);
  29         #endif
  30         // Long-Float allozieren,
  31         // Mantisse mit intDsize*len-DF_mant_len-1 Nullbits auffüllen:
  32         var Lfloat y = allocate_lfloat(len,exp+LF_exp_mid,sign);
  33         var uintD* ptr = arrayMSDptr(TheLfloat(y)->data,len);
  34         // erste k := ceiling(DF_mant_len+1,intDsize) Digits mit mant füllen:
  35         var const int shiftcount = (ceiling(DF_mant_len+1,intDsize)*intDsize-(DF_mant_len+1));
  36         #if (cl_word_size==64)
  37         mant = mant<<shiftcount;
  38         set_max64_Dptr(DF_mant_len+1,ptr,mant);
  39         #else
  40         manthi = (manthi<<shiftcount) | (mantlo>>(32-shiftcount));
  41         mantlo = mantlo<<shiftcount;
  42         set_max32_Dptr(DF_mant_len+1-32,ptr,manthi);
  43         set_32_Dptr(ptr mspop ceiling(DF_mant_len+1-32,intDsize),mantlo);
  44         #endif
  45         clear_loop_msp(ptr mspop ceiling(DF_mant_len+1,intDsize),len-ceiling(DF_mant_len+1,intDsize));
  46         return y;
  47 }