неожиданные результаты: microbenchmark

Question

Я всегда был немного обеспокоен отсутствием точности. Я вижу табличку с маркировкой system.time и rbenchmark (в том, что точность синхронизации может отсутствовать), и недавно увидел Hadley ссылку на пакет microbenchmark. Поэтому я решил дать ему вихрь, как показано ниже. Я провел mean против f <- function(x) {sum(x)/length(x)} и ожидал mean, чтобы сделать способ лучше, чем f, но результаты, как я их понимаю, не указывают, что это правда.

1. Я не понимаю результаты?
2. Действительно ли это быстрее, чем означает?
3. Микробиблиотека все еще находится в бета-фазе, и ее нужно гладить ?

Я запускаю R2.15 на машине с выигрышем 7 (так как микробиблиотека выполняет тайминги по-разному в зависимости от вашей ОС).

Результаты

Unit: microseconds 
    expr min  lq median  uq max 
1 f(x) 19.130 20.529 20.529 20.996 286.00 
2 mean(x) 28.927 29.860 30.327 30.327 672.31 

Кодекс

library(microbenchmark) 

x <- 1:10000 
f <- function(x) {sum(x)/length(x)} 
mean(x) 

res <- microbenchmark(
    mean(x), 
    f(x), 
times=1000L) 

print(res) 
boxplot(res) 

Answer 1

Я могу ошибаться, но это не кажется, что для меня удивительно. Перед тем, как mean.default может позвонить .Internal(mean(x)), необходимо проверить операторы 0 if, рассчитать длину x, а затем проверить еще одно сообщение if. И разница во времени это довольно мало.

Вызов .Internal(mean(x) непосредственно немного еще быстрее:

library(microbenchmark) 

x <- 1:10000 
f1 <- function(x) {sum(x)/length(x)} 
f2 <- function(x) {.Internal(mean(x))} 

res <- microbenchmark(
    mean(x), 
    f1(x), 
    f2(x), 
times=1000L) 

print(res) 

Unit: microseconds 
    expr min  lq median  uq  max 
1 f1(x) 32.195 32.4605 32.8850 33.4645 106.997 
2 f2(x) 21.840 22.0580 22.2015 22.6270 55.316 
3 mean(x) 35.393 35.9840 36.1860 36.4420 91.203 

Answer 2

Я думаю, вы обнаружите, что если вы увеличите до размера X с коэффициентом 10, вы увидите более стабильные результаты. Честно говоря, я был бы удивлен, если бы вы действительно смогли получить микросекундную точность времени на компьютере с многозадачной операционной системой.

Вы могли бы также рассмотреть следующие вопросы:

• Вы работаете на ноутбуке или компьютере, который имеет автоматическое масштабирование частоты процессора?
• Warmup?
• Привязка вашего процесса к одному ядру.