Add unit tests

- Add LUTs for a sample of irreducible polynomials from Gadiel Seroussi's table
- Add test_minimum_terms_from_database

tests/polys/luts/irreducible_polys_database.py

@@ -0,0 +1,15 @@
+"""
+A module containing LUTs for irreducible polynomials with min terms from irreducible_polys.db
+"""
+
+# LUT items are poly nonzero degrees and coefficients in degree-descending order
+
+# Gadiel Seroussi's table (1998)
+# LUT items obtained by randomly picking degrees and checking the PDF
+# sorted(numpy.random.default_rng(1337).integers(size=5, low=500, high=10_000, endpoint=True))
+
+IRREDUCIBLE_POLY_MIN_TERMS_2_2262 = [[2262, 57, 0], [1, 1, 1]]
+IRREDUCIBLE_POLY_MIN_TERMS_2_5632 = [[5632, 17, 15, 5, 0], [1, 1, 1, 1, 1]]
+IRREDUCIBLE_POLY_MIN_TERMS_2_5690 = [[5690, 1623, 0], [1, 1, 1]]
+IRREDUCIBLE_POLY_MIN_TERMS_2_7407 = [[7407, 27, 21, 17, 0], [1, 1, 1, 1, 1]]
+IRREDUCIBLE_POLY_MIN_TERMS_2_8842 = [[8842, 4143, 0], [1, 1, 1]]

tests/polys/test_irreducible_polys.py

@@ -1,6 +1,7 @@
 """
 A pytest module to test generating irreducible polynomials over finite fields.
 """
+import time
 import numpy as np
 import pytest
 
@@ -42,6 +43,14 @@
 )
 from .luts.irreducible_polys_25 import IRREDUCIBLE_POLYS_25_1, IRREDUCIBLE_POLYS_25_2
 
+from .luts.irreducible_polys_database import (
+ IRREDUCIBLE_POLY_MIN_TERMS_2_2262,
+ IRREDUCIBLE_POLY_MIN_TERMS_2_5632,
+ IRREDUCIBLE_POLY_MIN_TERMS_2_5690,
+ IRREDUCIBLE_POLY_MIN_TERMS_2_7407,
+ IRREDUCIBLE_POLY_MIN_TERMS_2_8842,
+)
+
 PARAMS = [
  (2, 1, IRREDUCIBLE_POLYS_2_1),
  (2, 2, IRREDUCIBLE_POLYS_2_2),
@@ -71,6 +80,14 @@
  (5**2, 2, IRREDUCIBLE_POLYS_25_2),
 ]
 
+PARAMS_DB = [
+ (2, 2262, IRREDUCIBLE_POLY_MIN_TERMS_2_2262),
+ (2, 5632, IRREDUCIBLE_POLY_MIN_TERMS_2_5632),
+ (2, 5690, IRREDUCIBLE_POLY_MIN_TERMS_2_5690),
+ (2, 7407, IRREDUCIBLE_POLY_MIN_TERMS_2_7407),
+ (2, 8842, IRREDUCIBLE_POLY_MIN_TERMS_2_8842),
+]
+
 
 def test_irreducible_poly_exceptions():
  with pytest.raises(TypeError):
@@ -160,6 +177,18 @@ def test_minimum_terms(order, degree, polys):
  assert f.coeffs.tolist() in min_term_polys
 
 
+@pytest.mark.parametrize("order,degree,polys", PARAMS_DB)
+def test_minimum_terms_from_database(order, degree, polys):
+ tick = time.time()
+ p = galois.irreducible_poly(order, degree, terms="min")
+ tock = time.time()
+ assert tock - tick < 1.0
+ db_degrees = p.nonzero_degrees.tolist()
+ db_coeffs = p.nonzero_coeffs.tolist()
+ exp_degrees, exp_coeffs = polys
+ assert db_degrees == exp_degrees and db_coeffs == exp_coeffs
+
+
 def test_large_degree():
  """
  See https://github.com/mhostetter/galois/issues/360.