我有使用OpenMPI.Since我有点新的使用的openmpi的概念编译此代码的问题,这将是巨大的,如果你有人可以给我一个提示的错误在这里。 编译的作品就好了,但如果我运行代码,我得到这个消息:
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.
NOTE: A common cause for this error is misspelling a mpirun command
line parameter option (remember that mpirun interprets the first
unrecognized command line token as the executable).
我使用的编译:
mpic++ matmult.cpp -o matmult
与运行它:
mpirun -n 2 matmult
...这里是使用的代码:
#include <stdio.h>
#include <stdlib.h>
#include <mpi.h>
#define MASTER 0
#define FROM_MASTER 1
#define FROM_WORKER 2
// ---------------------------------------------------------------------------
// allocate space for empty matrix A[row][col]
// access to matrix elements possible with:
// - A[row][col]
// - A[0][row*col]
float **alloc_mat(int row, int col)
{
float **A1, *A2;
A1 = (float **)calloc(row, sizeof(float *)); // pointer on rows
A2 = (float *)calloc(row*col, sizeof(float)); // all matrix elements
for (int i = 0; i < row; i++)
A1[i] = A2 + i*col;
return A1;
}
// ---------------------------------------------------------------------------
// random initialisation of matrix with values [0..9]
void init_mat(float **A, int row, int col)
{
for (int i = 0; i < row*col; i++)
A[0][i] = (float)(rand() % 10);
}
// ---------------------------------------------------------------------------
// DEBUG FUNCTION: printout of all matrix elements
void print_mat(float **A, int row, int col, char *tag)
{
int i, j;
printf("Matrix %s:\n", tag);
for (i = 0; i < row; i++)
{
for (j = 0; j < col; j++)
printf("%6.1f ", A[i][j]);
printf("\n");
}
}
// ---------------------------------------------------------------------------
int main(int argc, char *argv[]) {
int numtasks;
int taskid;
int numworkers;
int source;
int dest;
int mtype;
int rows;
int averow, extra, offset;
double starttime, endtime;
float **A, **B, **C; // matrices
int d1, d2, d3; // dimensions of matrices
int i, j, k, rc; // loop variables
MPI_Status status;
MPI_Init(&argc,&argv);
MPI_Comm_rank(MPI_COMM_WORLD,&taskid);
MPI_Comm_size(MPI_COMM_WORLD,&numtasks);
if (argc != 4) {
printf ("Matrix multiplication: C = A x B\n");
printf ("Usage: %s <NumRowA> <NumColA> <NumColB>\n", argv[0]);
return 0;
}
if (numtasks < 2 ) {
printf("Need at least two MPI tasks. Quitting...\n");
MPI_Abort(MPI_COMM_WORLD,rc);
exit(1);
}
/* read user input */
d1 = atoi(argv[1]); // rows of A and C d1
d2 = atoi(argv[2]); // cols of A and rows of B d2
d3 = atoi(argv[3]); // cols of B and C d3
printf("Matrix sizes C[%d][%d] = A[%d][%d] x B[%d][%d]\n", d1, d3, d1, d2, d2, d3);
/* prepare matrices */
A = alloc_mat(d1, d2);
init_mat(A, d1, d2);
B = alloc_mat(d2, d3);
init_mat(B, d2, d3);
C = alloc_mat(d1, d3);
/* Code für den Manager */
if (taskid == MASTER) {
/*printf("matrix multiplikation withMPI\n");
printf("initializing arrays ...\n");
for (i=0; i<d1; i++)
for (j=0; j<d2; j++)
A[i][j]=i+j;
for (i=0; i<d2; i++)
for (j=0; j<d3; j++)
B[i][j]=i*j;*/
/* Matrizen versenden */
averow = d1/numworkers;
extra = d1%numworkers;
offset = 0;
mtype = FROM_MASTER;
starttime=MPI_Wtime();
for (dest=1;dest<=numworkers;dest++) {
rows = (dest <= extra) ? averow+1 :averow;
printf("Sending %drows to task %doffset=%d\n",rows,dest,offset);
MPI_Send(&offset, 1, MPI_INT,dest,mtype, MPI_COMM_WORLD);
MPI_Send(&rows, 1, MPI_INT,dest,mtype, MPI_COMM_WORLD);
MPI_Send(&A[offset][0],rows*d2, MPI_DOUBLE,dest,mtype, MPI_COMM_WORLD);
MPI_Send(&B, d2*d3, MPI_DOUBLE,dest,mtype, MPI_COMM_WORLD);
offset =offset+rows;
}
/* Ergebnisse empfangen */
mtype = FROM_WORKER;
for (i=1; i<=numworkers; i++) {
source = i;
MPI_Recv(&offset, 1, MPI_INT,source,mtype, MPI_COMM_WORLD, &status);
MPI_Recv(&rows, 1, MPI_INT,source,mtype, MPI_COMM_WORLD, &status);
MPI_Recv(&C[offset][0],rows*d3,
MPI_DOUBLE,source,mtype,MPI_COMM_WORLD,&status);
printf("Received results from task %d\n",source);
}
endtime=MPI_Wtime();
printf("\nIt took %fseconds.\n",endtime-starttime);
}
/* Code für die Arbeiter */
if (taskid > MASTER) {
mtype = FROM_MASTER;
MPI_Recv(&offset, 1, MPI_INT, MASTER,mtype, MPI_COMM_WORLD, &status);
MPI_Recv(&d1, 1, MPI_INT, MASTER,mtype, MPI_COMM_WORLD, &status);
MPI_Recv(&A,rows*d2, MPI_DOUBLE, MASTER,mtype, MPI_COMM_WORLD, &status);
MPI_Recv(&B, d2*d3, MPI_DOUBLE, MASTER,mtype, MPI_COMM_WORLD, &status);
/* print user instruction */
// no initialisation of C, because it gets filled by matmult
/* serial version of matmult */
printf("Perform matrix multiplication...\n");
for (i = 0; i < d1; i++)
for (j = 0; j < d3; j++)
for (k = 0; k < d2; k++)
C[i][j] += A[i][k] * B[k][j];
mtype = FROM_WORKER;
MPI_Send(&offset, 1, MPI_INT, MASTER,mtype, MPI_COMM_WORLD);
MPI_Send(&d1, 1, MPI_INT, MASTER,mtype, MPI_COMM_WORLD);
MPI_Send(&C,rows*d3, MPI_DOUBLE, MASTER,mtype, MPI_COMM_WORLD);
}
MPI_Finalize();
/* test output
print_mat(A, d1, d2, "A");
print_mat(B, d2, d3, "B");
print_mat(C, d1, d3, "C"); */
printf ("\nDone.\n");
//return 0;
}
运行结果 mpirun matmult
(默认设置,单个进程):
的mpirun已因与PID 77202上节点juliuss-MBP-3不当离开退出到过程秩0。 有三个原因可能发生这种情况:
这个过程在退出之前没有叫“初始化”,但其他人的工作做。 这可能会导致作业无限期挂在等待所有进程调用“初始化”。 根据规则,如果一个进程调用“初始化”,则所有进程必须调用结束之前“初始化”。
这个过程被称为“初始化”,但退出而不调用“完成”。 通过规则,即所谓的“初始化”的所有进程必须调用“完成”退出之前,否则将被视为“异常终止”
这个过程被称为“MPI_Abort”或“orte_abort”和MCA参数orte_create_session_dirs设置为false。 在这种情况下,运行时无法检测中止呼叫是异常终止。 因此,您将收到的唯一的错误消息是这一个。 这可能导致在应用其他进程通过的mpirun发送的信号被终止(如这里报道的)。 您可以通过在mpirun命令行上指定-quiet避免此消息。