1 files changed, 125 insertions, 2 deletions
diff --git a/libj2_ji_rsh_to_j2i.c b/libj2_ji_rsh_to_j2i.c
index 128cec0..0982e91 100644
--- a/libj2_ji_rsh_to_j2i.c
+++ b/libj2_ji_rsh_to_j2i.c
@@ -8,7 +8,130 @@ extern inline void libj2_ji_rsh_to_j2i(intmax_t a, unsigned b, struct libj2_j2i
 
 #else
 
-CONST int main(void) { return 0; }
-/* TODO test libj2_ji_rsh_to_j2i, libj2_ji_rsh_to_j2i_underflow, libj2_ji_rsh_underflow_p */
+static const char *patterns[] = {
+	"0",
+	"1111111111111111111111111111111111111111111111111111111111111111",
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111",
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111",
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111"
+	"1111111111111111111111111111111111111111111111111111111111111111",
+	"1010101010101001010100100010101010101011000101101011101001000110",
+	"1110100100100010000010111111011101101001000010001000011110101010",
+	"11111111111100001001011",
+	"0010011011001",
+	"110010001011110100101001010011010111101",
+	NULL
+};
+
+
+static int
+set(struct libj2_j2i *a, const char *bits, size_t leading_signbits)
+{
+	size_t leading_ones = leading_signbits;
+	size_t i, n;
+	int punderflow = 0;
+	int nunderflow = 0;
+	int negative = 0;
+
+	a->high = 0;
+	a->low = 0;
+
+	n = strlen(bits);
+	if (n > LIBJ2_JU_BIT) {
+		bits = &bits[n - LIBJ2_JU_BIT];
+		n = LIBJ2_JU_BIT;
+	}
+
+	if (n == LIBJ2_JU_BIT)
+		negative = bits[0] == '1';
+
+	while (n && leading_signbits--) {
+		--n;
+		punderflow |= bits[n] == '1';
+		nunderflow |= bits[n] == '0';
+	}
+
+	for (i = 0; n && i < LIBJ2_JU_BIT; i++)
+		if (bits[--n] == '1')
+			a->low |= (uintmax_t)1 << i;
+
+	if (!negative)
+		return +punderflow;
+
+	a->high = UINTMAX_MAX;
+
+	for (i = LIBJ2_JU_BIT; i-- && leading_ones; leading_ones--)
+		a->low |= (uintmax_t)1 << i;
+
+	return -nunderflow;
+}
+
+
+static void
+check(const char *pattern)
+{
+	struct libj2_j2i a, r, expected;
+	unsigned i;
+	int underflow;
+	intmax_t v;
+
+	set(&a, pattern, 0);
+	EXPECT(a.high == 0 || a.high == UINTMAX_MAX);
+	if (a.high)
+		v = -(intmax_t)~a.low - 1;
+	else
+		v = (intmax_t)a.low;
+
+	for (i = 0; i < LIBJ2_J2I_BIT + 8U; i++) {
+		underflow = set(&expected, pattern, i);
+
+		r = (struct libj2_j2i){111, 222};
+		libj2_ji_rsh_to_j2i(v, i, &r);
+		EXPECT(libj2_j2i_eq_j2i(&r, &expected));
+
+		r = (struct libj2_j2i){111, 222};
+		EXPECT(libj2_ji_rsh_to_j2i_underflow(v, i, &r) == underflow);
+		EXPECT(libj2_j2i_eq_j2i(&r, &expected));
+
+		EXPECT(libj2_ji_rsh_underflow_p(v, i) == underflow);
+	}
+}
+
+
+int
+main(void)
+{
+	size_t i, j;
+	char pattern[LIBJ2_J2I_BIT + 1U];
+
+	srand((unsigned)time(NULL));
+
+	for (i = 0; patterns[i]; i++)
+		check(patterns[i]);
+
+	for (i = 0; i < 64; i++) {
+		for (j = 0; j < LIBJ2_JU_BIT; j++)
+			pattern[j] = '0' + (rand() < rand());
+		pattern[LIBJ2_JU_BIT] = '\0';
+
+		pattern[0] = '0';
+		check(pattern);
+
+		pattern[0] = '1';
+		check(pattern);
+	}
+
+	return 0;
+}
 
 #endif