Trisurf Monte Carlo simulator
Samo Penic
2013-11-30 bc4e1ee61725fdf7ec81a9583b7fa545f847e07b
commit | author | age
d7639a 1 #include<stdlib.h>
SP 2 #include "general.h"
3 #include "cell.h"
f74313 4 #include "frame.h"
d7639a 5 ts_bool centermass(ts_vesicle *vesicle){
f74313 6     ts_uint i, n=vesicle->vlist->n;
SP 7     ts_vertex **vtx=vesicle->vlist->vtx;
d7639a 8     vesicle->cm[0]=0;
SP 9     vesicle->cm[1]=0;
10     vesicle->cm[2]=0;
f74313 11     for(i=0;i<n;i++){
8f6a69 12         vesicle->cm[0]+=vtx[i]->x;
SP 13         vesicle->cm[1]+=vtx[i]->y;
14         vesicle->cm[2]+=vtx[i]->z; 
d7639a 15     } 
dac2e5 16     vesicle->cm[0]/=(ts_float)n;
SP 17     vesicle->cm[1]/=(ts_float)n;
18     vesicle->cm[2]/=(ts_float)n;
919555 19
SP 20     for(i=0;i<n;i++){
21         vtx[i]->x-=vesicle->cm[0];
22         vtx[i]->y-=vesicle->cm[1];
23         vtx[i]->z-=vesicle->cm[2]; 
24     } 
25
a63f17 26     vesicle->cm[0]=0;
SP 27     vesicle->cm[1]=0;
28     vesicle->cm[2]=0;
29
d7639a 30     return TS_SUCCESS;
SP 31 }
32
f74313 33 ts_bool cell_occupation(ts_vesicle *vesicle){
SP 34     ts_uint i,cellidx, n=vesicle->vlist->n;
d7639a 35     ts_double shift;
SP 36     ts_double dcell;
f74313 37     shift=(ts_double) vesicle->clist->ncmax[0]/2;
d7639a 38     dcell=1.0/(1.0 + vesicle->stepsize);
34d3de 39     //`fprintf(stderr, "Bil sem tu\n"); 
d7639a 40
f74313 41     cell_list_cell_occupation_clear(vesicle->clist);
SP 42     for(i=0;i<n;i++){
43     cellidx=vertex_self_avoidance(vesicle, vesicle->vlist->vtx[i]);
a63f17 44 //    already done in cell_add_vertex
SP 45 //    vesicle->vlist->vtx[i]->cell=vesicle->clist->cell[cellidx];
d7639a 46
f74313 47     cell_add_vertex(vesicle->clist->cell[cellidx],vesicle->vlist->vtx[i]);
d7639a 48
SP 49     }
dac2e5 50     //fprintf(stderr, "Bil sem tu\n"); 
41a035 51     if(dcell);
SP 52     if(shift);
d7639a 53     return TS_SUCCESS;
SP 54 }