From d5cea99d8ba1dd12e506a1691c9986cabacb7ad0 Mon Sep 17 00:00:00 2001
From: Lionel GUEZ <guez@lmd.ens.fr>
Date: Fri, 29 Apr 2022 21:36:55 +0200
Subject: [PATCH] Do not recompute Omega for each call
---
cost_function.py | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/cost_function.py b/cost_function.py
index c87c0ffa..d4454b27 100755
--- a/cost_function.py
+++ b/cost_function.py
@@ -26,6 +26,8 @@ import util_eddies
import bisect
import argparse
+Omega = 2 * math.pi / 86164.
+
def calculate_radii_rossby(list_eddies, e_overestim, handlers, array_d_init):
"""Compute average on list_eddies of Rossby number and radius of
maximum speed contour.
@@ -35,7 +37,6 @@ def calculate_radii_rossby(list_eddies, e_overestim, handlers, array_d_init):
radii = 0 # in m
rossby = 0
n_valid_Rossby = 0
- Omega = 2 * math.pi / 86164.
n_eddies = len(list_eddies)
for n in list_eddies:
--
GitLab