R/primesieve.R

Defines functions primesieve

Documented in primesieve

primesieve <- function(sieved, unsieved) {
    # sieved: sorted vector of sieved numbers
    # unsieved: sorted vector of unsieved numbers

    # cat("sieved", sieved, "\n")
    # cat("unsieved", unsieved, "\n")
    p <- unsieved[1]
    n <- unsieved[length(unsieved)]
    if (p^2 > n) {
        return(c(sieved, unsieved))
    } else {
        unsieved <- unsieved[unsieved %% p != 0]
        sieved <- c(sieved, p)
        return(primesieve(sieved, unsieved))
    }
}

Try the spuRs package in your browser

Any scripts or data that you put into this service are public.

spuRs documentation built on May 2, 2019, 12:44 p.m.