矩阵乘法MPI

Question

在大学时，我被要求一个实验室进行并行编程：并行化矩阵乘法。我写了一个类似于这个例子的算法，但是这个实现提供了一个断点。请告诉我，我能修好些什么？

#define _CRT_SECURE_NO_WARNINGS 
#include <cstdlib> 
#include <iostream> 
#include <math.h> 
#include "mpi.h" 
#include <stdio.h> 
const int N=4; 
int main(int argc, char *argv[]) { 
    int r,q,myid,numprocs; 
    int i0; 
    int *b,*c,*loc_a,*loc_c; 
    MPI_Init(&argc,&argv); 
    MPI_Comm_size(MPI_COMM_WORLD,&numprocs); 
    MPI_Comm_rank(MPI_COMM_WORLD,&myid); 
    MPI_Status status; 
    q=N/numprocs; 
    b=new int [N*N]; 
    c=new int [N*N]; 
    loc_c=new int[N*N]; 
    loc_a=new int[q]; 
    for(int i=0;i<N*N;i++)  { 
     c[i]=0;  loc_c[i]=0;  
    }  
    if(myid==0)  { 
     for (int j = 0; j<numprocs; j++) { 
      for (r = 0; r<q*N; r++) { 
       loc_a[r] = 1; 
      }   
      MPI_Send(&loc_a[0], q*N, MPI_INT, j, 0, MPI_COMM_WORLD); 
     }  
     for (int i = 0; i<N*N; i++) { 
      b[i] = 1; 
     }   
    }  
    MPI_Recv(&loc_a[0], q*N, MPI_INT, 0, 0, MPI_COMM_WORLD, &status); 
    for (r = 0; r<N; r++) { 
     MPI_Bcast(&b[r*N], N, MPI_INT, 0, MPI_COMM_WORLD); 
     i0 = myid*q; 
     for (int i = 0; i<q; i++) { 
      for (int j = 0; j<N; j++) { 
       loc_c[r*N + i0] += loc_a[i*N + j] * b[r*N + j]; 
      }  
      i0++; 
     }  
     MPI_Reduce(loc_c, c, N*N, MPI_INT, MPI_SUM, 0, MPI_COMM_WORLD); 
    }  
    if (myid == 0) { 
     FILE *f = fopen("result.txt", "w"); 
     for (int i = 0; i<N; i++) { 
      for (int j = 0; j<N; j++) { 
       fprintf(f, "%d\t", c[j*N + i]); 
      }  
      fprintf(f, "\n"); 
     }  
     fclose(f); 
    } 
return 0; 
MPI_Finalize(); 
} 

Answer 1

你得到该错误消息似乎非常简单：

要调用本身作为目标（myid == 0和j == 0）而不调用无阻塞第一RECV，这意味着你的进程上的进程阻塞MPI_SEND在继续之前将无限等待自己接收消息，导致死锁。

你有多种选择这里，这取决于你的目标：

1. 通话无阻塞的recv第一

2. 不使用MPI_SEND发送到它本身（检查是否my_id == j和不发如果为true）