Monte Carlo pi in R when the number of iteration is too large

calc_pi.R

calc_pi <- function(n, k) {
   stopifnot(n %% k == 0)
   r <- n %/% k
   s <- 0
   for (loop in seq_len(r)) {
      x <- runif(k)
      y <- runif(k)
      s <- s + sum(x * x + y * y <= 1)
   }
   4 * s / n
}

system.time(calc_pi(100000000L, 10000L))