mpi程序在c编译但不运行

Question

我是MPI编程的初学者。我试图执行一个代码，其中使用A new [i] [j] = Aold [i] [j] + Aold [i] [j + 1]从旧的矩阵形成新的矩阵，并且这发生在10倍。我写了下面的代码：

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


int main(int argc, char *argv[]) 

{ 
void calculate(); 
int taskid,numtasks,numworkers,rowmean,offset,t, 
BEGIN=1, 
msgtype, 
start, 
RTAG=3, 
LTAG=2, 
NX=20, 
NY=20, 
MASTER=0, 
NONE=0, 
DONE=4, 
TIMESTEP=10, 
source, 
ext, 
end, 
left, 
rows, 
right, 
dest; 
MPI_Status status; 
float u[2][NX][NY]; 

MPI_Init(&argc,&argv); 
MPI_Comm_rank(MPI_COMM_WORLD,&taskid); 
MPI_Comm_size(MPI_COMM_WORLD,&numtasks); 
numworkers=numtasks-1; 

rowmean=NX/numworkers; 
ext=NX%numworkers; 

if (taskid==MASTER) 
{ 
    offset=0; 
    for (int i=1;i<=numworkers;i++) 
    { 
     if (i==1) 
     { 
      left=NONE; 
     } 
     else 
     { 
      left=i-1; 
     } 
     if (i==numworkers) 
     { 
      right=NONE; 
     } 
     else 
     { 
      right=i+1; 
     } 
     dest=i; 
     rows=(i<=ext) ? rowmean+1:rowmean; 
     MPI_Send(&u[0][offset][0],rows*NY,MPI_FLOAT,dest,BEGIN,MPI_COMM_WORLD); 
     MPI_Send(&offset,1,MPI_INT,dest,BEGIN,MPI_COMM_WORLD); 
     MPI_Send(&rows,1,MPI_INT,dest,BEGIN,MPI_COMM_WORLD); 
     MPI_Send(&right,1,MPI_INT,dest,BEGIN,MPI_COMM_WORLD); 
     MPI_Send(&left,1,MPI_INT,dest,BEGIN,MPI_COMM_WORLD); 
     offset=offset+rows; 
    } 
    for(int i=1;i<=numworkers;i++) 
    { 
     source=i; 
     msgtype=DONE; 
     MPI_Recv(&u[1][offset][0],rows*NY,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
     MPI_Recv(&offset,1,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
     MPI_Recv(&rows,1,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
    } 
    MPI_Finalize(); 
    /*end of master code*/ 
} 

/************************************工人代码************************************ **************/

if (taskid!=MASTER) 
{ 
    /* initialization*/ 
    for (int i=0;i<2;i++) 
    { 
     for (int j=0;j<NX;j++) 
     { 
      for(int k=0;k<NY;k++) 
      { 
       u[i][j][k]=i+j+k; 
       //cout<<u[i][j][k]<<" "; 
      } 
      //cout<<"\n"; 
     } 
     //cout<<"\n\n"; 
    } 
    source=MASTER; 
    msgtype=BEGIN; 
    MPI_Recv(&u[0][offset][0],rows*NY,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
    MPI_Recv(&offset,1,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
    MPI_Recv(&rows,1,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
    MPI_Recv(&left,1,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
    MPI_Recv(&right,1,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
    int iz=0; 

    for (int t=1;t<=TIMESTEP;t++) 
    { 
     start=offset; 
     if (offset==0) 
     { 
      start=1; 
     } 
     end=offset+rows-1; 
     if (end==NY-1) 
     { 
      end=NY-2; 
     } 
     if (left!=NONE) 
     { 
      MPI_Send(&u[iz][offset][0],NY,MPI_FLOAT,left,RTAG,MPI_COMM_WORLD); 
      source=left; 
      msgtype=LTAG; 
      MPI_Recv(&u[iz][offset-1][0],NY,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
     } 
     if (right!=NONE) 
     { 
      MPI_Send(&u[iz][offset+rows-1][0],NY,MPI_FLOAT,right,LTAG,MPI_COMM_WORLD); 
      source=right; 
      msgtype=RTAG; 
      MPI_Recv(&u[iz][offset+rows][0],NY,MPI_FLOAT,source,msgtype,MPI_COMM_WORLD,&status); 
     } 
     iz=1-iz; 
     calculate(start,end,NY,&u[iz][offset][0],&u[1-iz][offset][0]); 
    } 
    MPI_Send(&u[iz][offset][0],rows*NY,MPI_FLOAT,MASTER,DONE,MPI_COMM_WORLD); 
    MPI_Send(&offset, 1, MPI_INT, MASTER, DONE, MPI_COMM_WORLD); 
    MPI_Send(&rows, 1, MPI_INT, MASTER, DONE, MPI_COMM_WORLD); 
    MPI_Finalize; 
} 
return 0; 

}

/* calculate function*/ 
void calculate(int start,int end,int NY,float *u1, float *u2) 
{ 
for(int i=start;i<=end;i++) 
{ 
    for (int j=0;j<NY;j++) 
    { 
     *(u2+i*NY+j)=*(u1+i*NY+j)+*(u1+(i+1)*NY+j); 
    } 
} 

}

代码编译细使用

mpicc mpi_matrix.c 

但使用

mpirun -n 4 mpi_matrix 

不运行提供了错误

mpirun was unable to find the specified executable file, and therefore 
did not launch the job. This error was first reported for process 
rank 0; it may have occurred for other processes as well. 

任何人都可以帮助吗？ 谢谢

编辑1：尝试@mko说什么之后。该程序编译并运行，但有一个错误。以下是错误

[shekhar-HP-Pavilion-Notebook:3296] *** An error occurred in MPI_Recv 
[shekhar-HP-Pavilion-Notebook:3296] *** reported by process [3817603073,3] 
[shekhar-HP-Pavilion-Notebook:3296] *** on communicator MPI_COMM_WORLD 
[shekhar-HP-Pavilion-Notebook:3296] *** MPI_ERR_TRUNCATE: message truncated 
[shekhar-HP-Pavilion-Notebook:3296] *** MPI_ERRORS_ARE_FATAL (processes in this communicator will now abort, 
[shekhar-HP-Pavilion-Notebook:3296] *** and potentially your MPI job) 

Answer 1

如果用mpicc编译你需要指定-o提供目标文件名。

试试这个

mpicc -o mpi_matrix mpi_matrix.c 
mpirun -np 4 ./mpi_matrix 

或者象这样开始

mpirun -np 4 ./a.out 

默认情况下，如果你不提供-o，mpicc将创建的a.out可执行文件。