1 b3f61791 2004-03-21 devnull #include "os.h"
2 b3f61791 2004-03-21 devnull #include <mp.h>
5 b3f61791 2004-03-21 devnull testcrt(mpint **p)
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;
12 b3f61791 2004-03-21 devnull fmtinstall('B', mpconv);
14 b3f61791 2004-03-21 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);
20 b3f61791 2004-03-21 devnull // do the precomputation for crt conversion
21 b3f61791 2004-03-21 devnull crt = crtpre(2, p);
23 b3f61791 2004-03-21 devnull // convert x to residues
24 b3f61791 2004-03-21 devnull res = crtin(crt, x);
26 b3f61791 2004-03-21 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);
36 b3f61791 2004-03-21 devnull main(void)
39 b3f61791 2004-03-21 devnull mpint *p[2];
40 b3f61791 2004-03-21 devnull long start;
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]);
51 b3f61791 2004-03-21 devnull print("%d secs with more\n", time(0)-start);
52 b3f61791 2004-03-21 devnull exits(0);