Found a sneaky bug (it's all a mess now)

3 files changed, 33 insertions, 0 deletions

diff --git a/test/103_cocsep_ttrep/00_all.in b/test/103_cocsep_ttrep/00_all.in
new file mode 100644
index 0000000..e69de29
--- /dev/null
+++ b/test/103_cocsep_ttrep/00_all.in
diff --git a/test/103_cocsep_ttrep/00_all.out b/test/103_cocsep_ttrep/00_all.out
new file mode 100644
index 0000000..e69de29
--- /dev/null
+++ b/test/103_cocsep_ttrep/00_all.out
diff --git a/test/103_cocsep_ttrep/cocsep_ttrep_tests.c b/test/103_cocsep_ttrep/cocsep_ttrep_tests.c
new file mode 100644
index 0000000..78f8014
--- /dev/null
+++ b/test/103_cocsep_ttrep/cocsep_ttrep_tests.c
@@ -0,0 +1,33 @@
+#include "../test.h"
+
+#define COCSEP_CLASSES 3393
+
+uint8_t inverse_trans(uint8_t);
+cube_t transform_corners(cube_t);
+int64_t coord_cocsep(cube_t);
+size_t gendata_cocsep(void *, uint64_t *, cube_t *);
+
+void run(void) {
+        uint8_t t, tinv;
+        uint32_t buf[300000], tt;
+        uint64_t selfsim[COCSEP_CLASSES];
+        int64_t i, j;
+        cube_t rep[COCSEP_CLASSES], c, d;
+
+        gendata_cocsep(buf, selfsim, rep);
+
+        for (i = 0; i < COCSEP_CLASSES; i++) {
+                c = rep[i];
+                for (t = 0; t < 48; t++) {
+                        tinv = inverse_trans(t);
+                        d = transform_corners(c);
+                        j = coord_cocsep(d);
+                        tt = (buf[j] & (0xFF << 8)) >> 8;
+                        if (tt != tinv)
+                                printf("cocsep %" PRId64 " <- %" PRId64 ": "
+                                    "expected t %" PRIu8 " (inverse of %"
+                                    PRIu8 "), got %" PRIu32 "\n",
+                                    i, j, tinv, t, tt);
+                }
+        }
+}