MPI C 中的一般运行时间测量

Question

目标是衡量运行时间与#of 进程。

我只是 MPI 的初学者，在某个地方卡住了。

我编写了一个 hello world 程序并想测试全局运行时。

我尝试使用屏障，以确保所有进程在测量系统时间之前终止，但我遇到了分段错误。

我的代码：

#include <mpi.h> 
#include <stdio.h> 
int main(int argc, char *argv[]) {
  double time1, time2;
  double duration=0.0000;
  int npes, myrank;
  time1 = clock();
  MPI_Init(&argc, &argv);
  MPI_Comm_size(MPI_COMM_WORLD, &npes);
  MPI_Comm_rank(MPI_COMM_WORLD, &myrank);
  printf("From process %d out of %d, Hello World!\n", myrank, npes);
  time2 = clock();
  if (time2-time1>duration) {
    duration = time2-time1;
  }
  duration = time2-time1;
  MPI_BARRIER(MPI_COMM_WORLD);
  printf("runtime is %f ", duration);
  MPI_Finalize();
  return 0;
}

帮我弄清楚为什么会出现分段错误？

Answer 1

我从您的代码中注意到的第一件事是您已经测量了MPI_Barrier 之前的时间，这意味着甚至可以在所有进程打印“hello world"”之前测量运行时间。为了确保正确性测量MPI_Barrier 之后的时间。

您可能还想使用MPI_Wtime() 来测量 MPI 进程中经过的时间。

您的代码只会在每台机器上打印运行时，要计算全局运行时，您必须使用MPI_Reduce。此函数将计算指定的操作（在本例中为 MAX）并将结果存储在根中。

你的代码应该是这样的：

#include <mpi.h> 
#include <stdio.h> 
int main(int argc, char *argv[]) {
   double time1, time2,duration,global;
   int npes, myrank;
   MPI_Init(&argc, &argv);
   time1 = MPI_Wtime();
   MPI_Comm_size(MPI_COMM_WORLD, &npes);
   MPI_Comm_rank(MPI_COMM_WORLD, &myrank);
   printf("From process %d out of %d, Hello World!\n", myrank, npes);
   MPI_Barrier(MPI_COMM_WORLD);
   time2 = MPI_Wtime();
   duration = time2 - time1;
   MPI_Reduce(&duration,&global,1,MPI_DOUBLE,MPI_MAX,0,MPI_COMM_WORLD);
   if(myrank == 0) {
       printf("Global runtime is %f\n",global);
   }
   printf("Runtime at %d is %f \n", myrank,duration);
   MPI_Finalize();
   return 0;
}