CコードでMPIフィルタを実装する方法は？

Question

私は下のフィルタコードのMPIを実装しようとしていますが、それをやるのが難しいです。これは、MPIの6つの機能を以下の間、私が試したものである

int A[100000][100000]; 
int B[100000][100000]; 

for (int i=1; i<(100000 - 1); i++) 
for (int i=1; j<(100000 - 1); j++) 
    B[i][j] = A[i-1][j] + A[i+1][j] + A[i][j-1] + A[i][j+1] - 4*A[i][j]; 

：：どのようにそれは？：

フィルタのコードを実行する必要があります

int myrank; /* Rank of process */ 
    int numprocs; /* Number of processes */ 
    int source; /* Rank of sender */ 
    int dest; /* Rank of receiver */ 

    char message[100]; /* Storage for the message */ 
    MPI_Status status; /* Return status for receive */ 
    MPI_Init(& argc, & argv); 
    MPI_Comm_size(MPI_COMM_WORLD, & numprocs); 
    MPI_Comm_rank(MPI_COMM_WORLD, & myrank); 

    if (myrank != 0) 
    { 
     dest = 0; 
     MPI_Send(message, strlen(message) + 1, 
      MPI_CHAR, dest, 15, MPI_COMM_WORLD); 
     } else { 
     for (source = 1; source < numprocs; source++) { 
      MPI_Recv(message, 100, MPI_CHAR, source, 
      15, MPI_COMM_WORLD, & status); 
     } 
     } 
     MPI_Finalize(); 

Answer 1

私はこのように行くと思います。まず第一に、私はこのコードを持っていました

int A[100000][100000]; 
int B[100000][100000]; 

は動的割り当てに置き換えられました。それぞれのプロセスごとにすべてのメモリが必要なわけではありません。

次に、配列Aを別のプロセスに送信します。行ごと。

データフレーム（行数）の「高さ」とは何か：

delta = (100000 - 2)/(numprocs-1);  // we don't count first and last row 
reminder = (100000 - 2) % (numprocs-1); // it might be that we need to give 
             // little bit more to calculate 
             // to one of the processes 

// we are starting from row with idx=1 (second row) and we want to finish when 
// we hit last row 
if(myrank == 0) { 
    for(int i=1; i < numprocs; i++) { 
    // +100000 - we need two more rows to calculate data 
    int how_many_bytes = delta * 100000 + 200000; 
    if(reminder != 0 && i == (numprocs-1)) { 
     how_many_bytes += reminder * 100000; 
    } 
    MPI_Send(&(A[(i-1)*delta][0]), how_many_bytes, MPI_INT, i, 0, 
       MPI_COMM_WORLD); 
    } 
} else { 
    // allocate memory for bytes 
    int *local_array = NULL; 
    int how_many_bytes = delta * 100000 + 200000; 
    if(reminder != 0 && i == (numprocs-1)) { 
    how_many_bytes += reminder * 100000; 
    } 
    local_array = malloc(how_many_bytes * sizeof(int)); 
    MPI_Status status; 

    MPI_Recv(
    local_array, 
    how_many_bytes, 
    MPI_INT, 
    0, 
    0, 
    MPI_COMM_WORLD, 
    &status); 
} 

// perform calculations for each and every slice 
// remembering that we always have on extra row on 
// top and one at the bottom 
// send data back to master (as above, but vice versa).