From c58cd2c3a13498cd200a6e708181c3427e68bd30 Mon Sep 17 00:00:00 2001
From: Samo Penic <samo.penic@gmail.com>
Date: Tue, 05 Nov 2013 18:38:11 +0000
Subject: [PATCH] Testing code without bondflip
---
src/timestep.c | 9 +++------
src/tape | 4 ++--
2 files changed, 5 insertions(+), 8 deletions(-)
diff --git a/src/tape b/src/tape
index 0351c83..695664d 100644
--- a/src/tape
+++ b/src/tape
@@ -19,9 +19,9 @@
#how many MC sweeps between subsequent records of states to disk
mcsweeps=100
#how many initial mcsweeps*inititer MC sweeps before recording to disk?
-inititer=100
+inititer=10
#how many records do you want on the disk iteration are there in a run?
-iterations=1000
+iterations=10
#shut up if we are using cluster!!!
diff --git a/src/timestep.c b/src/timestep.c
index 1ac1eb4..c25c1fb 100644
--- a/src/timestep.c
+++ b/src/timestep.c
@@ -31,7 +31,7 @@
ts_bool single_timestep(ts_vesicle *vesicle){
ts_bool retval;
ts_double rnvec[3];
- ts_uint i, b;
+ ts_uint i;// b;
for(i=0;i<vesicle->vlist->n;i++){
rnvec[0]=drand48();
rnvec[1]=drand48();
@@ -48,12 +48,8 @@
}
// ts_int cnt=0;
+/*
for(i=0;i<vesicle->vlist->n;i++){
-//why is rnvec needed in bondflip?
-/* rnvec[0]=drand48();
- rnvec[1]=drand48();
- rnvec[2]=drand48();
-*/
b=rand() % vesicle->blist->n;
//find a bond and return a pointer to a bond...
//call single_bondflip_timestep...
@@ -61,6 +57,7 @@
// if(retval==TS_SUCCESS) cnt++;
}
// printf("Bondflip success rate in one sweep: %d/%d=%e\n", cnt,vesicle->blist->n,(double)cnt/(double)vesicle->blist->n);
+*/
if(retval);
return TS_SUCCESS;
}
--
Gitblit v1.9.3