From d292edd0dbfd4aa2c546078e8db697a7f3501d50 Mon Sep 17 00:00:00 2001
From: Gamelin Alexis <alexis.gamelin@synchrotron-soleil.fr>
Date: Fri, 7 Jan 2022 15:22:03 +0100
Subject: [PATCH] [Fix] Coating class with frequency array
---
collective_effects/resistive_wall.py | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/collective_effects/resistive_wall.py b/collective_effects/resistive_wall.py
index eb5b198..1bb3854 100644
--- a/collective_effects/resistive_wall.py
+++ b/collective_effects/resistive_wall.py
@@ -304,7 +304,7 @@ class Coating(WakeField):
bracket = ( (np.sign(f) - 1j) * skin1 *
(alpha * tanh + 1) / (alpha + tanh) )
else:
- valid_approx = self.thickness / skin1
+ valid_approx = self.thickness / np.min(skin1)
if valid_approx < 0.01:
print("Approximation is not valid. Returning impedance anyway.")
bracket = ( (np.sign(f) - 1j) * skin2 - 2 * 1j * self.thickness *
@@ -352,7 +352,7 @@ class Coating(WakeField):
bracket = ( (1 - 1j*np.sign(f)) * skin1 *
(alpha * tanh + 1) / (alpha + tanh) )
else:
- valid_approx = self.thickness / skin1
+ valid_approx = self.thickness / np.min(skin1)
if valid_approx < 0.01:
print("Approximation is not valid. Returning impedance anyway.")
bracket = ( (1 - 1j*np.sign(f)) * skin2 - 2 * 1j * self.thickness
--
GitLab