pgen_jumpctx j;
pgen_gcdstepctx g;
unsigned nb = nbits/2;
- mp *x = MP_NEWSEC, *t = MP_NEW;
+ mp *x = MP_NEWSEC;
/* --- Generate @p@ --- */
g.r = mp_lsr(MP_NEW, bp->p, 1);
g.g = MP_NEW;
g.max = MP_ONE;
- t = mp_lsl(t, MP_ONE, nbits - 1);
- mp_div(&t, 0, t, bp->p);
- if (MP_CMP(x, <, t)) x = mp_leastcongruent(x, t, x, g.jp.m);
bp->q = pgen("q", MP_NEWSEC, x, event, ectx, n, pgen_gcdstep, &g,
rabin_iters(nb), pgen_test, &rb);
pfilt_destroy(&g.jp);
mp_drop(g.r);
mp_drop(g.g);
- mp_drop(t);
if (!bp->q) goto fail_q;
/* --- Compute @n@ --- */
bp->n = mp_mul(MP_NEW, bp->p, bp->q);
+ if (mp_bits(bp->n) != nbits) goto fail_n;
mp_drop(x);
return (PGEN_DONE);
/* --- Tidy up if things went wrong --- */
+fail_n:
+ mp_drop(bp->n);
fail_q:
mp_drop(bp->p);
fail_p: