我看了一些关于SO的其他问题,但似乎都没有解决类似的问题。
我有一个函数,它对数组进行排序(使用堆排序)并计算中值。堆排序例程直接取自Numerical Recipies。
我是calloc
和free
,它们是中值函数内的一个数组,但free
似乎没有在时间上释放空间。这里有一些代码来说明我的意思:
int calcMedian(int n1, int n2, int *dat)
{
int ii, npt;
int *inparr, retval;
npt = n2 - n1 + 1; /* Number of elements in array */
inparr = calloc(npt+1, sizeof(*inparr));
for(ii = n1; ii <= n2; ii++)
inparr[ii-n1+1] = dat[ii]; /* ii-n1+1 because heapsort function likes arrays to
start from 1 */
heapsortInt(npt, inparr); /* The error isn't here, function has been previously
debugged. Sorting is in-place.*/
if (npt % 2)
retval = inparr[(npt+1)/2];
else
retval = (inparr[npt/2]+inparr[npt/2+1])/2;
free(inparr);
return(retval);
}
函数heapsortInt
经过了相当彻底的调试,并在其他几个地方使用过,没有出现任何问题。现在我在循环中调用函数calcMedian
,如下所示:
for(ii = 0; ii < maxval; ii++) {
index = ii * maxpt;
med1 = calcMedian(index, index+npt1[ii]-1, data1+index);
med2 = calcMedian(index, index+npt2[ii]-1, data2+index);
}
其中相关变量定义如下:
int *data1, *data2;
int *npt1, *npt2;
data1 = calloc(maxval * maxpt, sizeof(*data1));
data2 = calloc(maxval * maxpt, sizeof(*data2));
npt1 = calloc(maxval, sizeof(*npt1));
npt2 = calloc(maxval, sizeof(*npt2));
所以我基本上是将一个大数组的不同部分传递到calcMedian
中,并返回必要的中值。
问题:calcMedian
在调用第二个函数时似乎崩溃了。我在valgrind上运行了它,它告诉我:
Invalid read of size 4
at 0x43F67E: calcMedian /* Line no. pointing to calloc in calcMedian */
by 0x4416C9: main /* Line no pointing to second call of calcMedian */
Address 0x128ffdc0 is 6,128 bytes inside a block of size 110,788 free'd
at 0x4A063F0: free
by 0x43F728: calcMedian /* Line no. pointing to free in calcMedian */
by 0x4416C9: main /* Line no pointing to first call of calcMedian */
这是free
的问题吗?我是不是free
和calloc
太频繁了?我不知道从哪里开始调试。任何帮助都将是美妙的!
免责声明:带有实际代码的计算机无法访问互联网。我在这里尽可能准确地复制了导致问题的代码。如果有任何丢失的分号等,那是我的错,它肯定不在原始代码中。
编辑:修复了一些转录错误。我会尽快得到原始代码,但从我现在看,这似乎很好。
问题在于对calcMedian的调用。
您要添加索引两次,一次是在调用中,然后是在calcMedian中。
应该是这样的:
med1 = calcMedian(index, index+npt1[ii]-1, data1);
med2 = calcMedian(index, index+npt2[ii]-1, data2);