如何运行这个编译的开放MPI程序(C)?

时间:2014-02-26 22:00:20

标签: c mpi

我正在尝试在以下网址上运行示例代码。我用“mpicc twoGroups.c”编译了程序,并尝试将其作为“./a.out”运行,但得到以下消息:必须指定MP_PROCS = 8.终止。

我的问题是你如何设置MP_PROCS = 8?

此处的小组和通讯例行示例https://computing.llnl.gov/tutorials/mpi/

#include "mpi.h"
#include <stdio.h>
#define NPROCS 8

main(int argc, char *argv[])  {
int        rank, new_rank, sendbuf, recvbuf, numtasks,
           ranks1[4]={0,1,2,3}, ranks2[4]={4,5,6,7};
MPI_Group  orig_group, new_group;
MPI_Comm   new_comm;

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

if (numtasks != NPROCS) {
  printf("Must specify MP_PROCS= %d. Terminating.\n",NPROCS);
  MPI_Finalize();
  exit(0);
  }

sendbuf = rank;

/* Extract the original group handle */
MPI_Comm_group(MPI_COMM_WORLD, &orig_group);

/* Divide tasks into two distinct groups based upon rank */
if (rank < NPROCS/2) {
  MPI_Group_incl(orig_group, NPROCS/2, ranks1, &new_group);
  }
else {
  MPI_Group_incl(orig_group, NPROCS/2, ranks2, &new_group);
  }

/* Create new new communicator and then perform collective communications */
MPI_Comm_create(MPI_COMM_WORLD, new_group, &new_comm);
MPI_Allreduce(&sendbuf, &recvbuf, 1, MPI_INT, MPI_SUM, new_comm);

MPI_Group_rank (new_group, &new_rank);
printf("rank= %d newrank= %d recvbuf= %d\n",rank,new_rank,recvbuf);

MPI_Finalize();
}

1 个答案:

答案 0 :(得分:2)

执行MPI程序时,需要使用适当的包装器。大多数时候它看起来像这样:

mpiexec -n <number_of_processes> <executable_name> <executable_args>

所以对于你的简单例子:

mpiexec -n 8 ./a.out

您还会看到使用mpirun代替mpiexec-np代替-n。两者在大多数情况下都很好。

如果您刚刚开始,那么确保您使用的是最新版本的MPI也是一个好主意,这样您就不会遇到旧的错误或奇怪的执行环境。 MPICH和Open MPI是两种最流行的实现。 MPICH刚刚发布3.1版here,而Open MPI版本1.7.4可用here。您通常也可以通过友好的社区包经理获得其中任何一个。