From 95df02c90a455c2e539e795acb1921b688e8bc66 Mon Sep 17 00:00:00 2001
From: Jaron Kent-Dobias <jaron@kent-dobias.com>
Date: Wed, 17 Feb 2021 15:46:00 +0100
Subject: Algorithmic tweaks to rope relaxation now suceeds in finding Stokes
 lines relatively quickly.

---
 dynamics.hpp | 9 ++-------
 1 file changed, 2 insertions(+), 7 deletions(-)

(limited to 'dynamics.hpp')

diff --git a/dynamics.hpp b/dynamics.hpp
index a27ccb4..22d590a 100644
--- a/dynamics.hpp
+++ b/dynamics.hpp
@@ -7,16 +7,11 @@
 #include "p-spin.hpp"
 #include "stereographic.hpp"
 
-template <class Scalar>
-Vector<Scalar> normalize(const Vector<Scalar>& z) {
-  return z * sqrt((double)z.size() / (Scalar)(z.transpose() * z));
-}
-
 class gradientDescentStallException: public std::exception {
   virtual const char* what() const throw() {
     return "Gradient descent stalled.";
   }
-} gradientDescentStall;
+};
 
 template <class Scalar, int p>
 std::tuple<double, Vector<Scalar>> gradientDescent(const Tensor<Scalar, p>& J, const Vector<Scalar>& z0, double ε, double γ0 = 1, double δγ = 2) {
@@ -41,7 +36,7 @@ std::tuple<double, Vector<Scalar>> gradientDescent(const Tensor<Scalar, p>& J, c
     }
 
     if (γ < 1e-50) {
-      throw gradientDescentStall;
+      throw gradientDescentStallException();
     }
   }
 
-- 
cgit v1.2.3-54-g00ecf