Sieve of Pritchard: Difference between revisions

void Compress(uint32_t w[], uint32_t length, uint32_t N, bool d[], uint32_t imaxto, uint32_t &w_end) {

/* Removes deleted values in w[0..imaxto], and if length < Nto=w_end, updates w_end, otherwise pads with zeros on right */

for (i=1; i <= imaxto; i++) {

if (!d[w[i]]) w[++j] = w[i];{

if (lengthto <== Nw_end) {

for (uint32_t k=j+1; k <= imaxto; k++) w[k] = 0;

Compress(w, length, N, d, (length < N ? w_end : imaxf), w_end);