From 4a1b277c5b9ad5369441d028e53f108adc8ef62f Mon Sep 17 00:00:00 2001 From: mihaf <miha.fosnaric@gmail.com> Date: Fri, 21 Feb 2014 08:36:20 +0000 Subject: [PATCH] Nothing special. --- src/frame.c | 45 +++++++++++++++++++++++++++++++++------------ 1 files changed, 33 insertions(+), 12 deletions(-) diff --git a/src/frame.c b/src/frame.c index a4a9d32..1ee270b 100644 --- a/src/frame.c +++ b/src/frame.c @@ -9,34 +9,55 @@ vesicle->cm[1]=0; vesicle->cm[2]=0; for(i=0;i<n;i++){ - vesicle->cm[0]+=vtx[i]->data->x; - vesicle->cm[1]+=vtx[i]->data->y; - vesicle->cm[2]+=vtx[i]->data->z; + vesicle->cm[0]+=vtx[i]->x; + vesicle->cm[1]+=vtx[i]->y; + vesicle->cm[2]+=vtx[i]->z; } vesicle->cm[0]/=(ts_float)n; vesicle->cm[1]/=(ts_float)n; vesicle->cm[2]/=(ts_float)n; + + for(i=0;i<n;i++){ + vtx[i]->x-=vesicle->cm[0]; + vtx[i]->y-=vesicle->cm[1]; + vtx[i]->z-=vesicle->cm[2]; + } + + vesicle->cm[0]=0; + vesicle->cm[1]=0; + vesicle->cm[2]=0; + return TS_SUCCESS; } ts_bool cell_occupation(ts_vesicle *vesicle){ - ts_uint i,cellidx, n=vesicle->vlist->n; - ts_double shift; - ts_double dcell; - shift=(ts_double) vesicle->clist->ncmax[0]/2; - dcell=1.0/(1.0 + vesicle->stepsize); + ts_uint i,j,cellidx, n=vesicle->vlist->n; + //ts_double shift; + //ts_double dcell; + //shift=(ts_double) vesicle->clist->ncmax[0]/2; + //dcell=1.0/(1.0 + vesicle->stepsize); //`fprintf(stderr, "Bil sem tu\n"); cell_list_cell_occupation_clear(vesicle->clist); for(i=0;i<n;i++){ cellidx=vertex_self_avoidance(vesicle, vesicle->vlist->vtx[i]); - vesicle->vlist->vtx[i]->data->cell=vesicle->clist->cell[cellidx]; +// already done in cell_add_vertex +// vesicle->vlist->vtx[i]->cell=vesicle->clist->cell[cellidx]; cell_add_vertex(vesicle->clist->cell[cellidx],vesicle->vlist->vtx[i]); - } + + for(i=0;i<vesicle->poly_list->n;i++){ + for(j=0;j<vesicle->poly_list->poly[i]->vlist->n;j++){ + cellidx=vertex_self_avoidance(vesicle, vesicle->poly_list->poly[i]->vlist->vtx[j]); + cell_add_vertex(vesicle->clist->cell[cellidx],vesicle->poly_list->poly[i]->vlist->vtx[j]); + } + } + + + //fprintf(stderr, "Bil sem tu\n"); - if(dcell); - if(shift); + //if(dcell); + //if(shift); return TS_SUCCESS; } -- Gitblit v1.9.3