From 958e0e46e248cc7fd408fc2c4f7eb119ed93b1e5 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@fe.uni-lj.si>
Date: Mon, 28 Apr 2014 13:40:04 +0000
Subject: [PATCH] Only reaction on successful acceptance of MC step needs to be programmed.
---
src/vertexmove.c | 6 +++---
1 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/src/vertexmove.c b/src/vertexmove.c
index 44ef63e..a27a776 100644
--- a/src/vertexmove.c
+++ b/src/vertexmove.c
@@ -23,7 +23,7 @@
ts_double delta_energy,oenergy,dvol=0.0;
ts_double costheta,sintheta,phi,r;
//This will hold all the information of vtx and its neighbours
- ts_vertex backupvtx[20], *constvol_vtx_moved, *constvol_vtx_backup;
+ ts_vertex backupvtx[20], *constvol_vtx_moved=NULL, *constvol_vtx_backup=NULL;
memcpy((void *)&backupvtx[0],(void *)vtx,sizeof(ts_vertex));
//Some stupid tests for debugging cell occupation!
@@ -157,7 +157,7 @@
for(i=0;i<vtx->tristar_no;i++) triangle_normal_vector(vtx->tristar[i]);
if(vesicle->tape->constvolswitch == 1){
- ts_bool constvolumerestore(constvol_vtx_backup);
+ constvolumerestore(constvol_vtx_moved,constvol_vtx_backup);
}
return TS_FAIL;
@@ -173,7 +173,7 @@
}
if(vesicle->tape->constvolswitch == 1){
- ts_bool constvolumeaccept(constvol_vtx_backup);
+ constvolumeaccept(constvol_vtx_moved,constvol_vtx_backup);
}
// if(oldcellidx);
//END MONTE CARLOOOOOOO
--
Gitblit v1.9.3