From: Dariusz Murakowski <murakdar@mit.edu>
Date: Tue, 21 Jul 2015 20:56:37 +0000 (-0400)
Subject: Add option to initialize state completely randomly.
X-Git-Url: http://src.murakowski.org/?a=commitdiff_plain;h=c47a14c59df60e8c4f16da8603bedb991f355877;p=VirEvoDyn.git

Add option to initialize state completely randomly.
---

diff --git a/mainQEE.cpp b/mainQEE.cpp
index 3cd263b..04676cd 100755
--- a/mainQEE.cpp
+++ b/mainQEE.cpp
@@ -134,6 +134,7 @@ int main(int argc, char *argv[]) {
         else if (strcmp(argv[i],"-v")==0)   { r.useVerbose=true;                                                            }
 
         else if (strcmp(argv[i],"-trackE")==0)   { r.trackE=true;                                                           }
+        else if (strcmp(argv[i],"-randInit")==0) { r.randomStart=true;                                                      }
         
         else printf("Unrecognized command! '%s'\n",argv[i]);
         
diff --git a/qEpitopeEval.cpp b/qEpitopeEval.cpp
index 4807def..ba6d3bc 100755
--- a/qEpitopeEval.cpp
+++ b/qEpitopeEval.cpp
@@ -76,6 +76,16 @@ void runEpitopeEval(RunParameters &r, std::vector<std::vector<int> > &eStates, s
         for (int i=0;i<N;i++) fscanf(startIn,"%d",&lattice[i]);
     
     }
+
+    else if (r.randomStart) {
+        srand((unsigned)time(0));
+        int numThreads=1;
+        unsigned long seed[numThreads];
+        for (int i=0;i<numThreads;i++) seed[i]=rand();
+        MT::MersenneTwist mt;
+        mt.init_genrand(seed[0]);
+        for (int i=0; i<N; i++) lattice[i] = mt.uniform_int_closed(0, (int)p[i].size());
+    }
     
     else { for (int i=0;i<N;i++) lattice[i]=(int) p[i].size(); }
     
diff --git a/qEpitopeEval.h b/qEpitopeEval.h
index a4d3e24..6ee3345 100755
--- a/qEpitopeEval.h
+++ b/qEpitopeEval.h
@@ -36,6 +36,7 @@ public:
     bool savePk;                    // If true, record P(k) mutations curve
 
     bool trackE;                    // if true, save relative contribution to energy from fields and couplings
+    bool randomStart;               // if true, set the initial state randomly
     
     RunParameters() {
         
@@ -53,6 +54,7 @@ public:
         savePk=false;
 
         trackE = false;
+        randomStart = false;
         
     }
     std::string getInfile()                   { return (directory+"/"+infile+".j");                   }