From 94b6ed7d9c591ae2df0957a25e73bd2e5c69f6e1 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Tue, 4 Jan 2022 17:42:30 +0100
Subject: Added option for spiking the interaction tensor to study minima below
 the threshold.

---
 pureStokes.cpp | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

(limited to 'pureStokes.cpp')

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;
-- 
cgit v1.2.3-54-g00ecf