huangcm
2025-02-28 b45e871a67cd1272e3da9ba5bd383f832b0f1824
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
//===----------------------------------------------------------------------===//
//
//                     The LLVM Compiler Infrastructure
//
// This file is dual licensed under the MIT and the University of Illinois Open
// Source Licenses. See LICENSE.TXT for details.
//
//===----------------------------------------------------------------------===//
 
// Define a hexfloat literal emulator since we can't depend on being able to
//   for hexfloat literals
 
// 0x10.F5p-10 == hexfloat<double>(0x10, 0xF5, -10)
 
#ifndef HEXFLOAT_H
#define HEXFLOAT_H
 
#include <cmath>
#include <climits>
 
template <class T>
class hexfloat
{
    T value_;
 
    static int CountLeadingZeros(unsigned long long n) {
        const std::size_t Digits = sizeof(unsigned long long) * CHAR_BIT;
        const unsigned long long TopBit = 1ull << (Digits - 1);
        if (n == 0) return Digits;
        int LeadingZeros = 0;
        while ((n & TopBit) == 0) {
            ++LeadingZeros;
            n <<= 1;
        }
        return LeadingZeros;
    }
 
public:
    hexfloat(long long m1, unsigned long long m0, int exp)
    {
        const std::size_t Digits = sizeof(unsigned long long) * CHAR_BIT;
        int s = m1 < 0 ? -1 : 1;
        int exp2 = -static_cast<int>(Digits - CountLeadingZeros(m0)/4*4);
        value_ = std::ldexp(m1 + s * std::ldexp(T(m0), exp2), exp);
    }
 
    operator T() const {return value_;}
};
 
#endif