gistfile1.txt

struct
{
    mp_limb_t p;
    mp_limb_t i;
    void * sieve;
}
n_primes_iter_t

static __inline__ void
n_primes_iter_init(n_primes_iter_t iter)
{
    iter->p = 2;
    iter->i = 0;
    iter->sieve = NULL;
}

static __inline__ mp_limb_t
n_primes_iter_next(n_primes_iter_t iter)
{
    iter->i++;

    if (iter->i < FLINT_MAX_PRIMES)
    {
        iter->p = flint_primes[iter->i];
    }
    else
    {
        /* generate 1M (or whatever) primes at a time */
        if (out_of_primes(iter->sieve))
            grow(iter->sieve);
        iter->p = next_prime(sieve);
    }

    return iter->p;
}