为什么微基准的第一次运行总是最慢的？答案

【问题标题】：Why is the first run of microbenchmark always the slowest?为什么微基准的第一次运行总是最慢的？
【发布时间】：2020-09-28 16:32:04
【问题描述】：

使用Microbenchmark 时，我注意到第一次执行总是比其他执行慢很多。这种效果在不同的机器和不同的功能上是相同的。这是否与图书馆有关，或者这是预期的某种预热？

library(microbenchmark)
X <- matrix(rnorm(100), nrow = 10)
microbenchmark(solve(X))$time
#[1] 82700 23700 18300 17700 19700 19100 16900 17500 17300 16600 16700 16700 18500 16900 17700 16900 17000 16200 17400 17000 16800 16600 17000 16700 16800 17100
#[27] 17300 17100 16800 17800 17400 18100 17400 18100 18000 16700 17400 17300 17000 16800 16400 17300 16700 16900 16900 16700 17200 17800 16600 17100 16800 17800
#[53] 17000 17200 17500 17200 17200 17300 17800 17600 17600 17200 16600 16700 16800 16600 16400 16500 17300 17600 16800 17600 16300 16800 17100 16500 16800 16700
#[79] 16300 16700 16300 16700 16800 16700 16400 17100 16400 17100 17000 18000 16600 16600 16600 16800 16700 16500 17600 19100 17400 16900

【问题讨论】：

这与预热时间有关，请参阅help('microbenchmark')，部分详细信息，参数control。如果增加预热迭代次数，第一次运行可能不是最慢的：microbenchmark(solve(X), control = list(warmup = 10))$time。
@RuiBarradas 你应该把它作为答案发布

标签： r microbenchmark

【解决方案1】：

和预热时间有关，见help('microbenchmark')，章节详情，参数控制：

控制列表可以包含以下条目：

订购
[省略]

热身
在实际基准测试之前执行的预热迭代次数。这些用于估计时序开销以及将处理器从可能处于的任何睡眠或空闲状态中启动。默认值为 2。

如果您增加预热迭代的次数，第一次运行可能不是最慢的，尽管很多次都是。

library(microbenchmark)

set.seed(2020)
X <- matrix(rnorm(100), nrow = 10)
times <- microbenchmark(solve(X), control = list(warmup = 10))$time
times
#  [1] 145229  72724  65333  65305 115715  63797 689113  72101  64830  66392
# [11]  65776  66619  65531  64765  65351  65605  65745  65106  64661  65790
# [21]  65435  64964  66138  65952  66893  65654  65585  75141  74666  69060
# [31]  72725  66650  65486  65894  66808  65381  66039  65959  64842  65029
# [41]  65673  66439  64394  70585  68899  73875  73180  67807  65891  65699
# [51]  64693  63679  65504  80190  66150  65048  64372  64842  65845  65144
# [61]  65543  65297  65485  64695  66580  64921  65453  64840  65559  65805
# [71]  64362  66098  65464  65227  64998  64007  65659  63919  64727  64796
# [81]  65231  64030  65871  65735  64217  65195  65181  65130  66015  63891
# [91]  63755  65274  65116  64573  64244  64214  64148  64457  65346  64228

现在看看order是第一个：

order(times, decreasing = TRUE)
#  [1]   7   1   5  54  28  29  46  47  31   2   8  44  30  45  48  25  35  32  12
# [20]  65  42  10  55  23  72  37  89  38  24  34  49  83  59  70  20  11  17  84
# [39]  50  41  77  26  16  27  69  61  13  53  33  63  73  67  21  36  15  99   3
# [58]   4  62  92  81  74  86  87  60  88  93  18  56  40  75  22  66  39  58  68
# [77]   9  80  14  79  64  51  19  94  98  43  57  71  95 100  85  96  97  82  76
# [96]  78  90   6  91  52

在这种情况下，最慢的是第七次，而不是第一次。

【讨论】：