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

标签 c segmentation-fault mpi

对于我的并行计算类(class)中的一个项目,我需要实现一个并行版本的生命游戏。

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

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

enter image description here

我们需要清理函数代码,使其符合我们的 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]);

关于c - 生命游戏实现中的段错误问题,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28934147/

相关文章:

c - 由于C文件中的一行代码而发生段错误,并且整个程序无法运行

c - 段错误解决方案

c++ - 如何使用 MPI 收集行和列分区矩阵?

c++ - 带有特殊字符的字符串的控制台输出

c - 返回一个 C 数组并在 main.c 中读取它

java - 在 JNI 下以 native 代码获取 java 缓冲区的最有效方法是什么?

C编程: Error in writing values of a particular column from a file to another file

c - 两个 MPI_Allreduce() 函数不起作用,给出 NULL Communicator 错误

fortran - 如何在保持兼容性的同时从 mpif.h 切换到 fortran 中的 mpi_f08?

c++ - 如何从 bmp 文件制作动画 png?