Trisurf Monte Carlo simulator
Samo Penic
2020-04-13 a7fb2b94260897ec914f0bbf31698e166263036c
src/cell.c
@@ -6,7 +6,7 @@
ts_cell_list  *init_cell_list(ts_uint ncmax1, ts_uint ncmax2, ts_uint ncmax3, ts_double stepsize){
    ts_uint i;
    ts_uint nocells=ncmax1*ncmax2*ncmax3;
    ts_cell_list *clist=(ts_cell_list *)malloc(sizeof(ts_cell_list));
    ts_cell_list *clist=(ts_cell_list *)calloc(1,sizeof(ts_cell_list));
    if(clist==NULL) fatal("Error while allocating memory for cell list!",100);
    clist->ncmax[0]=ncmax1;
@@ -15,8 +15,8 @@
    clist->cellno=nocells;
    clist->dcell=1.0/(1.0 + stepsize);
    clist->shift=(ts_double) clist->ncmax[0]/2;
    clist->cell=(ts_cell **)malloc(nocells*sizeof(ts_cell *));
    clist->max_occupancy=16; /* hard coded max occupancy? */
    clist->cell=(ts_cell **)calloc(nocells,sizeof(ts_cell *));
    if(clist->cell==NULL) fatal("Error while allocating memory for cell list! ncmax too large?",101);
    for(i=0;i<nocells;i++){
@@ -142,6 +142,7 @@
                cell_occupation=clist->cell[neigh_cidx]->nvertex;
          //      fprintf(stderr, "cell_occupation=%i\n",cell_occupation);
                if(cell_occupation>clist->max_occupancy){
         ts_fprintf(stderr,"max occupancy= %d, cell occupation= %d", clist->max_occupancy, cell_occupation);
                    fatal("Neighbouring cell occupation more than set max_occupancy value.",2000);
                }
// Now we check whether we didn't come close to some other vertices in the same