查找具有最大值 mpi 的进程

Question

我正在尝试使用 mpi 点对点通信在线程中找到最大值（也想添加最小值，但直到我弄清楚这一点） 我假设进程形成了一个环，从 0 开始发送到他的右边，即进程 1 我只希望具有最大值的进程打印结果，但不幸的是在某些情况下 2 进程打印出结果 谁能告诉我问题出在哪里？？？...

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

int main(int argc, char **argv)
{
    int rank, size, tag = 100;
    int rightrank, leftrank;
    int v, min, max, i, tmpmax, tmpmin;
    int *data;
    
    MPI_Init(&argc, &argv);
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    MPI_Comm_size(MPI_COMM_WORLD, &size);
    MPI_Status status;
    data = (int*) malloc(size * sizeof(int));
    
    rightrank = rank + 1;
    leftrank = rank - 1;
    
    if (rank == 0) {
        leftrank = size - 1;
    }
    if (rank == size - 1) {
        rightrank = 0;
    }
    
    for (i = 0; i < size; i++) {
        data[i] = rand() % 60;
    }
    v = data[rank];
    
//v=rank*10;
    printf("{rank %d , v = %d }\n", rank, v);
    
    if (rank == 0) {
        
        MPI_Send(&v, 1, MPI_INT, rightrank, tag, MPI_COMM_WORLD);
        
        MPI_Recv(&tmpmax, 1, MPI_INT, leftrank, tag, MPI_COMM_WORLD, &status);
    }
    else {
        
        MPI_Recv(&tmpmax, 1, MPI_INT, leftrank, tag, MPI_COMM_WORLD, &status);
        max = tmpmax;
        
        if (v > max) {
            max = v;
        }
        
        MPI_Send(&max, 1, MPI_INT, rightrank, tag, MPI_COMM_WORLD);
        
    }
    
    if (v == max) {
        printf("Me process of rank %d i have max value v = %d\n ", rank, v);
    }
    
    MPI_Finalize();
    return 0;
}

输出 ： as u can see rank 1 and 2 print ...

Answer 1

无论如何，非常感谢你，我已经弄清楚了...... 我的解决方案是让进程 0 循环并检查哪个索引具有最大值，如下所示：

#include <stdlib.h>
#include<unistd.h>
#include <mpi.h>

int main(int argc,char **argv)
{
int rank,size,tag=100;
int rightrank,leftrank;
int v, min, max,i,tmpmax,tmpmin;
int* data;

MPI_Init(&argc,&argv);
MPI_Comm_rank(MPI_COMM_WORLD,&rank);
MPI_Comm_size(MPI_COMM_WORLD,&size);
MPI_Status status;
data=(int*)malloc(size*sizeof(int));

rightrank=rank+1;
leftrank=rank-1;

if(rank==0){
leftrank=size-1;
}
if(rank==size-1){
rightrank=0;
}

for(i=0 ; i < size ; i ++){
    data[i]=rand()%80;
}
v=data[rank];

//v=rank*10;
printf("{rank %d , v = %d }
",rank,v);

if(rank==0){

    MPI_Send(&v,1,MPI_INT,rightrank,tag,MPI_COMM_WORLD);
    MPI_Send(&v,1,MPI_INT,rightrank,tag,MPI_COMM_WORLD);
    
    MPI_Recv(&tmpmax,1,MPI_INT,leftrank,tag,MPI_COMM_WORLD,&status);
    MPI_Recv(&tmpmin,1,MPI_INT,leftrank,tag,MPI_COMM_WORLD,&status);

}
else{
MPI_Recv(&tmpmax,1,MPI_INT,leftrank,tag,MPI_COMM_WORLD,&status);
max=tmpmax;
if(v>max){
max=v;
}
MPI_Send(&max,1,MPI_INT,rightrank,tag,MPI_COMM_WORLD);
sleep(0.1);
/////cal of min
MPI_Recv(&tmpmin,1,MPI_INT,leftrank,tag,MPI_COMM_WORLD,&status);
min=tmpmin;
if(v<min){
min=v;
}
MPI_Send(&min,1,MPI_INT,rightrank,tag,MPI_COMM_WORLD);
}

if(rank==0){
for(int i = 0 ; i < size ; i ++){
    if(tmpmax==data[i]){
    printf("Me Process of rank %d i have the max value ,v = %d
",i,tmpmax);
    }
}
}

if(rank==0){
for(int i = 0 ; i < size ; i ++){
    if(tmpmin==data[i]){
    printf("Me Process of rank %d i have the min value ,v = %d
",i,tmpmin);
    }
}
}

MPI_Finalize();
    return 0;
}