4 #include "base/cl_sysdep.h"
7 #include "cln/sfloat.h"
12 #include "float/sfloat/cl_SF.h"
16 const cl_SF fround (const cl_SF& x)
19 // x = 0.0 oder e<0 -> Ergebnis 0.0
20 // 0<=e<=16 -> letzte (17-e) Bits der Mantisse wegrunden,
21 // Exponent und Vorzeichen beibehalten.
23 var uintL uexp = SF_uexp(x); // e + SF_exp_mid
24 if (uexp < SF_exp_mid) // x = 0.0 oder e<0 ?
27 { if (uexp > SF_exp_mid+SF_mant_len) // e > 16 ?
30 if (uexp > SF_exp_mid+1) // e>1 ?
31 { var cl_uint bitmask = // Bitmaske: Bit 16-e gesetzt, alle anderen gelöscht
32 bit(SF_mant_len+SF_mant_shift + SF_exp_mid-uexp);
33 var cl_uint mask = // Bitmaske: Bits 15-e..0 gesetzt, alle anderen gelöscht
34 bitmask - bit(SF_mant_shift);
35 if ( ((x.word & bitmask) ==0) // Bit 16-e =0 -> abrunden
36 || ( ((x.word & mask) ==0) // Bit 16-e =1 und Bits 15-e..0 >0 -> aufrunden
37 // round-to-even, je nach Bit 17-e :
38 && ((x.word & (bitmask<<1)) ==0)
41 { mask |= bitmask; // Bitmaske: Bits 16-e..0 gesetzt, alle anderen gelöscht
42 return cl_SF_from_word(x.word & ~mask);
46 { return cl_SF_from_word(
47 (x.word | mask) // alle diese Bits 15-e..0 setzen (Bit 16-e schon gesetzt)
48 + bit(SF_mant_shift) // letzte Stelle erhöhen, dabei evtl. Exponenten incrementieren
52 elif (uexp == SF_exp_mid+1) // e=1 ?
53 // Wie bei 1 < e <= 16, nur daß Bit 17-e stets gesetzt ist.
54 { if ((x.word & bit(SF_mant_len+SF_mant_shift-1)) ==0) // Bit 16-e =0 -> abrunden
56 { return cl_SF_from_word(x.word & ~(bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))); }
59 { return cl_SF_from_word(
60 (x.word | (bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))) // alle diese Bits 16-e..0 setzen
61 + bit(SF_mant_shift) // letzte Stelle erhöhen, dabei evtl. Exponenten incrementieren
66 // Wie bei 1 < e <= 16, nur daß Bit 16-e stets gesetzt
67 // und Bit 17-e stets gelöscht ist.
68 { if ((x.word & (bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))) ==0)
69 // abrunden von +-0.5 zu 0.0
73 { return cl_SF_from_word(
74 (x.word | (bit(SF_mant_len+SF_mant_shift)-bit(SF_mant_shift))) // alle Bits 15-e..0 setzen
75 + bit(SF_mant_shift) // letzte Stelle erhöhen, dabei Exponenten incrementieren