https://github.com/cran/pracma
Revision 26e049d70b4a1c237987e260cba68f6a9413736c authored by Hans W. Borchers on 09 April 2019, 04:10:07 UTC, committed by cran-robot on 09 April 2019, 04:10:07 UTC
1 parent bf07673
Tip revision: 26e049d70b4a1c237987e260cba68f6a9413736c authored by Hans W. Borchers on 09 April 2019, 04:10:07 UTC
version 2.2.5
version 2.2.5
Tip revision: 26e049d
isprime.R
###
### ISPRIME.R Prime number property
###
isprime <- function(x) {
if (is.null(x) || length(x) == 0)
stop("Argument 'x' must be a nonempty vector or matrix.")
if (!is.numeric(x) || any(x < 0) || any(x != round(x)))
stop("All entries of 'x' must be nonnegative integers.")
n <- length(x)
X <- x[1:n]
L <- logical(n)
p <- primes(ceiling(sqrt(max(x))))
for (i in 1:n) {
L[i] <- all(X[i] %% p[p < X[i]] != 0) # all(rem(X[k], p[p < X[k]]))
}
L[X == 1 | X == 0] <- FALSE
N <- as.numeric(L)
dim(N) <- dim(x)
return(N)
}
Computing file changes ...