用mpi定义全局变量

时间:2013-05-21 14:30:52

标签: c++ c mpi

我在下面写了一个示例代码:

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

double x;

int main (int argc, char **argv) { 
   MPI_Init(&argc, &argv); 
   MPI_Comm_rank(MPI_COMM_WORLD, &rank); 
   MPI_Comm_size(MPI_COMM_WORLD, &size); 

   if (rank==0) x=10.1;

   MPI_Barrier(MPI_COMM_WORLD);
   printf("%f\n", x);

   MPI_Finalize(); 
   return 0; 
}

正如人们可能注意到的,这个程序实际上定义了一个名为x的全局变量,而第0个线程试图为它分配一些值。当我在具有4个内核的SMP(对称多处理)机器上运行该程序时,我得到以下结果:

10.1
0
0
0

更有趣的是,当我更改我的代码以便每个线程打印变量x的地址,即&amp; x时,它们都会打印相同的内容。

我的问题是SMP系统上的一些线程如何共享变量地址的相同值,而它们不共享相同的值?

我的第二个问题是我应该如何更改上面的代码,以便得到以下结果?

10.1
10.1
10.1
10.1

1 个答案:

答案 0 :(得分:5)

您可以使用广播:

MPI_Bcast(&x,1,MPI_DOUBLE,0,MPI_COMM_WORLD);

这会将进程0上的x值发送给所有其他进程。