command line flag for temperature convergence precision

1 files changed, 6 insertions, 2 deletions

diff --git a/hadamard.cpp b/hadamard.cpp
index 8eba184..fef5af2 100644
--- a/hadamard.cpp
+++ b/hadamard.cpp
@@ -57,13 +57,14 @@ int main(int argc, char* argv[]) {
   unsigned N = 16;
   unsigned k = 2;
   double ε = 0.01;
+  double ε2 = 0.01;
 
   unsigned M = 10;
   unsigned m = 1e4;
 
   int opt;
 
-  while ((opt = getopt(argc, argv, "k:b:c:n:t:N:M:e:m:")) != -1) {
+  while ((opt = getopt(argc, argv, "k:b:c:n:t:N:M:e:m:f:")) != -1) {
     switch (opt) {
     case 'k':
       k = atoi(optarg);
@@ -77,6 +78,9 @@ int main(int argc, char* argv[]) {
     case 'e':
       ε = atof(optarg);
       break;
+    case 'f':
+      ε2 = atof(optarg);
+      break;
     case 'n':
       m = (unsigned)atof(optarg);
       break;
@@ -112,7 +116,7 @@ int main(int argc, char* argv[]) {
   std::cout << "Beginning simulation of " << n << ".\n";
   std::cout << "Beginning " << n_tuning << " tuning tempering updates of " << M
             << " sweeps each.\n";
-  std::vector<double> f = p.tune(n_tuning, M, ε, 0.05);
+  std::vector<double> f = p.tune(n_tuning, M, ε, ε2);
   std::cout << "Finished tuning, beginning " << m << " measurement tempering updates of " << M
             << " sweeps each.\n";
   p.run(m, M);