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

3 files changed, 52 insertions, 0 deletions

diff --git a/test/110_coord_invcoord_h48/00_all.in b/test/110_coord_invcoord_h48/00_all.in
new file mode 100644
index 0000000..a2f9540
--- /dev/null
+++ b/test/110_coord_invcoord_h48/00_all.in
@@ -0,0 +1,8 @@
+UB0 DF0 DB0 UF0 UR0 UL0 DL0 DR0 FR0 FL0 BL0 BR0 UFR0 UBL0 DFL0 DBR0 UFL0 UBR0 DFR0 DBL0
+UF0 UB0 DB0 DF0 UR0 UL0 DL0 DR0 FR0 FL0 BL0 BR0 UFR0 UBL0 DFL0 DBR0 UFL0 UBR0 DFR0 DBL0
+UB0 UF0 DB0 DF0 UL0 UR0 DL0 DR0 FR0 FL0 BL0 BR0 UBL0 UFR0 DFL0 DBR0 UBR0 UFL0 DFR0 DBL0
+UL0 BL0 BR1 DL0 FR0 DF0 DB1 DR1 UB0 FL0 UF0 UR1 DFL0 UFR1 DBR1 UBR2 DBL2 DFR0 UFL1 UBL2
+UF0 UB0 DB0 DF0 FR0 UL0 DL0 BR0 DR0 FL0 BL0 UR0 DFR2 UBL0 DFL0 UBR2 UFL0 UFR1 DBR1 DBL0
+FL1 UB0 DB0 FR1 UR0 UL0 DL0 DR0 UF1 DF1 BL0 BR0 UFL1 UBL0 DFR1 DBR0 DFL2 UBR0 UFR2 DBL0
+FR0 BR0 BL0 FL0 UR0 UL0 DL0 DR0 DF0 UF0 UB0 DB0 UFR0 UBL0 DFL0 DBR0 UFL0 UBR0 DFR0 DBL0
+UR0 UL0 DL0 DR0 UB0 UF0 DF0 DB0 FR0 FL0 BL0 BR0 UFR0 UBL0 DFL0 DBR0 UFL0 UBR0 DFR0 DBL0
diff --git a/test/110_coord_invcoord_h48/00_all.out b/test/110_coord_invcoord_h48/00_all.out
new file mode 100644
index 0000000..e84be4b
--- /dev/null
+++ b/test/110_coord_invcoord_h48/00_all.out
@@ -0,0 +1,8 @@
+1 ok
+2 ok
+3 ok
+4 ok
+5 ok
+6 ok
+7 ok
+8 ok
diff --git a/test/110_coord_invcoord_h48/coord_invcoord_h48_tests.c b/test/110_coord_invcoord_h48/coord_invcoord_h48_tests.c
new file mode 100644
index 0000000..77dea81
--- /dev/null
+++ b/test/110_coord_invcoord_h48/coord_invcoord_h48_tests.c
@@ -0,0 +1,36 @@
+#include "../test.h"
+
+#define COCSEP_CLASSES 3393
+
+size_t gendata_cocsep(void *, uint64_t *, cube_t *);
+int64_t coord_h48(cube_t, const uint32_t *, uint8_t);
+cube_t invcoord_h48(int64_t, const cube_t *, uint8_t);
+cube_t transform(cube_t, uint8_t);
+
+void run(void) {
+        char str[STRLENMAX];
+        int i;
+        bool found;
+        uint8_t h, t;
+        uint32_t cocsepdata[300000];
+        uint64_t selfsim[COCSEP_CLASSES];
+        int64_t c, cc;
+        cube_t cube, invc, rep[COCSEP_CLASSES];
+
+        gendata_cocsep(cocsepdata, selfsim, rep);
+
+        i = 1;
+        h = 11;
+        while (fgets(str, STRLENMAX, stdin) != NULL) {
+                cube = readcube("H48", str);
+                c = coord_h48(cube, cocsepdata, h);
+                invc = invcoord_h48(c, rep, h);
+                for (t = 0, found = false; t < 48; t++) {
+                        cube = transform(invc, t);
+                        cc = coord_h48(cube, cocsepdata, h);
+                        found = found || cc == c;
+                }
+                printf("%d %s\n", i, found ? "ok" : "ERROR");
+                i++;
+        }
+}