CUDA：实现使用蒙特卡罗方法模拟π(附完整源码)

源代码大师

于 2024-03-05 00:08:50 发布

阅读量262

点赞数 6

CC 4.0 BY-SA版权

分类专栏： CUDA实战教程文章标签： CUDA

不予转载，严禁转载，违者必纠。

本文链接：https://blog.csdn.net/it_xiangqiang/article/details/136466531

CUDA实战教程专栏收录该内容

246 篇文章 ¥29.90 ¥99.00

订阅专栏

超级会员免费看

本文提供了一段使用CUDA在GPU上通过蒙特卡罗方法模拟π的代码示例。每个线程生成一个随机点，判断是否在单位圆内，通过统计并原子累加落在圆内的点数，最终估算π的值。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

CUDA：实现使用蒙特卡罗方法模拟π

下面是一个使用蒙特卡罗方法模拟π的示例代码：

#include <iostream>
#include <random>
#include <chrono>

__global__ void calculatePi(int *count, unsigned int seed) {
    int tid = blockIdx.x * blockDim.x + threadIdx.x;

    // Use XORShift32 algorithm for random number generation
    unsigned int x = seed + tid;
    x ^= (x << 13);
    x ^= (x >> 17);
    x ^= (x << 5);

    // Transform the random number to a point in the unit square
    float xPos = (float)(x & 0xFFFF) / 65535.0f;
    float yPos = (float)((x >> 16) & 0xFFFF) / 65535.0f;

    // Check if the point is inside the unit circle
    if (xPos * xPos + yPos * yPos <= 1.0f) {
        atomicAdd(count, 1);
    }
}

int main() {
    const int numPoints = 1000000;
    co