]> www.ginac.de Git - cln.git/blob - src/float/conv/cl_SF_to_float.cc
46087a6901dc0389af7824dba6c80a8691a0de4d
[cln.git] / src / float / conv / cl_SF_to_float.cc
1 // cl_SF_to_float().
2
3 // General includes.
4 #include "base/cl_sysdep.h"
5
6 // Specification.
7 #include "cln/sfloat.h"
8
9
10 // Implementation.
11
12 #include "float/sfloat/cl_SF.h"
13 #include "float/ffloat/cl_FF.h"
14
15 namespace cln {
16
17 float float_approx (const cl_SF& x)
18 {
19         // x entpacken:
20         var cl_signean sign;
21         var sintL exp;
22         var uint32 mant;
23         SF_decode(x, { return 0.0; }, sign=,exp=,mant=);
24         // Mantisse um 23-16=7 Bits nach links schieben:
25         union { ffloat eksplicit; float machine_float; } u;
26         if (((sintL)(SF_exp_high-SF_exp_mid) > (sintL)(FF_exp_high-FF_exp_mid))
27             && (exp > (sintL)(FF_exp_high-FF_exp_mid)))
28           { u.eksplicit = make_FF_word(sign,bit(FF_exp_len)-1,0); } // Infinity
29         else
30         if (((sintL)(SF_exp_low-SF_exp_mid) < (sintL)(FF_exp_low-FF_exp_mid))
31             && (exp < (sintL)(FF_exp_low-FF_exp_mid)))
32           { u.eksplicit = make_FF_word(sign,0,0); } // 0.0
33         else
34           { u.eksplicit = make_FF_word(sign,exp+FF_exp_mid,mant<<(FF_mant_len-SF_mant_len)); }
35         return u.machine_float;
36 }
37
38 }  // namespace cln