X-Git-Url: https://git.distorted.org.uk/u/mdw/catacomb/blobdiff_plain/ba6e6b64033b1f9de49feccb5c9cd438354481f7..0f00dc4c8eb47e67bc0f148c2dd109f73a451e0a:/math/t/mpbarrett

diff --git a/math/t/mpbarrett b/math/t/mpbarrett
new file mode 100644
index 0000000..c492ece
--- /dev/null
+++ b/math/t/mpbarrett
@@ -0,0 +1,133 @@
+# Test vectors for Barrett modular reduction
+
+mpbarrett-reduce {
+  17 11 6;
+  -17 11 5;
+
+  0x8ab316d0d1a2e88535cf77c1172881ead70d592c59e9c5fbc16e4b0c4dc49481
+  0x18ca3bf7ee3c6d7bab3f144b015ccc6c25472843d346b461
+  0x02c1815029b766b96ad4507dc1af8151307961c6d161d065;
+
+  0x8117d1663ee63341eb8faeff304549f0f8b32d587acc2fd5597ea6a31625881d
+  0xdc85df77dfb61876805623bcbed325b99d00c2cd65c252c879
+  0x395da02e8a6c66476467c4e04f328d8208cc411e3d1e96e14c;
+
+  0x63791966f2ad44a6df11bcc87c6b7c2400c74e69f7e3ca02fcac12b3bf56238b
+  0xa49e473b8f7539d89cdb002d73182558773eec10db93cc6049d8c5533e
+  0x65caf6833baa118b53c7bdc44a831605ca382b5993beead59f3971d13f;
+
+  0x9ca438db3e0f79305987292e8ec6174e6c313f7904ebb35a349a700e3ae63a37
+  0xb24c93d499c7073b8f7aac718c1f12da1a8fc8bccdd47b49
+  0x46393cb15e38cbbc8a85698151a113f28081b4c8f6ed232e;
+
+  0x8214fd17858a4a913015412b5331eb9654faeb5156a674b1e5f6550a68957146
+  0xc4f0ebaad6c0ee0111c57667ea8e0a254f3068f212949e20ededa89a7da6
+  0x3fde916ba21d19414d4316041420ca59d8b01aa2acf3f3ef106245c1915c;
+
+  0x367aa8f5ba9ac4e8e2ea198b8af2c3b3081deab392ffc05715783b245a62a6fa
+  0x72e2c37447f8bca34c4a39b130ea8e5c9a7d8b54564aa88ea773
+  0x08e8c03ebf398c63d71d8fd7ca4ece12367a8dde180ca650afb6;
+
+  0xae2d84438ac6643fc601c1634351aa75b284fecbbe5faf3a132be9dd1a326e6c
+  0xc33c890f030644d88cc65f8ccf99c625c9b9fa21d4eb153e52ef89df54130855
+  0xae2d84438ac6643fc601c1634351aa75b284fecbbe5faf3a132be9dd1a326e6c;
+
+  0x65901dcdad8dd0625d4d158f99b666fee10480d1df15e3bdac640584b9b746bc
+  0xd8a1d326fee87d55f39f15b5b2cfe71f5146083928
+  0x859c41164983547c03134b99530e25a0f874315964;
+
+  -0x65901dcdad8dd0625d4d158f99b666fee10480d1df15e3bdac640584b9b746bc
+  0xd8a1d326fee87d55f39f15b5b2cfe71f5146083928
+  0x53059210b56528d9f08bca1c5fc1c17e58d1d6dfc4;
+}
+
+mpbarrett-exp {
+  4325987397987458979875737589783
+  435365332435654643667
+  8745435676786567758678547
+  2439674515119108242643169132064;
+
+  8939489893434234331 1804289383 454353454354565 6139425926295484741;
+  8939489893434234331 1804289383 8939489893434234330 1;
+
+  # --- Negative exponents ---
+
+  0xfffffffdffffffffffffffffffffffff 0xfffffffdfffffffffffffffffffffffe 0 1;
+  8939489893434234331 1804289383 -8939035539979879765 6139425926295484741;
+
+  # --- DSA public key derivation ---
+
+  0xc9c7feaeaedb16505389c5582df1858d0fdb3eecfe61c230d612661bef8c1bc5
+  0x5cd41fc97d0db5322bab7d659354db2ed9f88e39d2c6fae9f29acab5a522131e
+  0x1234
+  0x51812af9600c89ffe0f73902eb09015c03b4e0fbf6ccf073931c12f9aad1fb47;
+
+  0xdde5808744e1cd37c88667e7033694b2513a7429f035f11c0bafc4dff2b96a672bd0a3ca16aba2ea526df00c8571106ba4a1d83eb62605fc9274ab70bef0a111cd070cca2d8b10edf042d6c44f863c36fabea8bb0d7340eb8c169da27a4b0ba2713c166152a0244235093391c5f71aee8c03dcaf2335a2e4689ccb27ba365ec7
+  0x65985e4c2d6027a8afdeb9b44cc619e1c4d46bde873e0d4b45325412a2f8365e51245324f888704295fe8233a6666624d9a4701172dbfcab5c9643e1caab79eb2a0c85284d1b858688b8f16804326321f53a723502a6d6ae08dcbffccf2187a799f6281c2478ef0faed5c5c80adeabc5ee435cff8b9ae0b603e47fb08d73b014
+  0x23a252f60bae4907a8ed5b6203e2b1da32848cd9
+  0x9720498d8ec1208585635faaf952c1204c37119acccc64ed7942867be24770e33db39ffcfa1194549ead8495a7918a20e15144e68125860ef4f8c1a3d771bad690938bdb2c8817e2b89a8fc615d067084a7a2f2f9280e15fb9ccebfe713584260d5ed30545b69745d7b22977bfd44d60d7c5e657aab1c79dc5cb33ff29ee9074;
+
+  # --- Quick RSA test ---
+
+  905609324890967090294090970600361		# This is p
+  3
+  905609324890967090294090970600360		# This is (p - 1)
+  1;						# Fermat test: p is prime
+
+  734589569806680985408670989082927		# This is q
+  5
+  734589569806680985408670989082926		# And this is (q - 1)
+  1;						# Fermat again: q is prime
+
+  # --- Encrypt a message ---
+  #
+  # The public and private exponents are from the GCD test.  The message
+  # is just obvious.  The modulus is the product of the two primes above.
+
+  665251164384574309450646977867045404520085938543622535546005136647
+  123456789012345678901234567890123456789012345678901234567890
+  5945908509680983480596809586040589085680968709809890671
+  25906467774034212974484417859588980567136610347807401817990462701;
+
+  # --- And decrypt it again ---
+
+  665251164384574309450646977867045404520085938543622535546005136647
+  25906467774034212974484417859588980567136610347807401817990462701
+  514778499400157641662814932021958856708417966520837469125919104431
+  123456789012345678901234567890123456789012345678901234567890;
+
+  # --- Strange bug ---
+  #
+  # Strange memory corruption bug causing mpbarrett_reduce to hang.
+  # Seems to have vanished on a clean rebuild so maybe this was just build
+  # fluff or something.  Leave it here to pick up regression.
+
+  30916300184138066757562815127455898754391418637651479989257882011612153138316483396290143656363307389355307404813205623380378799585065876990874918745873480056891738964426752
+  2 44329
+  10975270953088118598390528358454369816072092233138947812823085741431681835772475760228604878294514156761050649439713085882651077696501584961485141097470430535290622475239424;
+}
+
+# --- Simultaneous exponentiation ---
+
+mexp-1 {
+  4325987397987458979875737589783
+  435365332435654643667
+  8745435676786567758678547
+  2439674515119108242643169132064;
+}
+
+mexp-2 {
+  0x8df2a494492276aa3d25759bb06869cbeac0d83afb8d0cf7cbb8324f0d7882e5d0762fc5b7210eafc2e9adac32ab7aac49693dfbf83724c2ec0736ee31c80291
+  0x626d027839ea0a13413163a55b4cb500299d5522956cefcb3bff10f399ce2c2e71cb9de5fa24babf58e5b79521925c9cc42e9f6f464b088cc572af53e6d78802
+  0xbf655bd046f0b35ec791b004804afcbb8ef7d69d
+  0x19131871d75b1612a819f29d78d1b0d7346f7aa77bb62a859bfd6c5675da9d212d3a36ef1672ef660b8c7c255cc0ec74858fba33f44c06699630a76b030ee333
+  0x821a926312e97adeabcc8d082b5278978a2df4b0
+  0x2fc6cb9ac3be0eac3daf02eefb96fca3846708a28dd05730165fe50942f7f07edfef8e52fcb9369e3814aa24607e80475d0e61ad461d6b16b6cec5baae58946e;
+
+  0x8df2a494492276aa3d25759bb06869cbeac0d83afb8d0cf7cbb8324f0d7882e5d0762fc5b7210eafc2e9adac32ab7aac49693dfbf83724c2ec0736ee31c80291
+  0x626d027839ea0a13413163a55b4cb500299d5522956cefcb3bff10f399ce2c2e71cb9de5fa24babf58e5b79521925c9cc42e9f6f464b088cc572af53e6d78802
+  0xbf655bd046f0b35ec791b004804afcbb8ef7d69d
+  0x19131871d75b1612a819f29d78d1b0d7346f7aa77bb62a859bfd6c5675da9d212d3a36ef1672ef660b8c7c255cc0ec74858fba33f44c06699630a76b030ee333
+  -0x8df2a494492276aa3d25759bb06869cbeac0d83afb8d0cf7cbb8324f0d7882e5d0762fc5b7210eafc2e9adabb090e849367fc31d4c6a97bac0b4be56a79a0de0
+  0x2fc6cb9ac3be0eac3daf02eefb96fca3846708a28dd05730165fe50942f7f07edfef8e52fcb9369e3814aa24607e80475d0e61ad461d6b16b6cec5baae58946e;
+}