From 57a7520545134ab95f7bb0397dcbe991c906a3e9 Mon Sep 17 00:00:00 2001
From: Sebastiano Tronto <sebastiano@tronto.net>
Date: Wed, 6 Aug 2025 16:34:21 +0200
Subject: Added HTR solver

---
 src/arch/common.h | 39 +++++++++++++++++++++++++++++++++++++++
 1 file changed, 39 insertions(+)

(limited to 'src/arch/common.h')

diff --git a/src/arch/common.h b/src/arch/common.h
index d76b65f..f06f8d7 100644
--- a/src/arch/common.h
+++ b/src/arch/common.h
@@ -29,6 +29,8 @@ STATIC_INLINE uint64_t coord_cocsep(cube_t);
 STATIC_INLINE uint64_t coord_eo(cube_t);
 STATIC_INLINE uint64_t coord_esep(cube_t);
 STATIC_INLINE cube_t invcoord_esep(uint64_t);
+STATIC_INLINE uint64_t coord_epudsep(cube_t);
+STATIC_INLINE cube_t invcoord_epudsep(uint64_t);
 
 STATIC_INLINE bool is_eo_even(cube_t);
 
@@ -38,6 +40,8 @@ STATIC_INLINE void set_eo(cube_t [static 1], uint64_t);
 
 STATIC_INLINE void invcoord_esep_array(uint64_t, uint64_t, uint8_t[static 12]);
 STATIC_INLINE cube_t invcoord_eoesep(uint64_t);
+STATIC_INLINE uint64_t coord_epudsep_array(const uint8_t [8]);
+STATIC_INLINE void invcoord_epudsep_array(uint64_t, uint8_t [8]);
 
 STATIC_INLINE uint64_t coord_cp(cube_t);
 STATIC_INLINE cube_t invcoord_cp(uint64_t);
@@ -85,3 +89,38 @@ invcoord_eoesep(uint64_t i)
 
 	return c;
 }
+
+STATIC_INLINE uint64_t
+coord_epudsep_array(const uint8_t e[8])
+{
+	uint8_t i, k, is;
+	uint64_t ret;
+
+	ret = 0;
+	k = 4;
+	for (i = 0; i < 8; i++) {
+		is = (e[i] & UINT8_C(4)) >> UINT8_C(2);
+		ret += is * binomial[7-i][k];
+		k -= is;
+	}
+
+	return ret;
+}
+
+STATIC_INLINE void
+invcoord_epudsep_array(uint64_t c, uint8_t ret[8])
+{
+	uint8_t i, k, is, x, y;
+
+	k = 4;
+	x = 0;
+	y = 4;
+	for (i = 0; i < 8; i++) {
+		is = c >= binomial[7-i][k];
+		ret[i] = is*y + (1-is)*x;
+		y += is;
+		x += 1-is;
+		c -= is * binomial[7-i][k];
+		k -= is;
+	}
+}
-- 
cgit v1.3