From 203d0a83670c65fdc13c00c7246e2b5429d4fecb Mon Sep 17 00:00:00 2001
From: su530201 <olivier.kaufmann@umons.ac.be>
Date: Wed, 18 Oct 2023 23:45:31 +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 883412fb..c6ad2a89 100644
--- a/ohmpi/hardware_system.py
+++ b/ohmpi/hardware_system.py
@@ -512,11 +512,11 @@ class OhmPiHardware:
# ax.plot([0, vab_max], [0, vmn_lower_bound * vab_max / vab[k]], '-g', alpha=(k + 1) / n_steps)
# bounds on rab
print(f'rab_lb: {vab_list[k] / iab_upper_bound}')
- rab_lower_bound = np.min([0.1, np.abs(vab_list[k] / iab_upper_bound)])
- rab_upper_bound = np.min([0.1, np.abs(vab_list[k] / iab_lower_bound)])
+ rab_lower_bound = np.max([0.1, np.abs(vab_list[k] / iab_upper_bound)])
+ rab_upper_bound = np.max([0.1, np.abs(vab_list[k] / iab_lower_bound)])
# bounds on r
- r_lower_bound = np.min([0.01, np.abs(vmn_lower_bound / iab_upper_bound)])
- r_upper_bound = np.min([0.01, np.abs(vmn_upper_bound / iab_lower_bound)])
+ r_lower_bound = np.max([0.01, np.abs(vmn_lower_bound / iab_upper_bound)])
+ r_upper_bound = np.max([0.01, np.abs(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