Added option for spiking the interaction tensor to study minima below the threshold.

1 files changed, 6 insertions, 2 deletions

diff --git a/pureStokes.cpp b/pureStokes.cpp
index 448d5ef..b1ed1dc 100644
--- a/pureStokes.cpp
+++ b/pureStokes.cpp
@@ -16,10 +16,11 @@ int main(int argc, char* argv[]) {
   Real δz₀ = 1e-3;
   unsigned n = 10;
   bool useMinima = false;
+  bool useSpike = false;
 
   int opt;
 
-  while ((opt = getopt(argc, argv, "N:e:d:n:m")) != -1) {
+  while ((opt = getopt(argc, argv, "N:e:d:n:ms")) != -1) {
     switch (opt) {
     case 'N':
       N = (unsigned)atof(optarg);
@@ -36,6 +37,9 @@ int main(int argc, char* argv[]) {
     case 'm':
       useMinima = true;
       break;
+    case 's':
+      useSpike = true;
+      break;
     default:
       exit(1);
     }
@@ -45,7 +49,7 @@ int main(int argc, char* argv[]) {
 
   for (unsigned i = 0; i < n; i++) {
     auto tag = std::chrono::high_resolution_clock::now();
-    collectStokesData<3>(std::to_string(tag.time_since_epoch().count()), N, r.engine(), ε, δz₀, useMinima, 1.0);
+    collectStokesData<3>(std::to_string(tag.time_since_epoch().count()), N, r.engine(), ε, δz₀, useMinima, useSpike, 1.0);
   }
 
   return 0;