Commit | Line | Data |
062a4e99 |
1 | #!/usr/bin/perl -w |
2 | |
3 | use strict; |
4 | use Test; |
5 | |
6 | BEGIN |
7 | { |
8 | $| = 1; |
9 | chdir 't' if -d 't'; |
10 | unshift @INC, '../lib'; # for running manually |
11 | unshift @INC, '../blib/arch'; # for running manually |
4e99e077 |
12 | plan tests => 361; |
062a4e99 |
13 | } |
14 | |
15 | use Math::BigInt::FastCalc; |
16 | |
17 | my ($BASE_LEN, $AND_BITS, $XOR_BITS, $OR_BITS, $BASE_LEN_SMALL, $MAX_VAL) = |
18 | Math::BigInt::FastCalc->_base_len(); |
19 | |
20 | print "# BASE_LEN = $BASE_LEN\n"; |
21 | print "# MAX_VAL = $MAX_VAL\n"; |
22 | print "# AND_BITS = $AND_BITS\n"; |
23 | print "# XOR_BITS = $XOR_BITS\n"; |
24 | print "# IOR_BITS = $OR_BITS\n"; |
25 | |
26 | # testing of Math::BigInt::FastCalc |
27 | |
28 | my $C = 'Math::BigInt::FastCalc'; # pass classname to sub's |
29 | |
30 | # _new and _str |
31 | my $x = $C->_new("123"); my $y = $C->_new("321"); |
32 | ok (ref($x),'ARRAY'); ok ($C->_str($x),123); ok ($C->_str($y),321); |
33 | |
34 | ############################################################################### |
4e99e077 |
35 | # _new(0xffffffff) (the test is important for 32 bit Perls) |
36 | |
37 | my $ff = $C->_new(0xffffffff); |
38 | |
39 | ok ($C->_str($ff),"4294967295"); # must not be -1 |
40 | ok (scalar @{ $ff }, 2); # must be two parts |
41 | |
42 | ############################################################################### |
062a4e99 |
43 | # _add, _sub, _mul, _div |
44 | ok ($C->_str($C->_add($x,$y)),444); |
45 | ok ($C->_str($C->_sub($x,$y)),123); |
46 | ok ($C->_str($C->_mul($x,$y)),39483); |
47 | ok ($C->_str($C->_div($x,$y)),123); |
48 | |
49 | ############################################################################### |
50 | # check that mul/div doesn't change $y |
51 | # and returns the same reference, not something new |
52 | ok ($C->_str($C->_mul($x,$y)),39483); |
53 | ok ($C->_str($x),39483); ok ($C->_str($y),321); |
54 | |
55 | ok ($C->_str($C->_div($x,$y)),123); |
56 | ok ($C->_str($x),123); ok ($C->_str($y),321); |
57 | |
58 | $x = $C->_new("39483"); |
59 | my ($x1,$r1) = $C->_div($x,$y); |
60 | ok ("$x1","$x"); |
61 | $C->_inc($x1); |
62 | ok ("$x1","$x"); |
63 | ok ($C->_str($r1),'0'); |
64 | |
65 | $x = $C->_new("39483"); # reset |
66 | |
67 | ############################################################################### |
68 | my $z = $C->_new("2"); |
69 | ok ($C->_str($C->_add($x,$z)),39485); |
70 | my ($re,$rr) = $C->_div($x,$y); |
71 | |
72 | ok ($C->_str($re),123); ok ($C->_str($rr),2); |
73 | |
74 | # is_zero, _is_one, _one, _zero |
75 | ok ($C->_is_zero($x)||0,0); |
76 | ok ($C->_is_one($x)||0,0); |
77 | |
78 | ok ($C->_str($C->_zero()),"0"); |
79 | ok ($C->_str($C->_one()),"1"); |
80 | |
81 | # _two() and _ten() |
82 | ok ($C->_str($C->_two()),"2"); |
83 | ok ($C->_str($C->_ten()),"10"); |
84 | ok ($C->_is_ten($C->_two())||0,0); |
85 | ok ($C->_is_two($C->_two()),1); |
86 | ok ($C->_is_ten($C->_ten()),1); |
87 | ok ($C->_is_two($C->_ten())||0,0); |
88 | |
89 | ok ($C->_is_one($C->_one()),1); |
90 | ok ($C->_is_one($C->_two()) || 0,0); |
91 | ok ($C->_is_one($C->_ten()) || 0,0); |
92 | |
93 | ok ($C->_is_one($C->_zero()) || 0,0); |
94 | |
95 | ok ($C->_is_zero($C->_zero()),1); |
96 | |
97 | ok ($C->_is_zero($C->_one()) || 0,0); |
98 | |
99 | # is_odd, is_even |
100 | ok ($C->_is_odd($C->_one()),1); ok ($C->_is_odd($C->_zero())||0,0); |
101 | ok ($C->_is_even($C->_one()) || 0,0); ok ($C->_is_even($C->_zero()),1); |
102 | |
103 | # _len |
c039cd6f |
104 | for my $method (qw/_alen _len/) |
062a4e99 |
105 | { |
c039cd6f |
106 | $x = $C->_new("1"); ok ($C->$method($x),1); |
107 | $x = $C->_new("12"); ok ($C->$method($x),2); |
108 | $x = $C->_new("123"); ok ($C->$method($x),3); |
109 | $x = $C->_new("1234"); ok ($C->$method($x),4); |
110 | $x = $C->_new("12345"); ok ($C->$method($x),5); |
111 | $x = $C->_new("123456"); ok ($C->$method($x),6); |
112 | $x = $C->_new("1234567"); ok ($C->$method($x),7); |
113 | $x = $C->_new("12345678"); ok ($C->$method($x),8); |
114 | $x = $C->_new("123456789"); ok ($C->$method($x),9); |
115 | |
116 | $x = $C->_new("8"); ok ($C->$method($x),1); |
117 | $x = $C->_new("21"); ok ($C->$method($x),2); |
118 | $x = $C->_new("321"); ok ($C->$method($x),3); |
119 | $x = $C->_new("4321"); ok ($C->$method($x),4); |
120 | $x = $C->_new("54321"); ok ($C->$method($x),5); |
121 | $x = $C->_new("654321"); ok ($C->$method($x),6); |
122 | $x = $C->_new("7654321"); ok ($C->$method($x),7); |
123 | $x = $C->_new("87654321"); ok ($C->$method($x),8); |
124 | $x = $C->_new("987654321"); ok ($C->$method($x),9); |
125 | |
126 | $x = $C->_new("0"); ok ($C->$method($x),1); |
127 | $x = $C->_new("20"); ok ($C->$method($x),2); |
128 | $x = $C->_new("320"); ok ($C->$method($x),3); |
129 | $x = $C->_new("4320"); ok ($C->$method($x),4); |
130 | $x = $C->_new("54320"); ok ($C->$method($x),5); |
131 | $x = $C->_new("654320"); ok ($C->$method($x),6); |
132 | $x = $C->_new("7654320"); ok ($C->$method($x),7); |
133 | $x = $C->_new("87654320"); ok ($C->$method($x),8); |
134 | $x = $C->_new("987654320"); ok ($C->$method($x),9); |
135 | |
136 | for (my $i = 1; $i < 9; $i++) |
137 | { |
138 | my $a = "$i" . '0' x ($i-1); |
139 | $x = $C->_new($a); |
140 | print "# Tried len '$a'\n" unless ok ($C->_len($x),$i); |
141 | } |
062a4e99 |
142 | } |
143 | |
144 | # _digit |
145 | $x = $C->_new("123456789"); |
146 | ok ($C->_digit($x,0),9); |
147 | ok ($C->_digit($x,1),8); |
148 | ok ($C->_digit($x,2),7); |
149 | ok ($C->_digit($x,-1),1); |
150 | ok ($C->_digit($x,-2),2); |
151 | ok ($C->_digit($x,-3),3); |
152 | |
153 | # _copy |
154 | foreach (qw/ 1 12 123 1234 12345 123456 1234567 12345678 123456789/) |
155 | { |
156 | $x = $C->_new("$_"); |
157 | ok ($C->_str($C->_copy($x)),"$_"); |
158 | ok ($C->_str($x),"$_"); # did _copy destroy original x? |
159 | } |
160 | |
161 | # _zeros |
162 | $x = $C->_new("1256000000"); ok ($C->_zeros($x),6); |
163 | $x = $C->_new("152"); ok ($C->_zeros($x),0); |
164 | $x = $C->_new("123000"); ok ($C->_zeros($x),3); |
165 | $x = $C->_new("0"); ok ($C->_zeros($x),0); |
166 | |
167 | # _lsft, _rsft |
168 | $x = $C->_new("10"); $y = $C->_new("3"); |
169 | ok ($C->_str($C->_lsft($x,$y,10)),10000); |
170 | $x = $C->_new("20"); $y = $C->_new("3"); |
171 | ok ($C->_str($C->_lsft($x,$y,10)),20000); |
172 | |
173 | $x = $C->_new("128"); $y = $C->_new("4"); |
174 | ok ($C->_str($C->_lsft($x,$y,2)), 128 << 4); |
175 | |
176 | $x = $C->_new("1000"); $y = $C->_new("3"); |
177 | ok ($C->_str($C->_rsft($x,$y,10)),1); |
178 | $x = $C->_new("20000"); $y = $C->_new("3"); |
179 | ok ($C->_str($C->_rsft($x,$y,10)),20); |
180 | $x = $C->_new("256"); $y = $C->_new("4"); |
181 | ok ($C->_str($C->_rsft($x,$y,2)),256 >> 4); |
182 | |
183 | $x = $C->_new("6411906467305339182857313397200584952398"); |
184 | $y = $C->_new("45"); |
185 | ok ($C->_str($C->_rsft($x,$y,10)),0); |
186 | |
187 | # _acmp |
188 | $x = $C->_new("123456789"); |
189 | $y = $C->_new("987654321"); |
190 | ok ($C->_acmp($x,$y),-1); |
191 | ok ($C->_acmp($y,$x),1); |
192 | ok ($C->_acmp($x,$x),0); |
193 | ok ($C->_acmp($y,$y),0); |
194 | $x = $C->_new("12"); |
195 | $y = $C->_new("12"); |
196 | ok ($C->_acmp($x,$y),0); |
197 | $x = $C->_new("21"); |
198 | ok ($C->_acmp($x,$y),1); |
199 | ok ($C->_acmp($y,$x),-1); |
200 | $x = $C->_new("123456789"); |
201 | $y = $C->_new("1987654321"); |
202 | ok ($C->_acmp($x,$y),-1); |
203 | ok ($C->_acmp($y,$x),+1); |
204 | |
205 | $x = $C->_new("1234567890123456789"); |
206 | $y = $C->_new("987654321012345678"); |
207 | ok ($C->_acmp($x,$y),1); |
208 | ok ($C->_acmp($y,$x),-1); |
209 | ok ($C->_acmp($x,$x),0); |
210 | ok ($C->_acmp($y,$y),0); |
211 | |
212 | $x = $C->_new("1234"); |
213 | $y = $C->_new("987654321012345678"); |
214 | ok ($C->_acmp($x,$y),-1); |
215 | ok ($C->_acmp($y,$x),1); |
216 | ok ($C->_acmp($x,$x),0); |
217 | ok ($C->_acmp($y,$y),0); |
218 | |
219 | # _modinv |
220 | $x = $C->_new("8"); |
221 | $y = $C->_new("5033"); |
222 | my ($xmod,$sign) = $C->_modinv($x,$y); |
223 | ok ($C->_str($xmod),'629'); # -629 % 5033 == 4404 |
224 | ok ($sign, '-'); |
225 | |
226 | # _div |
227 | $x = $C->_new("3333"); $y = $C->_new("1111"); |
228 | ok ($C->_str(scalar $C->_div($x,$y)),3); |
229 | $x = $C->_new("33333"); $y = $C->_new("1111"); ($x,$y) = $C->_div($x,$y); |
230 | ok ($C->_str($x),30); ok ($C->_str($y),3); |
231 | $x = $C->_new("123"); $y = $C->_new("1111"); |
232 | ($x,$y) = $C->_div($x,$y); ok ($C->_str($x),0); ok ($C->_str($y),123); |
233 | |
234 | # _num |
235 | foreach (qw/1 12 123 1234 12345 1234567 12345678 123456789 1234567890/) |
236 | { |
237 | $x = $C->_new("$_"); |
238 | ok (ref($x)||'','ARRAY'); ok ($C->_str($x),"$_"); |
239 | $x = $C->_num($x); ok (ref($x)||'',''); ok ($x,$_); |
240 | } |
241 | |
242 | # _sqrt |
243 | $x = $C->_new("144"); ok ($C->_str($C->_sqrt($x)),'12'); |
244 | $x = $C->_new("144000000000000"); ok ($C->_str($C->_sqrt($x)),'12000000'); |
245 | |
246 | # _root |
247 | $x = $C->_new("81"); my $n = $C->_new("3"); # 4*4*4 = 64, 5*5*5 = 125 |
248 | ok ($C->_str($C->_root($x,$n)),'4'); # 4.xx => 4.0 |
249 | $x = $C->_new("81"); $n = $C->_new("4"); # 3*3*3*3 == 81 |
250 | ok ($C->_str($C->_root($x,$n)),'3'); |
251 | |
252 | # _pow (and _root) |
253 | $x = $C->_new("0"); $n = $C->_new("3"); # 0 ** y => 0 |
254 | ok ($C->_str($C->_pow($x,$n)), 0); |
255 | $x = $C->_new("3"); $n = $C->_new("0"); # x ** 0 => 1 |
256 | ok ($C->_str($C->_pow($x,$n)), 1); |
257 | $x = $C->_new("1"); $n = $C->_new("3"); # 1 ** y => 1 |
258 | ok ($C->_str($C->_pow($x,$n)), 1); |
259 | $x = $C->_new("5"); $n = $C->_new("1"); # x ** 1 => x |
260 | ok ($C->_str($C->_pow($x,$n)), 5); |
261 | |
262 | $x = $C->_new("81"); $n = $C->_new("3"); # 81 ** 3 == 531441 |
263 | ok ($C->_str($C->_pow($x,$n)),81 ** 3); |
264 | |
265 | ok ($C->_str($C->_root($x,$n)),81); |
266 | |
267 | $x = $C->_new("81"); |
268 | ok ($C->_str($C->_pow($x,$n)),81 ** 3); |
269 | ok ($C->_str($C->_pow($x,$n)),'150094635296999121'); # 531441 ** 3 == |
270 | |
271 | ok ($C->_str($C->_root($x,$n)),'531441'); |
272 | ok ($C->_str($C->_root($x,$n)),'81'); |
273 | |
274 | $x = $C->_new("81"); $n = $C->_new("14"); |
275 | ok ($C->_str($C->_pow($x,$n)),'523347633027360537213511521'); |
276 | ok ($C->_str($C->_root($x,$n)),'81'); |
277 | |
278 | $x = $C->_new("523347633027360537213511520"); |
279 | ok ($C->_str($C->_root($x,$n)),'80'); |
280 | |
281 | $x = $C->_new("523347633027360537213511522"); |
282 | ok ($C->_str($C->_root($x,$n)),'81'); |
283 | |
284 | my $res = [ qw/ 9 31 99 316 999 3162 9999/ ]; |
285 | |
286 | # 99 ** 2 = 9801, 999 ** 2 = 998001 etc |
287 | for my $i (2 .. 9) |
288 | { |
289 | $x = '9' x $i; $x = $C->_new($x); |
290 | $n = $C->_new("2"); |
291 | my $rc = '9' x ($i-1). '8' . '0' x ($i-1) . '1'; |
292 | print "# _pow( ", '9' x $i, ", 2) \n" unless |
293 | ok ($C->_str($C->_pow($x,$n)),$rc); |
294 | |
295 | if ($i <= 7) |
296 | { |
297 | $x = '9' x $i; $x = $C->_new($x); |
298 | $n = '9' x $i; $n = $C->_new($n); |
299 | print "# _root( ", '9' x $i, ", ", 9 x $i, ") \n" unless |
300 | ok ($C->_str($C->_root($x,$n)),'1'); |
301 | |
302 | $x = '9' x $i; $x = $C->_new($x); |
303 | $n = $C->_new("2"); |
304 | print "# _root( ", '9' x $i, ", ", 9 x $i, ") \n" unless |
305 | ok ($C->_str($C->_root($x,$n)), $res->[$i-2]); |
306 | } |
307 | } |
308 | |
309 | ############################################################################## |
310 | # _fac |
311 | $x = $C->_new("0"); ok ($C->_str($C->_fac($x)),'1'); |
312 | $x = $C->_new("1"); ok ($C->_str($C->_fac($x)),'1'); |
313 | $x = $C->_new("2"); ok ($C->_str($C->_fac($x)),'2'); |
314 | $x = $C->_new("3"); ok ($C->_str($C->_fac($x)),'6'); |
315 | $x = $C->_new("4"); ok ($C->_str($C->_fac($x)),'24'); |
316 | $x = $C->_new("5"); ok ($C->_str($C->_fac($x)),'120'); |
317 | $x = $C->_new("10"); ok ($C->_str($C->_fac($x)),'3628800'); |
318 | $x = $C->_new("11"); ok ($C->_str($C->_fac($x)),'39916800'); |
319 | $x = $C->_new("12"); ok ($C->_str($C->_fac($x)),'479001600'); |
320 | $x = $C->_new("13"); ok ($C->_str($C->_fac($x)),'6227020800'); |
321 | |
322 | # test that _fac modifes $x in place for small arguments |
323 | $x = $C->_new("3"); $C->_fac($x); ok ($C->_str($x),'6'); |
324 | $x = $C->_new("13"); $C->_fac($x); ok ($C->_str($x),'6227020800'); |
325 | |
326 | ############################################################################## |
327 | # _inc and _dec |
328 | foreach (qw/1 11 121 1231 12341 1234561 12345671 123456781 1234567891/) |
329 | { |
330 | $x = $C->_new("$_"); $C->_inc($x); |
331 | print "# \$x = ",$C->_str($x),"\n" |
332 | unless ok ($C->_str($x),substr($_,0,length($_)-1) . '2'); |
333 | $C->_dec($x); ok ($C->_str($x),$_); |
334 | } |
335 | foreach (qw/19 119 1219 12319 1234519 12345619 123456719 1234567819/) |
336 | { |
337 | $x = $C->_new("$_"); $C->_inc($x); |
338 | print "# \$x = ",$C->_str($x),"\n" |
339 | unless ok ($C->_str($x),substr($_,0,length($_)-2) . '20'); |
340 | $C->_dec($x); ok ($C->_str($x),$_); |
341 | } |
342 | foreach (qw/999 9999 99999 9999999 99999999 999999999 9999999999 99999999999/) |
343 | { |
344 | $x = $C->_new("$_"); $C->_inc($x); |
345 | print "# \$x = ",$C->_str($x),"\n" |
346 | unless ok ($C->_str($x), '1' . '0' x (length($_))); |
347 | $C->_dec($x); ok ($C->_str($x),$_); |
348 | } |
349 | |
350 | $x = $C->_new("1000"); $C->_inc($x); ok ($C->_str($x),'1001'); |
351 | $C->_dec($x); ok ($C->_str($x),'1000'); |
352 | |
353 | my $BL; |
354 | { |
355 | no strict 'refs'; |
356 | $BL = &{"$C"."::_base_len"}(); |
357 | } |
358 | |
359 | $x = '1' . '0' x $BL; |
360 | $z = '1' . '0' x ($BL-1); $z .= '1'; |
361 | $x = $C->_new($x); $C->_inc($x); ok ($C->_str($x),$z); |
362 | |
363 | $x = '1' . '0' x $BL; $z = '9' x $BL; |
364 | $x = $C->_new($x); $C->_dec($x); ok ($C->_str($x),$z); |
365 | |
366 | # should not happen: |
367 | # $x = $C->_new("-2"); $y = $C->_new("4"); ok ($C->_acmp($x,$y),-1); |
368 | |
369 | ############################################################################### |
370 | # _mod |
371 | $x = $C->_new("1000"); $y = $C->_new("3"); |
372 | ok ($C->_str(scalar $C->_mod($x,$y)),1); |
373 | $x = $C->_new("1000"); $y = $C->_new("2"); |
374 | ok ($C->_str(scalar $C->_mod($x,$y)),0); |
375 | |
376 | # _and, _or, _xor |
377 | $x = $C->_new("5"); $y = $C->_new("2"); |
378 | ok ($C->_str(scalar $C->_xor($x,$y)),7); |
379 | $x = $C->_new("5"); $y = $C->_new("2"); |
380 | ok ($C->_str(scalar $C->_or($x,$y)),7); |
381 | $x = $C->_new("5"); $y = $C->_new("3"); |
382 | ok ($C->_str(scalar $C->_and($x,$y)),1); |
383 | |
c039cd6f |
384 | # _from_hex, _from_bin, _from_oct |
062a4e99 |
385 | ok ($C->_str( $C->_from_hex("0xFf")),255); |
386 | ok ($C->_str( $C->_from_bin("0b10101011")),160+11); |
c039cd6f |
387 | ok ($C->_str( $C->_from_oct("0100")), 8*8); |
388 | ok ($C->_str( $C->_from_oct("01000")), 8*8*8); |
389 | ok ($C->_str( $C->_from_oct("010001")), 8*8*8*8+1); |
390 | ok ($C->_str( $C->_from_oct("010007")), 8*8*8*8+7); |
062a4e99 |
391 | |
c039cd6f |
392 | # _as_hex, _as_bin, as_oct |
062a4e99 |
393 | ok ($C->_str( $C->_from_hex( $C->_as_hex( $C->_new("128")))), 128); |
394 | ok ($C->_str( $C->_from_bin( $C->_as_bin( $C->_new("128")))), 128); |
c039cd6f |
395 | ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("128")))), 128); |
396 | |
397 | ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("123456")))), 123456); |
398 | ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("123456789")))), "123456789"); |
399 | ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("1234567890123")))), "1234567890123"); |
400 | |
401 | # _1ex |
402 | ok ($C->_str($C->_1ex(0)), "1"); |
403 | ok ($C->_str($C->_1ex(1)), "10"); |
404 | ok ($C->_str($C->_1ex(2)), "100"); |
405 | ok ($C->_str($C->_1ex(12)), "1000000000000"); |
406 | ok ($C->_str($C->_1ex(16)), "10000000000000000"); |
062a4e99 |
407 | |
408 | # _check |
409 | $x = $C->_new("123456789"); |
410 | ok ($C->_check($x),0); |
411 | ok ($C->_check(123),'123 is not a reference'); |
412 | |
413 | ############################################################################### |
414 | # __strip_zeros |
415 | |
416 | { |
417 | no strict 'refs'; |
418 | # correct empty arrays |
419 | $x = &{$C."::__strip_zeros"}([]); ok (@$x,1); ok ($x->[0],0); |
420 | # don't strip single elements |
421 | $x = &{$C."::__strip_zeros"}([0]); ok (@$x,1); ok ($x->[0],0); |
422 | $x = &{$C."::__strip_zeros"}([1]); ok (@$x,1); ok ($x->[0],1); |
423 | # don't strip non-zero elements |
424 | $x = &{$C."::__strip_zeros"}([0,1]); |
425 | ok (@$x,2); ok ($x->[0],0); ok ($x->[1],1); |
426 | $x = &{$C."::__strip_zeros"}([0,1,2]); |
427 | ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); |
428 | |
429 | # but strip leading zeros |
430 | $x = &{$C."::__strip_zeros"}([0,1,2,0]); |
431 | ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); |
432 | |
433 | $x = &{$C."::__strip_zeros"}([0,1,2,0,0]); |
434 | ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); |
435 | |
436 | $x = &{$C."::__strip_zeros"}([0,1,2,0,0,0]); |
437 | ok (@$x,3); ok ($x->[0],0); ok ($x->[1],1); ok ($x->[2],2); |
438 | |
439 | # collapse multiple zeros |
440 | $x = &{$C."::__strip_zeros"}([0,0,0,0]); |
441 | ok (@$x,1); ok ($x->[0],0); |
442 | } |
443 | |
444 | # done |
445 | |
446 | 1; |
447 | |