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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
|
/* See LICENSE file for copyright and license details. */
#include "common.h"
/* Disable warnings in <math.h> */
#if defined(__clang__)
# pragma clang diagnostic ignored "-Wdouble-promotion"
# pragma clang diagnostic ignored "-Wconversion"
#endif
USAGE("[-s]")
#define CONV(ITYPE, OTYPE, SOTYPE, EXPONENT, HA2EXPONENT, FRACTION)\
do {\
static int cache_i = 0;\
static ITYPE cache_in[] = {0, 0, 0, 0};\
static OTYPE cache_out[] = {0, 0, 0, 0};\
OTYPE signb, fraction, ret;\
SOTYPE exponent;\
ITYPE u, dexponent;\
if (host == cache_in[cache_i]) {\
ret = cache_out[cache_i++];\
cache_i &= 3;\
return ret;\
}\
cache_in[cache_i] = host;\
signb = (OTYPE)signbit(host);\
u = signb ? -host : host;\
if (g_isnan(host) || g_isinf(host)) {\
ret = ((((OTYPE)1 << EXPONENT) - (OTYPE)1) << FRACTION) | (OTYPE)!g_isinf(host);\
} else if (u == (ITYPE)0.0) {\
ret = 0;\
} else {\
dexponent = log2(u);\
exponent = (SOTYPE)dexponent;\
if (u == pow((ITYPE)2.0, (ITYPE)exponent)) {\
exponent += HA2EXPONENT;\
fraction = 0;\
} else {\
/* TODO subnormals are a bit rounded off */\
u *= pow((ITYPE)2.0, (ITYPE)(FRACTION + 1 - exponent));\
fraction = (OTYPE)u;\
while (fraction >= (OTYPE)2 << FRACTION) {\
fraction >>= 1;\
exponent += 1;\
}\
fraction &= ((OTYPE)1 << FRACTION) - 1;\
exponent += HA2EXPONENT - 1;\
}\
if (exponent < 1) {\
/* TODO subnormal result */\
exponent = 0;\
fraction = 0;\
} else if (exponent >= ((SOTYPE)1 << EXPONENT) - 1) { \
exponent = ((SOTYPE)1 << EXPONENT) - 1;\
fraction = 0;\
}\
ret = ((OTYPE)exponent << FRACTION) + fraction;\
}\
ret |= signb << (FRACTION + EXPONENT);\
cache_out[cache_i++] = ret;\
cache_i &= 3;\
return ret;\
} while (0)
#define PROCESS(ITYPE, OTYPE, BITS)\
do {\
size_t i, n;\
ITYPE *ibuf = (ITYPE *)(stream->buf);\
OTYPE *obuf = sizeof(ITYPE) == sizeof(OTYPE) ? (OTYPE *)(stream->buf)\
: alloca(sizeof(stream->buf) / sizeof(ITYPE) * sizeof(OTYPE));\
strict *= !USING_BINARY##BITS;\
if (!strict && sizeof(ITYPE) != sizeof(OTYPE))\
eprintf("-s is required on this machine\n");\
do {\
n = stream->ptr / sizeof(ITYPE);\
if (strict) {\
for (i = 0; i < n; i++)\
obuf[i] = htole##BITS(conv_##ITYPE(ibuf[i]));\
} else {\
for (i = 0; i < n; i++)\
obuf[i] = htole##BITS(*(OTYPE *)&ibuf[i]);\
}\
ewriteall(STDOUT_FILENO, obuf, n * sizeof(OTYPE), "<stdout>");\
n *= sizeof(ITYPE);\
memmove(stream->buf, stream->buf + n, stream->ptr -= n);\
} while (eread_stream(stream, SIZE_MAX));\
if (stream->ptr)\
eprintf("%s: incomplete frame\n", stream->file);\
} while (0)
static uint64_t conv_double(double host) {CONV(double, uint64_t, int64_t, 11, 1023, 52);}
static uint32_t conv_float (float host) {CONV(float, uint32_t, int32_t, 8, 127, 23);}
static void process_lf(struct stream *stream, int strict) {PROCESS(double, uint64_t, 64);}
static void process_f (struct stream *stream, int strict) {PROCESS(float, uint32_t, 32);}
int
main(int argc, char *argv[])
{
struct stream stream;
int strict = 1;
void (*process)(struct stream *stream, int strict);
ARGBEGIN {
case 's':
strict = 0;
break;
default:
usage();
} ARGEND;
if (argc)
usage();
eopen_stream(&stream, NULL);
SELECT_PROCESS_FUNCTION(&stream);
fprint_stream_head(stdout, &stream);
efflush(stdout, "<stdout>");
process(&stream, strict);
return 0;
}
|
