使用free时下一个大小无效



或者,复制" 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。

相关内容

  • 没有找到相关文章

最新更新