Trisurf Monte Carlo simulator
Samo Penic
2012-06-07 5bb11d36948f9ccc243711d4e9caf4b19b3c384b
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"
d7639a 8
SP 9 ts_bool single_timestep(ts_vesicle *vesicle){
10     ts_bool retval;
11     ts_double rnvec[3];
12     ts_uint i;
aec47d 13     for(i=0;i<vesicle->vlist->n;i++){
d7639a 14         rnvec[0]=drand48();
SP 15         rnvec[1]=drand48();
16         rnvec[2]=drand48();
aec47d 17         retval=single_verticle_timestep(vesicle,vesicle->vlist->vtx[i],rnvec);
d7639a 18     }
SP 19
aec47d 20     for(i=0;i<vesicle->blist->n;i++){
d7639a 21         rnvec[0]=drand48();
SP 22         rnvec[1]=drand48();
23         rnvec[2]=drand48();
24         //find a bond and return a pointer to a bond...
25         //call single_bondflip_timestep...
aec47d 26 //        retval=single_bondflip_timestep(vesicle,&vesicle->blist.bond[i],rnvec);
d7639a 27         
SP 28     } 
41a035 29     if(retval);
d7639a 30     return TS_SUCCESS;
SP 31 }
32
33
34