Trisurf Monte Carlo simulator
Samo Penic
2014-03-06 083e03adeaf985c9da1ed31d259ced9bb54c207c
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"
d7a113 9 #include "frame.h"
SP 10 #include "io.h"
fedf2b 11
d7a113 12 ts_bool run_simulation(ts_vesicle *vesicle, ts_uint mcsweeps, ts_uint inititer, ts_uint iterations){
SP 13     ts_uint i, j;
14
15     centermass(vesicle);
16     cell_occupation(vesicle);
17     ts_fprintf(stdout, "Starting simulation (first %d x %d MC sweeps will not be recorded on disk)\n", inititer, mcsweeps);
18     for(i=0;i<inititer+iterations;i++){
19         for(j=0;j<mcsweeps;j++){
20             single_timestep(vesicle);
21         }
22         centermass(vesicle);
23         cell_occupation(vesicle);
f8e6ba 24         ts_fprintf(stdout,"Done %d out of %d iterations (x %d MC sweeps).\n",i+1,inititer+iterations,mcsweeps);
083e03 25             dump_state(vesicle);
d7a113 26         if(i>inititer){
SP 27             write_vertex_xml_file(vesicle,i-inititer);
28         }
29     }
30     return TS_SUCCESS;
31 }
d7639a 32
SP 33 ts_bool single_timestep(ts_vesicle *vesicle){
34     ts_bool retval;
35     ts_double rnvec[3];
fedf2b 36     ts_uint i,j,b;
aec47d 37     for(i=0;i<vesicle->vlist->n;i++){
d7639a 38         rnvec[0]=drand48();
SP 39         rnvec[1]=drand48();
40         rnvec[2]=drand48();
aec47d 41         retval=single_verticle_timestep(vesicle,vesicle->vlist->vtx[i],rnvec);
d7639a 42     }
SP 43
b7c9b3 44 //    ts_int cnt=0;
fedf2b 45     for(i=0;i<3*vesicle->vlist->n;i++){
dcf17d 46 //why is rnvec needed in bondflip?
SP 47 /*        rnvec[0]=drand48();
48         rnvec[1]=drand48();
49         rnvec[2]=drand48();
50 */ 
142a67 51     b=rand() % vesicle->blist->n;
d7639a 52         //find a bond and return a pointer to a bond...
SP 53         //call single_bondflip_timestep...
142a67 54         retval=single_bondflip_timestep(vesicle,vesicle->blist->bond[b],rnvec);
b7c9b3 55 //    if(retval==TS_SUCCESS) cnt++;        
fedf2b 56     }
M 57
58     for(i=0;i<vesicle->poly_list->n;i++){
59     for(j=0;j<vesicle->poly_list->poly[i]->vlist->n;j++){
60         rnvec[0]=drand48();
61         rnvec[1]=drand48();
62         rnvec[2]=drand48();
63         retval=single_poly_vertex_move(vesicle,vesicle->poly_list->poly[i],vesicle->poly_list->poly[i]->vlist->vtx[j],rnvec);    
64     }
65
66     }
67  
68 //    printf("Bondflip success rate in one sweep: %d/%d=%e\n", cnt,3*vesicle->blist->n,(double)cnt/(double)vesicle->blist->n/3.0);
41a035 69     if(retval);
d7639a 70     return TS_SUCCESS;
SP 71 }
72
73
74