[u/mdw/catacomb] / tests / mp

# Test vectors for MP functions
#
# $Id: mp,v 1.15.2.2 2004/03/20 00:14:03 mdw Exp $

add {
  5 4 9; 5 -4 1; -5 4 -1; -5 -4 -9;
  0xffffffff 1 0x100000000;
}

sub {
  5 4 1; 5 -4 9; -5 4 -9; -5 -4 -1;
  4 5 -1; 4 -5 9; -4 5 -9; -4 -5 1;
}

mul {
  5 4 20; -5 4 -20; 5 -4 -20; -5 -4 20;
  0x10000 0x10000 0x100000000;
}

div {
  9 4 2 1; -9 4 -3 3; 9 -4 -3 -3; -9 -4 2 -1;
  -3 6277101735386680763835789423207666416083908700390324961279
    -1 6277101735386680763835789423207666416083908700390324961276;
  3131675836296406071791252329528905062261497366991742517193
    1110875761630725856340142297645383444629395595869672555585
    2 909924313034954359110967734238138173002706175252397406023;
  3131675836296406071791252329528905062261497366991742517193
    53
    59088223326347284373419855274130284193613157867768726739 26;
}

bin2c {
  and 5 3 1;
  or 5 3 7;
  xor 5 3 6;
  1111 0 0 -1;
  or 45 -7 -3;
  xor 0x343cd5 -0x6a49c -0x32984f;
}

lsr2c {
  -1 5 -1;
  1 5 0;
  -6 2 -2;
  5 0 5;
  -4 0 -4;
  7 2 1;
 -7 2 -2;
  -7 20 -1;
}

lsl2c {
  -1 5 -32;
  5 0 5;
  -4 0 -4;
  7 2 28;
  -7 2 -28;
  0xc0000000 1 0x180000000;
  -0xc0000000 1 -0x180000000;
  -1 32 -0x100000000;
}

setbit {
  0 40 0x10000000000;
  0x87348 40 0x10000087348;
  5 1 7;
  7 1 7;
  -3 1 -1;
}

clrbit {
  0x10000000000 40 0;
  0x87348 40 0x87348;
  5 1 5;
  7 1 5;
  -1 1 -3;
}

neg {
  0 0;
  15 -15;
  -15 15;
}

odd {
  1 0 1;
  2 1 1;
  4 2 1;
  12 2 3;
  0x10000000000000 52 1;
  0x10000000400000 22 0x40000001;
}

sqrt {
  0 0;
  1 1;
  4 2;
  9 3;
 16 4;
 99 9;
100 10;
101 10;
120 10;
121 11;

10106623487257186586 3179091613;
14565040310136678240 3816417208;
}

gcd {
  # --- Simple tests ---

  16 12 4 -11 15;
  12 16 4 -1 1;
  693 609 21 -7 8;
  4398082908043 90980984098081324 1 -32483863573352089 1570292150447;

  # --- Negative argument tests ---

  16 -12 4 -11 -15;
  -16 12 4 11 15;
  -12 -16 4 1 -1;
  -12 16 4 1 1;
  -693 609 21 7 8;
  693 -609 21 -7 -8;

  # --- Zero argument tests ---

  15 0 15 1 0;
  0 15 15 0 1;
  -5 0 5 -1 0;
  0 -5 5 0 -1;
  0 0 0 0 0;

  # --- Random number tests ---

  829561629303257626084392170900075 32498098450983560651904114638965
    5 -29340810037249902802634060204608 748967211613630574419802053172497;

  5509672937670943767152343650729669537671508
  398326674296699796695672966992514673531
  17
  -4158709420138833210339208344965073815
  57523460582278135926717203882531035926727;

  324098408098290809832490802984098208098324
  23430980840982340982098409823089098443
  1
  -4158709420138833210339208344965073815
  57523460582278135926717203882531035926727;

  # --- RSA test ---
  #
  # The first number is (p - 1)(q - 1) from `mpmont'.  The second is a
  # random number (it's actually prime, but that doesn't matter) which I
  # can use as an RSA encryption exponent.  The last is the partner
  # decryption exponent, produced using the extended GCD algorithm.

  665251164384574309450646977867043764321191240895546832784045453360
  5945908509680983480596809586040589085680968709809890671
  1
  -4601007896041464028712478963832994007038251361995647370
  514778499400157641662814932021958856708417966520837469125919104431;

  # --- Misery ---
  #
  # Some bugs discovered during RSA testing.

  100000423751500546004561515884626739136961367515520675987004088469753859696407139054406989735113827981148062449057870561788973142250811838720214530386151198455545176591384352343648452329042764530196327665219224050630680827543991306749402959935685172017409062967157813233001567797128414009962262840951763040181
  44895767034162990997987303986882660674722497505237491649296190658571471979065889234144353811843706629535512848235473808330181517421970135930320187227697512315919757806204341545022714991717913006031724818461724742069401359454784533576615919680949125073761586043027941204059690093447093117249681641020785611986
  1
  -44146175664861261172356293340716833133750232401287328189797639296698679436925232375473973898100363205157703913050824405116878299310008848005045714833814493992539429428295945643439440068026313232881493081836812480325977761600303456915493177366981470223898994906470419007730670657168179659899713837827764669213
  98330790743257232930640417364963717704786040860302439189781385170246412183980882564239377268174203679366339563908361674571088519452885615348465535190260914996055274486493192655677181637142116473172979503236297658204730543049175626205461452256333155750566288282331419748434569978343545573401114593095927172889;

  44895767034162990997987303986882660674722497505237491649296190658571471979065889234144353811843706629535512848235473808330181517421970135930320187227697512315919757806204341545022714991717913006031724818461724742069401359454784533576615919680949125073761586043027941204059690093447093117249681641020785611986
  100000423751500546004561515884626739136961367515520675987004088469753859696407139054406989735113827981148062449057870561788973142250811838720214530386151198455545176591384352343648452329042764530196327665219224050630680827543991306749402959935685172017409062967157813233001567797128414009962262840951763040181
  1
  -1669633008243313073921098519663021432175326655218236797222703299507447512426256490167612466939624301781722885149508887217884622797926223371748995195890283459489902104891159687971270691900648057023348161982926392425950284494815680543941507679352016266842774684826393484566997818784868436561148247855835867292
  749591369301729825631010646165827540972265103950163459498551361872792542140656858670379913743343424377808935184649403213303218111961287925274472393883018323380328377908395901583274923691599773150231736624912261743423597854481076661122742313967654849862591136557522196329019436278913457349967803193020942773;

  # --- Some other bugs ---

  19504439280113284806725522136967618725661733412699408177537810327183285842670
  1
  1
  0
  1;
}

jacobi {
  4 5 1;
  6 7 -1;
  15 27 0;
  2132498039840981 98729378979237498798347932749951 1;
}

modsqrt {
  1 3 1;
  4 5 3;
  9775592058107450692 13391974640168007623 3264570455655810730;
  8155671698868891620 10189552848261357803 2073812183305821596;
  3248339460720824413 8976233780911635437 1220523478429582717;
  3447751741648956439 10155704720805654949 2812971608818169892;
  1453601744816463433 3095659104519735473 1260511572497628526;
  3366261317119810224 3756232416311497601 610261287187759737;
  3869491397135339653 5762828162167967567 2974328005712882420;
  660864223630638896 1729533840094059799 1058197842375219723;
}

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