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