From 12f15f49cd8cc4ab9c809700e8cb88a0efe198d8 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 17 Feb 2021 16:14:33 +0100
Subject: Rearranged some functions among files, and wrote the normalize
 function to take generic Eigen expressions.

---
 dynamics.hpp | 6 ++----
 1 file changed, 2 insertions(+), 4 deletions(-)

(limited to 'dynamics.hpp')

diff --git a/dynamics.hpp b/dynamics.hpp
index 22d590a..d421d13 100644
--- a/dynamics.hpp
+++ b/dynamics.hpp
@@ -21,8 +21,7 @@ std::tuple<double, Vector<Scalar>> gradientDescent(const Tensor<Scalar, p>& J, c
   auto [W, dW] = WdW(J, z);
 
   while (W > ε) {
-    Vector<Scalar> zNewTmp = z - γ * dW.conjugate();
-    Vector<Scalar> zNew = normalize(zNewTmp);
+    Vector<Scalar> zNew = normalize(z - γ * dW.conjugate());
 
     auto [WNew, dWNew] = WdW(J, zNew);
 
@@ -102,8 +101,7 @@ std::tuple<double, Vector<Scalar>> metropolis(const Tensor<Scalar, p>& J, const
   std::uniform_real_distribution<double> D(0, 1);
 
   for (unsigned i = 0; i < N; i++) {
-    Vector<Scalar> zNewTmp = z + γ * randomVector<Scalar>(z.size(), d, r);
-    Vector<Scalar> zNew = normalize(zNewTmp);
+    Vector<Scalar> zNew = normalize(z + γ * randomVector<Scalar>(z.size(), d, r));
 
     double ENew = energy(J, zNew);
 
-- 
cgit v1.2.3-54-g00ecf