Catch errors in the J-perturbing loop and save entire vectors, not just their norm difference.

1 files changed, 9 insertions, 4 deletions

diff --git a/langevin.cpp b/langevin.cpp
index acac6c7..a4c9cf3 100644
--- a/langevin.cpp
+++ b/langevin.cpp
@@ -1,3 +1,4 @@
+#include <exception>
 #include <getopt.h>
 
 #include <eigen3/Eigen/LU>
@@ -131,7 +132,7 @@ int main(int argc, char* argv[]) {
 
   // simulation parameters
   double ε = 1e-4;
-  double εJ = 5e-2;
+  double εJ = 1e-2;
   double δ = 1e-2;   // threshold for determining saddle
   double Δ = 1e-3;
   double γ = 1e-2;   // step size
@@ -219,10 +220,14 @@ int main(int argc, char* argv[]) {
 
     iterateOver<Scalar, p>(Jp, perturbJ);
 
-    Vector zSaddleNew = findSaddle(Jp, zSaddle, ε);
-    Vector zSaddlePrevNew = findSaddle(Jp, zSaddlePrev, ε);
+    try {
+      Vector zSaddleNew = findSaddle(Jp, zSaddle, ε);
+      Vector zSaddlePrevNew = findSaddle(Jp, zSaddlePrev, ε);
 
-    std::cout << (zSaddleNew - zSaddlePrevNew).norm() << std::endl;
+      std::cout << zSaddleNew.transpose() << " " << zSaddlePrevNew.transpose() << std::endl;
+    } catch (std::exception& e) {
+      std::cerr << "Couldn't find a saddle with new couplings, skipping." << std::endl;
+    }
   }
 
   return 0;