mpi编程的计算π的程序

以下是使用MPI编写的估算π值的程序：

#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;
    double pi, sum = 0.0, x;
    const int n = 100000000;

    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);

    for (i = rank; i < n; i += size) {
        x = (i + 0.5) / n;
        sum += f(x);
    }

    MPI_Reduce(&sum, &pi, 1, MPI_DOUBLE, MPI_SUM, 0, MPI_COMM_WORLD);

    if (rank == 0) {
        pi /= n;
        printf("pi = %.16f\n", pi);
    }

    MPI_Finalize();

    return 0;
}

该程序使用MPI实现了一个简单的求π值的算法。每个进程都计算一部分小矩形的面积，然后将它们的和相加，最后通过MPI_Reduce函数将结果发送给rank为0的进程，该进程计算最终的π值并输出结果。

要编译该程序，可以使用以下命令：

mpicc -o pi pi.c

然后可以使用以下命令运行程序：

mpirun -np 4 ./pi

其中“-np 4”指定使用4个进程运行程序。