MPI（異步）循環迭代

我有一個類似於下面的程序。在下面的代碼中，所有進程都知道所有其他進程的當前迭代步驟。但是，我很好奇，如果有一種方法可以在沒有諸如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; 
}

來源

2011-07-13 dnanto

MPI_Put不是一個集體呼叫。 – powerrox

啊，我認爲我與MPI_Win_fence混淆...... – dnanto

在使用RMA調用之前，您無法避免使用MPI_Win_fence。它的一個要求。您可以通過將柵欄移出外部for循環來優化代碼。現在的方式是，NOPRECEDE＆NUCCEDED在當前位置不合適。 – powerrox

不需要您使用MPI_Win_fence來同步MPI中的RMA呼叫。也可以使用MPI_Lock（及其變體）來執行所謂的「被動目標」模式。這種模式不太需要同步，因爲它不需要「主動目標」模式所做的所有集體呼叫。

這裏提供完整的答案太複雜了，但您可以在網上找到許多有關該主題的論文。完整的MPI-3標準可用here。看看第11章的片面交流。

來源

2013-08-21 13:09:31

MPI（異步）循環迭代

回答

相關問題