Added CPEPE coordinate in preparation of full DRFIN solver

3 files changed, 33 insertions, 0 deletions

diff --git a/test/082_invcoord_epe/00_all.in b/test/082_invcoord_epe/00_all.in
new file mode 100644
index 0000000..e69de29
--- /dev/null
+++ b/test/082_invcoord_epe/00_all.in
diff --git a/test/082_invcoord_epe/00_all.out b/test/082_invcoord_epe/00_all.out
new file mode 100644
index 0000000..d6109db
--- /dev/null
+++ b/test/082_invcoord_epe/00_all.out
@@ -0,0 +1 @@
+All good
diff --git a/test/082_invcoord_epe/invcoord_epe_tests.c b/test/082_invcoord_epe/invcoord_epe_tests.c
new file mode 100644
index 0000000..69c3c50
--- /dev/null
+++ b/test/082_invcoord_epe/invcoord_epe_tests.c
@@ -0,0 +1,32 @@
+#include "../test.h"
+
+#define FACT_4 24
+
+uint64_t coord_epe(cube_t);
+cube_t invcoord_epe(uint64_t);
+
+void run(void) {
+        oriented_cube_t cube;
+        uint64_t coord, coord2;
+
+        cube.orientation = 0;
+
+        /* Test all possible values for CP coordinate */
+        for (coord = 0; coord < FACT_4; coord++) {
+                cube.cube = invcoord_epe(coord);
+
+                if (!isconsistent(cube)) {
+                        printf("Not consistent\n");
+                        return;
+                }
+
+                coord2 = coord_epe(cube.cube);
+                if (coord != coord2) {
+                        printf("Error: invcoord of %" PRIu64
+                            " returns %" PRIu64 "\n", coord, coord2);
+                        return;
+                }
+        }
+
+        printf("All good\n");
+}