Add option to initialize state completely randomly.

diff --git a/mainQEE.cpp b/mainQEE.cpp
index 3cd263b..04676cd 100755 (executable)
--- 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;                                                            }\r
 \r
         else if (strcmp(argv[i],"-trackE")==0)   { r.trackE=true;                                                           }\r
+        else if (strcmp(argv[i],"-randInit")==0) { r.randomStart=true;                                                      }\r
         \r
         else printf("Unrecognized command! '%s'\n",argv[i]);\r
         \r

diff --git a/qEpitopeEval.cpp b/qEpitopeEval.cpp
index 4807def..ba6d3bc 100755 (executable)
--- 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]);\r
     \r
     }\r
+\r
+    else if (r.randomStart) {\r
+        srand((unsigned)time(0));\r
+        int numThreads=1;\r
+        unsigned long seed[numThreads];\r
+        for (int i=0;i<numThreads;i++) seed[i]=rand();\r
+        MT::MersenneTwist mt;\r
+        mt.init_genrand(seed[0]);\r
+        for (int i=0; i<N; i++) lattice[i] = mt.uniform_int_closed(0, (int)p[i].size());\r
+    }\r
     \r
     else { for (int i=0;i<N;i++) lattice[i]=(int) p[i].size(); }\r
     \r

diff --git a/qEpitopeEval.h b/qEpitopeEval.h
index a4d3e24..6ee3345 100755 (executable)
--- a/qEpitopeEval.h
+++ b/qEpitopeEval.h
@@ -36,6 +36,7 @@ public:
     bool savePk;                    // If true, record P(k) mutations curve\r
 \r
     bool trackE;                    // if true, save relative contribution to energy from fields and couplings\r
+    bool randomStart;               // if true, set the initial state randomly\r
     \r
     RunParameters() {\r
         \r
@@ -53,6 +54,7 @@ public:
         savePk=false;\r
 \r
         trackE = false;\r
+        randomStart = false;\r
         \r
     }\r
     std::string getInfile()                   { return (directory+"/"+infile+".j");                   }\r