2011-07-13 126 views
0

我有一个类似于下面的程序。在下面的代码中,所有进程都知道所有其他进程的当前迭代步骤。但是,我很好奇,如果有一种方法可以在没有诸如MPI_PUT,等集体调用的情况下执行此操作,特别是在每个进程以不同的速率迭代的情况下。MPI(异步)循环迭代

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


/* OUTPUT 
* $ mpirun -np 4 so00.exe 
* @[0]: p |  0  1  2  3 
* @[0]: p |  4  5  6  7 
* @[0]: p |  8  9  10  11 
* @[0]: p |  12  13  14  15 
* @[0]: p |  16  17  18  19 
* @[1]: p |  0  1  2  3 
* @[1]: p |  4  5  6  7 
* @[1]: p |  8  9  10  11 
* @[1]: p |  12  13  14  15 
* @[1]: p |  16  17  18  19 
* @[2]: p |  0  1  2  3 
* @[2]: p |  4  5  6  7 
* @[2]: p |  8  9  10  11 
* @[2]: p |  12  13  14  15 
* @[2]: p |  16  17  18  19 
* @[3]: p |  0  1  2  3 
* @[3]: p |  4  5  6  7 
* @[3]: p |  8  9  10  11 
* @[3]: p |  12  13  14  15 
* @[3]: p |  16  17  18  19 
*/ 
int main(int argc, char *argv[]) 
{ 
    int i, n, rank, np; 

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

    //////////////////////////////////////////////////////////////////////////// 

    int *pos; 
    MPI_Win win; 
    MPI_Alloc_mem(sizeof(int)*np, MPI_INFO_NULL, &pos); 
    MPI_Win_create(pos, np, sizeof(int), MPI_INFO_NULL, MPI_COMM_WORLD, &win); 

    for (i=rank; i<(np*5); i+=np) 
    { 
     MPI_Win_fence(MPI_MODE_NOPRECEDE, win); 
     for (n = 0; n < np; n++) 
     { 
      MPI_Put(&i, 1, MPI_INT, n, rank, 1, MPI_INT, win); 
     } 
     MPI_Win_fence((MPI_MODE_NOSTORE | MPI_MODE_NOSUCCEED), win); 
     printf("@[%d]: p | ", rank); 
     for (n = 0; n < np; n++) printf("\t%d", pos[n]); 
     printf("\n"); 
    } 

    //////////////////////////////////////////////////////////////////////////// 

    MPI_Win_free(&win); 
    MPI_Free_mem(pos); 

    MPI_Finalize(); 

    return EXIT_SUCCESS; 
} 
+0

MPI_Put不是一个集体呼叫。 – powerrox

+0

啊,我认为我与MPI_Win_fence混淆...... – dnanto

+0

在使用RMA调用之前,您无法避免使用MPI_Win_fence。它的一个要求。您可以通过将栅栏移出外部for循环来优化代码。现在的方式是,NOPRECEDE&NUCCEDED在当前位置不合适。 – powerrox

回答

0

不需要您使用MPI_Win_fence来同步MPI中的RMA呼叫。也可以使用MPI_Lock(及其变体)来执行所谓的“被动目标”模式。这种模式不太需要同步,因为它不需要“主动目标”模式所做的所有集体呼叫。

这里提供完整的答案太复杂了,但您可以在网上找到许多有关该主题的论文。完整的MPI-3标准可用here。看看第11章的片面交流。