当我在程序中添加此函数时,会得到此错误:"分割故障"。我发现它可能与指针有关。
这是我的功能:
void deBoor(double* D, double* t, double x, int k, int n, double* R){
int i,j,l;
double ret[3];
double* P=(double*)malloc(n*sizeof(double));
for(i=0; i<n; ++i) P[i]=D[i];
j=k;
while(x<t[j]) ++j;
while(x>=t[j+1]) ++j;
for(l=1; l<=k; ++l){
for(i=j; i>=j-k+l; --i)
{
P[i]=(x-t[i])/(t[i+k+1-l]-t[i])*P[i]+(t[i+k+1-l]-x)/(t[i+k+1-l]-t[i])*P[i-1];
}
}
R[0]=P[j];
free(P);
}
这是功能的呼吁:
deBoor (kontrolne_tocke, prosirena_particija, xx, k, duljina, R);
代码相当不安全。有些问题,不仅与安全有关,而且其中一些肯定会导致您的崩溃:
- 输入参数应该是
const
指数,为了清楚起见。 - 未检查
malloc()
的返回值,因此,如果分配失败,则您是敬酒。 - 不要在c。 中施放
- 不要使用循环复制,请使用
memcpy()
。 - 无法保证
j
不会超出界限。 -
ret
数组从未使用过。
malloc()
的返回值