-
-/*----- Data structures ---------------------------------------------------*/
-
-typedef struct gcdctx {
- mp *q;
- mp *r;
- pfilt p;
-} gcdctx;
-
-/*----- Custom stepper ----------------------------------------------------*/
-
-static int gcdstep(int rq, pgen_event *ev, void *p)
-{
- gcdctx *g = p;
- int rc = PGEN_ABORT;
- mp *z = MP_NEW;
-
- switch (rq) {
- case PGEN_BEGIN: {
- mp *p = ev->m = mp_split(ev->m);
- p->v[0] |= 3;
- g->q = mp_lsr(MP_NEW, p, 1);
- rc = pfilt_create(&g->p, p);
- } break;
- case PGEN_TRY:
- g->q = mp_add(g->q, g->q, MP_FOUR);
- rc = pfilt_step(&g->p, 4);
- break;
- case PGEN_DONE:
- pfilt_destroy(&g->p);
- mp_drop(g->q);
- return (PGEN_DONE);
- }
- for (;;) {
- if (rc != PGEN_FAIL) {
- mp_gcd(&z, 0, 0, g->r, g->q);
- if (MP_CMP(z, !=, MP_ONE))
- rc = PGEN_FAIL;
- }
- if (rc != PGEN_FAIL)
- break;
- g->q = mp_add(g->q, g->q, MP_FOUR);
- rc = pfilt_step(&g->p, 4);
- }
-
- mp_drop(z);
- mp_drop(ev->m);
- ev->m = MP_COPY(g->p.m);
- return (rc);
-}