MPI_BCastからデータを受信する方法

Question

重力の干渉が問題です。システムにN個の粒子があり、M個のプロセスがある。 現在のプロセスの個々のブロックにパーティクル（すべて11個のパラメータがあります）の新しい位置を計算し、その後、すべての他のプロセスに新しいデータをブロードキャストします。

double * particles; 
... 
int startForProcess = numberOfParticlesPerThread * 11 * rank; 
int endForProcess = startForProcess + numberOfParticlesPerThread * 11; 
calculateNewPosition(particles, startForProcess, endForProcess, nParticles, rank); 
MPI_Barrier(MPI_COMM_WORLD); 
MPI_Bcast(particles+startForProcess, endForProcess - startForProcess, MPI_DOUBLE, rank, MPI_COMM_WORLD); 

残念ながら、各スレッドで私はこのスレッドで行われた変更を見ることができます。 はここに私のコードです。プロセス間に通信はありません。

私に間違っていると教えてください。

Answer 1

MPI_Bcastは、ルートからコミュニケータの残りのランクにデータをブロードキャストします。 MPI_Bcastコール内のすべてのプロセスは、同じプロセスからデータをすべて受信しているため、同じルートパラメータを使用する必要があります。

あなたが探している機能はMPI_Allgatherだと思います。更新されたデータがすべてのプロセス間で収集され、伝達されます。

コードは次のようなものになります：余談として

double * particles; 
double * my_particles; 
... 
int startForProcess = numberOfParticlesPerThread * 11 * rank; 
int endForProcess = startForProcess + numberOfParticlesPerThread * 11; 
calculateNewPosition(particles, startForProcess, endForProcess, nParticles, rank); 

memcpy(my_particles, particles + startForProcess * sizeof(double), 
     (endForProcess - startForProcess) * sizeof(double)); 
MPI_Allgather(my_particles, endForProcess - startForProcess, MPI_DOUBLE, 
       particles, endForProcess - startForProcess, MPI_DOUBLE, MPI_COMM_WORLD); 

をもちろん、あなたが他のいくつかの理由で同期が必要な場合を除き、あなたは、集団的呼び出しの前MPI_Barrierは必要ありません。