[发明专利]一种基于众核处理器的多重积分计算方法

权利要求书

1.一种基于众核处理器的多重积分计算方法, 其特征在于利用随机数生成以及蒙特卡洛方法中各打靶点间数据无关性的特点，将随机数的生成部分和打靶点函数值计算这两部分计算密集的块交由GPU来完成，将对积分精度以及收敛性判断的操作交由CPU来完成，从而有效的加速多重数值积分的求解，程序架构和组织流程为：

1）设置多重积分的计算精度、第一次计算所需要计算的被积函数值个数n及精度所需新增被积函数值个数δn；

2）调用GPU设备端函数生成N*M个随机数， m为多重积分的维度，将结果存放在设备端显存的线性数组Array中；

3）划分和生成GPU线程空间；

4）调用设备端函数（kernel函数）计算多重被积函数值并归约，为避免访存冲突，计算第n个被积函数值时的第m个靶点从线性数组Array[M*N+n]中选取；

5）将计算完成的多重积分值拷回主内存，并在CPU端判断结果是否已达到所要求精度，以决定是否需要进一步计算；

6）完成计算。

2.根据权利要求1所述的方法，其特征在于，存储在显存中的随机数，先将其扩展到函数的变量空间，然后在GPU的各核上分别计算被积函数的值，并对每个block内的线程归约。

3.根据权利要求2所述的方法，其特征在于积分精度以及收敛性判断是将计算完成的积分结果返回CPU端，并由CPU判断是否满足所要求的计算精度以决定是否生成新的靶点以进一步计算。