Введение в практику разработки параллельных программ в стандарте MPI. Баканов В.М - 55 стр.

UptoLike

ВУЗ: 

МИРЭА | Москва

Составители: 

Рубрика: 

- 55 -
printf("Here is the result matrix:\n");
for (i=0; i<NRA; i++)
{
printf("\n");
for (j=0; j<NCB; j++)
printf("%6.2f ", c[i][j]);
}
printf ("\n");
}
/**************************** worker task ************************************/
if (taskid > MASTER)
{
MPI_Recv(&offset, 1, MPI_INT, MASTER, FROM_MASTER, M_C_W, &status);
MPI_Recv(&rows, 1, MPI_INT, MASTER, FROM_MASTER, M_C_W, &status);
MPI_Recv(&a, rows*NCA, MPI_DOUBLE, MASTER, FROM_MASTER,
M_C_W, status);
MPI_Recv(&b, NCA*NCB, MPI_DOUBLE, MASTER, FROM_MASTER,
M_C_W, &status);
for (k=0; k<NCB; k++)
for (i=0; i<rows; i++)
{
c[i][k] = 0.0;
for (j=0; j<NCA; j++)
c[i][k] += a[i][j] * b[j][k];
}
MPI_Send(&a_offset, 1, MPI_INT, MASTER, FROM_WORKER, M_C_W);
MPI_Send(&a_rows, 1, MPI_INT, MASTER, FROM_WORKER, M_C_W);
MPI_Send(&c, a_rows*NCB, MPI_DOUBLE, MASTER, FROM_WORKER,
M_C_W);
}
MPI_Finalize();
} // end of MM_MPI_0.C program
Cтроки матрицы [A] распределяются между numworkers процессов по av-
rows штук; причем каждому процессу dest=1
÷
numworkers посредством вызова
функции MPI_Send(&a[offset][0], rows*NCA, MPI_DOUBLE, dest, FROM_MASTER,
MPI_COMM_WORLD) пересылаются rows строк (rows=averow, если NRA делится
на numworkers нацело и rows=averow+1 в противном случае, offset – номер
первой пересылаемой строки). Матрица [B] пересылается посредством
MPI_Send(&b, NCA*NCB, MPI_DOUBLE, dest, FROM_MASTER,
MPI_COMM_WORLD) всем numworkers процессам целиком.
После получения процессами данных выполняется умножение посредст-
вом внешнего цикла по k=1 ÷ NCA (т.е. по столбцам [A] и строкам [B] и внут-
реннего цикла по i=1 ÷rows (т.е. переданной части строк [A] и [C]):
for (k=0; k<NCB; k++)
for (i=0; i<rows; i++) // …на узле присутствует только rows строк матрицы [A] 
     printf("Here is the result matrix:\n");
     for (i=0; i MASTER)
  {
     MPI_Recv(&offset, 1, MPI_INT, MASTER, FROM_MASTER, M_C_W, &status);
     MPI_Recv(&rows, 1, MPI_INT, MASTER, FROM_MASTER, M_C_W, &status);
     MPI_Recv(&a, rows*NCA, MPI_DOUBLE, MASTER, FROM_MASTER,
                                                               M_C_W, status);
     MPI_Recv(&b, NCA*NCB, MPI_DOUBLE, MASTER, FROM_MASTER,
                                                               M_C_W, &status);

     for (k=0; k

Страницы