chiark - git - mdw - catacomb/blob - tests/mp

   1 # Test vectors for MP functions
   2 #
   3 # $Id: mp,v 1.13 2002/10/19 17:56:50 mdw Exp $
   4
   5 add {
   6   5 4 9; 5 -4 1; -5 4 -1; -5 -4 -9;
   7   0xffffffff 1 0x100000000;
   8 }
   9
  10 sub {
  11   5 4 1; 5 -4 9; -5 4 -9; -5 -4 -1;
  12   4 5 -1; 4 -5 9; -4 5 -9; -4 -5 1;
  13 }
  14
  15 mul {
  16   5 4 20; -5 4 -20; 5 -4 -20; -5 -4 20;
  17   0x10000 0x10000 0x100000000;
  18 }
  19
  20 div {
  21   9 4 2 1; -9 4 -3 3; 9 -4 -3 -3; -9 -4 2 -1;
  22 }
  23
  24 bin2c {
  25   and 5 3 1;
  26   or 5 3 7;
  27   xor 5 3 6;
  28   1111 0 0 -1;
  29   or 45 -7 -3;
  30   xor 0x343cd5 -0x6a49c -0x32984f;
  31 }
  32
  33 lsr2c {
  34   -1 5 -1;
  35   1 5 0;
  36   -6 2 -2;
  37 }
  38
  39 setbit {
  40   0 40 0x10000000000;
  41   0x87348 40 0x10000087348;
  42   5 1 7;
  43   7 1 7;
  44   -3 1 -1;
  45 }
  46
  47 clrbit {
  48   0x10000000000 40 0;
  49   0x87348 40 0x87348;
  50   5 1 5;
  51   7 1 5;
  52   -1 1 -3;
  53 }
  54
  55 neg {
  56   0 0;
  57   15 -15;
  58   -15 15;
  59 }
  60
  61 odd {
  62   1 0 1;
  63   2 1 1;
  64   4 2 1;
  65   12 2 3;
  66   0x10000000000000 52 1;
  67   0x10000000400000 22 0x40000001;
  68 }
  69
  70 sqrt {
  71   0 0;
  72   1 1;
  73   4 2;
  74   9 3;
  75  16 4;
  76  99 9;
  77 100 10;
  78 101 10;
  79 120 10;
  80 121 11;
  81
  82 10106623487257186586 3179091613;
  83 14565040310136678240 3816417208;
  84 }
  85
  86 gcd {
  87   # --- Simple tests ---
  88
  89   16 12 4 -11 15;
  90   12 16 4 -1 1;
  91   693 609 21 -7 8;
  92   4398082908043 90980984098081324 1 -32483863573352089 1570292150447;
  93
  94   # --- Negative argument tests ---
  95
  96   16 -12 4 -11 -15;
  97   -16 12 4 11 15;
  98   -12 -16 4 1 -1;
  99   -12 16 4 1 1;
 100   -693 609 21 7 8;
 101   693 -609 21 -7 -8;
 102
 103   # --- Zero argument tests ---
 104
 105   15 0 15 1 0;
 106   0 15 15 0 1;
 107   -5 0 5 -1 0;
 108   0 -5 5 0 -1;
 109   0 0 0 0 0;
 110
 111   # --- Random number tests ---
 112
 113   829561629303257626084392170900075 32498098450983560651904114638965
 114     5 -29340810037249902802634060204608 748967211613630574419802053172497;
 115
 116   5509672937670943767152343650729669537671508
 117   398326674296699796695672966992514673531
 118   17
 119   -4158709420138833210339208344965073815
 120   57523460582278135926717203882531035926727;
 121
 122   324098408098290809832490802984098208098324
 123   23430980840982340982098409823089098443
 124   1
 125   -4158709420138833210339208344965073815
 126   57523460582278135926717203882531035926727;
 127
 128   # --- RSA test ---
 129   #
 130   # The first number is (p - 1)(q - 1) from `mpmont'.  The second is a
 131   # random number (it's actually prime, but that doesn't matter) which I
 132   # can use as an RSA encryption exponent.  The last is the partner
 133   # decryption exponent, produced using the extended GCD algorithm.
 134
 135   665251164384574309450646977867043764321191240895546832784045453360
 136   5945908509680983480596809586040589085680968709809890671
 137   1
 138   -4601007896041464028712478963832994007038251361995647370
 139   514778499400157641662814932021958856708417966520837469125919104431;
 140
 141   # --- Misery ---
 142   #
 143   # Some bugs discovered during RSA testing.
 144
 145   100000423751500546004561515884626739136961367515520675987004088469753859696407139054406989735113827981148062449057870561788973142250811838720214530386151198455545176591384352343648452329042764530196327665219224050630680827543991306749402959935685172017409062967157813233001567797128414009962262840951763040181
 146   44895767034162990997987303986882660674722497505237491649296190658571471979065889234144353811843706629535512848235473808330181517421970135930320187227697512315919757806204341545022714991717913006031724818461724742069401359454784533576615919680949125073761586043027941204059690093447093117249681641020785611986
 147   1
 148   -44146175664861261172356293340716833133750232401287328189797639296698679436925232375473973898100363205157703913050824405116878299310008848005045714833814493992539429428295945643439440068026313232881493081836812480325977761600303456915493177366981470223898994906470419007730670657168179659899713837827764669213
 149   98330790743257232930640417364963717704786040860302439189781385170246412183980882564239377268174203679366339563908361674571088519452885615348465535190260914996055274486493192655677181637142116473172979503236297658204730543049175626205461452256333155750566288282331419748434569978343545573401114593095927172889;
 150
 151   44895767034162990997987303986882660674722497505237491649296190658571471979065889234144353811843706629535512848235473808330181517421970135930320187227697512315919757806204341545022714991717913006031724818461724742069401359454784533576615919680949125073761586043027941204059690093447093117249681641020785611986
 152   100000423751500546004561515884626739136961367515520675987004088469753859696407139054406989735113827981148062449057870561788973142250811838720214530386151198455545176591384352343648452329042764530196327665219224050630680827543991306749402959935685172017409062967157813233001567797128414009962262840951763040181
 153   1
 154   -1669633008243313073921098519663021432175326655218236797222703299507447512426256490167612466939624301781722885149508887217884622797926223371748995195890283459489902104891159687971270691900648057023348161982926392425950284494815680543941507679352016266842774684826393484566997818784868436561148247855835867292
 155   749591369301729825631010646165827540972265103950163459498551361872792542140656858670379913743343424377808935184649403213303218111961287925274472393883018323380328377908395901583274923691599773150231736624912261743423597854481076661122742313967654849862591136557522196329019436278913457349967803193020942773;
 156
 157   # --- Some other bugs ---
 158
 159   19504439280113284806725522136967618725661733412699408177537810327183285842670
 160   1
 161   1
 162   0
 163   1;
 164 }
 165
 166 jacobi {
 167   4 5 1;
 168   6 7 -1;
 169   15 27 0;
 170   2132498039840981 98729378979237498798347932749951 1;
 171 }
 172
 173 modsqrt {
 174   1 3 1;
 175   4 5 3;
 176   9775592058107450692 13391974640168007623 3264570455655810730;
 177   8155671698868891620 10189552848261357803 2073812183305821596;
 178   3248339460720824413 8976233780911635437 1220523478429582717;
 179   3447751741648956439 10155704720805654949 2812971608818169892;
 180   1453601744816463433 3095659104519735473 1260511572497628526;
 181   3366261317119810224 3756232416311497601 610261287187759737;
 182   3869491397135339653 5762828162167967567 2974328005712882420;
 183   660864223630638896 1729533840094059799 1058197842375219723;
 184 }
 185
 186 factorial {
 187   0 1;
 188   1 1;
 189   2 2;
 190   3 6;
 191   4 24;
 192   5 120;
 193   30 265252859812191058636308480000000;
 194   100 93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000;
 195   500
 196    1220136825991110068701238785423046926253574342803192842192413588385845373153881997605496447502203281863013616477148203584163378722078177200480785205159329285477907571939330603772960859086270429174547882424912726344305670173270769461062802310452644218878789465754777149863494367781037644274033827365397471386477878495438489595537537990423241061271326984327745715546309977202781014561081188373709531016356324432987029563896628911658974769572087926928871281780070265174507768410719624390394322536422605234945850129918571501248706961568141625359056693423813008856249246891564126775654481886506593847951775360894005745238940335798476363944905313062323749066445048824665075946735862074637925184200459369692981022263971952597190945217823331756934581508552332820762820023402626907898342451712006207714640979456116127629145951237229913340169552363850942885592018727433795173014586357570828355780158735432768888680120399882384702151467605445407663535984174430480128938313896881639487469658817504506926365338175055478128640000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000;
 197 }