nissy-classic

commit 374153308ceae1a594408ec5d6e770e28e8867cd
parent 0304953af37005488c40f81bde0eae8d802e0a53
Author: Sebastiano Tronto <sebastiano.tronto@gmail.com>
Date:   Fri, 12 Nov 2021 11:25:36 +0100

Fixed anti-index for eofbepos. It did not compute a value for eorl, which
causes problems when using this coordinate combined with symmetries (e.g.
in symcoord khuge, drud and similar). This actually reverts a change that I
made before the first commit of v2.

Diffstat:
A
nissy.o
 | 
0

M
src/coord.c
 | 
29
++++++++++++++++++++++++++++-

2 files changed, 28 insertions(+), 1 deletion(-)
diff --git a/nissy.o b/nissy.o
Binary files differ.
diff --git a/src/coord.c b/src/coord.c
@@ -153,7 +153,13 @@ coord_drud_eofb = {
 static Cube
 antindex_eofb(uint64_t ind)
 {
-	return (Cube){ .eofb = ind, .eorl = ind, .eoud = ind };
+	Cube ret = {0};
+
+	ret.eofb = ind;
+	ret.eorl = ind;
+	ret.eoud = ind;
+
+	return ret;
 }
 
 static Cube
@@ -161,8 +167,29 @@ antindex_eofbepos(uint64_t ind)
 {
 	Cube ret = {0};
 
+	/* We need eorl for sym16 coordinate */
+	static int initialized = false;
+	static uint64_t eorl_aux[POW2TO11][BINOM12ON4];
+	static int eo_aux[12], ep_aux[12];
+	unsigned int i, j, k;
+
+	if (!initialized) {
+		for (i = 0; i < POW2TO11; i++) {
+			for (j = 0; j < BINOM12ON4; j++) {
+				int_to_sum_zero_array(i, 2, 12, eo_aux);
+				index_to_subset(j, 12, 4, ep_aux);
+				for (k = 0; k < 12; k++)
+					if (ep_aux[k])
+						eo_aux[k] = 1 - eo_aux[k];
+			}
+		}
+
+		initialized = true;
+	}
+
 	ret.eofb = ind % POW2TO11;
 	ret.epose = (ind / POW2TO11) * 24;
+	ret.eorl = eorl_aux[ret.eofb][ret.epose/24];
 
 	return ret;
 }