summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--langevin.cpp2
1 files changed, 1 insertions, 1 deletions
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};