From fb75f0d923ba89296fc49a20d5296309a0ce5b1e Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Tue, 29 Dec 2020 13:43:52 +0100
Subject: Fixed mistake in constraint step.

---
 langevin.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/langevin.cpp b/langevin.cpp
index 73cf4c0..5287cd7 100644
--- a/langevin.cpp
+++ b/langevin.cpp
@@ -126,7 +126,7 @@ std::tuple<double, Vector, std::complex<double>, Vector> WdW(const Tensor& J, co
 
   Vector dWdz = grad.contract(conjHess, ip) - (pow(p, 2) / 2.0 * ((double)z.size() - zz(0))) * zc;
   Scalar dWdε = (-(double)p) * grad.contract(zc, ip);
-  Vector dεdz = (1 / (double)z.size()) * ((1 - 1/(double)p) * grad - ((double)p * ε) * z - (1 /(double)p) * z.contract(hess, ip));
+  Vector dεdz = (1 / (double)z.size()) * ((1 - 1/(double)p) * grad + ((double)p * ε) * z - (1 /(double)p) * z.contract(hess, ip));
 
   return {W, dWdz, dWdε(0), dεdz};
 
-- 
cgit v1.2.3-70-g09d2