Question

：）

关于mpi程序我遇到了一个令人不安的问题。这个想法是：每个进程（从属）都将数据发送给master，以便计算mandelbrot分形。

首先，每个奴隶发送点和它的工作人员。然后他们发了一条线就行了！

但是现在，我试图让他们发送一个行块（让我们假设有5行，所以是一个子矩阵）。

我的想法是将这五行分成一行。主人收到第一个＆＃34;新＆＃34;但是对其他人来说并不是O_o。我很不安。

我收到其他人（＆gt; 1）：信号分割故障信号代码：地址未映射在地址上失败

请帮帮我！因为很长一段时间，我一直在寻找：（

Ps：我是法国人（这就是为什么我的英语不好）



  //the whole table to be used in a master

  //int table[NX*NY];
  //int count =0;

  if (rank == 0) {

    int res;
    int line[MAXY+MAXY+1];
    int block[5*(MAXY+MAXY+1)];
    int count = 0;    

    /* Begin User Program  - the master */

   //MPI_Recv(&line, MAXY+MAXY+1, MPI_INT,MPI_ANY_SOURCE, DATATAG, MPI_COMM_WORLD, &status);
   MPI_Recv(&block, 5*(MAXY+MAXY+1), MPI_INT,MPI_ANY_SOURCE, DATATAG, MPI_COMM_WORLD, &status);
   printf("sizeof of datablock received is = %d \n",sizeof(block)/sizeof(block[0]));
   recvd = status.MPI_SOURCE;
   printf("i have received blockdata from %d \n",recvd);
   /* remplissage du case */

   for(i = -MAXX; i <= MAXX; i++) { 
      for(j = -MAXY; j <= MAXY; j++) { 
        cases[i + MAXX][j + MAXY] = block[count%(MAXY+MAXY+1)];
        //printf("j'ai fait un bloc[count], pas credible\n"); 
        count++;
        }
    } 

    dump_ppm("mandel.ppm", cases);
    printf("Fini.\n");
  }

  else {

    /* On est l'un des fils */
    /* for the block;let's suppose each son send 5 rows*/ 
    double x, y;
    int i, j, res, rc, rank,count;
    //int line[MAXY + MAXY + 1];
    int block[5*(MAXY+MAXY+1)];
    count = 0;
    MPI_Comm_rank(MPI_COMM_WORLD, &rank);
    for(i = -MAXX; i <= MAXX; i++) {
      for(j = -MAXY; j <= MAXY; j++) {
        x = 2 * i / (double)MAXX;
        y = 1.5 * j / (double)MAXY;
        res = mandel(x, y);
        //line[j+MAXY] = res;
        block[count] =res;
        if (count % (5*(MAXY+MAXY+1)) == 0){ //we send each five rows
            MPI_Send(&block,5*(MAXY+MAXY+1), MPI_INT, 0, DATATAG, MPI_COMM_WORLD);
            printf("me slave %d, have sent datablock to master\n",rank);
            printf("sizeof of datablock sent is = %d\n",sizeof(block)/sizeof(block[0]));
            }
        count++;
        }
    //MPI_Send(&line, MAXY+MAXY+1 , MPI_INT, 0, DATATAG, MPI_COMM_WORLD);
    }
  }

  MPI_Finalize();
  return 0;
}

Answer 1

函数MPI_Recv()需要接收数据的缓冲区的地址。 MPI_Send()也是如此。由于int block[5*(MAXY+MAXY+1)]是一个数组，block指向数组block[0]的第一项：这是所需的地址。另一方面，&block指向block：它类似于指向int指针的指针。但&block的值不是数组第一项的地址！

因此，你能尝试一下：

int block[5*(MAXY+MAXY+1)]
...
MPI_Send(block,5*(MAXY+MAXY+1), MPI_INT, 0, DATATAG, MPI_COMM_WORLD);
...
MPI_Recv(block, 5*(MAXY+MAXY+1), MPI_INT,MPI_ANY_SOURCE, DATATAG, MPI_COMM_WORLD, &status);

相当于：

int block[5*(MAXY+MAXY+1)]
...
MPI_Send(&block[0],5*(MAXY+MAXY+1), MPI_INT, 0, DATATAG, MPI_COMM_WORLD);
...
MPI_Recv(&block[0], 5*(MAXY+MAXY+1), MPI_INT,MPI_ANY_SOURCE, DATATAG, MPI_COMM_WORLD, &status);

如果要发送单个整数int a怎么办？ a（&a）的地址可以提供给MPI_Send()，正如专门针对MPI_Send()的许多示例中所执行的那样：

int a=42;
MPI_Send(&a,1, MPI_INT, 0, DATATAG, MPI_COMM_WORLD);

最后，请确保MPI_Send()的调用次数与MPI_Recv()的次数相同。实际上，在您发布的代码中，MPI_Recv()仅由根进程调用一次，而每个非根进程将向根发送消息。因此，该程序将适用于2个进程，如果使用更多进程或使用单个进程，则可能会失败。

尝试使用Mpi recv时出错：信号分段故障

1 个答案: