ref: 2bbc75bddc6f2a07056ff017108e35f14061041b
dir: /libmp/crttest.c/
#include "os.h" #include <mp.h> void testcrt(mpint **p) { CRTpre *crt; CRTres *res; mpint *m, *x, *y; fmtinstall('B', mpfmt); // get a modulus and a test number m = mpnew(1024+160); mpmul(p[0], p[1], m); x = mpnew(1024+160); mpadd(m, mpone, x); // do the precomputation for crt conversion crt = crtpre(2, p); // convert x to residues res = crtin(crt, x); // convert back y = mpnew(1024+160); crtout(crt, res, y); print("x %B\ny %B\n", x, y); mpfree(m); mpfree(x); mpfree(y); } void main(void) { int i; mpint *p[2]; long start; start = time(0); for(i = 0; i < 10; i++){ p[0] = mpnew(1024); p[1] = mpnew(1024); DSAprimes(p[0], p[1], nil); testcrt(p); mpfree(p[0]); mpfree(p[1]); } print("%ld secs with more\n", time(0)-start); exits(0); }