MPI:反复并行化缓冲区

Question

假设我有一个大的输入文件。

假设这个文件有我想并行处理的项目。

std::vector<std::string> items(100000,"");
for(int i = 0; i < 1000000; i++)
items[i] = pop_item(file);

接下来，我想通过与MPI:并行处理这些项目来加快处理速度

std::vector<MyObj> processed_items(100000); // pseudo-code, i handle the memory mallocing
int size; rank;
MPI_INIT();
MPI_Comm_size(MPI_COMM_WORLD,&size);
MPI_Comm_rank(MPI_COMM_WORLD,&rank);
for(i = rank; i < 100000; i += size)
processed_items[i] = process_item(items[i]);
MPI_FINALIZE();

好的，太好了，很管用。

现在，我想在一段时间内反复做一遍：

while(!done){
done = fill_items(&items, file); 
MPI_INIT();
...;
MPI_FINALIZE();
print_items(&processed_items);
}

然而，我失败了，返回"error:mpi_init called after mpi finalize invoked">

---

我在MPI中处理此问题的预期方式是什么？

Answer 1

MPI_INIT((和MPI_FINALIZE每个程序只能调用一次，正如您的错误所暗示的那样。这个半年前的老答案概述了如何让MPI并行运行程序的某些部分：

int main(int argc, char *argv[]) {
MPI_Init(&argc, &argv);  
MPI_Comm_size(MPI_COMM_WORLD,&numprocs);  
MPI_Comm_rank(MPI_COMM_WORLD,&myid);
if (myid == 0) { // Do the serial part on a single MPI thread
printf("Performing serial computation on cpu %dn", myid);
PreParallelWork();
}
ParallelWork();  // Every MPI thread will run the parallel work
if (myid == 0) { // Do the final serial part on a single MPI thread
printf("Performing the final serial computation on cpu %dn", myid);
PostParallelWork();
}
MPI_Finalize();  
return 0;  
}