From 1094eaa434d791d3a511f853df9d6d66d169b2ab Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 6 Jan 2021 11:16:20 +0100
Subject: Lots of work, and fixed bug stemming from choice of matrix algorithm.

---
 langevin.cpp | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

(limited to 'langevin.cpp')

diff --git a/langevin.cpp b/langevin.cpp
index 3bc12ed..05b436e 100644
--- a/langevin.cpp
+++ b/langevin.cpp
@@ -7,6 +7,8 @@
 #include "pcg-cpp/include/pcg_random.hpp"
 #include "randutils/randutils.hpp"
 
+#include <eigen3/Eigen/LU>
+
 #include "complex_normal.hpp"
 #include "p-spin.hpp"
 #include "stereographic.hpp"
@@ -57,7 +59,8 @@ Vector findSaddle(const Tensor& J, const Vector& z0, double γ0, double δ, doub
   while (W > ε) {
     double γ = pow(r.variate<double, std::normal_distribution>(0, γ0), 2);
 
-    Vector ζNew = ζ - ddH.ldlt().solve(dH);
+    Vector dζ = ddH.partialPivLu().solve(dH);
+    Vector ζNew = ζ - dζ;
 
     double WNew;
     Vector dWNew;
-- 
cgit v1.2.3-54-g00ecf