并行计算mpi [ PI ]

最新推荐文章于 2020-10-24 16:33:29 发布

原创

最新推荐文章于 2020-10-24 16:33:29 发布 · 5.7k 阅读

7 ·

CC 4.0 BY-SA版权

文章标签：

#intervals #include

该博客通过MPI（Message Passing Interface）展示了并行计算PI的示例代码。首先，通过MPI_Bcast广播变量n，然后根据每个进程的ID进行计算，使用MPI_Reduce将所有进程的局部结果累加到根进程，从而得出全局PI值。计算完成后，根进程输出PI的值和误差，并记录执行时间。实验中使用了12个进程进行计算。

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

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

int main( int argc, char *argv[])
{
    int done = 0, n, myid, numprocs, i;
    double PI25DT = 3.141592653589793238462643;
    double mypi, pi, h, sum, x;
    double startwtime = 0.0, endwtime;
    int namelen;
    char processor_name[MPI_MAX_PROCESSOR_NAME];

    MPI_Init(&argc,&argv);
    MPI_Comm_size(MPI_COMM_WORLD,&numprocs);
    MPI_Comm_rank(MPI_COMM_WORLD,&myid);
    MPI_Get_processor_name(processor_name,&namelen);

    fprintf(stderr,"Process %d on %s\n",