]> www.ginac.de Git - cln.git/blob - src/float/conv/cl_FF_to_LF.cc
Use paths relative the `src' directory in the #include directives.
[cln.git] / src / float / conv / cl_FF_to_LF.cc
1 // cl_FF_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/ffloat/cl_FF.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_FF_to_LF (const cl_FF& x, uintC len)
20 {
21         // x entpacken:
22         var cl_signean sign;
23         var sintL exp;
24         #if (intDsize==64)
25         var uint64 mant;
26         #else
27         var uint32 mant;
28         #endif
29         FF_decode(x, { return encode_LF0(len); }, sign=,exp=(sintL),mant=);
30         // Long-Float allozieren,
31         // Mantisse mit intDsize*len-FF_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(FF_mant_len+1,intDsize) Digits mit mant füllen:
35         mant = mant << (ceiling(FF_mant_len+1,intDsize)*intDsize-(FF_mant_len+1));
36         #if (intDsize==64)
37         set_max64_Dptr(FF_mant_len+1,ptr,mant);
38         #else
39         set_max32_Dptr(FF_mant_len+1,ptr,mant);
40         #endif
41         clear_loop_msp(ptr mspop ceiling(FF_mant_len+1,intDsize),len-ceiling(FF_mant_len+1,intDsize));
42         return y;
43 }
44
45 }  // namespace cln