Julia 笔记（2）

关于多核运算

 

先在当前文档下建立mc_pi.jl文件

function mcpi5(num_try)

　　local sum_pi=0;

　　for i = 1:num_try;

　　　　sum_pi += (rand()^2+rand()^2)<=1 ? 1 : 0;

　　end

　　return (sum_pi*4/num);

end

 

插入以上代码

 

而后在命令行

julia -p 3 打开具有3个worker的julia REPL

 

 

@everywhere include("mc_pi.jl")

将函数导入到每个worker中，因为当先对用户可见的worker只要worker1.

如果在当前REPL中编写函数则只会对worker1可见。

使用remotecall()函数时，会显示

ERROR: On worker 2:
UndefVarError: #mcpi5 not defined

 

分别使用了1,2,3个worker去运行mcpi5

 

num=100000000

julia> @time begin;
r1=remotecall(mcpi5,1,num*5);
r2=remotecall(mcpi5,2,num*5);
print(fetch(r1));
print(fetch(r2));
end

15.70699224 15.7074298 6.795770 seconds (173 allocations: 6.719 KiB)

 

julia> @time begin;
r1=remotecall(mcpi5,1,num*5);
r2=remotecall(mcpi5,2,num*5);
r3=remotecall(mcpi5,3,num*5);
print(fetch(r1));
print(fetch(r2));
print(fetch(r3));
end

15.70741836 15.70805016 15.70769144 13.473580 seconds (288 allocations: 9.984 KiB)

julia> @time begin;
r1=remotecall(mcpi5,1,num*5);
print(fetch(r1));
end

15.70778692 5.999748 seconds (66 allocations: 2.578 KiB)

　

结果可以发现使用1个2个worker的运行时间几乎一样。说明代码时并行在2个cpu核心中的。

而由于本计算机只有2个cpu核心，当运行到第3个worker时，由于没有额外的cpu核心空闲，worker3属于挂起状态。当worker1,2结束工作后，worker3才开始工作。所以时间刚好翻倍。此外从任务管理器中也可以观察到，前6秒core1,2均饱和工作。6-12秒core1饱和工作，core2空闲。

结论

1.remotecall()确实是多核运算，而非单核多线程运算。

2.系统会自动安排worker工作在空闲的core上。因此worker数与计算机core数最好一致，或略大于才能保证最高的cpu利用效率。