import matplotlib matplotlib.use('TkAgg') import matplotlib.pyplot as plt from ohmpi.utils import change_config change_config('../configs/config_mb_2023_4_mux_2023.py', verbose=False) import importlib import os import time import logging from ohmpi.config import HARDWARE_CONFIG import pandas as pd stand_alone = False part_of_hardware_system = False within_ohmpi = True # Stand alone if stand_alone: ctl_module = importlib.import_module(f'ohmpi.hardware_components.{HARDWARE_CONFIG["ctl"].pop("model")}') pwr_module = importlib.import_module(f'ohmpi.hardware_components.{HARDWARE_CONFIG["pwr"].pop("model")}') tx_module = importlib.import_module(f'ohmpi.hardware_components.{HARDWARE_CONFIG["tx"].pop("model")}') rx_module = importlib.import_module(f'ohmpi.hardware_components.{HARDWARE_CONFIG["rx"].pop("model")}') ctl = ctl_module.Ctl() HARDWARE_CONFIG['tx'].update({'ctl': ctl, 'exec_logger': ctl.exec_logger, 'soh_logger': ctl.soh_logger}) HARDWARE_CONFIG['rx'].update({'ctl': ctl, 'exec_logger': ctl.exec_logger, 'soh_logger': ctl.soh_logger}) HARDWARE_CONFIG['tx'].update({'connection': HARDWARE_CONFIG['tx'].pop('connection', ctl.interfaces[ HARDWARE_CONFIG['tx'].pop( 'interface_name', 'i2c')])}) HARDWARE_CONFIG['rx'].update({'connection': HARDWARE_CONFIG['rx'].pop('connection', ctl.interfaces[ HARDWARE_CONFIG['rx'].pop( 'interface_name', 'i2c')])}) HARDWARE_CONFIG['pwr'].update({'connection': HARDWARE_CONFIG['pwr'].pop('connection', ctl.interfaces[ HARDWARE_CONFIG['pwr'].pop( 'interface_name', None)])}) rx = rx_module.Rx(**HARDWARE_CONFIG['rx']) tx = tx_module.Tx(**HARDWARE_CONFIG['tx']) pwr = pwr_module.Pwr(**HARDWARE_CONFIG['pwr']) role = 'A' mux_id = f'mux_{role}' mux_boards = [] mux_module = importlib.import_module( f'ohmpi.hardware_components.{HARDWARE_CONFIG["mux"]["boards"][mux_id].pop("model")}') MUX_CONFIG = HARDWARE_CONFIG['mux']['boards'][mux_id] MUX_CONFIG.update({'ctl': ctl, 'connection': MUX_CONFIG.pop('connection', ctl.interfaces[ MUX_CONFIG.pop('interface_name', 'i2c_ext')]), 'exec_logger': ctl.exec_logger, 'soh_logger': ctl.soh_logger}) MUX_CONFIG.update({'id': mux_id}) mux = mux_module.Mux(**MUX_CONFIG) mux.reset() mux.test({role: [i for i in range(1, 65)]}, activation_time=.5) mux.reset() # mux as part of a OhmPiHardware system if part_of_hardware_system: from ohmpi.hardware_system import OhmPiHardware print('Starting test of as part of an OhmPiHardware system.') # mux_id = 'mux_03' k = OhmPiHardware() k.exec_logger.setLevel(logging.DEBUG) # Test mux switching k.reset_mux() # k.switch_mux(electrodes=[1, 4, 2, 3], roles=['A', 'B', 'M', 'N'], state='on') # time.sleep(1.) # k.switch_mux(electrodes=[1, 4, 2, 3], roles=['A', 'B', 'M', 'N'], state='off') # k.mux_boards[mux_id].test(activation_time=.4) k.test_mux() k.reset_mux() if within_ohmpi: from ohmpi.ohmpi import OhmPi # from ohmpi.plots import plot_exec_log print('Starting test with OhmPi.') k = OhmPi() # k.get_data() k.load_sequence(os.path.join(os.path.dirname(__file__), '../sequences/test_circuit_1423.txt')) k.reset_mux() k.run_multiple_sequences(sequence_delay=20, nb_meas=3) # k.run_sequence(injection_duration=0.2) # k.rs_check(tx_volt=4) # k.test_mux(mux_id=None, activation_time=0.2) # k._hw.switch_mux([A, B, M, N], state='on') # k._hw.vab_square_wave(12.,1., cycles=2) # k._hw.switch_mux([A, B, M, N], state='off') # k._hw.calibrate_rx_bias() # electrodes 1 4 2 3 should be connected to a reference circuit # k._hw.rx._bias = -1.38 # print(f'Resistance: {k._hw.last_rho :.2f} ohm, dev. {k._hw.last_dev:.2f} %, rx bias: {k._hw.rx._bias:.2f} mV') # k._hw._plot_readings() # k._hw.switch_mux([A, B, M, N], state='on') # k._hw.vab_square_wave(12., cycle_duration=10., cycles=3) # k._hw.switch_mux([A, B, M, N], state='off') # print(f'OhmPiHardware Resistance: {k._hw.last_rho :.2f} ohm, dev. {k._hw.last_dev:.2f} %, rx bias: {k._hw.rx._bias:.2f} mV') # k._hw._plot_readings() # A, B, M, N = (16, 13, 15, 14) # A, B, M, N = (1, 4, 2, 3) # d = k.run_measurement([A, B, M, N], injection_duration=0.5, nb_stack=2, duty_cycle=0.5) # print(d) # k._hw._plot_readings() # print(f'OhmPiHardware: Resistance: {k._hw.last_resistance() :.2f} ohm, dev. {k._hw.last_dev():.2f} %, sp: {k._hw.sp:.2f} mV, rx bias: {k._hw.rx._bias:.2f} mV') # print(f'OhmPi: Resistance: {d["R [ohm]"] :.2f} ohm, dev. {d["R_std [%]"]:.2f} %, rx bias: {k._hw.rx._bias:.2f} mV') # k._hw._plot_readings(save_fig=False) # k.run_sequence(nb_stack=2, injection_duration=0.2, duty_cycle=0.5) # ddir = os.path.join(os.path.dirname(__file__), '../data/') # print(sorted(os.listdir(ddir))) # df = pd.read_csv(os.path.join(ddir, sorted(os.listdir(ddir))[-3])) # fig, ax = plt.subplots() # df['R [Ohm]'].plot(ax=ax) # fig.show(block=True) # plot_exec_log('ohmpi/logs/exec.log') change_config('../configs/config_default.py', verbose=False)