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()))")

您是否即将以演出告罄？