Eu estava sample
comparando a função em R e comparando-a com igraph:sample_seq
um resultado estranho.
Quando corro algo como:
library(microbenchmark)
library(igraph)
set.seed(1234)
N <- 55^4
M <- 500
(mbm <- microbenchmark(v1 = {sample(N,M)},
v2 = {igraph::sample_seq(1,N,M)}, times=50))
Eu recebo um resultado como este:
Unit: microseconds
expr min lq mean median uq max neval
v1 21551.475 22655.996 26966.22166 23748.2555 28340.974 47566.237 50
v2 32.873 37.952 82.85238 81.7675 96.141 358.277 50
Mas quando eu corro, por exemplo,
set.seed(1234)
N <- 100^4
M <- 500
(mbm <- microbenchmark(v1 = {sample(N,M)},
v2 = {igraph::sample_seq(1,N,M)}, times=50))
Eu obtenho um resultado muito mais rápido para sample
:
Unit: microseconds
expr min lq mean median uq max neval
v1 52.165 55.636 64.70412 58.2395 78.636 88.120 50
v2 39.174 43.504 62.09600 53.5715 73.253 176.419 50
Parece que quando N
é uma potência de 10 (ou algum outro número especial?), sample
É muito mais rápido que outros menores N
que não são potências de 10. Esse é o comportamento esperado ou estou perdendo alguma coisa?