DR coordinate solver

4 files changed, 61 insertions, 1 deletions

diff --git a/test/121_coorddata_dr/00_all.in b/test/121_coorddata_dr/00_all.in
new file mode 100644
index 0000000..e69de29
--- /dev/null
+++ b/test/121_coorddata_dr/00_all.in
diff --git a/test/121_coorddata_dr/00_all.out b/test/121_coorddata_dr/00_all.out
new file mode 100644
index 0000000..d6109db
--- /dev/null
+++ b/test/121_coorddata_dr/00_all.out
@@ -0,0 +1 @@
+All good
diff --git a/test/121_coorddata_dr/coorddata_dr.c b/test/121_coorddata_dr/coorddata_dr.c
new file mode 100644
index 0000000..92e7c2a
--- /dev/null
+++ b/test/121_coorddata_dr/coorddata_dr.c
@@ -0,0 +1,59 @@
+#include "../test.h"
+
+#define POW_3_7 2187
+#define BOUND   (64430 * POW_3_7 / 100)
+#define TGROUP  UINT64_C(4278190335)
+
+cube_t transform(cube_t, uint8_t);
+uint64_t coordinate_dr_coord(cube_t, const void *);
+cube_t coordinate_dr_cube(uint64_t, const void *);
+uint64_t coordinate_dr_gendata(void *);
+
+void run(void) {
+        bool found;
+        uint64_t t;
+        char str[STRLENMAX];
+        void *data;
+        size_t size;
+        cube_t cube;
+        uint64_t coord, coord2;
+
+        size = coordinate_dr_gendata(NULL);
+        data = malloc(size);
+        coordinate_dr_gendata(data);
+
+        /* Test all possible values for CO coordinate */
+        for (coord = 0; coord < BOUND; coord++) {
+                cube = coordinate_dr_cube(coord, data);
+
+                if (!isconsistent(cube)) {
+                        printf("Error: invcoord of %" PRId64
+                            " is not consistent\n", coord);
+                        goto cleanup;
+                }
+
+                for (t = 0, found = false; t < 48; t++) {
+                        if (!((UINT64_C(1) << t) & TGROUP))
+                                continue;
+
+                        coord2 = coordinate_dr_coord(transform(cube, t), data);
+                        if (coord == coord2) {
+                                found = true;
+                                break;
+                        }
+                }
+
+                if (!found) {
+                        printf("Error: invcoord of %" PRId64 " returns %"
+                            PRId64 " with cube:\n", coord, coord2);
+                        writecube("H48", cube, STRLENMAX, str);
+                        printf("%s\n", str);
+                        goto cleanup;
+                }
+        }
+
+        printf("All good\n");
+
+cleanup:
+        free(data);
+}
diff --git a/test/test.h b/test/test.h
index e753977..615402b 100644
--- a/test/test.h
+++ b/test/test.h
@@ -26,7 +26,7 @@ bool isconsistent(cube_t);
 bool issolvable(cube_t);
 bool issolved(cube_t);
 cube_t readcube(char *, char *);
-int64_t writecube(char *, cube_t, size_t n, char [n]);
+int64_t writecube(const char *, cube_t, size_t n, char [n]);
 
 /* Test function to be implemented by all tests */
 void run(void);