X-Git-Url: https://git.distorted.org.uk/u/mdw/catacomb/blobdiff_plain/043613346cc5a0a99ad54781212bf7fd2c62c415..d92b3cee9432cfb30b85acd42ed4b4f8dc48b4d7:/limlee.c diff --git a/limlee.c b/limlee.c index a25c034..dd03866 100644 --- a/limlee.c +++ b/limlee.c @@ -1,6 +1,6 @@ /* -*-c-*- * - * $Id: limlee.c,v 1.1 2000/07/09 21:30:58 mdw Exp $ + * $Id: limlee.c,v 1.3 2000/07/29 09:58:32 mdw Exp $ * * Generate Lim-Lee primes * @@ -30,6 +30,13 @@ /*----- Revision history --------------------------------------------------* * * $Log: limlee.c,v $ + * Revision 1.3 2000/07/29 09:58:32 mdw + * (limlee): Bug fix. Old versions didn't set the filter step if @ql@ was + * an exact divisor of @pl@. + * + * Revision 1.2 2000/07/26 18:00:00 mdw + * No footer line! + * * Revision 1.1 2000/07/09 21:30:58 mdw * Lim-Lee prime generation. * @@ -176,10 +183,10 @@ mp *limlee(const char *name, mp *d, mp *newp, if (oev && oev(PGEN_BEGIN, &ev, oec) == PGEN_ABORT) goto fail; + pf.step = 2; if (qql) { dstr_putf(&dn, "%s [+]", name); qq = mprand(d, qql, r, 1); - pf.step = 2; qq = pgen(dn.buf, qq, qq, iev, iec, 0, pgen_filter, &pf, rabin_iters(qql), pgen_test, &rb); } @@ -298,7 +305,7 @@ done: { *nf = nn; *f = vv = xmalloc(nn * sizeof(mp *)); } - + for (i = 0; i < mm; i++) { if (c[i] && vv) *vv++ = v[i]; @@ -330,3 +337,4 @@ fail: return (0); } +/*----- That's all, folks -------------------------------------------------*/