]> www.ginac.de Git - cln.git/blob - src/float/conv/cl_FF_to_LF.cc
Initial revision
[cln.git] / src / float / conv / cl_FF_to_LF.cc
1 // cl_FF_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_FF.h"
13 #include "cl_LF.h"
14 #include "cl_LF_impl.h"
15 #include "cl_DS.h"
16
17 const cl_LF cl_FF_to_LF (const cl_FF& x, uintC len)
18 {
19         // x entpacken:
20         var cl_signean sign;
21         var sintL exp;
22         #if (intDsize==64)
23         var uint64 mant;
24         #else
25         var uint32 mant;
26         #endif
27         FF_decode(x, { return encode_LF0(len); }, sign=,exp=(sintL),mant=);
28         // Long-Float allozieren,
29         // Mantisse mit intDsize*len-FF_mant_len-1 Nullbits auffüllen:
30         var Lfloat y = allocate_lfloat(len,exp+LF_exp_mid,sign);
31         var uintD* ptr = arrayMSDptr(TheLfloat(y)->data,len);
32         // erste k := ceiling(FF_mant_len+1,intDsize) Digits mit mant füllen:
33         mant = mant << (ceiling(FF_mant_len+1,intDsize)*intDsize-(FF_mant_len+1));
34         #if (intDsize==64)
35         set_max64_Dptr(FF_mant_len+1,ptr,mant);
36         #else
37         set_max32_Dptr(FF_mant_len+1,ptr,mant);
38         #endif
39         clear_loop_msp(ptr mspop ceiling(FF_mant_len+1,intDsize),len-ceiling(FF_mant_len+1,intDsize));
40         return y;
41 }