nsqdriver 0.12.17__cp313-cp313-manylinux2014_x86_64.manylinux_2_17_x86_64.whl

nsqdriver/NS_MCI.py

@@ -0,0 +1,597 @@
+import enum
+import traceback
+import pickle
+import socket
+import atexit
+import struct
+import warnings
+import xmlrpc.client
+from xmlrpc.client import Transport
+from multiprocessing import shared_memory
+from functools import wraps
+from enum import Enum
+
+import numpy as np
+
+try:
+    import waveforms
+    HAS_WAVEFORMS = True
+except ImportError as e:
+    HAS_WAVEFORMS = False
+
+try:
+    from .common import BaseDriver, Quantity, get_coef
+except ImportError as e:
+    class BaseDriver:
+        def __init__(self, addr, timeout, **kw):
+            self.addr = addr
+            self.timeout = timeout
+
+
+    class Quantity(object):
+        def __init__(self, name: str, value=None, ch: int = 1, unit: str = ''):
+            self.name = name
+            self.default = dict(value=value, ch=ch, unit=unit)
+
+
+    def get_coef(*args):
+        return '', '', '', ''
+
+DEBUG_PRINT = False
+
+
+def print_debug(*args, **kwargs):
+    if DEBUG_PRINT:
+        print(*args, **kwargs)
+
+
+@atexit.register
+def global_system_exit():
+    """
+    关闭共享内存
+
+    :return:
+    """
+    SHARED_DEVICE_MEM.close()
+    SHARED_DEVICE_MEM.unlink()
+
+
+def solve_rpc_exception(func):
+    @wraps(func)
+    def wrapper(*args, **kwargs):
+        try:
+            return func(*args, **kwargs)
+        except xmlrpc.client.Fault as e:
+            pack = RPCFaultPack.from_fault(e)
+            print(f'************{args[0].__class__}: {args[0].addr}************')
+            print(f'*-*-*-*-*-*-*-远程函数报错: {pack.type}::{pack.name} -*-*-*-*-*-*-*-*')
+            print(pack.message)
+            print(f'*-*-*-*-*-*-*-远程函数报错: {pack.type}::{pack.name} -*-*-*-*-*-*-*-*')
+        except TimeoutError as e:
+            print(f'************{args[0].__class__}: {args[0].addr}************')
+            print(f'Driver报错: 无法连接到 {args[0].addr}\n请检查网络、设备开机状态、设备是否正在程序更新等')
+            print(f'{e}')
+            print('*****************************')
+            raise e
+        except socket.timeout as e:
+            print(f'************{args[0].__class__}: {args[0].addr}************')
+            print(f'Driver报错: 获取数据超时 {e}')
+            print('*****************************')
+    return wrapper
+
+
+class DataMode(str, Enum):
+    DIRECT = 'direct'
+    NORMAL = 'normal'
+    RDMA = 'rdma'
+
+
+class Driver(BaseDriver):
+    CHs = list(range(1, 25))
+    segment = ('ns', '111|112|113|114|115')
+
+    quants = [
+        Quantity('ReInit', value={}, ch=1),  # set, 设备重新初始化
+        Quantity('Instruction', value=None, ch=1),  # set   参数化波形指令队列配置
+        # 采集运行参数
+        Quantity('Shot', value=1024, ch=1),  # set,运行次数
+        Quantity('PointNumber', value=16384, unit='point'),  # set/get,AD采样点数
+        Quantity('TriggerDelay', value=0, ch=1, unit='s'),  # set/get,AD采样延时
+        Quantity('FrequencyList', value=[], ch=1, unit='Hz'),  # set/get,解调频率列表，list，单位Hz
+        Quantity('PhaseList', value=[], ch=1, unit='Hz'),  # set/get,解调频率列表，list，单位Hz
+        Quantity('Coefficient', value=None, ch=1),
+        Quantity('DemodulationParam', value=None, ch=1),
+        Quantity('CaptureMode'),
+        Quantity('StartCapture'),  # set,开启采集（执行前复位）
+        Quantity('TraceIQ', ch=1),  # get,获取原始时域数据
+        # 返回：array(shot, point)
+        Quantity('IQ', ch=1),  # get,获取解调后数据,默认复数返回
+        # 系统参数，宏定义修改，open时下发
+        # 复数返回：array(shot,frequency)
+        # 实数返回：array(IQ,shot,frequency)
+
+        # 任意波形发生器
+        Quantity('Waveform', value=np.array([]), ch=1),  # set/get,下发原始波形数据
+        Quantity('Delay', value=0, ch=1),  # set/get,播放延时
+        Quantity('KeepAmp', value=0),  # set, 电平是否维持在波形最后一个值, 0：波形播放完成后归0，1：保持波形最后一个值，2:保持波形第一个值
+        Quantity('Biasing', value=0, ch=1),  # set, 播放延迟
+        Quantity('LinSpace', value=[0, 30e-6, 1000], ch=1),  # set/get, np.linspace函数，用于生成timeline
+        Quantity('Output', value=True, ch=1),  # set/get,播放通道开关设置
+        Quantity('GenWave', value=None, ch=1),  # set/get, 设备接收waveform对象，根据waveform对象直接生成波形
+        # set/get, 设备接收IQ分离的waveform对象列表，根据waveform对象列表直接生成波形
+        Quantity('GenWaveIQ', value=None, ch=1),
+        Quantity('MultiGenWave', value={1: np.ndarray([])}),  # 多通道波形同时下发
+        Quantity('EnableWaveCache', value=False),  # 是否开启waveform缓存
+        Quantity('PushWaveCache'),  # 使waveform缓存中的波形数据生效
+        # 混频相关配置
+        Quantity('EnableDAMixer', value=False, ch=1),  # DA通道混频模式开关
+        Quantity('MixingWave',),  # 修改完混频相关参数后，运行混频器
+        Quantity('DAIQRate', value=1e9, ch=1),  # 基带信号采样率
+        Quantity('DALOFreq', value=100e6, ch=1),  # 中频信号频率
+        Quantity('DALOPhase', value=0, ch=1),  # 基带信号相位，弧度制
+        Quantity('DASideband', value='lower', ch=1),  # 混频后取的边带
+        Quantity('DAWindow', value=None, ch=1),
+        # 基带信号升采样率时所使用的窗函数，默认不使用任何窗，
+        # 可选：None、boxcar、triang、blackman、hamming、hann、bartlett、flattop、parzen、bohman、blackmanharris、nuttall、
+        # barthann、cosine、exponential、tukey、taylor
+
+        # 内触发
+        Quantity('GenerateTrig', value=1e7, unit='ns'),  # set/get,触发周期单位ns，触发数量=shot
+        Quantity('UpdateFirmware', value='', ch=1),  # qsync固件更新
+    ]
+
+    SystemParameter = {
+        'ExcLimit': -1,   # 设备traceback打印信息的回溯深度
+        'Warning': True,  # 开启后台警告上报
+        'MixMode': 2,  # Mix模式，1：第一奈奎斯特去； 2：第二奈奎斯特区
+        'PLLFreq': 100e6,  # 参考时钟频率, 单位为Hz
+        'RefClock': 'out',  # 参考时钟选择： ‘out’：外参考时钟；‘in’：内参考时钟
+        'ADrate': 4e9,  # AD采样率，单位Hz
+        'DArate': 6e9,  # DA采样率，单位Hz
+        'KeepAmp': 0,  # DA波形发射完毕后，保持最后一个值
+        'Delay': 0,  # 配置DA的原生Delay
+        'Shot':1024, # 配置实验重复次数
+        'SegmentSampling': [],  # 配置分段采样
+        # 'EnableDAMixer': False,  # 默认关闭DA混频器
+        # 'DASideband': 'lower',  # da混频默认使用下边带
+    }
+
+    def __init__(self, addr: str = '', timeout: float = 10.0, **kw):
+        super().__init__(addr, timeout, **kw)
+        self.fast_rpc = None
+        self.handle = None
+        self.model = 'NS_MCI'  # 默认为设备名字
+        self.srate = 6e9
+        self.data_mode = DataMode.NORMAL
+        self.device_online = True
+        self.has_start_capture = False
+        self.backend_version = (0, 0, 0)
+        print_debug(f'Driver: 实例化成功{addr}')
+
+    @solve_rpc_exception
+    def open(self, **kw):
+        """
+        输入IP打开设备，配置默认超时时间为5秒
+        打开设备时配置RFSoC采样时钟，采样时钟以参数定义
+        """
+        tran = Transport(use_builtin_types=True)
+        tran.accept_gzip_encoding = False
+        self.handle = xmlrpc.client.ServerProxy(f'http://{self.addr}:10801', transport=tran, allow_none=True,
+                                                use_builtin_types=True)
+        self.fast_rpc = FastRPC(self.addr, self.timeout)
+        debug = kw.get('debug', False)
+        if debug:
+            return
+        print_debug('Driver: RPC句柄建立成功')
+        # 检查目标设备是否在位设备
+        result = self.init_device(**kw)
+        if result:
+            # 在共享内存中注册本设备的ip
+            SHARED_DEVICE_MEM.ip = self.addr
+            self._show_system_status()
+            print(f'{self.addr}开启成功')
+        elif result is None:
+            ...
+        else:
+            print(self.handle.get_all_status())
+
+    def init_device(self, **kw):
+        try:
+            self._connect(self.addr, timeout=1).close()
+        except:
+            print(f'Driver: {self.addr}不在位，请检查网络重新open')
+            self.device_online = False
+            return
+        # 此时会连接rfsoc的指令接收tcp server
+        result = True
+        result &= self.handle.start_command(self.timeout)
+        print_debug('Driver: 后端rfs网络连接成功')
+
+        # 配置系统初始值
+        system_parameter = kw.get('system_parameter', {})
+        values = self.SystemParameter.copy()
+        values.update(system_parameter)
+        for name, value in values.items():
+            if value is not None:
+                res = self.fast_rpc.rpc_set(name, value, 0, False)
+                if isinstance(res, str):
+                    print("*********后台提示*********\n", res)
+        print_debug('Driver: 后端参数配置成功')
+
+        # 系统开启前必须进行过一次初始化
+        result &= self.__init_system()
+        print_debug('Driver: 后端RF配置成功')
+        result &= self.fast_rpc.rpc_set('Reset')
+        print_debug('Driver: 后端缓存配置成功')
+        return result
+
+    @solve_rpc_exception
+    def close(self, **kw):
+        """
+        关闭设备
+        """
+        ...
+
+    def write(self, name: str, value, **kw):
+        channel = kw.get('ch', 1)
+        if name in {'Coefficient'}:
+            data, f_list, numberOfPoints, phases = get_coef(value, 4e9)
+            judg_param = []
+            for qu in value['wlist']:
+                judg_param.append([qu['phi'], qu['threshold']])
+            self.set('JudgParam', judg_param, channel)
+            self.set('DemodulationParam', data, channel)
+        else:
+            return self.set(name, value, channel)
+
+    def read(self, name: str, **kw):
+        channel = kw.get('ch', 1)
+        result = self.get(name, channel)
+        return result
+
+    @solve_rpc_exception
+    def set(self, name, value=0, channel=1):
+        """
+        设置设备属性
+
+        :param name: 属性名
+                "Waveform"| "Amplitude" | "Offset"| "Output"
+        :param value: 属性值
+                "Waveform" --> 1d np.ndarray & -1 <= value <= 1
+                "Amplitude"| "Offset"| "Phase" --> float    dB | dB | °
+                “PRFNum”   采用内部PRF时，可以通过这个参数控制开启后prf的数量
+                "Output" --> bool
+        :param channel：通道号
+        """
+        if not self.device_online:
+            return
+        print_debug(f'Driver: set操作被调用{name}')
+        if name in {'UpdateFirmware'}:
+            self.update_firmware(value)
+            return
+        elif name in {'ReInit'}:
+            if not isinstance(value, dict):
+                value = {}
+            self.init_device(system_parameter=value)
+            return
+        if name in {'ProgramSrc', 'ProgramOUTSrc', 'ProgramINSrc'}:
+            source = [value.top_src, value.func_args, value.func_kwargs]
+            value = source
+        value = RPCValueParser.dump(value)
+        func = self.fast_rpc.rpc_set
+        name = 'GenWave' if name == 'Waveform' and value[0] == RPCValueParser.dump_tag_waveform else name
+
+        if self.has_start_capture and name == 'StartCapture':
+            return
+        self.has_start_capture = name == 'StartCapture'
+
+        if not func(name, value, channel):
+            ...
+            # raise xmlrpc.client.Fault(400, f'指令{name}执行失败, 请重新open板卡')
+
+    @solve_rpc_exception
+    def get(self, name, channel=1, value=0):
+        """
+        查询设备属性，获取数据
+
+        """
+        if not self.device_online:
+            return
+        print_debug(f'Driver: get操作被调用{name}')
+        self.has_start_capture = False
+        if name in {'Trace', 'TraceIQ'} and self.backend_version >= (2, 8, 0):
+            data_mode = self.get('DataMode', channel)
+            data_mode = DataMode(data_mode) if data_mode else DataMode.NORMAL
+            if data_mode == DataMode.DIRECT:
+                return self._get_trace_data(channel)
+
+        func = self.fast_rpc.rpc_get
+        tmp = func(name, channel, value)
+        tmp = RPCValueParser.load(tmp)
+        return tmp
+
+    def update_firmware(self, file_path, target='all', _id=0):
+        if not self.device_online:
+            return
+        import os
+        if not os.path.exists(file_path):
+            raise ValueError(f'文件路径: {file_path} 不存在')
+        with open(file_path, 'rb') as fp:
+            result = self.handle.update_rfs_firmware(fp.read(), target, _id)
+        if result:
+            print(f'设备{self.addr} 固件更新成功，3s后设备自动重启')
+        else:
+            print(f'设备{self.addr} 固件更新失败')
+        return result
+
+    def __init_system(self):
+        version = self.get('Version')
+        if isinstance(version, dict):
+            for key in version:
+                key: str
+                if key.endswith('version'):
+                    version = version[key][1:].split('-')[0]
+                    version = tuple(int(v) for v in version.split('.'))
+                    break
+                else:
+                    version = (0, 0, 0)
+        self.backend_version = version
+        if version >= (1, 2, 7):
+            self.set('InitSystem')
+            return True
+        else:
+            return self.__exec_command('初始化')
+
+    def __exec_command(self, button_name: str,
+                       need_feedback=True, file_name=None, check_feedback=True,
+                       callback=lambda *args: True, wait: int = 0):
+        flag = self.handle.execute_command(button_name, need_feedback, file_name, check_feedback)
+        if not flag:
+            print(f'设备{self.addr}-指令{button_name}执行失败')
+        return flag
+
+    def _get_trace_data(self, channel):
+        res = self.get('ShmPath', channel)
+        if not res:
+            return
+        _path, chnl, chnl_num, shots, points = res
+        shm = shared_memory.SharedMemory(name=_path, create=False)
+        data = np.frombuffer(shm.buf, dtype=np.int16).reshape((chnl_num, shots, points))
+        return data[chnl], data
+
+    def _connect(self, addr=None, port=5001, timeout=None):
+        """
+        获取到指定ip的tcp连接
+
+        :param addr:
+        :param port:
+        :return:
+        """
+        timeout = self.timeout if timeout is None else timeout
+        addr = self.addr if addr is None else addr
+        sock = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
+        sock.settimeout(timeout+1)
+        sock.connect((addr, port))
+        return sock
+
+    def _show_system_status(self):
+        keys = ['device_type', 'backend_version', 'ad_num', 'da_num', 'cpu_temp', 'memory_use']
+        if self.backend_version >= (2, 1, 0):
+            status: dict = self.get("Status")
+        else:
+            status: dict = self.handle.get_all_status(False)
+        _string = [f'*********设备{self.addr}开启成功*********']
+        for key in keys:
+            _string.append(f'{key}: {status.get(key, "nan")}')
+        print('\n'.join(_string))
+        if self.backend_version >= (2, 0, 1):
+            print('available chnl: ')
+            print(self.get('ChnlInfo'))
+
+
+class RPCFaultPack:
+    split_flag = '|+*+-*-+*+|*-*|'
+
+    class FaultType(str, enum.Enum):
+        EXCEPTION = 'exception'
+        WARNING = 'warning'
+        INFO = 'info'
+
+    def __init__(self, _type: str, name: str, message: str):
+        self.type: RPCFaultPack.FaultType = self.FaultType(_type)
+        self.name = name
+        self.message = message
+
+    @classmethod
+    def from_string(cls, string: str):
+        self = cls(cls.FaultType.INFO, '', '')
+        self.string = string
+        return self
+
+    @classmethod
+    def from_fault(cls, fault: xmlrpc.client.Fault):
+        return cls.from_string(fault.faultString)
+
+    @classmethod
+    def from_exception(cls, exc: Exception, limit=0):
+        message = traceback.format_exception(exc, limit=limit)
+        return cls(cls.FaultType.EXCEPTION, exc.__class__.__name__, ''.join(message))
+
+    @classmethod
+    def from_warning(cls, warns: "list[warnings.WarningMessage]"):
+        warn = [str(warn.category.__name__) for warn in warns]
+        warn = str(warn[0]) if len(warn) == 1 else str(warn)
+        messages = [str(warn.message) for warn in warns]
+        return cls(cls.FaultType.WARNING, warn, '\n'.join(messages))
+
+    @property
+    def string(self):
+        return self.split_flag.join([self.type, self.name, self.message])
+
+    @string.setter
+    def string(self, value: str):
+        _p = value.split(self.split_flag)
+        if len(_p) == 1:
+            _p = (self.FaultType.INFO, 'Default', _p[0])
+        elif len(_p) == 2:
+            _p = (self.FaultType.INFO, _p[0], _p[1])
+        else:
+            _p = _p[:3]
+        self.type, self.name, self.message = _p
+
+    def __str__(self):
+        return self.string
+
+
+class FastRPC:
+    def __init__(self, address, timeout):
+        self.addr = address
+        self.timeout = timeout
+
+    def connect(self):
+        sock = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
+        sock.connect((self.addr, 10800))
+        sock.settimeout(self.timeout+1)
+        return sock
+
+    def rpc_set(self, *param):
+        sock = self.connect()
+        a = pickle.dumps(param)
+        head = struct.pack('=IIII', *[0x5F5F5F5F, 0x32000001, 0, 16 + len(a)])
+        sock.sendall(head)
+        sock.sendall(a)
+        head = struct.unpack("=IIIII", sock.recv(20))
+        # print(head)
+        data = sock.recv(head[3] - 20)
+        data = pickle.loads(data)
+        if head[4]:
+            print(data)
+            return False
+        sock.close()
+        return data
+
+    def rpc_get(self, *param):
+        sock = self.connect()
+        a = pickle.dumps(param)
+        head = struct.pack('=IIII', *[0x5F5F5F5F, 0x32000002, 0, 16 + len(a)])
+        sock.sendall(head)
+        sock.sendall(a)
+        head = struct.unpack("=IIIII", sock.recv(20))
+        length = head[3] - 20
+        data = []
+        while length > 0:
+            _data = sock.recv(length)
+            length -= len(_data)
+            data.append(_data)
+        data = pickle.loads(b''.join(data))
+        if head[4]:
+            raise xmlrpc.client.Fault(400, data)
+        sock.close()
+        return data
+
+    def debug_param(self, *param):
+        sock = self.connect()
+        a = pickle.dumps(param)
+        head = struct.pack('=IIII', *[0x5F5F5F5F, 0x32000003, 0, 16 + len(a)])
+        sock.sendall(head)
+        sock.sendall(a)
+        head = struct.unpack("=IIIII", sock.recv(20))
+        length = head[3] - 20
+        data = []
+        while length > 0:
+            _data = sock.recv(length)
+            length -= len(_data)
+            data.append(_data)
+        data = pickle.loads(b''.join(data))
+        if head[4]:
+            print(data)
+            return False
+        sock.close()
+        return data
+
+
+class RPCValueParser:
+    """
+    rpc调用格式化工具集
+
+    """
+    dump_tag_waveform = 'waveform.Waveforms'
+    dump_tag_ndarray = 'numpy.ndarray'
+    dump_tag_complex = 'complex'
+
+    @staticmethod
+    def dump(value):
+        if isinstance(value, np.ndarray):
+            value = [RPCValueParser.dump_tag_ndarray, value.tobytes(), str(value.dtype), value.shape]
+        elif isinstance(value, float):
+            value = float(value)
+        elif isinstance(value, int):
+            value = int(value)
+        elif isinstance(value, np.uint):
+            value = int(value)
+        elif isinstance(value, complex):
+            value = [RPCValueParser.dump_tag_complex, value.real, value.imag]
+        elif HAS_WAVEFORMS and isinstance(value, waveforms.Waveform):
+            value = [RPCValueParser.dump_tag_waveform, pickle.dumps(value)]
+        elif isinstance(value, (list, tuple)):
+            value = [RPCValueParser.dump(_v) for _v in value]
+
+        return value
+
+    @staticmethod
+    def load(value):
+        if isinstance(value, list) and len(value) >= 2:
+            if value[0] == RPCValueParser.dump_tag_ndarray:
+                data = np.frombuffer(value[1], dtype=value[2])
+                value = data.reshape(value[3])
+            # elif: value[0] == 'numpy.float'
+            elif value[0] == RPCValueParser.dump_tag_waveform:
+                value = pickle.loads(value[1])
+            elif value[0] == RPCValueParser.dump_tag_complex:
+                value = complex(value[1], value[2])
+            else:
+                value = [RPCValueParser.load(_v) for _v in value]
+        elif isinstance(value, (list, tuple)):
+            value = [RPCValueParser.load(_v) for _v in value]
+        return value
+
+class InfoSharedList:
+    memory_name = 'NS_DEVICE_MEMORY'
+    memory_size = 1024 ** 2
+    head_size = 32  # memory中前head_size的长度为头部预留信息
+    def __init__(self):
+        try:
+            self._memory = shared_memory.SharedMemory(name=self.memory_name, create=True, size=self.memory_size)
+        except FileExistsError:
+            self._memory = shared_memory.SharedMemory(name=self.memory_name, create=False, size=self.memory_size)
+
+    def clear_ip(self):
+        _exit_pkl = pickle.dumps(self.ip)
+        self._memory.buf[self.head_size:len(_exit_pkl) + self.head_size] = b'\x00' * len(_exit_pkl)
+
+    @property
+    def ip(self):
+        try:
+            return pickle.loads(self._memory.buf[self.head_size:])
+        except pickle.UnpicklingError:
+            return []
+
+    @ip.setter
+    def ip(self, value):
+        ips = self.ip
+        ips.append(value)
+        ips = list(set(ips))
+        _pkl = pickle.dumps(ips)
+        self._memory.buf[self.head_size:len(_pkl) + self.head_size] = _pkl
+
+    def close(self):
+        self._memory.close()
+
+    def unlink(self):
+        try:
+            self._memory.unlink()
+        except FileNotFoundError as e:
+            pass
+
+
+SHARED_DEVICE_MEM = InfoSharedList()