1 b3f61791 2004-03-21 devnull #include "os.h"
 2 b3f61791 2004-03-21 devnull #include <mp.h>
 3 b3f61791 2004-03-21 devnull 
 4 b3f61791 2004-03-21 devnull void
 5 b3f61791 2004-03-21 devnull testcrt(mpint **p)
 6 b3f61791 2004-03-21 devnull {
 7 b3f61791 2004-03-21 devnull 	CRTpre *crt;
 8 b3f61791 2004-03-21 devnull 	CRTres *res;
 9 b3f61791 2004-03-21 devnull 	mpint *m, *x, *y;
10 b3f61791 2004-03-21 devnull 	int i;
11 b3f61791 2004-03-21 devnull 
12 b3f61791 2004-03-21 devnull 	fmtinstall('B', mpconv);
13 b3f61791 2004-03-21 devnull 
14 cbeb0b26 2006-04-01 devnull 	/* get a modulus and a test number */
15 b3f61791 2004-03-21 devnull 	m = mpnew(1024+160);
16 b3f61791 2004-03-21 devnull 	mpmul(p[0], p[1], m);
17 b3f61791 2004-03-21 devnull 	x = mpnew(1024+160);
18 b3f61791 2004-03-21 devnull 	mpadd(m, mpone, x);
19 b3f61791 2004-03-21 devnull 
20 cbeb0b26 2006-04-01 devnull 	/* do the precomputation for crt conversion */
21 b3f61791 2004-03-21 devnull 	crt = crtpre(2, p);
22 b3f61791 2004-03-21 devnull 
23 cbeb0b26 2006-04-01 devnull 	/* convert x to residues */
24 b3f61791 2004-03-21 devnull 	res = crtin(crt, x);
25 b3f61791 2004-03-21 devnull 
26 cbeb0b26 2006-04-01 devnull 	/* convert back */
27 b3f61791 2004-03-21 devnull 	y = mpnew(1024+160);
28 b3f61791 2004-03-21 devnull 	crtout(crt, res, y);
29 b3f61791 2004-03-21 devnull 	print("x %B\ny %B\n", x, y);
30 b3f61791 2004-03-21 devnull 	mpfree(m);
31 b3f61791 2004-03-21 devnull 	mpfree(x);
32 b3f61791 2004-03-21 devnull 	mpfree(y);
33 b3f61791 2004-03-21 devnull }
34 b3f61791 2004-03-21 devnull 
35 b3f61791 2004-03-21 devnull void
36 b3f61791 2004-03-21 devnull main(void)
37 b3f61791 2004-03-21 devnull {
38 b3f61791 2004-03-21 devnull 	int i;
39 b3f61791 2004-03-21 devnull 	mpint *p[2];
40 b3f61791 2004-03-21 devnull 	long start;
41 b3f61791 2004-03-21 devnull 
42 b3f61791 2004-03-21 devnull 	start = time(0);
43 b3f61791 2004-03-21 devnull 	for(i = 0; i < 10; i++){
44 b3f61791 2004-03-21 devnull 		p[0] = mpnew(1024);
45 b3f61791 2004-03-21 devnull 		p[1] = mpnew(1024);
46 b3f61791 2004-03-21 devnull 		DSAprimes(p[0], p[1], nil);
47 b3f61791 2004-03-21 devnull 		testcrt(p);
48 b3f61791 2004-03-21 devnull 		mpfree(p[0]);
49 b3f61791 2004-03-21 devnull 		mpfree(p[1]);
50 b3f61791 2004-03-21 devnull 	}
51 b3f61791 2004-03-21 devnull 	print("%d secs with more\n", time(0)-start);
52 b3f61791 2004-03-21 devnull 	exits(0);
53 b3f61791 2004-03-21 devnull }