我正在编写一个C代码来解决欧拉方程。我的代码在群集上运行良好,但在我的 PC 上则不然。似乎是 malloc() 的问题。它无法分配请求的内存并失败。
如何使其工作?它与碎片整理有关吗?但系统设置显示(0% 碎片整理)。
只需在此处包含 malloc() 代码的一部分。
double **u, **rho_u, **rho,
int Size = 1000;
u = (double**)malloc(Size*sizeof(double*));
for(i=0;i<=Size;i++)
u[i] = (double*)malloc(Size*sizeof(double));
rho_u = (double**)malloc(Size*sizeof(double*));
for(i=0;i<=Size;i++)
rho_u[i] = (double*)malloc(Size*sizeof(double));
你可能在这里破坏了你的堆:
for(i=0;i<=Size;i++)
u[i] = (double*)malloc(Size*sizeof(double));
您分配了 1001 个指针,但只分配了 1000 个指针。正确版本:
for(i=0;i<Size;i++)
u[i] = (double*)malloc(Size*sizeof(double));
第二个循环相同。
仔细阅读malloc
的文档。它可能会失败,当它失败时,malloc
返回NULL
(失败原因由您经常使用 perror
显示的errno
给出)。
所以你应该测试malloc
失败。典型的代码至少是:
u = (double**)malloc(Size*sizeof(double*));
if (u==NULL) { perror ("malloc u"); exit(EXIT_FAILURE); };
同样,为了你的rho_u
和每一个rho_u[i]
某些操作系统可能会提供内存过载。这是我不喜欢的功能。
考虑完全初始化每个内存区域。在有效内存区域(或有效地址)之外使用内存是未定义的行为(您的程序有一个,Ctx 的答案注意到了这一点)。害怕。
我还建议使用valgrind。这是一个非常方便的工具,可以寻找与内存相关的错误,并且可以检测到您的错误。
观察:
- 避免类型转换
malloc()
,阅读 我是否转换 malloc 的结果? - 检查
malloc()
的返回值并进行正确的错误处理。 - 将循环条件从
i<=Size
更改为i<Size
,因为它会导致缓冲区溢出,因为早期的内存仅分配给Size
行而不是size+1
行。
试试这个版本:
int Size = 1000;
double **u = malloc(Size * sizeof(*u)); /* typecasting is not needed */
if(u == NULL) {
/* @TODO error handling */
}
for(i=0;i<Size;i++) { /* loop should rotate only 1000 times not 1001 times */
u[i] = malloc(Size * sizeof(**u));
if(u[i] == NULL) {
/* @TODO error handling */
}
}
同样适用于rho_u
和rho
.