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
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
// Template header for 64 and 32 bit floating point numbers
// This file expects `FLOATX_H__BITS` to be defined before including:
//
// #define FLOATX_H__BITS 64
// #include "floatX.h"
//
#include <stdbool.h>
#include <stdint.h>
#include "datatypes.h"
#include "stdlib.h"
#include "types.h"
#include "util.h"
#ifndef FLOATX_H__BITS
#define FLOATX_H__BITS 64
#endif
#if FLOATX_H__BITS == 64
#define FLOAT_T double
#define OPT_T double
#define NAMESPACED(x) Float64$##x
#define TYPE_STR "Float64"
#define SUFFIXED(x) x
#elif FLOATX_H__BITS == 32
#define FLOAT_T float
#define OPT_T float
#define NAMESPACED(x) Float32$##x
#define TYPE_STR "Float32"
#define SUFFIXED(x) x##f
#else
#error "Unsupported bit width for Float"
#endif
Text_t NAMESPACED(as_text)(const void *x, bool colorize, const TypeInfo_t *type);
Text_t NAMESPACED(value_as_text)(FLOAT_T x);
PUREFUNC int32_t NAMESPACED(compare)(const void *x, const void *y, const TypeInfo_t *type);
PUREFUNC bool NAMESPACED(equal)(const void *x, const void *y, const TypeInfo_t *type);
CONSTFUNC bool NAMESPACED(near)(FLOAT_T a, FLOAT_T b, FLOAT_T ratio, FLOAT_T absolute);
Text_t NAMESPACED(percent)(FLOAT_T x, FLOAT_T precision);
FLOAT_T CONSTFUNC NAMESPACED(with_precision)(FLOAT_T n, FLOAT_T precision);
FLOAT_T NAMESPACED(mod)(FLOAT_T n, FLOAT_T modulus);
FLOAT_T NAMESPACED(mod1)(FLOAT_T n, FLOAT_T modulus);
CONSTFUNC bool NAMESPACED(isinf)(FLOAT_T n);
CONSTFUNC bool NAMESPACED(finite)(FLOAT_T n);
CONSTFUNC bool NAMESPACED(isnan)(FLOAT_T n);
bool NAMESPACED(is_none)(const void *n, const TypeInfo_t *info);
FLOAT_T NAMESPACED(nan)(Text_t tag);
CONSTFUNC FLOAT_T NAMESPACED(mix)(FLOAT_T amount, FLOAT_T x, FLOAT_T y);
OPT_T NAMESPACED(parse)(Text_t text, Text_t *remainder);
CONSTFUNC bool NAMESPACED(is_between)(const FLOAT_T x, const FLOAT_T low, const FLOAT_T high);
CONSTFUNC FLOAT_T NAMESPACED(clamped)(FLOAT_T x, FLOAT_T low, FLOAT_T high);
#if FLOATX_H__BITS == 64
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_float32)(float n) { return (FLOAT_T)n; }
#elif FLOATX_H__BITS == 32
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_float64)(double n) { return (FLOAT_T)n; }
#endif
#ifdef __GNUC__
#pragma GCC diagnostic push
#pragma GCC diagnostic ignored "-Wfloat-equal"
#endif
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_int)(Int_t i, bool truncate) {
if likely (i.small & 0x1) {
FLOAT_T ret = (FLOAT_T)(i.small >> 2);
if unlikely (!truncate && (int64_t)ret != (i.small >> 2))
fail("Could not convert integer to " TYPE_STR " without losing precision: ", i.small >> 2);
return ret;
} else {
FLOAT_T ret = mpz_get_d(i.big);
if (!truncate) {
mpz_t roundtrip;
mpz_init_set_d(roundtrip, (double)ret);
if unlikely (mpz_cmp(i.big, roundtrip) != 0)
fail("Could not convert integer to " TYPE_STR " without losing precision: ", i);
}
return ret;
}
}
#ifdef __GNUC__
#pragma GCC diagnostic pop
#endif
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_int64)(Int64_t i, bool truncate) {
FLOAT_T n = (FLOAT_T)i;
if unlikely (!truncate && (Int64_t)n != i)
fail("Could not convert integer to " TYPE_STR " without losing precision: ", i);
return n;
}
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_int32)(Int32_t i) { return (FLOAT_T)i; }
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_int16)(Int16_t i) { return (FLOAT_T)i; }
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_int8)(Int8_t i) { return (FLOAT_T)i; }
MACROLIKE CONSTFUNC FLOAT_T NAMESPACED(from_byte)(Byte_t i) { return (FLOAT_T)i; }
extern const TypeInfo_t NAMESPACED(info);
#undef FLOAT_T
#undef OPT_T
#undef NAMESPACED
#undef TYPE_STR
#undef SUFFIXED
#undef FLOATX_H__BITS
|
