Julia 从标准正态分布中绘制

Julia drawing from standard normal distribution

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

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

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

可能是分配变量 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()))")

您是否在演出中快要用完内存了?