蒙特卡罗算法并行计算

因为蒙特卡罗算法求圆周率这个问题没有数据相关性，完全可以用并行计算来模拟之。

首先回顾下蒙特卡罗算法求圆周率，扔了几颗“豆子”，就是说需要执行多少次的随机数产生。比如我们扔了10000次，那么我们要创新获取10000次随机数；我们容忍不了，完全没有数据相关性的事情，完全可以一起做的事情，何必呢！拉出10000人的队伍，没人同时扔一颗不就完事儿了。对的，我们用gpu吧。

开10000个线程，一起来做这件事情！

算法如下，这是《并行计算》里面给出的一个任务分配，很好就用了。

npoints = 10000

circle_count = 0

p = number of tasks

num = npoints/p

find out if I am MASTER or WORKER

do j = 1,num

generate 2 random numbers between 0 and 1

xcoordinate = random1

ycoordinate = random2

if (xcoordinate, ycoordinate) inside circle

then circle_count = circle_count + 1

end do 

if I am MASTER

receive from WORKERS their circle_counts

compute PI (use MASTER and WORKER calculations)

else if I am WORKER

send to MASTER circle_count

endif

 

cuda程序很简单的，自己试着写写。