From 3de2893d882d2c407375146376fbf9bb54f33e8b Mon Sep 17 00:00:00 2001 From: Samo Penic <samo.penic@fe.uni-lj.si> Date: Wed, 30 Apr 2014 13:44:12 +0000 Subject: [PATCH] Constant volume works. CAVEAT: in centermass I had to recalculate all normals of triangles again! This could influence behaviour of the system. --- src/timestep.c | 10 ++++++++++ 1 files changed, 10 insertions(+), 0 deletions(-) diff --git a/src/timestep.c b/src/timestep.c index 7fd321c..a57687a 100644 --- a/src/timestep.c +++ b/src/timestep.c @@ -32,11 +32,16 @@ for(i=start_iteration;i<inititer+iterations;i++){ vmsr=0.0; bfsr=0.0; +/* vesicle_volume(vesicle); + fprintf(stderr,"Volume before TS=%1.16e\n", vesicle->volume); */ for(j=0;j<mcsweeps;j++){ single_timestep(vesicle, &vmsrt, &bfsrt); vmsr+=vmsrt; bfsr+=bfsrt; } +/* + vesicle_volume(vesicle); + fprintf(stderr,"Volume after TS=%1.16e\n", vesicle->volume); */ vmsr/=(ts_double)mcsweeps; bfsr/=(ts_double)mcsweeps; centermass(vesicle); @@ -87,6 +92,8 @@ } ts_bool single_timestep(ts_vesicle *vesicle,ts_double *vmsr, ts_double *bfsr){ +// vesicle_volume(vesicle); +// fprintf(stderr,"Volume before TS=%1.16e\n", vesicle->volume); ts_bool retval; ts_double rnvec[3]; ts_uint i,j, b; @@ -105,6 +112,7 @@ //find a bond and return a pointer to a bond... //call single_bondflip_timestep... retval=single_bondflip_timestep(vesicle,vesicle->blist->bond[b],rnvec); + // b++; retval=TS_FAIL; if(retval==TS_SUCCESS) bfsrcnt++; } @@ -131,6 +139,8 @@ // printf("Bondflip success rate in one sweep: %d/%d=%e\n", cnt,3*vesicle->blist->n,(double)cnt/(double)vesicle->blist->n/3.0); *vmsr=(ts_double)vmsrcnt/(ts_double)vesicle->vlist->n; *bfsr=(ts_double)bfsrcnt/(ts_double)vesicle->vlist->n/3.0; +// vesicle_volume(vesicle); +// fprintf(stderr,"Volume after TS=%1.16e\n", vesicle->volume); return TS_SUCCESS; } -- Gitblit v1.9.3