performance - Julia 从标准正态分布中绘制

Question

我需要从标准正态分布中提取 53000000 个观测值。我当前的代码需要很长时间才能在 Julia 中运行(实际上，它已经运行了过去二十分钟)，我想知道我是否可以做些什么来加快它的速度。这是我尝试过的:

using Distributions
d = Normal()
shock = rand(d, 1, 53000000)

当我在 REPL 中执行代码时(我在 Juno/Atom 中工作)，代码会立即工作，但是当我逐步使用调试器时，此时会滞后(从标准法线绘制)。所以我认为调试器可能是这里的真正罪魁祸首。

Answer 1

可能是变量 shock 的分配使用了 1/2 gig 的内存有时会在加载调试器时导致交换。

尝试在调试器中运行它以查看:

using Distributions, Base.Sys

println("Free memory is $(Int(Sys.free_memory()))")
d = Normal()
shock = rand(d, 1, 53000000)
println("shock uses $(sizeof(shock)) bytes.")
println("Free memory is $(Int(Sys.free_memory()))")

你在演出中接近内存不足了吗？