如何计算π的求和值？OpenMP实现代码出错求助

问题根源分析 #

• 整数除法导致采样点错误：代码中i/N和(i+1)/N为整数除法（i和N均为int类型），当i < N时结果恒为0，仅i=N-1时(i+1)/N等于1。这使得绝大多数情况下调用f(0)，完全偏离梯形法需要的细分区间采样逻辑，最终结果必然错误。
• 未启用OpenMP并行机制：你提到要使用OpenMP计算，但当前代码无任何OpenMP指令，无法实现并行加速。

修正后的代码 #

#include <chrono>
#include <iostream>
#include <omp.h> // 引入OpenMP头文件

double f(double x)
{
    return 4.0 / (1.0 + x * x); // 使用double字面量，避免整数除法截断
}

int main()
{
    const int N = 60;
    double pi = 0.0;
    const double s = 1.0 / N; // 明确步长s，符合你给出的公式定义

    // OpenMP并行循环，通过归约操作安全累加pi，避免线程竞争
    #pragma omp parallel for reduction(+:pi)
    for (int i = 0; i < N; ++i)
    {
        double x_i = static_cast<double>(i) / N; // 强制类型转换，执行浮点数除法
        double x_i1 = static_cast<double>(i + 1) / N;
        // 梯形法面积公式：(f(x_i) + f(x_i1)) * 步长 / 2
        pi += (f(x_i) + f(x_i1)) * s / 2.0;
    }

    printf("Pi = %f\n", pi);
    return 0;
}

关键修正说明 #

1. 修复除法类型：将i/N改为static_cast<double>(i)/N，确保执行浮点数除法，生成正确的区间采样点。
2. 优化函数计算：f函数中使用4.0和1.0代替整数字面量，避免返回值因整数除法被错误截断。
3. 启用OpenMP并行：添加#pragma omp parallel for reduction(+:pi)，实现循环并行化，同时通过归约操作保证多线程下pi的累加安全。
4. 贴合公式逻辑：明确定义步长s = 1.0/N，让代码逻辑与你给出的求和公式一致。

内容的提问来源于stack exchange，提问作者IsLearning