我正在编写一个矩阵乘法程序,为了节省代码空间,我创建了一个函数,该函数make
接受一个指向双精度指针并动态分配给定大小的二维数组。数组a
和b
被初始化为 [-2, 2] 中的随机值。但是,当我尝试在 gcc 上运行代码时,出现了 seg 错误。
我通过gdb运行代码,当我尝试设置b[0][0]
为随机数时出现错误。当我尝试b[0][0]
在 gdb 中打印时,出现错误:
无法访问地址 0xbfebfdc3f5d77f80 处的内存
但是,我实际上可以b[0][0]
在此之前访问。我可以在分配后打印数组而不会出错。出于某种原因,似乎总是数组b
引起了问题。
我感觉这是一个指针问题,但我看不到在哪里。我不是一个没有经验的程序员,但我花了 2 天的时间试图找出为什么这个错误不断出现。任何人都可以解释一下吗?
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
void make(double ***a, int m, int n)
{
int i;
double *tmp;
(*a) = (double **) malloc(m*sizeof(double *));
tmp = (double *) calloc(m*n,sizeof(double));
if(tmp == NULL)
{
printf("Error with malloc.\n");
exit(EXIT_FAILURE);
}
for(i=0;i<m;i++)
(*a)[i] = &tmp[i*n];
free(tmp);
if(a == NULL)
{
fprintf(stderr, "Error with the matrix, dimensions: %d, %d. \n", m, n);
exit(EXIT_FAILURE);
}
}
int main()
{
int i, j;
int l, m, n;
double **a, **b, **c;
l = m = n = 8;
srand48(time(NULL));
make(&a, l, m);
make(&b, m, n);
make(&c, l, n);
for(i=0; i<l; i++)
for(j=0; j<m; j++)
a[i][j] = (drand48()*4.0) - 2.0;
for(i=0; i<m; i++)
for(j=0; j<n; j++)
b[i][j] = (drand48()*4.0) - 2.0;
free(a);
free(b);
free(c);
return 0;
}