From 7e8a688588e153ed7fd01307c2d9c56566f38c6c Mon Sep 17 00:00:00 2001
From: mhostetter <matthostetter@gmail.com>
Date: Mon, 7 Nov 2022 14:41:51 -0500
Subject: [PATCH] Make `lagrange_poly()` unit tests more robust

---
 scripts/generate_field_test_vectors.py        |  19 +++++++
 tests/polys/conftest.py                       |  10 ++++
 .../polys/data/GF(109987^4)/lagrange_poly.pkl | Bin 0 -> 1329 bytes
 tests/polys/data/GF(2)/lagrange_poly.pkl      | Bin 0 -> 439 bytes
 .../data/GF(2147483647)/lagrange_poly.pkl     | Bin 0 -> 967 bytes
 tests/polys/data/GF(2^100)/lagrange_poly.pkl  | Bin 0 -> 1693 bytes
 tests/polys/data/GF(2^2)/lagrange_poly.pkl    | Bin 0 -> 546 bytes
 tests/polys/data/GF(2^3)/lagrange_poly.pkl    | Bin 0 -> 769 bytes
 tests/polys/data/GF(2^32)/lagrange_poly.pkl   | Bin 0 -> 989 bytes
 tests/polys/data/GF(2^8)/lagrange_poly.pkl    | Bin 0 -> 877 bytes
 .../data/GF(2^8, 283, 19)/lagrange_poly.pkl   | Bin 0 -> 877 bytes
 tests/polys/data/GF(31)/lagrange_poly.pkl     | Bin 0 -> 877 bytes
 tests/polys/data/GF(3191)/lagrange_poly.pkl   | Bin 0 -> 906 bytes
 .../lagrange_poly.pkl                         | Bin 0 -> 1306 bytes
 tests/polys/data/GF(5)/lagrange_poly.pkl      | Bin 0 -> 607 bytes
 tests/polys/data/GF(7)/lagrange_poly.pkl      | Bin 0 -> 715 bytes
 tests/polys/data/GF(7^3)/lagrange_poly.pkl    | Bin 0 -> 887 bytes
 .../data/GF(7^3, 643, 244)/lagrange_poly.pkl  | Bin 0 -> 887 bytes
 tests/polys/test_lagrange_poly.py             |  49 ++++--------------
 19 files changed, 38 insertions(+), 40 deletions(-)
 create mode 100644 tests/polys/data/GF(109987^4)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2147483647)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2^100)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2^2)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2^3)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2^32)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2^8)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(2^8, 283, 19)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(31)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(3191)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(36893488147419103183)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(5)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(7)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(7^3)/lagrange_poly.pkl
 create mode 100644 tests/polys/data/GF(7^3, 643, 244)/lagrange_poly.pkl

diff --git a/scripts/generate_field_test_vectors.py b/scripts/generate_field_test_vectors.py
index 9b3305a50..020dd40fa 100644
--- a/scripts/generate_field_test_vectors.py
+++ b/scripts/generate_field_test_vectors.py
@@ -964,6 +964,25 @@ def make_luts(field, sub_folder, seed, sparse=False):
     d = {"X": X, "Y": Y, "Z": Z}
     save_pickle(d, folder, "crt.pkl")
 
+    set_seed(seed + 406)
+    X = [0,]*3  # The x values
+    Y = [0,]*3  # The y values
+    Z = [0,]*3  # The lagrange polynomial
+    for i in range(3):
+        N = min(10, order)
+        if dtype == object:
+            x = randint_matrix(0, order, (N,))
+        else:
+            x = np.array(random.sample(range(0, order), N), dtype=dtype)
+        y = randint_matrix(0, order, (N,))
+        points = [(F(x[i]), F(y[i])) for i in range(N)]
+        z = RING.lagrange_polynomial(points)
+        X[i] = x
+        Y[i] = y
+        Z[i] = poly_to_list(z)
+    d = {"X": X, "Y": Y, "Z": Z}
+    save_pickle(d, folder, "lagrange_poly.pkl")
+
     ###############################################################################
     # Special polynomials
     ###############################################################################
diff --git a/tests/polys/conftest.py b/tests/polys/conftest.py
index b402fe58d..62ecec357 100644
--- a/tests/polys/conftest.py
+++ b/tests/polys/conftest.py
@@ -208,6 +208,16 @@ def poly_modular_power(field_folder):
     return d
 
 
+@pytest.fixture(scope="session")
+def poly_lagrange_poly(field_folder):
+    GF, d = read_pickle(field_folder, "lagrange_poly.pkl")
+    d["GF"] = GF
+    d["X"] = [GF(X) for X in d["X"]]
+    d["Y"] = [GF(Y) for Y in d["Y"]]
+    d["Z"] = [galois.Poly(p, field=GF) for p in d["Z"]]
+    return d
+
+
 ###############################################################################
 # Fixtures for special polynomials
 ###############################################################################
diff --git a/tests/polys/data/GF(109987^4)/lagrange_poly.pkl b/tests/polys/data/GF(109987^4)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..466dc1d64eb9bd16fc2d99ac0e45d24cfb920648
GIT binary patch
literal 1329
zcmah|X;4#F6wV5H1kq9`%YcCzRw>|AOGszPu$3j1APh@@h>$!Di!6~;p;8fu$Rh}t
zM@OIqDN<x{kPJ#g5gmqtRa!)Xlwwf^7?9N^I!tBh3G_#QmHXqqGv|Hh`@Vb5y*ZR}
zTPk5|W#jg8Vk9nrk8z6hWLl=PFjay&pOhwNB*aN1ahbSWKURVYQ&Uc5NTk9HT!G7}
zTR>c{n<CnVC~z)8hI<p^@d7*qx91YMdNN$BBd3ZoGSg6P35mzT3omZ}CCiCt@t~lf
zz)e_N8<&%T#}|SJ46&q|rP<k;&k`}jw4zbm-o=}wVu*!|#9|*<FQ#C~2ku2|&46kR
zLbP5T@k8{;j1GpRJ0#&DZY6XQhBVAy3gc#)KZS5Q+S{=raXb&OKmTO>iT{F2DP(Q^
zTXEX^Lt-ils$!~`rZq3t`>%DwZ8N>uqC(SM=zB7TJnz0?Ib1l<1XN7a=}c9Y<stz?
zBABA_-IIOEx)`$SZF<|`J>SGYcuer2VDW36Y6?aRe5}krG^1o-NYC%-R~{kygK!$}
zcapGM0qhz$W0&l<XSt3?En!D5QwvB^H@PJiit^fgSL#QN%;7BYtfxkO=Y*cX8skqt
zC>VMY4q?mk*zEPce=`8HQJl_a?FnH$kZmu!prD6uG(o#Azp=pF_R%i@;C}wy@ro;Z
zXMlI&kvYG1DC!MaYX+4Q`PysozYV%N)&CG(op1qm0pIwgL!^NT^vZ8~SUUf%LjtX!
zcXTSQ#f0Smd;W$@?^B&Gj)3dqL4I5p4_z0KS39{8Y)20Zf#Jp8a<e?1DI2o(x{q%a
zb*J-zY)&-QSlTK7e&=b0pek{oGOrY<_T2ns&Y@}PM{r<sro!N}T&DJt)VWvYGD%Vv
z0EW5iYhU#@tbm!-MD?BXl%-&>s93gsw_CH>4%pMZ&n7wd<iP+?M)6E4ts0u4ETN5L
zs&ZH55_cLGtu&HnzzW8EF@zP6_ToY0bw40wh#XJpayPY&nfqfbl)Eep1=T;RY^*<@
z9dLncLAq$l!!<|^s&yvSfxcZ#C3M5F_;tAc$EGX*9C&?@#1Vy&$QZ4~NqKZ&_;Sp@
z!-~=lYwHR5FcRZ(Z;pmTWmb0R4aTv1ROsVrquKS9<5CUSDC2EpIp-eb0ADnwC~~f_
z9fa&PQu6ur2MQFldt8G9ZJHU{yQaQBQj_x}+zvv&vgE~*+l|#g);G@)c;^XL3Fww$
zEni^omd{cHNaI`is+$}=4Ugl9&}(GBAt@KQBkovRRm2NA0kBa;kK`|neFQQLKh^VI
zHl-K{jhfuf^<Uz~0ANlM|8{@$2~Tim_9eS<!_lx7LYe7=mvV`gLfTejithg9Q>)_u
zGl{C%Ocv_Z=>j2dTrgO4*1!w!sjd5(n-{J=Bw{pGhsNTiZ-6_7xW=x0x71bG2US)l
r@OpjC+7MJV4>h!w__!BB_$v8ZV}4RWG!ST1vtt}e^RHwSl{)_eWAYO|

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2)/lagrange_poly.pkl b/tests/polys/data/GF(2)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..48bdeee12bb358b94f46e6ae44b2c81de4d0604c
GIT binary patch
literal 439
zcmZo*nYxCN0Ss!VX!I~fOo;_jqIspc1(kZq`9-OExurQJnTbV3iIr1&c;bsvlk@Y6
zONvU9OQuYo(!&Z?IHiX@F9o7x@)U1|)+x@6NmJUV1WnQKX7pxiosz-U!<tf3S&$0S
z!jx$-rLzMf%iP0gGsVx(&+9)BfC+Dgk|{}@&H@l8F)%>g#xMmaoWYvG0Wv*<3)yTT
ss;ExwVT=TMov<A+ml5Y8(o`e+7U-fVkc(o0A({v!K;g-dnp&y{0Haxb761SM

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2147483647)/lagrange_poly.pkl b/tests/polys/data/GF(2147483647)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..e9409a886ecb91d17f97a7652aac3ed26890eae2
GIT binary patch
literal 967
zcmaLUe@IhN6bJClZ5q=gO%n+fo!X)tv6M_CwHfGnIZa74ojSJ{v*tFRo4L$1bfox*
za5p9r1chj>l0r#K@dqf26e7?{bOa)pM$|tM1=-HEKl&#hkA1k?x%YhUS*!TYxJW>s
zW+LZ>B}C&R++?e@*cCc6$0}-VMr)0h<Ft0d2UT*c&TOi;ayFfnI0-MIg@g|<>32$;
zgb_H1N?1*9l0qV9gpAZd4F0@CZ?#)k+~QZ0N;*4seiiY;bfQ+PVcYl64@O`m)t#z>
zo+Y^oQ#nNLAhf(c>MA<FU2zGGkB-uy5Bi^;MF)q952Gm=kcG+)sv1z^kJEHN@f~ZB
z81`T~gY@5Qp}uZylijpxHAhT`aaoy1Pw7zl#o#IZIRSC#a$u_t9S$iAMnBa(rYUD#
zMRr{C6EEcG?%hTQ?iFE`lW(LMU%zCEas8W&9YMQhJ`2(LL*_KJci54DrYv<{MDug<
z%h2V0$(bm2V|igI{&oMq<L&qEK&*~fJwWf!=hf1QhPv8mcAa#B-sSxD8w#uo+nNkO
z?`JDoP`__#E4r05%Avd8(fjjZC|1p`WsA{*Q;Yji=Uame{jerWM?LOG5_DXaWJGgk
z&s3osbI+5|GFcc+vwJsv_?rx@Ua)j)(XN8MG-K(?t1Gx(HZX{wuxf!$nrSSg-6ZBn
z5q_V!(r^OJ8i}Q8|4#G*PSp3p#xJOcb>s08C>q*|hq|#g4NNbtOTpVTl?u%x8Vg)|
zmItui5^06G5OXw(yRwz82=KT!6QI%KX@JD;rjrmTuoXl5F>f#!zRbkJ-CoWF=e%DN
zp`x!+KzH?g)WpGb@^~_2Nhiv|Y!~Ok*l3C#Ry;@jVKZ%tgTBFA8SrF483W(nC8}B0
GruYL1K_@)`

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2^100)/lagrange_poly.pkl b/tests/polys/data/GF(2^100)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..fe51f9965a47efeee640d280d91cfa4e6a0520f0
GIT binary patch
literal 1693
zcmai!dpuNW7{_Or8N)j6%A!RiTTv<HQn@X1kKC=$m_{<wxQtuICX5`D!l4^->4QRP
zm6C~e3MFNgOBAU@(nYS-(CmZ#WB=ODALo4D&+|Ud`~5x7@4SZzIaws!;={&txOgwj
z4?c<v7M&TZM>!Zy)uXdQBlib}hX=-DTuJ|MD&-&}A~KvsiNvxoE^%=q#uZ}(Ep5rh
z>~SfWC4MiKin(Gs_IP`V6fBs;B?d*tGN}tGBF=_b+R>$`C>I}(IXOAm3E^3M*_%dU
zd((gcMY0{jtUN+UA@$m44@#lPBPZnx=c_#)C<)4^#Zcr}XQC~0f=I(OHLppc$hc~T
z{P!*D8560gEkB`%;doR25YD`Mw)8ny8b#LlTSXmcXeE@r2zdgf`bXb6$cF8zn5!U$
zL)53_K=7^EFMBCJI}So+bdY2I=j_<FEy+d(P>OoGie8hqtS~SBafcWc=z@vC<OPwz
z691Jbw*(zTB{X&*QS6Rd`de`nVSG|eL{szGcpFVIGK$!JGfDdUnY5Ngi+d7q*z6ux
zNV$Hdz1fvFC4nNpTh!Rya|<3gv~z2&2#ScR=A_~^#@@@6Dn<}d#HC>3#dxC09&UhP
z&jR$zTO!JLjYdEF`E7ux(RD^kk#2N2kG#49N^z1}&3nI&9UGOYUnBBEXkmko_9}{l
zlM0R!0|%Q<B)9%@47)(#mAfM-l2=vy`eI!{ZJ}ID4+%wV_n;%ng0`+aUlJ7*Ub@q=
ztCo5j2_Y&y6h)DGafNktm)@+4G+y40M-eTx+pR6p=F9J>v*V$)b7isCW%5qnx%#0S
zVQ59gyy8~+SF7`dL$e?$^$+$K^F1+z#`bC=EkIMbc<%*2{%^;NiiVs-Csro72<#Mq
z78^b3p_VqBQMu`~0`NSRQ!^oZS6aS!8<Yo2<=oPwz}RcpmIO~B0Y#3d9#T<vJQ=0P
z>DdPJu)t}T_};W4<d5VCfOcv}5ngGSTbcKYMFYGU{-~#$RfHX3WBfGC@@ONKt|gP^
z{OGT|C$O1qDE-MmYS6o+ivR7$HhFhtQ)@%=yE$BdQZ$UquJ*;Q6lIN)8|<qppn;RX
zsgb|ZXPM6DjKMsTGYM;S<8nl5G^=mJNzqT{;#E&(M)|qi7eMrL@vw(zoOtnM;kq|~
zHJx!)Z!TTQwqjHNIcRZG#18##|5Gm$w%4n`Qt`5zZWhVA<>IZTqy_zPxgS4BX7V36
z8n)wqXcJ;8)zh3KzS{PbY%?U&K9i32WLisEaso`jf%^>)Me$_))<2Ux5e6D)E}AT4
zQm4Euq-Ov&r)xh=Y;zDbE!f=O4bZ4Tl|zP`_Oexuv2~DJ?$yp@+CL`{=RMy6Wz@{}
z60AHR*5dA))6g9QXA|GPf(l-CH_I7vSsq%wHyF$!sZi^7<X@ZkEZAi6)`%cKZ|`d`
zG1l1JG!N|%v%K$z5eD+rW7=r2Y|B`SPbRunNM$Y41bZGoLYg_!r~l>I(Bu~Q)i;~n
zl)2u3ogejd7Xt|&J$gTAq@Z{u;YMe`=umU7vB|XTV+gMbKvCp$soB}g(;;g+#qi(`
zKK&JkU>2aYp5;k~+_Qu>3iCnG#eTkVln8Z^tn+!V%}Gt)M)@PJbs*t8Z8TY$cX+Je
zF8x<X=w@PZLd8kuWSRCpAPd~NqWx3CKwEuiMKg3{Y-;{}ibiLUS8N*&TF|b-IN<0}
zaYd9P_64FV9eu*PEa+Q9?IgOu<PX^#qaAWL-$u*}8BmI6QY*v<8&K_KpALb@4NX%T
z=aA~=+}<!6X^E`0hFAAnrEbn-Z+{s4;lmUbzmaB5cO?qXfKbx|y#piN4*J70t4nZj
z4;jo0_MhtNuYZ}921{Rom9mZ*b4L8VyeW`bp3oOKju2xG`%yAru)DppjqfS7d?j16
zydlaYi;m@8-lm$FOw<R>kM?m^@Y7o>N600iaFzJoDRmt5&o^+rfJRGEsVu#J0Krt%
AfdBvi

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2^2)/lagrange_poly.pkl b/tests/polys/data/GF(2^2)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..f5f153a9ff477cd457858a5819fa0ef40a50e866
GIT binary patch
literal 546
zcma)(u?oU45QdYqRS_paTpgV{JGi+S=-{M+NLvFHw3V0+9R!!6LN?Fm)ue4e3Kf%Z
zcaZz;{!6a*Ye-PBE-;pO2^-}Gae6#uZkTZ99@B%^`<(k3$ljK-Fp1BCr=fr+NT)~y
z*^YN*$`dHcVTJ>^!2-q<Q#*&KC7qqfPD~?YKY_=+ytO2rg6Dbjw_1gvNI>wIjVev0
zGG%cUS(OKs7CBKzcZ<5;z4~8toigQq)24@^>smC;Bwp#KYSvJvO%+r9KMi&2zX=T0
b<f{qRdIEJQQ5#Xh${4F(Ln@}kXQ}%E7J`1l

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2^3)/lagrange_poly.pkl b/tests/polys/data/GF(2^3)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..5e8c9558ebc6f6e7f5e161a10e9e759ef444704d
GIT binary patch
literal 769
zcmaJ;J!``-5S3)h4y99}{~?p-4jls7BJj|uEtHZt0tXr=xc1<o&?U41ZMpu{uI<Uz
zgHWg?+t2s)?w#}=e|;;_uQ!@Hc~4Ia*HydScUiHkZMJP+n$5DVmmRsw=h_y#>b0rc
zqM;Kx-6xVutM!m_A|nnomn-_9JDM74CI>1Lr`Julx4uMeZs_<qL?b6}DbMr853l}X
zL_@3Ne3lNHL=sR33|fT@b{sGY2q2~@bvf~B%8PR?SImwqRDlKT=!Py(0afGd$O0>>
z!3K78`>V#YTPGj<Rs25#*ujQbV52M4APyX#0oBlmD%h9>h#62f&W>!bVistHjX1g@
p9!}`dpAbJLG1>@|@|OHY*83v%#cXT-Qzq$ss`q+Z@>bTi*)OR4h0_24

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2^32)/lagrange_poly.pkl b/tests/polys/data/GF(2^32)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..2a1491ba57b5ceae8c20079284dc5f8977944692
GIT binary patch
literal 989
zcmaLVdq`7J7y$5lwOb)(q*)<qVFXKQDPbmQ%tRMwx_m5rkdNsuPE)tmrp|oj1NAbE
z+?j?BH8U$kQ%KVk=d?dk15pcKDX1m%M<w|H(H{G*{n0;pxZEG#z27<Ccka1)`*aiL
z1Us@ZM?;3=C=<4+sua1J8!uDP+!R%^GEuBhh}BqQ8Ka=%WvO?R3RS!kYq5r66S3w%
zs$@4siv>gx_9f$RF%H2T0V%L3!qWX3N}^QDX_#V`cnO!5?7lVEkX~3Q6r$gz!QKRf
z633PL291C09)hl(&(-vUGd(Xwpybq6Ie4TsiGA+T9IyiM)Z(fT^gmvf1(vHjC&3nL
z+b$5BmX%V{gV3GS;@{q^P1X2;{#KBsi_X)7=XtYLU~oaM5iI;^mk;t++A_fOUq`jz
zh})q7u%SI87@T)siUjwrcD2Bn?bS23_BLeOt1Q{oGUd0DAdcJSg@Wxnj?-ZMf+PS;
zIL+Gxx7R)og8ElUJ|I-1HDm-F-T(I(v$%N<vik<iBfv>J!&~r1XC)WZKU$9h=`Qgc
zIM`Sh3-bBnY<DQ=t_Q?)ZNUb(_a6N|g(pH*KWUr?mtQ0}gN0-D2t3pKjROwFh?-cF
zhbKU-Ze#(JWU7CFbvdEU;NI1AUf04P%b*{10l#j_w?JcrV>jp?CScDmd+yD5h$|}8
zJ>boI_XEMvw16~l=jjLbulxH&BjJmhE{5W*4jTv*ZM^8jP)@}TI|SO!uW~^mPRVuj
z&So?nrCYZ7G1QkV?@9vMJuW=SP=-ePS%PL~i)3UO!BEGfm)A_V;u8%7vJqa9p@;rc
zT!vbDa$%W3W)qVf)ZslQVyM;`{LGYi*s=OOy8^25>7!7ot%E&siz>-rsH)ti2LxK^
h%_LDsM~)axZ-mRx>lU>YYCP*{jX3Mot~9OU{sB{WLSz5{

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2^8)/lagrange_poly.pkl b/tests/polys/data/GF(2^8)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..3e56b877e08549c7d0f13cb111cfc3b80db5411b
GIT binary patch
literal 877
zcmaLVJ4?e*6bJCcwkov<Dg_7KJJjl?h+8)gl%h*3i1ZOip)X=W?NAh?Y9&kQDSiMQ
z1W~%V2re!nD2ND(gARfYf-Z#$K9ckwaIqoemveLOJ@=6{XFcu^tgFO6U5xOefj6sV
z^2I<Rr$~W}mR3_SMTr$z_be$=B9~oO6)mB1jdhom$a+sU*-oi331tq7ao*%P_K`^L
zGRxh%E2$RqlG&vzHN};!_E(24&TuFcn*A{@>p_Cb@k%gUZ_i11gbd%H6RyJsJb{z&
z0J1Tx?2>!U*~#v|a~m|OO(W|X9>Es;hW2P9h*w|<p2N=7(#VR!4QN;0MXbRHe1_vN
z2hD1dE=J9x`+to%W#I<B9^ebSgJF0NU!gsVRm7dw_=GG8enNXqUc`QQ0SmAX+B@92
z+A*?|@DPr|CftYi-}8m|7JfkN=z_VU#!WIcC`477Fy5*(Pup}UneA89PfdEK0*z6I
bJVyJV)jmL{bVe;2B0ufW9_49RlC;1tF5tG7

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(2^8, 283, 19)/lagrange_poly.pkl b/tests/polys/data/GF(2^8, 283, 19)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..74a74ca37a703d0d23471aefc1bebdc88b8f6f21
GIT binary patch
literal 877
zcmaLVJxIe)5C`zYS}jEo3a$<gQYmzDa1aGM)q^6UU=^wT2&UREVlq@|sY0oRP{E#?
zZcgIhC<qGTtY8NzLfsq$LD9iM&?Nl_T<jy{mwUP7?!BbceGqd9)+n%77Z-ThU?8KV
zb9sL}t4RK|mQs^ZMTzEF?~5o>Je%226)moEopq;`$a-%k(N3u|31tq5F|P0|dr2gB
zndKhcnNagN$$X?MImXr9_FacAj&mp!n*K3{bs|CKST!&+)b2?*hZ=l=kMI&+z<XGS
z3CPB4WvA>iyOZ62d-E98mXI|C&)^ZXt8F5lfJOKQRoJ;2l7t2>TJQndy?-E{hqv$v
zUPH5*q>Ewm)BV539b}KdW4H|ium<no4XnaE?7YS~WNpEI7=%su1yA57oP<}fb2T5b
z_F)*>YmOowh0pL1dZ7JxVMiCu9W@T}l143BAqNdojY_mi4auBOQ;J@xX~e_yP7&Ip
YqBTEA-{hh-3eq~AQiTfSlO)am3j(&bQvd(}

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(31)/lagrange_poly.pkl b/tests/polys/data/GF(31)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..3de3960304f9996eba66de722dee6fdf362d4fa2
GIT binary patch
literal 877
zcma))%TK~U5XQSedC5b-qJntl#My&4ZcajsCp9sKN0XAsBegwnAn`(!WN(arukM!l
z5-t*(roYVYbiVm!_apN?N;}E3=B`$c{9@qE<Jof68!tk?H;dn*=`ak3E7qmg&>t`6
z?@<_!Bi^xglSJ0J`DCB6V{$g!S0n!92kw%h;)a8)b|=wl>DymY(;MD??RV2!-SS{C
zxcf0r@*yYU(YAkI-1p>M2pdA_)j5~=L|75lg)L!0$fj)I26@|EP(19dZq!PW)fA$7
zRpOfPQiyJ{5<|5u@&D93$wsfJ29Kw*uY`Hwnea+z)qJg<?9u)0u_M`-C<yP^M^#ty
zp`Mc%6F=?|p3p0LM>TkbD)KQ|oEvyt4Yvi8gBm6XPi64{HB`sB9(pqqeYO)dFQq9%
u6>=#>9ebn6qpEM$8&sw&mCU+F1@o<vO4hH<dWTxnq_$z*{8RJ&xc3V!N{h(=

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(3191)/lagrange_poly.pkl b/tests/polys/data/GF(3191)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..72ab4815ae6ef375f5d4592c420996c071326136
GIT binary patch
literal 906
zcma)3&npCB7~XGw?5^3pP;yc@$ze~DgXG6WrgxKzA0;t6wPY7#rtCo|vTBm}oXF1u
z;y{!OImk_1lu&ZPAK>7~%zlq@A=5NZ-}^q_^S;kB6Ux%K9LPOmyxe8&yvu{aw38aj
z#Ee19jHR3byRXNxdNSONc3Y+~m>#h$$FTV#chww`yP<SKP!_o!nC6YFm(TE4Uam7;
zo94-&t0wHskm+BM`l|VCR&-^T)$n*c-u&r3xsM*OdGBmvi;*Xm+6@U@(FE3F*(dZ(
zMbtl;z_T!Txc7`R_hgf5GVE`c)PC+(;7yx^`Xty_P(77|+6Jc1W9o@02W~LMiHO2x
zP~hLw`k`6~m3pY`V~aOcH14rY8Tywv4@X763tRfrOqX@|zwURB>zG!A>F2op0QEi8
zkKv&b{5sTEkZ|^29=G7-JIOgDw+0Vy=(j=r00}ODw}Zq>;q~uCAEEXRPj|p6Bv}F#
z0ye<IEj-=<SATf(67BR8^<FwuXjGz$FgYPQ3(~ei*|7gLWO|L5KBGdSVvQDpWNUOU
qQGSGq@+iURMWrJ($JmkRTGN}TN~4rawTupA`jBZ^Ar>`FC-w#Rle)qH

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(36893488147419103183)/lagrange_poly.pkl b/tests/polys/data/GF(36893488147419103183)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..6ced2b1877b43f01ebfbe2e91ac22014953c11c7
GIT binary patch
literal 1306
zcmaizdoWyg6vy{sxmRO728maKC`Lu8%1|?vh9t(Ot3nY}SeC8qZjkj@R=v6*zigha
zNR73qR6SM`G?~a#tHNZ1N>OwqWvHO_+^+HLv8FTqtGR#N`<Zjj?|i@K_xs%(jbdFY
zMfGHn_F^oQ{0u%uycBMtz$xZTGRujZ!r^}&otzvkAjP{*C9`7A@Xqp+Q)2jJ5h<ps
z3Q4gRFLtM-h@?|Qq&LPS3&<0sJsqQKi%7Pnm>SC$B(k=bsQJ5-=L>hT>S9btG8hcM
zEqGKfx)+~ho(B(z=m=D_XytEJVu;`w=BMnDTnhsXtTWBexgE9)=+|1X5P4b@2{@#Y
zkznhQbpg=gdWevgbB>Q8yqOc6xGy1158;hP<y#+(mhVNxMGMaB4&|}C6hx?7aa~i(
zzF)>zpvk7P@onpD?f)7#*s;>pSvY=9#@2)O6INF{%w6YO0QZJihs56scLdDp_cs1u
zveg?<%kD5eGw+8`&~Lqbq{Q=jl?x&yb<N#vI`az1bE>m1r;o2c1Oq=M4O_=2^xHzz
zMl>S(R<rEXo4vGcC?+Dtrv8RS^G5w;!G6P!03QbFX1$Z@7y@X&X`El{d##;<@H|&X
zEuM_F8xie&rTPt}{o5P#wSk&fRORv%<jbdSImdL22@r8tSuyiqGnfJS!%fQxnmU8g
zDl058c5?el{O|U>N|U}U=o6R0$hB;D`AF-HcM)zgv*`9SjY~tgLyLw_$)0D25dM1Z
zVy+aEf(YH%ce3j3{Y!w;a-kyij@wbdfzjF6)y8TO?2XB#YTHn5=L+<69iG!*O75Gy
z@tE2LM!LND=^*ITgB%j23&XvDYaPWj?OoNgFz;?&wwz;l-wP2`A??7cikd=D`tuV#
zKbp&f08qN2se1{d6{I4fVk<+&)o<U<2Q13Eq>r6P0PB~3yU@MTRY`f16?{ItC&aow
z(-+}s5{IaXqt87NA?OzTQqRoGgGNd-yLyp_{CDuAZrSf{y}?sMc$B%4nC?Cpgow^W
zxv}P*7x$o2K9sSPZ=72MXq7f#Ru{hC6K2xy{i~@jjcfAn8HH`1k?OB`S*I6K^60}s
zn8rEIXrL+kM=gYBM)0d^a{Tup!dtYd7?|+532FPZfUyStug4G`d8J4p^0y2`xM|YI
z4B>%vGelq#v#WzIlycC;WQ*~$dheG5vTt*XkGLtGP*^Ils6ldWmK`?)nHAh8>h$xf
zZ@}Emv&~j<&QB&GD)ez(dD=6BAZlLv$9uJFwND{h<Fd}I*WvWS^ggBH#~ZGZu?Ux|
zNfX9%4NMU}=so7xwzg)12nL#VU+ylk#8h&)ow6m!rgssp^U_y7q490ZGH5oL^wyzI
znZs0=Xf=HF;PHFTuuywu;mZa5yI9bTl&7~zCH^M?Pps2S4gK7-p%f*{7Ow;~>LY@0
Zwi?mRaVdd((-$iZ7SBRN6c#JR=^rKiBa{FD

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(5)/lagrange_poly.pkl b/tests/polys/data/GF(5)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..8884420dbd89ab4793f80fb17df29783d8d1de52
GIT binary patch
literal 607
zcmb7>KMTSz5XF=JQ4uG>@6f5UgPWUy4&tDKh}J+0S|z4K2f;;D$maQZHEDYw0YTG}
z*Q57de%GG8=fol9Izgc56c$=~S-v|IVZ0YS+~r%Du7p?>pgK#z<9&9NB9A3JK>1}O
zsAjgQS{}g21%`ADR~SRUDDy9nG?c%Q#eo}(llI|$tKPbb4j_u6(Mzv#V?@IGJ`CK7
ziPUVwt7{ZmKY(se{G?@0P1>Jlc6EyFMi-`I+HWbg+h#H;T+x~Ngx{I$6lDC53CYog
micF}#O@?k{v9+FTZYI;W{;Mu?%+9egwO3Xit1gH0JbVMlY=Kq)

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(7)/lagrange_poly.pkl b/tests/polys/data/GF(7)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..b9657c8eec1057f513a0f7c415905f5b028a43ca
GIT binary patch
literal 715
zcma)2v1-FG5S8r6DWy}P|B%UZhHlv^@X)C(l#)0C2NEZ!^5CJ+CA0x;x%^N5td%T1
z2!@0Rb5A<od-v|X2Zv+Ay7NUN%U|f-!9(4wHf>(6waQn`(k!N0Pg}BgA6k{``qSv9
zG_)rxx<azUdLBylB-xfGd`3U?Oe4wVa7$HY#oV+T<tfQxOuOIkH?{ncilTV>bEgZE
z8Jg`T<LHSEf;W~$<Pk+gKpHRP=#{8$yo>7g)Wf|~1q*6@Hq`swFRQ{MJe=#&SNE#Q
z@>jp|^Pb@q2+u&|f$)BvXRv_%888cGysU}|!GcGy(FJeR#ZA2VO}I}cQpilElDU|v
UkV#zDy?pF4X}G%@pHWTz589N3X#fBK

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(7^3)/lagrange_poly.pkl b/tests/polys/data/GF(7^3)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..11d458cc536231c6a3f29bdc932f12e9a2eb8945
GIT binary patch
literal 887
zcmaKqyGsK>5XNs4qai4W5X7j3Xb~TrAc~3#f}rBCNl+0Lf_ZolAurAg7*VkB3D-GR
z!Nyk9&cfbO{0}T_M6j?F8@<bY2R3Fo?zi90%+BoW9@JfAy-IcLbJ!HC9Jdh6mvY6j
znl9*?nk!|Et)#9e%WO7o=vumv-!}A8+Te3GeN{&`8}gfu%Q=&B#4{qrC%njE5)?S%
zOugybG|ENI4)JV_@!6?k^_pUmBaz6$k2R`4QVdR=&BWqPP09n<4t|0O@CKX~&Ljr?
z1ulVXZJY6B{C0Jjz~9<htZX~*n!_$0=wq-OjDjoRAy@&Q!P?o}n#SQhg{&d)5$m(i
z-_XO*LFjAH&Ze1S&Az(-&)ACS0YvqHZihReZ^1pV1AGS)VC^#!vZ?TPvo}Fc!rlu$
zVf_-sn^bn8YiD!sYXZBSAgT`xVBPJs`-<jZcRTCe06MyEchvgmMUD!J%1$98g6^rA
v{L(LENNBe4LoN$?r7oJMH@cIGpo;W3rZyQ9G)PyZ(2!IGjSy46rj^uRM-rz6

literal 0
HcmV?d00001

diff --git a/tests/polys/data/GF(7^3, 643, 244)/lagrange_poly.pkl b/tests/polys/data/GF(7^3, 643, 244)/lagrange_poly.pkl
new file mode 100644
index 0000000000000000000000000000000000000000..b5a2c84d75552959551aae6d6f977e6ee5266316
GIT binary patch
literal 887
zcmaKq%SyvQ6oxasDW#$aDyX0$NEceQAP5!|ih?~%5ryJ~BE1A!Y10yuMHhk;i>+kW
z>E4wK@fCalbt8g;JD2te1Q(jiA8=!bAz%JEC+9yilXG9E<dy7kCc*~W7ikkCie9eP
zBDsnxN6LChD`r$RQxit$Se0`X<y2GkoF>}B@Y{*P=u`4eN?VZB5OJ0jmtt3hiBYg2
z3IW5P*J@SS`ovRQ6|JTd^%^WDlF4M^#~k*9Bu!*n@tqmwG<D@bdIKZiBH~Zbd!PnB
zf>}_QWefg7ztvqK__w!V)9o2q2~_3aEY|Db7x)3jz$SPC_IA5QmJcp~w27B4;sIzK
z`T}(Kl!CV0WP=@8SNH!Oqo_6su7TgEl7h}b2f<}90D6$udygB)cHLxXw~G+;JeUT@
zP<;W~{SSJ(b&&1u$#s3fx_fVT&|7es0w=(C#BrjB)`^;vZ!>d}LihAZgFMdo0OP}q
tQnDGOFm)+Pb-LlxjMn)aqX#;&Bd<Kf_z<HCt?(sAqZFWVx|L-;@(bm*r#%1w

literal 0
HcmV?d00001

diff --git a/tests/polys/test_lagrange_poly.py b/tests/polys/test_lagrange_poly.py
index 60d658c47..2bca924fa 100644
--- a/tests/polys/test_lagrange_poly.py
+++ b/tests/polys/test_lagrange_poly.py
@@ -36,43 +36,12 @@ def test_exceptions():
         galois.lagrange_poly(GF([0, 1, 2, 0]), y)
 
 
-def test_binary_field():
-    GF = galois.GF(2)
-    x = GF([0, 1])
-    y = GF([0, 0])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("0", field=GF)
-
-    GF = galois.GF(2)
-    x = GF([0, 1])
-    y = GF([0, 1])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("x", field=GF)
-
-    GF = galois.GF(2)
-    x = GF([0, 1])
-    y = GF([1, 0])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("x + 1", field=GF)
-
-    x = GF([0, 1])
-    y = GF([1, 1])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("1", field=GF)
-
-
-def test_prime_field():
-    GF = galois.GF(251)
-    x = GF([10, 117, 142, 120, 163, 13, 37, 67, 135, 55])
-    y = GF([197, 191, 31, 24, 65, 222, 224, 58, 47, 178])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("207*x^9 + 222*x^8 + 227*x^7 + 88*x^6 + 6*x^5 + 52*x^4 + 82*x^3 + 198*x^2 + 27*x + 78", field=GF)
-
-
-def test_binary_extension_field():
-    GF = galois.GF(2**8)
-    x = GF([141, 144, 248, 59, 208, 32, 254, 46, 230, 15])
-    y = GF([250, 27, 81, 177, 62, 208, 221, 76, 182, 6])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("104*x^9 + 130*x^8 + 92*x^7 + 199*x^6 + 64*x^5 + 211*x^4 + 130*x^3 + 38*x^2 + 140*x + 114", field=GF)
-
-
-def test_prime_extension_field():
-    GF = galois.GF(3**5)
-    x = GF([114, 151, 235, 198, 129, 192, 73, 184, 186, 78])
-    y = GF([152, 50, 232, 129, 212, 226, 152, 26, 148, 239])
-    assert galois.lagrange_poly(x, y) == galois.Poly.Str("42*x^9 + 200*x^8 + 109*x^7 + 82*x^6 + x^5 + 241*x^4 + 65*x^3 + 153*x^2 + 143*x + 88", field=GF)
+def test_lagrange_poly(poly_lagrange_poly):
+    GF, X, Y, Z = poly_lagrange_poly["GF"], poly_lagrange_poly["X"], poly_lagrange_poly["Y"], poly_lagrange_poly["Z"]
+    for i in range(len(X)):
+        x = X[i]
+        y = Y[i]
+        poly = galois.lagrange_poly(x, y)
+        assert type(poly) == galois.Poly
+        assert poly.field == GF
+        assert poly == Z[i]