✅ Lyapunov exp of Lorentz sys + cleanup

vortexlab-uclouvain · Oct 25, 2020 · a58c46f · a58c46f
1 parent f41da7a
commit a58c46f
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 1 deletion.
diff --git a/multiflap/lorentz_main.py b/multiflap/lorentz_main.py
@@ -27,11 +27,13 @@
 
 x = [10., 10., 3.6]
 
+# Creating the object containing the eqs
 mymodel = Lorentz(a=10, b=28, c=8/3)
 
 t_f = 100
 n = 10000
 
+# calling the LyapunovExponents class
 lambda_t = LyapunovExponents(x, n, t_f, mymodel).get_lyapunov_exponent()
 
 print(lambda_t)
diff --git a/multiflap/ms_package/lyapunov_exponents.py b/multiflap/ms_package/lyapunov_exponents.py
@@ -83,10 +83,11 @@ def get_lyapunov_exponent(self):
             fx_pert = integration_pert[-1, :]
             d_j = fx_pert - fx
             norm_d_j = np.linalg.norm(d_j)
+
+            # Re-scaling the flow separation
             lambda_local[i] = np.log(norm_d_j/norm_d) 
             x = fx
             x_pert = x + (d_j*(norm_d/norm_d_j))
 
         lambda_t = np.sum(lambda_local)/(self.n*self.delta_time)
-        print(lambda_t)
         return lambda_t