+ * Returns: Minimum recommended amount of data to test.
+ *
+ * Use: Initializes a Maurer testing context. Note that @l@ values
+ * larger than 16 are not supported, and values less than 6 can
+ * give bizarre results.
+ */
+
+unsigned long maurer_init(maurer_ctx *m, unsigned l)
+{
+ size_t i;
+
+ assert(((void)"(maurer_init): chunks larger than 16 bits not supported",
+ 0 < l && l <= 16));
+
+ m->l = l;
+ m->a = 0;
+ m->b = 0;
+ m->n = 0;
+ m->x = 0;
+
+ m->t = xmalloc(sizeof(unsigned long) << l);
+ for (i = 0; i < (1 << l); i++)
+ m->t[i] = 0;
+ return (((1000ul << l) * l + 7)/8);
+}
+
+/* --- @bits@ --- *
+ *
+ * Arguments: @maurer_ctx *m@ = pointer to testing context
+ * @const octet **p@ = address of a buffer pointer
+ * @const octet *q@ = limit of the buffer pointer
+ * @unsigned *x@ = address to store next chunk
+ *
+ * Returns: Nonzero if some more bits arrived.
+ *
+ * Use: Reads some bits from a buffer.
+ */
+
+static int bits(maurer_ctx *m, const octet **p, const octet *q, unsigned *x)
+{
+ while (m->b < m->l) {
+ if (*p >= q)
+ return (0);
+ m->a |= U8(*(*p)++) << m->b;
+ m->b += 8;
+ }
+ *x = m->a & ((1 << m->l) - 1);
+ m->a >>= m->l;
+ m->b -= m->l;
+ m->n++;
+ return (1);
+}
+
+/* --- @maurer_test@ --- *
+ *
+ * Arguments: @maurer_ctx *m@ = pointer to a Maurer testing context
+ * @const void *buf@ = pointer to a buffer of data
+ * @size_t sz@ = size of the buffer
+ *
+ * Returns: ---
+ *
+ * Use: Scans a buffer of data and updates the testing context.
+ */
+
+void maurer_test(maurer_ctx *m, const void *buf, size_t sz)
+{
+ const octet *p = buf, *l = p + sz;
+ unsigned long q = 10 << m->l;
+ unsigned x;
+
+ /* --- Initialize the table --- */
+
+ while (m->n < q) {
+ if (!bits(m, &p, l, &x))
+ return;
+ m->t[x] = m->n;
+ }
+
+ /* --- Update the statistic --- */
+
+ while (bits(m, &p, l, &x)) {
+ m->x += log(m->n - m->t[x]);
+ m->t[x] = m->n;
+ }
+}
+
+/* --- @maurer_done@ --- *
+ *
+ * Arguments: @maurer_ctx *m@ = pointer to a Maurer testing context
+ *