]> git.tdb.fi Git - libs/datafile.git/blob - source/binfloat.h
Use <cstdint> in BinFloat
[libs/datafile.git] / source / binfloat.h
1 #ifndef MSP_DATAFILE_BINFLOAT_H_
2 #define MSP_DATAFILE_BINFLOAT_H_
3
4 #include <climits>
5 #include <cstdint>
6 #include <cstring>
7
8 namespace Msp {
9 namespace DataFile {
10
11 /**
12 Facilitates splitting floating-point numbers into parts and putting them back
13 together.  Supports arbitary sizes up to 64 bits.  The 16, 32 and 64 bit
14 formats exactly match those defined by ISO/IEC 60559:2011.
15
16 The exponent is stored in an unbiased form.  The mantissa is stored with the
17 integer part included, aligned to the high bits of a 64-bit integer.
18 */
19 struct BinFloat
20 {
21         struct Bits
22         {
23                 unsigned exponent;
24                 unsigned mantissa;
25
26                 Bits(unsigned);
27         };
28
29         template<typename T>
30         struct MatchingInt;
31
32         template<typename T>
33         union Conversion;
34
35         bool sign;
36         bool infinity;
37         int exponent;
38         std::uint64_t mantissa;
39
40         static BinFloat explode(std::uint64_t, const Bits &);
41
42         template<typename T>
43         static BinFloat explode_iec559(T v);
44
45         std::uint64_t compose(const Bits &);
46
47         template<typename T>
48         T compose_iec559();
49 };
50
51 template<>
52 struct BinFloat::MatchingInt<float> { typedef std::uint32_t Type; };
53
54 template<>
55 struct BinFloat::MatchingInt<double> { typedef std::uint64_t Type; };
56
57 template<typename T>
58 union BinFloat::Conversion
59 {
60         T f;
61         typename MatchingInt<T>::Type i;
62 };
63
64 template<typename T>
65 inline BinFloat BinFloat::explode_iec559(T v)
66 {
67         Conversion<T> c;
68         c.f = v;
69         return explode(c.i, sizeof(T)*CHAR_BIT);
70 }
71
72 template<typename T>
73 inline T BinFloat::compose_iec559()
74 {
75         Conversion<T> c;
76         c.i = compose(sizeof(T)*CHAR_BIT);
77         return c.f;
78 }
79
80 } // namespace DataFile
81 } // namespace Msp
82
83 #endif