From 9e138efea54627582de16f48aaae42d3b74cd0de Mon Sep 17 00:00:00 2001
From: su530201 <olivier.kaufmann@umons.ac.be>
Date: Wed, 18 Oct 2023 23:35:28 +0200
Subject: [PATCH] Fixes bug in _compute_tx_volt
---
ohmpi/hardware_system.py | 8 ++++----
1 file changed, 4 insertions(+), 4 deletions(-)
diff --git a/ohmpi/hardware_system.py b/ohmpi/hardware_system.py
index 3e73ed06..c8cc39bd 100644
--- a/ohmpi/hardware_system.py
+++ b/ohmpi/hardware_system.py
@@ -511,11 +511,11 @@ class OhmPiHardware:
# ax.plot([0, vab_max], [0, vmn_upper_bound * vab_max / vab[k]], '-r', alpha=(k + 1) / n_steps)
# ax.plot([0, vab_max], [0, vmn_lower_bound * vab_max / vab[k]], '-g', alpha=(k + 1) / n_steps)
# bounds on rab
- rab_lower_bound = vab_list[k] / iab_upper_bound
- rab_upper_bound = vab_list[k] / iab_lower_bound
+ rab_lower_bound = np.min([0.1, vab_list[k] / iab_upper_bound])
+ rab_upper_bound = np.min([0.1, vab_list[k] / iab_lower_bound])
# bounds on r
- r_lower_bound = vmn_lower_bound / iab_upper_bound
- r_upper_bound = vmn_upper_bound / iab_lower_bound
+ r_lower_bound = np.min([0.01, vmn_lower_bound / iab_upper_bound])
+ r_upper_bound = np.min([0.01, vmn_upper_bound / iab_lower_bound])
# conditions for vab update
cond_vmn_max = rab_lower_bound / r_upper_bound * vmn_max
cond_p_max = np.sqrt(p_max * rab_lower_bound)
--
GitLab