From a549778bfece819cf0b599168c571e400e6eb67d Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Wed, 10 Sep 2014 07:03:15 +0000
Subject: [PATCH] Added calculation of epsV.
---
src/timestep.c | 3 ++-
1 files changed, 2 insertions(+), 1 deletions(-)
diff --git a/src/timestep.c b/src/timestep.c
index 26f42ce..98b47dc 100644
--- a/src/timestep.c
+++ b/src/timestep.c
@@ -46,7 +46,8 @@
cell_occupation(vesicle);
vesicle_volume(vesicle); //needed for constant volume at this moment
V0=vesicle->volume;
- epsvol=V0*0.003e-2; //TODO! Follow Miha's derivation for exact formula;
+ epsvol=4.0*sqrt(2.0*M_PI)/pow(3.0,3.0/4.0)*V0/pow(vesicle->tlist->n,3.0/2.0);
+ // fprintf(stderr, "DVol=%1.16f (%1.16f), V0=%1.16f\n", epsvol,0.003e-2*V0,V0);
if(start_iteration<inititer) ts_fprintf(stdout, "Starting simulation (first %d x %d MC sweeps will not be recorded on disk)\n", inititer, mcsweeps);
for(i=start_iteration;i<inititer+iterations;i++){
vmsr=0.0;
--
Gitblit v1.9.3