New test cases for square root functions.
[u/mdw/catacomb] / tests / mp
1 # Test vectors for MP functions
2 #
3 # $Id: mp,v 1.5 2000/06/22 19:01:14 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 odd {
25 1 0 1;
26 2 1 1;
27 4 2 1;
28 12 2 3;
29 0x10000000000000 52 1;
30 0x10000000400000 22 0x40000001;
31 }
32
33 sqrt {
34 0 0;
35 1 1;
36 4 2;
37 9 3;
38 16 4;
39 99 9;
40 100 10;
41 101 10;
42 120 10;
43 121 11;
44
45 10106623487257186586 3179091613;
46 14565040310136678240 3816417208;
47 }
48
49 gcd {
50
51 # --- Simple tests ---
52
53 16 12 4 -2 3;
54 12 16 4 -1 1;
55 693 609 21 -181 206;
56 4398082908043 90980984098081324 1 -32483863573352089 1570292150447;
57
58 # --- Negative argument tests ---
59
60 16 -12 4 -2 -3;
61 -16 12 4 2 3;
62 -12 -16 4 1 -1;
63 -12 16 4 1 1;
64 -693 609 21 181 206;
65 693 -609 21 -181 -206;
66
67 # --- Zero argument tests ---
68
69 15 0 15 1 0;
70 0 15 15 0 1;
71 -5 0 5 -1 0;
72 0 -5 5 0 -1;
73 0 0 0 0 0;
74
75 # --- Random number tests ---
76
77 829561629303257626084392170900075 32498098450983560651904114638965
78 5 -22841190347053190672253237276815 583054885752979049202923618992482;
79
80 5509672937670943767152343650729669537671508
81 398326674296699796695672966992514673531
82 17
83 -191606556147997561067126486929677861359
84 2650310725368604614586643627755316700713319;
85
86 324098408098290809832490802984098208098324
87 23430980840982340982098409823089098443
88 1
89 -4158709420138833210339208344965073815
90 57523460582278135926717203882531035926727;
91
92 # --- RSA test ---
93 #
94 # The first number is (p - 1)(q - 1) from `mpmont'. The second is a
95 # random number (it's actually prime, but that doesn't matter) which I
96 # can use as an RSA encryption exponent. The last is the partner
97 # decryption exponent, produced using the extended GCD algorithm.
98
99 665251164384574309450646977867043764321191240895546832784045453360
100 5945908509680983480596809586040589085680968709809890671
101 1
102 -4601007896041464028712478963832994007038251361995647370
103 514778499400157641662814932021958856708417966520837469125919104431;
104
105 # --- Misery ---
106 #
107 # Some bugs discovered during RSA testing.
108
109 100000423751500546004561515884626739136961367515520675987004088469753859696407139054406989735113827981148062449057870561788973142250811838720214530386151198455545176591384352343648452329042764530196327665219224050630680827543991306749402959935685172017409062967157813233001567797128414009962262840951763040181
110 44895767034162990997987303986882660674722497505237491649296190658571471979065889234144353811843706629535512848235473808330181517421970135930320187227697512315919757806204341545022714991717913006031724818461724742069401359454784533576615919680949125073761586043027941204059690093447093117249681641020785611986
111 1
112 -44146175664861261172356293340716833133750232401287328189797639296698679436925232375473973898100363205157703913050824405116878299310008848005045714833814493992539429428295945643439440068026313232881493081836812480325977761600303456915493177366981470223898994906470419007730670657168179659899713837827764669213
113 98330790743257232930640417364963717704786040860302439189781385170246412183980882564239377268174203679366339563908361674571088519452885615348465535190260914996055274486493192655677181637142116473172979503236297658204730543049175626205461452256333155750566288282331419748434569978343545573401114593095927172889;
114
115 44895767034162990997987303986882660674722497505237491649296190658571471979065889234144353811843706629535512848235473808330181517421970135930320187227697512315919757806204341545022714991717913006031724818461724742069401359454784533576615919680949125073761586043027941204059690093447093117249681641020785611986
116 100000423751500546004561515884626739136961367515520675987004088469753859696407139054406989735113827981148062449057870561788973142250811838720214530386151198455545176591384352343648452329042764530196327665219224050630680827543991306749402959935685172017409062967157813233001567797128414009962262840951763040181
117 1
118 -1669633008243313073921098519663021432175326655218236797222703299507447512426256490167612466939624301781722885149508887217884622797926223371748995195890283459489902104891159687971270691900648057023348161982926392425950284494815680543941507679352016266842774684826393484566997818784868436561148247855835867292
119 749591369301729825631010646165827540972265103950163459498551361872792542140656858670379913743343424377808935184649403213303218111961287925274472393883018323380328377908395901583274923691599773150231736624912261743423597854481076661122742313967654849862591136557522196329019436278913457349967803193020942773;
120
121 # --- Some other bugs ---
122
123 19504439280113284806725522136967618725661733412699408177537810327183285842670
124 1
125 1
126 0
127 1;
128 }
129
130 jacobi {
131 4 5 1;
132 6 7 -1;
133 15 27 0;
134 2132498039840981 98729378979237498798347932749951 1;
135 }
136
137 modsqrt {
138 4 5 3;
139 9775592058107450692 13391974640168007623 3264570455655810730;
140 8155671698868891620 10189552848261357803 2073812183305821596;
141 3248339460720824413 8976233780911635437 1220523478429582717;
142 3447751741648956439 10155704720805654949 2812971608818169892;
143 1453601744816463433 3095659104519735473 1260511572497628526;
144 3366261317119810224 3756232416311497601 610261287187759737;
145 3869491397135339653 5762828162167967567 2974328005712882420;
146 660864223630638896 1729533840094059799 1058197842375219723;
147 }