尝试使用Mpi recv时发生错误：信号分段错误

我对mpi程序有一个令人困扰的问题。这个想法是：每个进程（从）发送数据给主，以计算mandelbrot分形。

首先，每个奴隶发送点和它的工人。然后，他们发送线路，它的工作！

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

我的想法是将这五行写成一行。主人收到第一个“新”线，但不为其他人O_o。我感到不安。

我收到了别人（> 1）：信号分割故障信号代码：ADRESS没有映射在您的地址

没有，请帮助我！因为它是一个很长一段时间，i'have一直在寻找:(

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; 
}

来源

2016-10-22 yann Dan

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个进程，并且如果使用更多进程或使用单个进程，该程序可能会失败。

来源

2016-10-23 11:21:48 francis

尝试使用Mpi recv时发生错误：信号分段错误

回答

相关问题