180ee5cbfSDavid du Colombier #include "os.h"
280ee5cbfSDavid du Colombier #include <mp.h>
380ee5cbfSDavid du Colombier #include <libsec.h>
480ee5cbfSDavid du Colombier
580ee5cbfSDavid du Colombier mpint*
egencrypt(EGpub * pub,mpint * in,mpint * out)680ee5cbfSDavid du Colombier egencrypt(EGpub *pub, mpint *in, mpint *out)
780ee5cbfSDavid du Colombier {
8*9a747e4fSDavid du Colombier mpint *m, *k, *gamma, *delta, *pm1;
980ee5cbfSDavid du Colombier mpint *p = pub->p, *alpha = pub->alpha;
10*9a747e4fSDavid du Colombier int plen = mpsignif(p);
11*9a747e4fSDavid du Colombier int shift = ((plen+Dbits)/Dbits)*Dbits;
1280ee5cbfSDavid du Colombier // in libcrypt version, (int)(LENGTH(pub->p)*sizeof(NumType)*CHARBITS);
1380ee5cbfSDavid du Colombier
1480ee5cbfSDavid du Colombier if(out == nil)
1580ee5cbfSDavid du Colombier out = mpnew(0);
16*9a747e4fSDavid du Colombier pm1 = mpnew(0);
1780ee5cbfSDavid du Colombier m = mpnew(0);
1880ee5cbfSDavid du Colombier gamma = mpnew(0);
1980ee5cbfSDavid du Colombier delta = mpnew(0);
2080ee5cbfSDavid du Colombier mpmod(in, p, m);
21*9a747e4fSDavid du Colombier while(1){
22*9a747e4fSDavid du Colombier k = mprand(plen, genrandom, nil);
23*9a747e4fSDavid du Colombier if((mpcmp(mpone, k) <= 0) && (mpcmp(k, pm1) < 0))
24*9a747e4fSDavid du Colombier break;
25*9a747e4fSDavid du Colombier }
2680ee5cbfSDavid du Colombier mpexp(alpha, k, p, gamma);
2780ee5cbfSDavid du Colombier mpexp(pub->key, k, p, delta);
2880ee5cbfSDavid du Colombier mpmul(m, delta, delta);
2980ee5cbfSDavid du Colombier mpmod(delta, p, delta);
3080ee5cbfSDavid du Colombier mpleft(gamma, shift, out);
3180ee5cbfSDavid du Colombier mpadd(delta, out, out);
32*9a747e4fSDavid du Colombier mpfree(pm1);
3380ee5cbfSDavid du Colombier mpfree(m);
3480ee5cbfSDavid du Colombier mpfree(k);
3580ee5cbfSDavid du Colombier mpfree(gamma);
3680ee5cbfSDavid du Colombier mpfree(delta);
3780ee5cbfSDavid du Colombier return out;
3880ee5cbfSDavid du Colombier }
39