diff options
| -rw-r--r-- | langevin.cpp | 2 |
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}; |