MPI_Barrier和递归

我想使用MPI_Barrier（OpenMPI）来强制所有进程处于递归调用的相同深度。MPI_Barrier和递归

这是代码

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

using namespace std; 

void recursive_function(int,int,int); 

int main(int argc, char *argv[]) { 
    int rank, size; 

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

    recursive_function(0,3,rank); 

    MPI_Finalize(); 
} 

void recursive_function(int depth, int limit, int rank) { 
    printf("depth: %d, processor %d\n", depth, rank); 

    MPI_Barrier(MPI_COMM_WORLD); 
    if(depth == limit) return; 
    else recursive_function(depth+1,limit,rank); 
}

我得到的是（与正在运行的mpirun -np 2垒）

depth: 0, processor 0 
depth: 1, processor 0 
depth: 2, processor 0 
depth: 3, processor 0 
depth: 0, processor 1 
depth: 1, processor 1 
depth: 2, processor 1 
depth: 3, processor 1

虽然我希望像

depth: 0, processor 0 
depth: 0, processor 1 
depth: 1, processor 0 
depth: 1, processor 1 
depth: 2, processor 1 
depth: 2, processor 0 
depth: 3, processor 1 
depth: 3, processor 0

来源

2011-06-30 mambro

有不能保证MPI进程之间的stdout以任何方式排序。

也就是说，您需要让进程进行通信以证明它们都处于相同的递归深度。例如。在屏障之后，每个进程！= 0发送一条消息给等级为0的打印内容。

来源

2011-06-30 17:07:45 janneb

+1，但也许值得一提的是代码看起来是正确的，因为没有进程可以下降到下一级递归，直到所有进程完成当前深度。只是stdout的无序性会扰乱屏幕上的结果。 – suszterpatt

谢谢你的帮助。问题只是标准输出。代码实际上是正确的。 – mambro

MPI_Barrier和递归

回答

相关问题