1 files changed, 128 insertions, 0 deletions
diff --git a/src/blind-to-portable.c b/src/blind-to-portable.c
new file mode 100644
index 0000000..9cb0448
--- /dev/null
+++ b/src/blind-to-portable.c
@@ -0,0 +1,128 @@
+/* See LICENSE file for copyright and license details. */
+#include "stream.h"
+#include "util.h"
+
+#include <alloca.h>
+#include <math.h>
+#include <string.h>
+
+USAGE("[-s]")
+
+#define USING_BINARY32 0
+#define USING_BINARY64 0
+
+#define CONV(ITYPE, OTYPE, SOTYPE, EXPONENT, HA2EXPONENT, FRACTION, SUFFIX)\
+	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;\
+		}\
+		signb = signbit(host);\
+		u = signb ? -host : host;\
+		if (isnan(host) || !isfinite(host)) {\
+			ret = ((((OTYPE)1 << EXPONENT) - 1) << FRACTION) | !isinf(host);\
+		} else if (u == (ITYPE)0.0) {\
+			ret = 0;\
+		} else {\
+			dexponent = log2##SUFFIX(u);\
+			exponent = (SOTYPE)dexponent;\
+			if (u == pow##SUFFIX(2.0, (ITYPE)exponent)) {\
+				exponent += HA2EXPONENT;\
+				fraction = 0;\
+			} else {\
+				/* TODO subnormals are a bit rounded off */\
+				u *= pow##SUFFIX(2.0, (ITYPE)(FRACTION + 1 - exponent));\
+				fraction = 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) - (SOTYPE)1) { \
+				exponent = ((SOTYPE)1 << EXPONENT) - (SOTYPE)1;\
+				fraction = 0;\
+			}\
+			ret = (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, f);}
+
+static void process_xyza (struct stream *stream, int strict) {PROCESS(double, uint64_t, 64);}
+static void process_xyzaf(struct stream *stream, int strict) {PROCESS(float, uint32_t, 32);}
+
+int
+main(int argc, char *argv[])
+{
+	struct stream stream;
+	int strict = 0;
+	void (*process)(struct stream *stream, int strict);
+
+	ARGBEGIN {
+	case 's':
+		strict = 1;
+		break;
+	default:
+		usage();
+	} ARGEND;
+	if (argc)
+		usage();
+
+	eopen_stream(&stream, NULL);
+
+	if (!strcmp(stream.pixfmt, "xyza"))
+		process = process_xyza;
+	else if (!strcmp(stream.pixfmt, "xyza f"))
+		process = process_xyzaf;
+	else
+		eprintf("pixel format %s is not supported\n", stream.pixfmt);
+
+	fprint_stream_head(stdout, &stream);
+	efflush(stdout, "<stdout>");
+	process(&stream, strict);
+	return 0;
+}