c语言梯形算法MPI并行数值积分求解pi

下面是一个使用MPI实现的梯形算法数值积分求解pi的C语言代码：

#include <stdio.h>
#include <mpi.h>

double f(double x) {
    return 4.0 / (1.0 + x * x);
}

int main(int argc, char* argv[]) {
    int rank, size, i, n = 100000000;
    double h, x, sum = 0.0, pi;
    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    h = 1.0 / (double)n;
    for (i = rank + 1; i <= n; i += size) {
        x = h * ((double)i - 0.5);
        sum += f(x);
    }
    MPI_Reduce(&sum, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);
    if (rank == 0) {
        pi *= h;
        printf("pi = %lf\n", pi);
    }
    MPI_Finalize();
    return 0;
}

这个代码中，每个进程都计算一部分梯形面积的和，最后使用MPI_Reduce函数将所有进程的和累加起来，得到最终的积分值pi。通过调整n的值可以控制计算精度，n越大，则计算精度越高。