shepherd-core 2024.8.1__tar.gz → 2024.8.2__tar.gz

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: shepherd_core
-Version: 2024.8.1
+Version: 2024.8.2
 Summary: Programming- and CLI-Interface for the h5-dataformat of the Shepherd-Testbed
 Author-email: Ingmar Splitt <ingmar.splitt@tu-dresden.de>
 Maintainer-email: Ingmar Splitt <ingmar.splitt@tu-dresden.de>

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/examples/eenv_generator.py

@@ -35,7 +35,7 @@ for _v, _c in product(voltages_V, currents_A):
             voltage_vector = np.linspace(_v, _v, int(file.samplerate_sps * duration_s))
             current_vector = np.linspace(_c, _c, int(file.samplerate_sps * duration_s))
 
-            for idx in trange(repetitions, desc="generate"):
+            for idx in trange(repetitions, desc="generate", leave=False):
                 timestamps = idx * duration_s + timestamp_vector
                 file.append_iv_data_si(timestamps, voltage_vector, current_vector)
 

shepherd_core-2024.8.2/examples/simulate_vharvester.py

@@ -0,0 +1,74 @@
+"""Demonstrate behavior of Virtual Harvester Algorithms.
+
+- simulation is based on ivsamples derived from a solar-isc-voc-recording during a jogging-trip
+- harvesting is done by various algorithms and preconfigured virtual harvesters
+- results are printed on console (harvested energy)
+
+Output:
+E_out = 0.000 mWs -> cv20
+E_out = 17.165 mWs -> cv10
+E_out = 17.427 mWs -> mppt_voc
+E_out = 17.242 mWs -> mppt_bq_solar
+E_out = 13.998 mWs -> mppt_bq_thermoelectric
+E_out = 15.202 mWs -> mppt_po
+E_out = 17.811 mWs -> mppt_opt
+"""
+
+from pathlib import Path
+
+from shepherd_core import Reader
+from shepherd_core.data_models import VirtualHarvesterConfig
+from shepherd_core.vsource import simulate_harvester
+from shepherd_data import ivonne
+
+# config simulation
+sim_duration = 32
+file_ivonne = Path(__file__).parent.parent.parent / "shepherd_data/examples/jogging_10m.iv"
+file_ivcurve = Path(__file__).parent / "jogging_ivcurve.h5"
+
+hrv_list = [
+    "cv20",
+    # ⤷ fails due to lower solar voltage
+    "cv10",
+    "mppt_voc",
+    "mppt_bq_solar",
+    # ⤷ bq needs 16 s to start -> bad performance for this demo
+    "mppt_bq_thermoelectric",
+    # ⤷ thermoelectric setpoint -> bad performance for solar
+    "mppt_po",
+    "mppt_opt",
+]
+
+save_files: bool = False
+
+# convert IVonne to IVCurve
+if not file_ivcurve.exists():
+    with ivonne.Reader(file_ivonne) as db:
+        db.convert_2_ivcurves(file_ivcurve, duration_s=sim_duration)
+
+# Input Statistics
+with Reader(file_ivcurve, verbose=False) as file:
+    window_size = file.get_window_samples()
+    I_in_max = 0.0
+    for _t, _v, _i in file.read_buffers():
+        I_in_max = max(I_in_max, _i.max())
+    print(
+        f"Input-file: \n"
+        f"\tE_in = {file.energy() * 1e3:.3f} mWs (not representative)\n"
+        f"\tI_in_max = {I_in_max * 1e3:.3f} mA\n"
+        f"\twindow_size = {window_size} n\n",
+    )
+
+# Simulation
+for hrv_name in hrv_list:
+    file_output = (
+        file_ivcurve.with_name(file_ivcurve.stem + "_" + hrv_name + file_ivcurve.suffix)
+        if save_files
+        else None
+    )
+    E_out_Ws = simulate_harvester(
+        config=VirtualHarvesterConfig(name=hrv_name),
+        path_input=file_ivcurve,
+        path_output=file_output,
+    )
+    print(f"E_out = {E_out_Ws * 1e3:.3f} mWs -> {hrv_name}")

shepherd_core-2024.8.2/examples/simulate_vsource.py

@@ -0,0 +1,55 @@
+"""Demonstrate behavior of Virtual Source Algorithms.
+
+The emulation recreates an observer-cape, the virtual Source and a virtual target
+- input = hdf5-file with a harvest-recording
+- output = hdf5-file
+- config is currently hardcoded, but it could be an emulation-task
+- target is currently a simple resistor
+
+The output file can be analyzed and plotted with shepherds tool suite.
+
+Output:
+E_out = 0.007 mWs -> direct
+E_out = 0.033 mWs -> diode+capacitor
+E_out = 0.033 mWs -> diode+resistor+capacitor
+E_out = 14.964 mWs -> BQ25504
+E_out = 15.042 mWs -> BQ25504s
+E_out = 13.909 mWs -> BQ25570
+E_out = 14.073 mWs -> BQ25570s
+"""
+
+from pathlib import Path
+
+from shepherd_core.data_models import VirtualSourceConfig
+from shepherd_core.vsource import ConstantCurrentTarget
+from shepherd_core.vsource import simulate_source
+
+# config simulation
+file_input = Path(__file__).parent / "jogging_ivcurve.h5"
+
+src_list = [
+    "direct",
+    "diode+capacitor",
+    "diode+resistor+capacitor",
+    "BQ25504",
+    "BQ25504s",
+    "BQ25570",
+    "BQ25570s",
+]
+tgt = ConstantCurrentTarget(i_active_A=1e-3, i_sleep_A=200e-9)
+save_files = False
+
+for src_name in src_list:
+    file_output = (
+        file_input.with_name(file_input.stem + "_emu_" + src_name + file_input.suffix)
+        if save_files
+        else None
+    )
+
+    e_out_Ws = simulate_source(
+        config=VirtualSourceConfig(name=src_name),
+        target=tgt,
+        path_input=file_input,
+        path_output=file_output,
+    )
+    print(f"E_out = {e_out_Ws * 1e3:.3f} mWs -> {src_name}")

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/shepherd_core/data_models/content/virtual_source_fixture.yaml

@@ -15,7 +15,7 @@
     interval_startup_delay_drain_ms: 0
 
     harvester:
-      name: mppt_opt # harvester only active if input is "ivcurves"
+      name: mppt_opt # harvester only active if input is "ivcurve"
 
     V_input_max_mV: 10000
     I_input_max_mA: 4200

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/shepherd_core/version.py

@@ -1,3 +1,3 @@
 """Separated string avoids circular imports."""
 
-version: str = "2024.8.1"
+version: str = "2024.8.2"

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/shepherd_core/vsource/__init__.py

@@ -1,13 +1,23 @@
 """Simulation model of the virtual source."""
 
+from .target_model import ConstantCurrentTarget
+from .target_model import ConstantPowerTarget
+from .target_model import ResistiveTarget
 from .virtual_converter_model import PruCalibration
 from .virtual_converter_model import VirtualConverterModel
 from .virtual_harvester_model import VirtualHarvesterModel
+from .virtual_harvester_simulation import simulate_harvester
 from .virtual_source_model import VirtualSourceModel
+from .virtual_source_simulation import simulate_source
 
 __all__ = [
     "PruCalibration",
     "VirtualConverterModel",
     "VirtualHarvesterModel",
     "VirtualSourceModel",
+    "simulate_harvester",
+    "simulate_source",
+    "ResistiveTarget",
+    "ConstantCurrentTarget",
+    "ConstantPowerTarget",
 ]

shepherd_core-2024.8.2/shepherd_core/vsource/target_model.py

@@ -0,0 +1,62 @@
+"""Virtual targets with different characteristics.
+
+TODO: add more targets
+  - diode
+  - constant power
+  - constant current
+  - msp430 (const I)
+  - nRF (constant power due to regulator)
+  - riotee
+"""
+
+from abc import ABC
+from abc import abstractmethod
+
+
+class TargetABC(ABC):
+    """Abstract base class for all targets."""
+
+    @abstractmethod
+    def step(self, voltage_uV: int, *, pwr_good: bool) -> float:
+        """Calculate one time step and return drawn current in nA."""
+
+
+class ResistiveTarget(TargetABC):
+    """Predictable target for matching the real world."""
+
+    def __init__(self, resistance_Ohm: float, *, controlled: bool = False) -> None:
+        if resistance_Ohm <= 1e-3:
+            raise ValueError("Resistance must be greater than 1 mOhm.")
+        self.r_kOhm = 1e-3 * resistance_Ohm
+        self.ctrl = controlled
+
+    def step(self, voltage_uV: int, *, pwr_good: bool) -> float:
+        if pwr_good or not self.ctrl:
+            return voltage_uV / self.r_kOhm  # = nA
+        return 0
+
+
+class ConstantCurrentTarget(TargetABC):
+    """Recreate simple MCU without integrated regulator."""
+
+    def __init__(self, i_active_A: float, i_sleep_A: float = 0) -> None:
+        if i_active_A <= 0 or i_sleep_A <= 0:
+            raise ValueError("Current must be greater than 0.")
+        self.i_active_nA = 1e9 * i_active_A
+        self.i_sleep_nA = 1e9 * i_sleep_A
+
+    def step(self, voltage_uV: int, *, pwr_good: bool) -> float:  # noqa: ARG002
+        return self.i_active_nA if pwr_good else self.i_sleep_nA
+
+
+class ConstantPowerTarget(TargetABC):
+    """Recreate MCU with integrated regulator."""
+
+    def __init__(self, p_active_W: float, p_sleep_W: float = 0) -> None:
+        if p_active_W <= 0 or p_sleep_W <= 0:
+            raise ValueError("Power must be greater than 0.")
+        self.p_active_fW = 1e15 * p_active_W
+        self.p_sleep_fW = 1e15 * p_sleep_W
+
+    def step(self, voltage_uV: int, *, pwr_good: bool) -> float:
+        return (self.p_active_fW if pwr_good else self.p_sleep_fW) / voltage_uV  # = nA

shepherd_core-2024.8.2/shepherd_core/vsource/virtual_harvester_simulation.py

@@ -0,0 +1,71 @@
+"""Simulate behavior of virtual harvester algorithms.
+
+The simulation recreates an observer-cape and the virtual harvester
+- input = hdf5-file with an ivcurve
+- output = optional as hdf5-file
+
+The output file can be analyzed and plotted with shepherds tool suite.
+"""
+
+from contextlib import ExitStack
+from pathlib import Path
+from typing import Optional
+
+from tqdm import tqdm
+
+from .. import CalibrationHarvester
+from .. import Reader
+from .. import Writer
+from ..data_models.content.virtual_harvester import HarvesterPRUConfig
+from ..data_models.content.virtual_harvester import VirtualHarvesterConfig
+from .virtual_harvester_model import VirtualHarvesterModel
+
+
+def simulate_harvester(
+    config: VirtualHarvesterConfig, path_input: Path, path_output: Optional[Path] = None
+) -> float:
+    """Simulate behavior of virtual harvester algorithms.
+
+    Fn return the harvested energy.
+    """
+    stack = ExitStack()
+    file_inp = Reader(path_input, verbose=False)
+    stack.enter_context(file_inp)
+    cal_inp = file_inp.get_calibration_data()
+
+    if path_output:
+        cal_hrv = CalibrationHarvester()
+        file_out = Writer(
+            path_output, cal_data=cal_hrv, mode="harvester", verbose=False, force_overwrite=True
+        )
+        stack.enter_context(file_out)
+        cal_out = file_out.get_calibration_data()
+        file_out.store_hostname("hrv_sim_" + config.name)
+
+    hrv_pru = HarvesterPRUConfig.from_vhrv(
+        config,
+        for_emu=True,
+        dtype_in=file_inp.get_datatype(),
+        window_size=file_inp.get_window_samples(),
+    )
+    hrv = VirtualHarvesterModel(hrv_pru)
+    e_out_Ws = 0.0
+
+    for _t, v_inp, i_inp in tqdm(
+        file_inp.read_buffers(is_raw=True), total=file_inp.buffers_n, desc="Buffers", leave=False
+    ):
+        v_uV = cal_inp.voltage.raw_to_si(v_inp) * 1e6
+        i_nA = cal_inp.current.raw_to_si(i_inp) * 1e9
+        length = min(v_uV.size, i_nA.size)
+        for _n in range(length):
+            v_uV[_n], i_nA[_n] = hrv.ivcurve_sample(
+                _voltage_uV=int(v_uV[_n]), _current_nA=int(i_nA[_n])
+            )
+        e_out_Ws += (v_uV * i_nA).sum() * 1e-15 * file_inp.sample_interval_s
+        if path_output:
+            v_out = cal_out.voltage.si_to_raw(v_uV / 1e6)
+            i_out = cal_out.current.si_to_raw(i_nA / 1e9)
+            file_out.append_iv_data_raw(_t, v_out, i_out)
+
+    stack.close()
+    return e_out_Ws

shepherd_core-2024.8.2/shepherd_core/vsource/virtual_source_simulation.py

@@ -0,0 +1,78 @@
+"""Simulate behavior of virtual source algorithms.
+
+The simulation recreates an observer-cape, the virtual Source and a virtual target
+- input = hdf5-file with a harvest-recording
+- output = optional as hdf5-file
+
+The output file can be analyzed and plotted with shepherds tool suite.
+"""
+
+from contextlib import ExitStack
+from pathlib import Path
+from typing import Optional
+
+from tqdm import tqdm
+
+from .. import CalibrationEmulator
+from .. import Reader
+from .. import Writer
+from ..data_models import VirtualSourceConfig
+from . import VirtualSourceModel
+from .target_model import TargetABC
+
+
+def simulate_source(
+    config: VirtualSourceConfig,
+    target: TargetABC,
+    path_input: Path,
+    path_output: Optional[Path] = None,
+) -> float:
+    """Simulate behavior of virtual source algorithms.
+
+    FN returns the consumed energy of the target.
+    """
+    stack = ExitStack()
+    file_inp = Reader(path_input, verbose=False)
+    stack.enter_context(file_inp)
+    cal_emu = CalibrationEmulator()
+    cal_inp = file_inp.get_calibration_data()
+
+    if path_output:
+        file_out = Writer(
+            path_output, cal_data=cal_emu, mode="emulator", verbose=False, force_overwrite=True
+        )
+        stack.enter_context(file_out)
+        file_out.store_hostname("emu_sim_" + config.name)
+        file_out.store_config(config.model_dump())
+        cal_out = file_out.get_calibration_data()
+
+    src = VirtualSourceModel(
+        config, cal_emu, log_intermediate=False, window_size=file_inp.get_window_samples()
+    )
+    i_out_nA = 0
+    e_out_Ws = 0.0
+
+    for _t, v_inp, i_inp in tqdm(
+        file_inp.read_buffers(is_raw=True), total=file_inp.buffers_n, desc="Buffers", leave=False
+    ):
+        v_uV = 1e6 * cal_inp.voltage.raw_to_si(v_inp)
+        i_nA = 1e9 * cal_inp.current.raw_to_si(i_inp)
+
+        for _n in range(len(_t)):
+            v_uV[_n] = src.iterate_sampling(
+                V_inp_uV=int(v_uV[_n]),
+                I_inp_nA=int(i_nA[_n]),
+                I_out_nA=i_out_nA,
+            )
+            i_out_nA = target.step(int(v_uV[_n]), pwr_good=src.cnv.get_power_good())
+            i_nA[_n] = i_out_nA
+            # TODO: src.cnv.get_I_mod_out_nA() has more internal drains
+
+        e_out_Ws += (v_uV * i_nA).sum() * 1e-15 * file_inp.sample_interval_s
+        if path_output:
+            v_out = cal_out.voltage.si_to_raw(1e-6 * v_uV)
+            i_out = cal_out.current.si_to_raw(1e-9 * i_nA)
+            file_out.append_iv_data_raw(_t, v_out, i_out)
+
+    stack.close()
+    return e_out_Ws

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/shepherd_core.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: shepherd_core
-Version: 2024.8.1
+Version: 2024.8.2
 Summary: Programming- and CLI-Interface for the h5-dataformat of the Shepherd-Testbed
 Author-email: Ingmar Splitt <ingmar.splitt@tu-dresden.de>
 Maintainer-email: Ingmar Splitt <ingmar.splitt@tu-dresden.de>

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/shepherd_core.egg-info/SOURCES.txt

@@ -8,11 +8,11 @@ examples/experiment_models.py
 examples/firmware_model.py
 examples/firmware_modification.py
 examples/inventory.py
+examples/simulate_vharvester.py
+examples/simulate_vsource.py
 examples/uart_decode_waveform.py
 examples/uart_raw2.csv
-examples/vharvester_simulation.py
-examples/vsource_emulation.py
-examples/vsource_simulation.py
+examples/vsource_debug_sim.py
 shepherd_core/__init__.py
 shepherd_core/calibration_hw_def.py
 shepherd_core/commons.py
@@ -88,9 +88,12 @@ shepherd_core/testbed_client/client_web.py
 shepherd_core/testbed_client/fixtures.py
 shepherd_core/testbed_client/user_model.py
 shepherd_core/vsource/__init__.py
+shepherd_core/vsource/target_model.py
 shepherd_core/vsource/virtual_converter_model.py
 shepherd_core/vsource/virtual_harvester_model.py
+shepherd_core/vsource/virtual_harvester_simulation.py
 shepherd_core/vsource/virtual_source_model.py
+shepherd_core/vsource/virtual_source_simulation.py
 tests/__init__.py
 tests/conftest.py
 tests/test_cal_hw.py

{shepherd_core-2024.8.1 → shepherd_core-2024.8.2}/tests/test_examples.py

@@ -18,8 +18,9 @@ examples = [
     "experiment_models.py",
     "inventory.py",
     "uart_decode_waveform.py",
-    "vharvester_simulation.py",
-    "vsource_simulation.py",
+    "simulate_vharvester.py",
+    "simulate_vsource.py",
+    "vsource_debug_sim.py",
 ]
 
 

shepherd_core-2024.8.1/examples/vharvester_simulation.py

@@ -1,96 +0,0 @@
-"""Demonstrate behavior of Virtual Harvester Algorithms.
-
-- simulation is based on ivsamples derived from a solar-isc-voc-recording during a jogging-trip
-- harvesting is done by various algorithms and preconfigured virtual harvesters
-- results are printed on console (harvested energy)
-
-"""
-
-from contextlib import ExitStack
-from pathlib import Path
-
-from shepherd_core import CalibrationHarvester
-from shepherd_core import Reader
-from shepherd_core import Writer
-from shepherd_core.data_models import VirtualHarvesterConfig
-from shepherd_core.data_models.content.virtual_harvester import HarvesterPRUConfig
-from shepherd_core.vsource import VirtualHarvesterModel
-from shepherd_data import ivonne
-
-# config simulation
-sim_duration = 32
-file_ivonne = Path(__file__).parent.parent.parent / "shepherd_data/examples/jogging_10m.iv"
-file_ivcurve = Path(__file__).parent / "jogging_ivcurve.h5"
-
-hrv_list = [
-    "cv20",
-    # ⤷ fails due to lower solar voltage
-    "cv10",
-    "mppt_voc",
-    "mppt_bq_solar",
-    # ⤷ bq needs 16 s to start -> bad performance for this demo
-    "mppt_bq_thermoelectric",
-    # ⤷ thermoelectric setpoint -> bad performance for solar
-    "mppt_po",
-    "mppt_opt",
-]
-
-save_files: bool = True
-
-# convert IVonne to IVCurve
-if not file_ivcurve.exists():
-    with ivonne.Reader(file_ivonne) as db:
-        db.convert_2_ivcurves(file_ivcurve, duration_s=sim_duration)
-
-# Input Statistics
-with Reader(file_ivcurve, verbose=False) as file:
-    window_size = file.get_window_samples()
-    I_in_max = 0.0
-    for _t, _v, _i in file.read_buffers():
-        I_in_max = max(I_in_max, _i.max())
-    print(
-        f"Input-file: \n"
-        f"\tE_in = {file.energy() * 1e3:.3f} mWs (not representative)\n"
-        f"\tI_in_max = {I_in_max * 1e3:.3f} mA\n"
-        f"\twindow_size = {window_size} n\n",
-    )
-
-# Simulation
-for hrv_name in hrv_list:
-    E_out_Ws = 0.0
-    _cal = CalibrationHarvester()
-    if save_files:
-        stack = ExitStack()
-        file_output = file_ivcurve.with_name(
-            file_ivcurve.stem + "_" + hrv_name + file_ivcurve.suffix
-        )
-        fh_output = Writer(
-            file_output, cal_data=_cal, mode="harvester", verbose=False, force_overwrite=True
-        )
-        stack.enter_context(fh_output)
-        fh_output.store_hostname("hrv_sim_" + hrv_name)
-
-    with Reader(file_ivcurve, verbose=False) as file:
-        hrv_config = VirtualHarvesterConfig(name=hrv_name)
-        hrv_pru = HarvesterPRUConfig.from_vhrv(
-            hrv_config,
-            for_emu=True,
-            dtype_in=file.get_datatype(),
-            window_size=4 * file.get_window_samples(),
-        )
-        hrv = VirtualHarvesterModel(hrv_pru)
-        for _t, _v, _i in file.read_buffers():
-            # TODO: _t should be handed to new file without conversions
-            length = max(_v.size, _i.size)
-            for _n in range(length):
-                _v[_n], _i[_n] = hrv.ivcurve_sample(
-                    _voltage_uV=_v[_n] * 10**6, _current_nA=_i[_n] * 10**9
-                )
-            E_out_Ws += (_v * _i).sum() * 1e-15 * file.sample_interval_s
-            if save_files:
-                fh_output.append_iv_data_si(_t, _v / 1e6, _i / 1e9)
-        if save_files:
-            stack.close()
-            with Reader(file_output) as fh_output:
-                fh_output.save_metadata()
-        print(f"E_out = {E_out_Ws * 1e3:.3f} mWs -> {hrv_name}")

shepherd_core-2024.8.1/examples/vsource_emulation.py

@@ -1,86 +0,0 @@
-"""Demonstrate behavior of Virtual Source Algorithms.
-
-The emulation recreates an observer-cape, the virtual Source and a virtual target
-- input = hdf5-file with a harvest-recording
-- output = hdf5-file
-- config is currently hardcoded, but it could be an emulation-task
-- target is currently a simple resistor
-
-The output file can be analyzed and plotted with shepherds tool suite.
-"""
-# TODO: `shepherd-data emulate config.yaml`
-
-from pathlib import Path
-
-import numpy as np
-from tqdm import tqdm
-
-from shepherd_core import CalibrationEmulator
-from shepherd_core import Reader
-from shepherd_core import Writer
-from shepherd_core.data_models import VirtualHarvesterConfig
-from shepherd_core.data_models import VirtualSourceConfig
-from shepherd_core.vsource import VirtualSourceModel
-
-# config simulation
-file_input = Path(__file__).parent / "jogging_ivcurve.h5"
-
-src_list = ["BQ25504"]
-
-I_mcu_sleep_A = 3e-3
-I_mcu_active_A = 3e-3
-R_Ohm = 1000
-
-for vs_name in src_list:
-    file_output = file_input.with_name(file_input.stem + "_emu_" + vs_name + file_input.suffix)
-
-    cal_emu = CalibrationEmulator()
-    src_config = VirtualSourceConfig(
-        inherit_from=vs_name,
-        V_intermediate_init_mV=3000,
-        harvester=VirtualHarvesterConfig(name="mppt_bq_solar"),
-        C_intermediate_uF=50,
-    )
-
-    with Reader(file_input, verbose=False) as f_inp, Writer(
-        file_output, cal_data=cal_emu, mode="emulator", verbose=False
-    ) as f_out:
-        window_size = f_inp.get_window_samples()
-        f_out.store_hostname("emu_sim_" + vs_name)
-        f_out.store_config(src_config.model_dump())
-        src = VirtualSourceModel(
-            src_config, cal_emu, log_intermediate=False, window_size=window_size
-        )
-
-        I_out_nA = 0
-
-        for _t, _V_inp, _I_inp in tqdm(f_inp.read_buffers(), total=f_inp.buffers_n):
-            V_out = np.empty(_V_inp.shape)
-            I_out = np.empty(_I_inp.shape)
-
-            for _iter in range(len(_t)):
-                V_out_uV = src.iterate_sampling(
-                    V_inp_uV=_V_inp[_iter] * 10**6,
-                    I_inp_nA=_I_inp[_iter] * 10**9,
-                    I_out_nA=I_out_nA,
-                )
-                I_out_nA = 1e3 * V_out_uV / R_Ohm
-
-                V_out[_iter] = V_out_uV / 1e6
-                I_out[_iter] = I_out_nA / 1e9
-
-                # TODO: src.cnv.get_I_mod_out_nA() has more internal drains
-
-            f_out.append_iv_data_si(_t, V_out, I_out)
-
-            # listen to power-good signal
-            """
-            if src.cnv.get_power_good():
-                I_out_nA = int(I_mcu_active_A * 10 ** 9)
-                N_good += 1
-            else:
-                I_out_nA = int(I_mcu_sleep_A * 10 ** 9)
-            """
-
-    with Reader(file_output, verbose=False) as f_out:
-        f_out.save_metadata()