私はMPIに非常に慣れていません。MPI_Sendrecv
onp
プロセスを使用して次のプログラムを実行しようとしています。n
各プロセスには、 double 型の任意の値を持つ長さの 4 つの配列 (A、B、C、D) があります。今、これが私がやりたいことですfor i=0,...,p-1
:
プロセス i は、プロセス (i+1) で B を A に送信します%p
プロセス i は、プロセス (i-1) で C を D に送信します%p
プロセス i は (i-1)%p から A に B を受け取ります
プロセス i は (i+1)%p から D に C を受け取ります
次のコードがあります。
#include<stdio.h>
#include<mpi.h>
#include<math.h>
main(int argc, char *argv[]){
/*Initialize the MPI environment */
/******************************************/
MPI_Init(&argc, &argv);
/*variables that will be used for sending*/
/*and receiving messages. */
/*****************************************/
int n=3;
int p;
int myRank;
int comm;
double A[n];
double C[n];
double B[n];
double D[n];
/*Get the number of processors and their */
/*rank. */
/******************************************/
comm = MPI_COMM_WORLD;
MPI_Comm_size(comm, &p);
MPI_Comm_rank(comm, &myRank);
/*Initialize the arrays with the given */
/*initial values. */
/******************************************/
for(int i=0; i<n; i++){
A[i]=0; D[i]=0;
B[i]=-myRank; C[i]=myRank;
}
int left = myRank-1; int right = myRank+1;
if(left<0){
left = p-1;
}
if(right==p){
right = 0;
}
MPI_Sendrecv(&B,n,MPI_DOUBLE,right,1231,
&A,n,MPI_DOUBLE,left,1231,
comm,MPI_STATUS_IGNORE);
MPI_Sendrecv(&C,n,MPI_DOUBLE,left,1232,
&D,n,MPI_DOUBLE,right,1232,
comm,MPI_STATUS_IGNORE);
/*Close the MPI environment. */
/******************************************/
MPI_Finalize();
/*Print the updated values. */
for(int i=0; i<p; i++){
printf("\n");
for(int j=0; j<n; j++){
if(myRank==i){
printf("rank= %d, Updated Values = A[%d]= %f, B[%d]= %f,
C[%d]= %f, D[%d]= %f\n", myRank,j,A[j],j,B[j],j,C[j],j,D[j]);
}
}
}
}
ただし、このコードを実行するたびに、次の結果が得られます。
vincent% mpicc -o hmw2Prob1 hmw2Prob1.c
vincent% mpirun -np 4 hmw2Prob1
2: signal 11 received, exiting..
2: Signal sent from unknown source.
3: signal 11 received, exiting..
3: Signal sent from unknown source.
0: signal 11 received, exiting..
0: Signal sent from unknown source.
1: signal 11 received, exiting..
1: Signal sent from unknown source.
vincent%
ここでいくつかの投稿を読み、そこでいくつかのアイデアを試しましたが、これまでのところ何も機能していません. あなたが提供できるどんな助けにも感謝します。