src/float/conv/cl_DF_to_LF.cc

   1 // cl_DF_to_LF().
   2
   3 // General includes.
   4 #include "base/cl_sysdep.h"
   5
   6 // Specification.
   7 #include "float/cl_F.h"
   8
   9
  10 // Implementation.
  11
  12 #include "float/dfloat/cl_DF.h"
  13 #include "float/lfloat/cl_LF.h"
  14 #include "float/lfloat/cl_LF_impl.h"
  15 #include "base/digitseq/cl_DS.h"
  16
  17 namespace cln {
  18
  19 const cl_LF cl_DF_to_LF (const cl_DF& x, uintC len)
  20 {
  21         // x entpacken:
  22         var cl_signean sign;
  23         var sintL exp;
  24         #if (cl_word_size==64)
  25         var uint64 mant;
  26         DF_decode(x, { return encode_LF0(len); }, sign=,exp=(sintL),mant=);
  27         #else
  28         var uint32 manthi;
  29         var uint32 mantlo;
  30         DF_decode2(x, { return encode_LF0(len); }, sign=,exp=(sintL),manthi=,mantlo=);
  31         #endif
  32         // Long-Float allozieren,
  33         // Mantisse mit intDsize*len-DF_mant_len-1 Nullbits auffüllen:
  34         var Lfloat y = allocate_lfloat(len,exp+LF_exp_mid,sign);
  35         var uintD* ptr = arrayMSDptr(TheLfloat(y)->data,len);
  36         // erste k := ceiling(DF_mant_len+1,intDsize) Digits mit mant füllen:
  37         var const int shiftcount = (ceiling(DF_mant_len+1,intDsize)*intDsize-(DF_mant_len+1));
  38         #if (cl_word_size==64)
  39         mant = mant<<shiftcount;
  40         set_max64_Dptr(DF_mant_len+1,ptr,mant);
  41         #else
  42         manthi = (manthi<<shiftcount) | (mantlo>>(32-shiftcount));
  43         mantlo = mantlo<<shiftcount;
  44         set_max32_Dptr(DF_mant_len+1-32,ptr,manthi);
  45         set_32_Dptr(ptr mspop ceiling(DF_mant_len+1-32,intDsize),mantlo);
  46         #endif
  47         clear_loop_msp(ptr mspop ceiling(DF_mant_len+1,intDsize),len-ceiling(DF_mant_len+1,intDsize));
  48         return y;
  49 }
  50
  51 }  // namespace cln