生命游戏实现中的分段错误问题
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]);
我的并行计算项目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]);