comptimes

test <- function(n){
  no_cores <- detectCores()- 1
  cl <- makeCluster(no_cores)
  registerDoParallel(cl)
  tree = sim.tree(seed=1)$tree.extant
  pars = c(0.8,0.1,40)
  trees <- foreach(i = 1:n, combine = list) %dopar% dmea02::rec.tree(tree=tree, pars=pars)
  stopCluster(cl)
}

system.time(test(10))

##    user  system elapsed 
##   0.776   0.152   3.280

system.time(test(100))

##    user  system elapsed 
##   0.556   0.184  13.100

#system.time(test(1000))
#system.time(test(10000))