生命游戏实现中的分段错误问题

Segmentation Fault Issues in Game of Life Implementation

我的并行计算项目class,我需要实现一个并行版本的生命游戏。

我正在使用教科书作者编写的函数 "read_row_stripped_matrix"。此函数从包含矩阵中的行数、矩阵中的列数和矩阵中的数据的文件中读取输入。

该函数通过分配一个名为 "storage" 的一维数组来设置二维矩阵,该数组包含矩阵的所有数据。二维矩阵的每一行都指向它在存储中的第一个元素,如下图所示:

我们需要清理函数代码,使其符合我们的 C 风格指南。所以我清理了一些东西,这样它会更具可读性。

我 运行 现在遇到的问题是将矩阵中的每一行指向存储中的第一个元素。我在连接这些指针时遇到分段错误,特别是在函数的这一部分:

   /* Dynamically allocate matrix. Allow double subscripting
      through 'a'. */

   *storage = my_malloc (id, local_rows * *n * sizeof(int));       
   *subs    = my_malloc (id, local_rows * PTR_SIZE);

   for (i = 0; i < local_rows; i++) {
      *subs[i]=&(*storage[i * *n]);
   }

让我感到困惑的是,我很确定我已经为数组分配了足够的内存。在我正在测试的示例中,*m 和 *n 等于 5,并且 local_rows 等于 5。所以我分配 25*sizeof(int) 用于存储,这应该足以容纳 5x5 矩阵的所有元素.

这是 my_malloc 函数,它为特定处理器分配内存:

/*
 *   Function 'my_malloc' is called when a process wants
 *   to allocate some space from the heap. If the memory
 *   allocation fails, the process prints an error message
 *   and then aborts execution of the program.
 */

void* my_malloc (
   int id,     /* IN - Process rank */
   int bytes)  /* IN - Bytes to allocate */
{
   void *buffer;
   if ((buffer = malloc ((size_t) bytes)) == NULL) {
      printf ("Error: Malloc failed for process %d\n", id);
      fflush (stdout);
      MPI_Abort (MPI_COMM_WORLD, MALLOC_ERROR);
   }
   return buffer;
}

老实说,我发现指针令人困惑,如果问题很明显,请原谅我。我做这个的时间比我应该做的要长,所以我的脑子可能炸了。

如果您需要更多代码,请随时询问。

首先,你这样做:

*subs    = my_malloc (id, local_rows * PTR_SIZE);

然后,你这样做:

*subs[i]=&(*storage[i * *n]);

很确定这就是您的问题,就在那里。在我看来它真的应该是:

(*subs)[i]=&(*storage[i * *n]);