From 1121faa13a2038facad22073f0fc610903d98691 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo@CAE-linux.(none)>
Date: Fri, 05 Sep 2014 20:18:05 +0000
Subject: [PATCH] First variant of constant volume constrant (a new one proposed by Miha after reading his article). It seems to work, however there are still some things to be done, such as Miha's derivation of the epsvol (0.003% is used at the moment) and solving the problem of additional global variables.
---
src/timestep.c | 2 ++
1 files changed, 2 insertions(+), 0 deletions(-)
diff --git a/src/timestep.c b/src/timestep.c
index f2fb478..26f42ce 100644
--- a/src/timestep.c
+++ b/src/timestep.c
@@ -45,6 +45,8 @@
centermass(vesicle);
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;
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