或者,复制" Facing an error - glibc detected free invalid next size (fast).
我有一个叫做bond的结构体,定义如下:
typedef struct{
int type1;
int type2;
int id_1;
int id_2;
float dist;
} bond;
我在嵌套循环中分配这些结构体的数组,并定期释放它。然而,由于某种原因,我得到一个免费的()无效的下一个大小错误。代码如下:
while(i<(C.num_type_A+C.num_type_B)){//1a
bond_arr=(bond*)malloc(100*sizeof(bond));
a=-N1;
while(a<=N1){//2a
b=-N2;
while(b<=N2){//3a
c=-N3;
while(c<=N3){//4a
j=0;
while(j<(C.num_type_A+C.num_type_B)){//5a
//if the same atom do nothing
if((i==j)&&(a==0)&&(b==0)&&(c==0)){//6a
j=j+1;
}//6b
else{//7a
//calculate bond length
bondlength=calc_dist(a,b,c,i,j,C);
if(bondlength<=cutoff){//8a
//store bond
temp_bond.type1=((i+1)<=C.num_type_A);
temp_bond.type2=((j+1)<=C.num_type_B);
temp_bond.id_1=i;
temp_bond.id_2=j;
temp_bond.dist=bondlength;
bond_arr[n]=temp_bond;
n=n+1;
//if out of memory allocate twice as much
if(n==(nmax-1)){//9a
printf("beginn");
temp_ptr=realloc(bond_arr,sizeof(bond)*2*nmax);
printf("endn");
if(temp_ptr==NULL){
printf("Memory allocation failedn");
}
else{
bond_arr=temp_ptr;
}
nmax=2*nmax;
}//9b
}//8b
}//7b
j=j+1;
}//5b
c=c+1;
}//4b
b=b+1;
}//3b
a=a+1;
}//2b
//sort bonds and update LI index
sort_bonds(bond_arr,n);
f=update_LI(bond_arr,LIcurr,n);
LIcurr=f;
printf("%dn",n);
free(bond_arr);
n=0;
i=i+1;
}//1b
看起来realloc逻辑可能是罪魁祸首。初始分配大小为100。然后在n
达到nmax-1
时生长。nmax
的初始值没有显示(我看到)。即使它从100开始,它也不会在循环的顶部重置。因此,如果n
增长超过100并引起realloc,则nmax将加倍,并且不再匹配原始大小100。