Trisurf Monte Carlo simulator
Samo Penic
2012-06-14 30ee9c409ed7c743c693f1e068c1d7393949207a
commit | author | age
d7639a 1 #include<stdlib.h>
SP 2 #include<stdio.h>
aec47d 3 #include<math.h>
SP 4 //#include "io.h"
5 #include "general.h"
6 #include "timestep.h"
7 #include "vertexmove.h"
30ee9c 8 #include "bondflip.h"
d7639a 9
SP 10 ts_bool single_timestep(ts_vesicle *vesicle){
11     ts_bool retval;
12     ts_double rnvec[3];
13     ts_uint i;
aec47d 14     for(i=0;i<vesicle->vlist->n;i++){
d7639a 15         rnvec[0]=drand48();
SP 16         rnvec[1]=drand48();
17         rnvec[2]=drand48();
aec47d 18         retval=single_verticle_timestep(vesicle,vesicle->vlist->vtx[i],rnvec);
d7639a 19     }
SP 20
aec47d 21     for(i=0;i<vesicle->blist->n;i++){
d7639a 22         rnvec[0]=drand48();
SP 23         rnvec[1]=drand48();
24         rnvec[2]=drand48();
25         //find a bond and return a pointer to a bond...
26         //call single_bondflip_timestep...
30ee9c 27         retval=single_bondflip_timestep(vesicle,vesicle->blist->bond[i],rnvec);
d7639a 28         
SP 29     } 
41a035 30     if(retval);
d7639a 31     return TS_SUCCESS;
SP 32 }
33
34
35