threads::shared 1.21
[p5sagit/p5-mst-13.2.git] / ext / Math / BigInt / FastCalc / t / bigintfc.t
CommitLineData
062a4e99 1#!/usr/bin/perl -w
2
3use strict;
4use Test;
5
6BEGIN
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
15use Math::BigInt::FastCalc;
16
17my ($BASE_LEN, $AND_BITS, $XOR_BITS, $OR_BITS, $BASE_LEN_SMALL, $MAX_VAL) =
18 Math::BigInt::FastCalc->_base_len();
19
20print "# BASE_LEN = $BASE_LEN\n";
21print "# MAX_VAL = $MAX_VAL\n";
22print "# AND_BITS = $AND_BITS\n";
23print "# XOR_BITS = $XOR_BITS\n";
24print "# IOR_BITS = $OR_BITS\n";
25
26# testing of Math::BigInt::FastCalc
27
28my $C = 'Math::BigInt::FastCalc'; # pass classname to sub's
29
30# _new and _str
31my $x = $C->_new("123"); my $y = $C->_new("321");
32ok (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
37my $ff = $C->_new(0xffffffff);
38
39ok ($C->_str($ff),"4294967295"); # must not be -1
40ok (scalar @{ $ff }, 2); # must be two parts
41
42###############################################################################
062a4e99 43# _add, _sub, _mul, _div
44ok ($C->_str($C->_add($x,$y)),444);
45ok ($C->_str($C->_sub($x,$y)),123);
46ok ($C->_str($C->_mul($x,$y)),39483);
47ok ($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
52ok ($C->_str($C->_mul($x,$y)),39483);
53ok ($C->_str($x),39483); ok ($C->_str($y),321);
54
55ok ($C->_str($C->_div($x,$y)),123);
56ok ($C->_str($x),123); ok ($C->_str($y),321);
57
58$x = $C->_new("39483");
59my ($x1,$r1) = $C->_div($x,$y);
60ok ("$x1","$x");
61$C->_inc($x1);
62ok ("$x1","$x");
63ok ($C->_str($r1),'0');
64
65$x = $C->_new("39483"); # reset
66
67###############################################################################
68my $z = $C->_new("2");
69ok ($C->_str($C->_add($x,$z)),39485);
70my ($re,$rr) = $C->_div($x,$y);
71
72ok ($C->_str($re),123); ok ($C->_str($rr),2);
73
74# is_zero, _is_one, _one, _zero
75ok ($C->_is_zero($x)||0,0);
76ok ($C->_is_one($x)||0,0);
77
78ok ($C->_str($C->_zero()),"0");
79ok ($C->_str($C->_one()),"1");
80
81# _two() and _ten()
82ok ($C->_str($C->_two()),"2");
83ok ($C->_str($C->_ten()),"10");
84ok ($C->_is_ten($C->_two())||0,0);
85ok ($C->_is_two($C->_two()),1);
86ok ($C->_is_ten($C->_ten()),1);
87ok ($C->_is_two($C->_ten())||0,0);
88
89ok ($C->_is_one($C->_one()),1);
90ok ($C->_is_one($C->_two()) || 0,0);
91ok ($C->_is_one($C->_ten()) || 0,0);
92
93ok ($C->_is_one($C->_zero()) || 0,0);
94
95ok ($C->_is_zero($C->_zero()),1);
96
97ok ($C->_is_zero($C->_one()) || 0,0);
98
99# is_odd, is_even
100ok ($C->_is_odd($C->_one()),1); ok ($C->_is_odd($C->_zero())||0,0);
101ok ($C->_is_even($C->_one()) || 0,0); ok ($C->_is_even($C->_zero()),1);
102
103# _len
c039cd6f 104for 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");
146ok ($C->_digit($x,0),9);
147ok ($C->_digit($x,1),8);
148ok ($C->_digit($x,2),7);
149ok ($C->_digit($x,-1),1);
150ok ($C->_digit($x,-2),2);
151ok ($C->_digit($x,-3),3);
152
153# _copy
154foreach (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");
169ok ($C->_str($C->_lsft($x,$y,10)),10000);
170$x = $C->_new("20"); $y = $C->_new("3");
171ok ($C->_str($C->_lsft($x,$y,10)),20000);
172
173$x = $C->_new("128"); $y = $C->_new("4");
174ok ($C->_str($C->_lsft($x,$y,2)), 128 << 4);
175
176$x = $C->_new("1000"); $y = $C->_new("3");
177ok ($C->_str($C->_rsft($x,$y,10)),1);
178$x = $C->_new("20000"); $y = $C->_new("3");
179ok ($C->_str($C->_rsft($x,$y,10)),20);
180$x = $C->_new("256"); $y = $C->_new("4");
181ok ($C->_str($C->_rsft($x,$y,2)),256 >> 4);
182
183$x = $C->_new("6411906467305339182857313397200584952398");
184$y = $C->_new("45");
185ok ($C->_str($C->_rsft($x,$y,10)),0);
186
187# _acmp
188$x = $C->_new("123456789");
189$y = $C->_new("987654321");
190ok ($C->_acmp($x,$y),-1);
191ok ($C->_acmp($y,$x),1);
192ok ($C->_acmp($x,$x),0);
193ok ($C->_acmp($y,$y),0);
194$x = $C->_new("12");
195$y = $C->_new("12");
196ok ($C->_acmp($x,$y),0);
197$x = $C->_new("21");
198ok ($C->_acmp($x,$y),1);
199ok ($C->_acmp($y,$x),-1);
200$x = $C->_new("123456789");
201$y = $C->_new("1987654321");
202ok ($C->_acmp($x,$y),-1);
203ok ($C->_acmp($y,$x),+1);
204
205$x = $C->_new("1234567890123456789");
206$y = $C->_new("987654321012345678");
207ok ($C->_acmp($x,$y),1);
208ok ($C->_acmp($y,$x),-1);
209ok ($C->_acmp($x,$x),0);
210ok ($C->_acmp($y,$y),0);
211
212$x = $C->_new("1234");
213$y = $C->_new("987654321012345678");
214ok ($C->_acmp($x,$y),-1);
215ok ($C->_acmp($y,$x),1);
216ok ($C->_acmp($x,$x),0);
217ok ($C->_acmp($y,$y),0);
218
219# _modinv
220$x = $C->_new("8");
221$y = $C->_new("5033");
222my ($xmod,$sign) = $C->_modinv($x,$y);
223ok ($C->_str($xmod),'629'); # -629 % 5033 == 4404
224ok ($sign, '-');
225
226# _div
227$x = $C->_new("3333"); $y = $C->_new("1111");
228ok ($C->_str(scalar $C->_div($x,$y)),3);
229$x = $C->_new("33333"); $y = $C->_new("1111"); ($x,$y) = $C->_div($x,$y);
230ok ($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
235foreach (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
248ok ($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
250ok ($C->_str($C->_root($x,$n)),'3');
251
252# _pow (and _root)
253$x = $C->_new("0"); $n = $C->_new("3"); # 0 ** y => 0
254ok ($C->_str($C->_pow($x,$n)), 0);
255$x = $C->_new("3"); $n = $C->_new("0"); # x ** 0 => 1
256ok ($C->_str($C->_pow($x,$n)), 1);
257$x = $C->_new("1"); $n = $C->_new("3"); # 1 ** y => 1
258ok ($C->_str($C->_pow($x,$n)), 1);
259$x = $C->_new("5"); $n = $C->_new("1"); # x ** 1 => x
260ok ($C->_str($C->_pow($x,$n)), 5);
261
262$x = $C->_new("81"); $n = $C->_new("3"); # 81 ** 3 == 531441
263ok ($C->_str($C->_pow($x,$n)),81 ** 3);
264
265ok ($C->_str($C->_root($x,$n)),81);
266
267$x = $C->_new("81");
268ok ($C->_str($C->_pow($x,$n)),81 ** 3);
269ok ($C->_str($C->_pow($x,$n)),'150094635296999121'); # 531441 ** 3 ==
270
271ok ($C->_str($C->_root($x,$n)),'531441');
272ok ($C->_str($C->_root($x,$n)),'81');
273
274$x = $C->_new("81"); $n = $C->_new("14");
275ok ($C->_str($C->_pow($x,$n)),'523347633027360537213511521');
276ok ($C->_str($C->_root($x,$n)),'81');
277
278$x = $C->_new("523347633027360537213511520");
279ok ($C->_str($C->_root($x,$n)),'80');
280
281$x = $C->_new("523347633027360537213511522");
282ok ($C->_str($C->_root($x,$n)),'81');
283
284my $res = [ qw/ 9 31 99 316 999 3162 9999/ ];
285
286# 99 ** 2 = 9801, 999 ** 2 = 998001 etc
287for 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
328foreach (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 }
335foreach (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 }
342foreach (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
353my $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");
372ok ($C->_str(scalar $C->_mod($x,$y)),1);
373$x = $C->_new("1000"); $y = $C->_new("2");
374ok ($C->_str(scalar $C->_mod($x,$y)),0);
375
376# _and, _or, _xor
377$x = $C->_new("5"); $y = $C->_new("2");
378ok ($C->_str(scalar $C->_xor($x,$y)),7);
379$x = $C->_new("5"); $y = $C->_new("2");
380ok ($C->_str(scalar $C->_or($x,$y)),7);
381$x = $C->_new("5"); $y = $C->_new("3");
382ok ($C->_str(scalar $C->_and($x,$y)),1);
383
c039cd6f 384# _from_hex, _from_bin, _from_oct
062a4e99 385ok ($C->_str( $C->_from_hex("0xFf")),255);
386ok ($C->_str( $C->_from_bin("0b10101011")),160+11);
c039cd6f 387ok ($C->_str( $C->_from_oct("0100")), 8*8);
388ok ($C->_str( $C->_from_oct("01000")), 8*8*8);
389ok ($C->_str( $C->_from_oct("010001")), 8*8*8*8+1);
390ok ($C->_str( $C->_from_oct("010007")), 8*8*8*8+7);
062a4e99 391
c039cd6f 392# _as_hex, _as_bin, as_oct
062a4e99 393ok ($C->_str( $C->_from_hex( $C->_as_hex( $C->_new("128")))), 128);
394ok ($C->_str( $C->_from_bin( $C->_as_bin( $C->_new("128")))), 128);
c039cd6f 395ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("128")))), 128);
396
397ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("123456")))), 123456);
398ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("123456789")))), "123456789");
399ok ($C->_str( $C->_from_oct( $C->_as_oct( $C->_new("1234567890123")))), "1234567890123");
400
401# _1ex
402ok ($C->_str($C->_1ex(0)), "1");
403ok ($C->_str($C->_1ex(1)), "10");
404ok ($C->_str($C->_1ex(2)), "100");
405ok ($C->_str($C->_1ex(12)), "1000000000000");
406ok ($C->_str($C->_1ex(16)), "10000000000000000");
062a4e99 407
408# _check
409$x = $C->_new("123456789");
410ok ($C->_check($x),0);
411ok ($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
4461;
447