emsutil 0.2.2__tar.gz → 0.4.0__tar.gz

{emsutil-0.2.2 → emsutil-0.4.0}/.bumpversion.toml

@@ -1,5 +1,5 @@
 [tool.bumpversion]
-current_version = "0.2.2"
+current_version = "0.4.0"
 parse = "(?P<major>\\d+)\\.(?P<minor>\\d+)\\.(?P<patch>\\d+)"
 serialize = ["{major}.{minor}.{patch}"]
 search = "{current_version}"

{emsutil-0.2.2 → emsutil-0.4.0}/PKG-INFO

@@ -1,12 +1,14 @@
 Metadata-Version: 2.4
 Name: emsutil
-Version: 0.2.2
+Version: 0.4.0
 Summary: Common utilities for Emerge projects EMerge, Optycal and Heavi
 Project-URL: Homepage, https://github.com/FennisRobert/emsutil
 Project-URL: Issues, https://github.com/FennisRobert/emsutil/issues
 Requires-Python: >=3.10
 Requires-Dist: loguru>=0.7.3
 Requires-Dist: matplotlib>=3.8.0
+Requires-Dist: msgpack-numpy>=0.4.8
+Requires-Dist: msgpack>=1.1.2
 Requires-Dist: numpy<2.3,>=1.24
 Requires-Dist: pyvista>=0.46.4
 Requires-Dist: scipy>=1.14.0

emsutil-0.4.0/emerge-install-mumps

@@ -0,0 +1,32 @@
+#!/usr/bin/env bash
+set -e
+
+# Resolve Python environment
+PYTHON="$(command -v python)"
+
+if [ -z "$PYTHON" ]; then
+  echo "ERROR: No active Python environment found."
+  exit 1
+fi
+
+echo "Using Python:"
+echo "  $PYTHON"
+
+# Location of bundled installer
+INSTALLER_ROOT="$HOME/.emerge/installer"
+
+if [ ! -f "$INSTALLER_ROOT/install.sh" ]; then
+  echo "ERROR: MUMPS installer not found."
+  echo "Expected:"
+  echo "  $INSTALLER_ROOT/install.sh"
+  exit 1
+fi
+
+# Run installer
+bash "$INSTALLER_ROOT/install.sh"
+
+# Install Python package into *this* environment
+"$PYTHON" -m pip install "$HOME/.emerge/mumps_py"
+
+echo
+echo "MUMPS successfully installed for this Python environment."

{emsutil-0.2.2 → emsutil-0.4.0}/pyproject.toml

@@ -4,13 +4,15 @@ allow-direct-references = true
 
 [project]
 name = "emsutil"
-version = "0.2.2"
+version = "0.4.0"
 description = "Common utilities for Emerge projects EMerge, Optycal and Heavi"
 readme = "README.md"
 requires-python = ">=3.10"
 dependencies = [
     "loguru>=0.7.3",
     "matplotlib>=3.8.0",
+    "msgpack>=1.1.2",
+    "msgpack-numpy>=0.4.8",
     "numpy>=1.24, <2.3",
     "pyvista>=0.46.4",
     "scipy>=1.14.0",

{emsutil-0.2.2 → emsutil-0.4.0}/src/emsutil/__init__.py

@@ -4,4 +4,5 @@ from .emdata import FarFieldComponent, EHFieldFF, EHField
 from .material import Material, MatProperty, FreqCoordDependent, FreqDependent, CoordDependent
 from .lib import EISO, EOMNI, EPS0, AIR, PEC, COPPER
 from .pyvista.display_settings import EMergeTheme
-from . import themes
+from . import themes
+from .file import Saveable, save_object, load_object

{emsutil-0.2.2 → emsutil-0.4.0}/src/emsutil/emdata.py

@@ -3,10 +3,13 @@ from dataclasses import dataclass, field
 from typing import Literal, Callable
 from .const import Z0, EPS0
 from .lib import EISO, EOMNI
+from .file import Saveable
 
 @dataclass
-class FarFieldComponent:
+class FarFieldComponent(Saveable):
     F: np.ndarray
+    _th: np.ndarray
+    _ph: np.ndarray
 
     @property
     def x(self) -> np.ndarray:
@@ -22,16 +25,16 @@ class FarFieldComponent:
     
     @property
     def theta(self) -> np.ndarray:
-        thx = np.cos(self.theta)*np.cos(self.phi)
-        thy = np.cos(self.theta)*np.sin(self.phi)
-        thz = -np.sin(self.theta)
+        thx = np.cos(self._th)*np.cos(self._ph)
+        thy = np.cos(self._th)*np.sin(self._ph)
+        thz = -np.sin(self._th)
         return thx*self.F[0,:] + thy*self.F[1,:] + thz*self.F[2,:]
     
     @property
     def phi(self) -> np.ndarray:
-        phx = -np.sin(self.phi)
-        phy = np.cos(self.phi)
-        phz = np.zeros_like(self.theta)
+        phx = -np.sin(self._ph)
+        phy = np.cos(self._ph)
+        phz = np.zeros_like(self._th)
         return phx*self.F[0,:] + phy*self.F[1,:] + phz*self.F[2,:]
     
     @property
@@ -55,7 +58,7 @@ class FarFieldComponent:
 
 
 @dataclass
-class FieldPlotData:
+class FieldPlotData(Saveable):
     x: np.ndarray
     y: np.ndarray
     z: np.ndarray
@@ -80,6 +83,16 @@ class FieldPlotData:
         """
         return (self.x, self.y, self.z, self.F)
 
+    @property
+    def xyzvec(self) -> tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]:
+        """Returns the X, Y, Z, F arrays.
+
+        Returns:
+            tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray]: The X, Y, Z, F arrays
+        """
+        return (self.x, self.y, self.z, self.vx, self.vy, self.vz)
+
+
     @property
     def xyzftri(self) -> tuple[np.ndarray, np.ndarray, np.ndarray, np.ndarray, np.ndarray]:
         """Returns the X, Y, Z, F + triangulation arrays.
@@ -115,13 +128,14 @@ class FieldPlotData:
             yield self.F
             
 @dataclass
-class EHFieldFF:
+class EHFieldFF(Saveable):
     _E: np.ndarray
     _H: np.ndarray
     theta: np.ndarray
     phi: np.ndarray
     Ptot: float
     ang: np.ndarray | None = field(default=None)
+    freq: float | None = field(default=None)
     
     def total_radiated_power_integral(
         self,
@@ -183,11 +197,11 @@ class EHFieldFF:
     
     @property
     def E(self) -> np.ndarray:
-        return FarFieldComponent(self._E)
+        return FarFieldComponent(self._E, self.theta, self.phi)
     
     @property
     def H(self) -> np.ndarray:
-        return FarFieldComponent(self._H)
+        return FarFieldComponent(self._H, self.theta, self.phi)
     
     @property
     def Ex(self) -> np.ndarray:
@@ -244,16 +258,16 @@ class EHFieldFF:
     @property
     def gain(self, kind: Literal['iso','omni'] = 'iso') -> FarFieldComponent:
         if kind=='iso':
-            return FarFieldComponent(self._E/EISO)
+            return FarFieldComponent(self._E/EISO, self.theta, self.phi)
         else:
-            return FarFieldComponent(self._E/EOMNI)
+            return FarFieldComponent(self._E/EOMNI, self.theta, self.phi)
         
     @property
     def dir(self, kind: Literal['iso','omni'] = 'iso') -> FarFieldComponent:
         if kind=='iso':
-            return FarFieldComponent(self._E/(EISO*(self.Ptot)**0.5))
+            return FarFieldComponent(self._E/(EISO*(self.Ptot)**0.5), self.theta, self.phi)
         else:
-            return FarFieldComponent(self._E/(EOMNI*(self.Ptot)**0.5))
+            return FarFieldComponent(self._E/(EOMNI*(self.Ptot)**0.5), self.theta, self.phi)
     
     @property
     def normE(self) -> np.ndarray:
@@ -306,7 +320,7 @@ class EHFieldFF:
         return FieldPlotData(x=xs, y=ys, z=zs, F=F)
 
 @dataclass
-class EHField:
+class EHField(Saveable):
     _E: np.ndarray
     _H: np.ndarray
     x: np.ndarray
@@ -626,6 +640,7 @@ class EHField:
         """
         Fx, Fy, Fz = getattr(self, field)
         
+        
         if metric=='real':
             Fx, Fy, Fz = Fx.real, Fy.real, Fz.real
         elif metric=='imag':

emsutil-0.4.0/src/emsutil/file.py

@@ -0,0 +1,192 @@
+from __future__ import annotations
+import numpy as np
+from scipy.sparse import csc_matrix, csr_matrix, isspmatrix
+from numbers import Number
+import msgpack
+import msgpack_numpy as m
+
+m.patch()
+
+
+class Saveable:
+    save_fields: list[str] = None
+    skip_fields: list[str] = None
+    _registry = dict()
+    
+    def __init_subclass__(cls, **kwargs):
+        super().__init_subclass__(**kwargs)
+        cls._registry[cls.__name__] = cls
+        
+    @classmethod
+    def is_saveable(cls, obj) -> bool:
+        if not isinstance(obj, dict):
+            return False
+        if '__SAVEABLE__' in obj:
+            return True
+    
+    def pack(self) -> dict:
+        if self.skip_fields is None:
+            skip = []
+        else:
+            skip = self.skip_fields
+            
+        if self.save_fields is None:
+            data = dict()
+            data['__CLASS__'] = self.__class__.__name__
+            data['__SAVEABLE__'] = True
+            data['__FIELDS__'] = {key: value for key, value in self.__dict__.items()
+                                  if '__' not in key and key not in skip}
+            return data
+        else:
+            data = dict()
+            data['__CLASS__'] = self.__class__.__name__
+            data['__SAVEABLE__'] = True
+            data['__FIELDS__'] = {key: value for key, value in self.__dict__.items() if key in self.save_fields and key not in skip}
+            return data
+    
+    @staticmethod 
+    def unpack(data: dict, collection: dict) -> Saveable:
+        cls = Saveable._registry[data['__CLASS__']]
+        fields = data['__FIELDS__']
+        obj = cls.__new__(cls)
+        for field, value in fields.items():
+            setattr(obj, field, _unpack_data(value, collection))
+        return obj
+    
+def _is_set(data) -> bool:
+    if not isinstance(data, dict):
+        return False
+    if '__ISSET__' in data:
+        return True
+    
+def _pack_set(data) -> dict:
+    return {
+        '__ISSET__': True,
+        '__DATA__': list(data)
+    }
+
+def _is_type(data) -> bool:
+    if not isinstance(data, dict):
+        return False
+    if '__CLASS_TYPE__' in data:
+        return True
+
+def _pack_type(data) -> dict:
+    return {
+        '__CLASS_TYPE__': True,
+        '__CLASS_NAME__': f"{data.__module__}.{data.__qualname__}"
+    }
+
+def _is_saved(data) -> bool:
+    if not isinstance(data, dict):
+        return False
+    if '__SAVED_OBJ__' in data:
+        return True
+
+def _pack_saved(data) -> dict:
+    return {
+        '__SAVED_OBJ__': True,
+        '__OBJ_ID__': id(data)
+    }
+
+def _unpack_saved(data: dict) -> int:
+    return data['__OBJ_ID__']
+
+def _unpack_type(data: dict):
+    class_name = data['__CLASS_NAME__']
+    module_name, _, cls_name = class_name.rpartition('.')
+    module = __import__(module_name, fromlist=[cls_name])
+    cls = getattr(module, cls_name)
+    return cls
+
+def _unpack_set(data: dict) -> set:
+    return set(data['__DATA__'])
+
+def _pack_object(obj, _obj_collector: dict = None) -> dict:
+    if obj is None:
+        return obj
+    if isinstance(obj, (str, Number, np.ndarray, bool)):
+        return obj
+    elif isinstance(obj, set):
+        return _pack_set(obj)
+    elif isspmatrix(obj):
+        return obj
+    elif isinstance(obj, list):
+        return [_pack_object(item, _obj_collector) for item in obj]
+    elif isinstance(obj, tuple):
+        return tuple(_pack_object(item, _obj_collector) for item in obj)
+    elif isinstance(obj, dict):
+        return {key: _pack_object(value, _obj_collector) for key, value in obj.items()}
+    elif isinstance(obj, Saveable):
+        if id(obj) in _obj_collector:
+            return _obj_collector[id(obj)]['reference']
+        data = _pack_saved(obj)
+        _obj_collector[data['__OBJ_ID__']] = {'object': _pack_object(obj.pack(), _obj_collector), 'reference': data}
+        return data
+    elif isinstance(obj, type):
+        return _pack_type(obj)
+    else:
+        print(f'Trying to unpack {obj} of type {type(obj)}')
+        raise TypeError(f"Cannot pack object of type {type(obj)}")
+    
+def _unpack_data(data, _obj_collection):
+    if _is_saved(data):
+        obj_id = _unpack_saved(data)
+        saved_data = _obj_collection[obj_id]['object']
+        return _unpack_data(saved_data, _obj_collection)
+    if Saveable.is_saveable(data):
+        return Saveable.unpack(data, _obj_collection)
+    elif data is None:
+        return None
+    elif _is_type(data):
+        return _unpack_type(data)
+    elif _is_set(data):
+        return _unpack_set(data)
+    elif isinstance(data, np.ndarray):
+        return data.copy()
+    elif isinstance(data, (str, Number, np.ndarray, bool)):
+        return data
+    elif isinstance(data, list):
+        return [_unpack_data(item, _obj_collection) for item in data]
+    elif isinstance(data, tuple):
+        return tuple(_unpack_data(item, _obj_collection) for item in data)
+    elif isinstance(data, dict):
+        return {key: _unpack_data(value, _obj_collection) for key, value in data.items()}
+    elif isspmatrix(data):
+        return data
+    else:
+        raise TypeError(f"Cannot unpack data of type {type(data)}")
+
+def save_object(filename: str, data: dict):
+    """Saves a class to a file using msgpack serialization.
+
+    Args:
+        filename (str): The filename for the final object
+        data (dict): The data to store
+    """
+    obj_collector = dict()
+    dataset = _pack_object(data, obj_collector)
+    save_dict = {
+        'data': dataset,
+        'objects':  obj_collector
+    }
+    
+    with open(filename, 'wb') as f:
+        f.write(msgpack.packb(save_dict, use_bin_type=True))
+
+def load_object(filename: str) -> dict:
+    """Loads a class from a file using msgpack serialization.
+
+    Args:
+        filename (str): The filename of the dataset to load
+
+    Returns:
+        dict: The dict() object with the loaded data.
+    """
+    with open(filename, 'rb') as f:
+        loaded = msgpack.unpackb(f.read(), raw=False, strict_map_key=False, use_list=False)
+    
+    data = loaded['data']
+    obj_collection = loaded['objects']
+    return _unpack_data(data, obj_collection)
+

emsutil-0.4.0/src/emsutil/inexport/ffdata.py

@@ -0,0 +1,304 @@
+import numpy as np
+from .strutil import arry_to_line, matrix_to_lines, arry_to_fwl
+from ..emdata import EHFieldFF
+from dataclasses import dataclass
+from typing import List
+from loguru import logger
+import os
+
+@dataclass
+class Column:
+    name: str
+    values: np.ndarray
+
+    @staticmethod
+    def from_array(name: str, data: np.ndarray) -> List["Column"]:
+        """
+        Returns:
+        - [Column] if data is real
+        - [Column(real), Column(imag)] if data is complex
+        """
+        data = np.asarray(data)
+
+        if np.iscomplexobj(data):
+            return [
+                Column(f"{name}[re]", data.real),
+                Column(f"{name}[im]", data.imag),
+            ]
+        else:
+            return [Column(name, data)]
+
+
+def to_lines(
+    columns: List[Column],
+    header_prefix: str = "$ ",
+    separator: str = "  ",
+    fmt: str = "{:g}",
+) -> List[str]:
+    """
+    Convert columns to equally spaced text lines.
+    """
+
+    # Convert all values to strings first
+    str_cols = [[fmt.format(v) for v in col.values] for col in columns]
+    headers = [col.name for col in columns]
+    headers[0] = header_prefix + headers[0]
+    for i in range(len(headers)-1):
+        headers[i] = headers[i] + ';'
+    # Compute column widths (including separator!)
+    widths = []
+    for h, col in zip(headers, str_cols):
+        widths.append(max(len(h), max(len(v) for v in col)))
+
+    def format_row(items, sep: str):
+        return sep.join(item.ljust(w) for item, w in zip(items, widths))
+
+    lines = []
+
+    # Header
+    lines.append(format_row(headers, separator))
+
+    # Data rows
+    nrows = len(columns[0].values)
+    for i in range(nrows):
+        row = [col[i] for col in str_cols]
+        lines.append(format_row(row, separator))
+
+    return lines
+
+class _fieldGetter:
+    
+    def __init__(self, fields: list[EHFieldFF]):
+        self.names: list[str] = []
+        self.data: list = fields
+    
+    @property
+    def name(self) -> str:
+        return '.'.join([n.capitalize() for n in self.names])
+        
+    def __getattr__(self, name: str) -> EHFieldFF:
+        self.names.append(name)
+        self.data = [getattr(item, name) for item in self.data]
+        return self
+    
+class FarFieldExporter:
+    """A class to export far-field data to a text file.
+    
+    """
+    def __init__(self, filename: str, fields: list[EHFieldFF], precision: int = 4, todeg: bool = True):
+        self.filename: str = filename
+        self.fields: list[EHFieldFF] = fields
+        self.precision: int = precision
+        self.columns: list[_fieldGetter] = []
+        self.todeg: bool = todeg
+    
+    def addcol(self) -> EHFieldFF:
+        """Returns a field getter to specify which field component to export.
+        
+        Example:
+        >>> exporter = FarFieldExporter('data.txt', field_data)
+        >>> exporter.addcol().Ex
+        >>> exporter.addcol().gain.theta
+        >>> exporter.write()
+
+        Returns:
+            EHFieldFF(_fieldGetter): Field getter to specify field component.
+        """
+        getter = _fieldGetter(self.fields)
+        self.columns.append(getter)
+        return getter
+    
+    def write(self) -> None:
+        """Saves the far-field data to a text file.
+        """
+        base = self.fields[0]
+        
+        thetas = base.theta
+        phis = base.phi
+        
+        theta_lin = thetas[0,:].flatten()
+        phi_lin = phis[:,0].flatten()
+        p = self.precision
+        mult = 1.0
+        unit = 'rad'
+        if self.todeg:
+            mult = 180.0 / np.pi
+            theta_lin = theta_lin * mult
+            phi_lin = phi_lin * mult
+            thetas = thetas * mult
+            phis = phis * mult
+            unit = 'deg'
+        
+        datalines = []
+        datalines.append(f'% Theta({unit})')
+        datalines.append(' '.join([f'{x:.{p}g}' for x in theta_lin]))
+        datalines.append(f'% Phi({unit})')
+        datalines.append(' '.join([f'{x:.{p}g}' for x in phi_lin]))
+        datalines.append('% Frequencies(Hz)')
+        datalines.append(' '.join([f'{field.freq:.{p}g}' for field in self.fields]))
+        datalines.append('')
+        thcol = Column(f'Theta({unit})', thetas.flatten())
+        phcol = Column(f'Phi({unit})', phis.flatten())
+        
+        for i, field in enumerate(self.fields):
+            freq = field.freq
+            cols = []
+            datalines.append(f'# {freq:.{p}g} Hz')
+            for col in self.columns:
+                cols.extend(Column.from_array(col.name, col.data[i].flatten()))
+            
+            datalines.extend(to_lines([thcol, phcol] + cols, fmt=f"{{:.{p}g}}"))
+            datalines.append('')
+        
+        #replace any extension of filename with.emff
+        filename = os.path.splitext(self.filename)[0] + '.emff'
+        self.filename = filename
+        with open(self.filename, 'w') as f:
+            f.write('\n'.join(datalines))
+        logger.info(f'Wrote far-field data to {self.filename}')
+        
+def import_farfield(filename: str, degrees: bool = True) -> dict[str, np.ndarray]:
+    """Reads emerge farfiel datafiles created by FarFieldExporter.
+    They can be recognized by the .emff extension.
+
+    Args:
+        filename (str): The path to the far-field data file.
+        degrees (bool): If True, theta and phi angles are converted to degrees.
+
+    Returns:
+        dict[str, np.ndarray]: The imported far-field data structured as:
+            {
+                "theta": np.ndarray,  # Theta grid (nTheta, nPhi)
+                "phi": np.ndarray,    # Phi grid (nTheta, nPhi)
+                "freq": np.ndarray,   # Frequencies (nF,)
+                "Ex": np.ndarray,     # Electric field component Ex (nF, nTheta, nPhi)
+                "Ey": np.ndarray,     # Electric field component Ey (nF, nTheta, nPhi)
+                ...
+    """
+    with open(filename, "r") as f:
+        lines = f.readlines()
+
+    # ---------------- helpers ----------------
+    def parse_vector(header):
+        i = lines.index(header) + 1
+        return np.array(lines[i].split(), float)
+
+    def strip_unit(name):
+        if "(" in name and ")" in name:
+            base, unit = name.split("(")
+            return base.strip(), unit.rstrip(")")
+        return name.strip(), None
+
+    def to_rad(values, unit):
+        if unit == "deg":
+            return np.deg2rad(values)
+        return values
+
+    # ---------------- global grids ----------------
+    theta_name, theta_unit = strip_unit("% Theta(deg)")
+    phi_name,   phi_unit   = strip_unit("% Phi(deg)")
+
+    # Find actual headers (deg or rad)
+    theta_header = next(l for l in lines if l.startswith("% Theta"))
+    phi_header   = next(l for l in lines if l.startswith("% Phi"))
+
+    theta_vals = parse_vector(theta_header)
+    phi_vals   = parse_vector(phi_header)
+    freqs      = parse_vector("% Frequencies(Hz)\n")
+
+    theta_unit = theta_header.split("(")[1].split(")")[0]
+    phi_unit   = phi_header.split("(")[1].split(")")[0]
+
+    theta_vals = to_rad(theta_vals, theta_unit)
+    phi_vals   = to_rad(phi_vals, phi_unit)
+
+    nTheta = len(theta_vals)
+    nPhi   = len(phi_vals)
+    nF     = len(freqs)
+
+    theta_grid, phi_grid = np.meshgrid(theta_vals, phi_vals, indexing="ij")
+
+    result = {
+        "theta": theta_grid,
+        "phi": phi_grid,
+        "freq": freqs,
+    }
+
+    # ---------------- frequency blocks ----------------
+    field_buffers = {}
+    column_map = {}
+    header_info = {}
+    fidx = -1
+    header_cols = None
+    data_rows = []
+
+    for line in lines:
+        line = line.strip()
+        if not line:
+            continue
+
+        # Frequency block
+        if line.startswith("#"):
+            fidx += 1
+            data_rows = []
+            continue
+
+        # Header line
+        if line.startswith("$"):
+            header_cols = []
+            column_map.clear()
+            header_info.clear()
+
+            raw_cols = line[1:].split(";")
+            for i, raw in enumerate(raw_cols):
+                name, unit = strip_unit(raw)
+                header_cols.append(name)
+                header_info[name] = unit
+
+                if "[" in name and "]" in name:
+                    base, comp = name.split("[")
+                    comp = comp.rstrip("]")
+                    column_map.setdefault(base, {})[comp] = i
+                else:
+                    column_map[name] = i
+            continue
+
+        if line.startswith("%"):
+            continue
+
+        parts = line.split()
+        if header_cols is None or len(parts) < len(header_cols):
+            continue
+
+        data_rows.append([float(p) for p in parts])
+
+        if len(data_rows) == nTheta * nPhi:
+            data = np.array(data_rows)
+
+            for field, comps in column_map.items():
+                unit = header_info.get(field)
+
+                if isinstance(comps, dict):  # complex
+                    re = data[:, comps["re"]]
+                    im = data[:, comps["im"]]
+                    arr = re + 1j * im
+                else:  # real
+                    arr = data[:, comps]
+
+                if unit in ("deg", "rad"):
+                    arr = to_rad(arr, unit)
+
+                arr = arr.reshape(nTheta, nPhi)
+                field_buffers.setdefault(field, []).append(arr)
+
+    
+    # Stack frequencies
+    for field, blocks in field_buffers.items():
+        result[field] = np.stack(blocks, axis=0)
+    
+    if degrees:
+        result["theta"] = np.rad2deg(result["theta"])
+        result["phi"] = np.rad2deg(result["phi"])
+        result["Theta"] = np.rad2deg(result["Theta"])
+        result["Phi"] = np.rad2deg(result["Phi"])
+    return result

{emsutil-0.2.2 → emsutil-0.4.0}/src/emsutil/material.py

@@ -19,6 +19,7 @@ from __future__ import annotations
 import numpy as np
 from typing import Callable
 import inspect
+from .file import Saveable
 
 C0 = 299792458
 
@@ -45,15 +46,16 @@ def _to_mat(value: float | complex | int | np.ndarray) -> np.ndarray:
     else:
         return ValueError(f'Trying to parse {value} as a material property tensor but it cant be identified as scalar, vector or matrix')
 
-class MatProperty:
+class MatProperty(Saveable):
     _freq_dependent: bool = False
     _coord_dependent: bool = False
     _pickle_exclude = {"_func","_fmax"}
+    skip_fields = ("_func","_fmax")
     """The MatProperty class is an interface for EMerge to deal with frequency and coordinate dependent material properties
     """
     
     def __init__(self, value: float | complex | int | np.ndarray):
-        self.value: np.ndarray = _to_mat(value)
+        self._value: np.ndarray = _to_mat(value)
         
         self._apply_to: np.ndarray = np.array([], dtype=np.int64)
         self._x: np.ndarray = np.array([], dtype=np.float64)
@@ -69,11 +71,15 @@ class MatProperty:
         self._z = np.concatenate([self._z, z])
         
     def __call__(self, f: float, data: np.ndarray) -> np.ndarray:
-        data[:,:,self._apply_to] = np.repeat(self.value[:,:,np.newaxis], self._apply_to.shape[0], axis=2)
+        data[:,:,self._apply_to] = np.repeat(self._value[:,:,np.newaxis], self._apply_to.shape[0], axis=2)
         return data
 
+    @property
+    def value(self) -> float:
+        return self._value[0,0]
+    
     def scalar(self, f: float):
-        return self.value[0,0]
+        return self._value[0,0]
 
     def reset(self) -> None:
         self._apply_to: np.ndarray = np.array([], dtype=np.int64)
@@ -93,10 +99,10 @@ class MatProperty:
         for k in self._pickle_exclude:
             setattr(self, k, None)
     
-class FreqDependent(MatProperty):
+class FreqDependent(MatProperty, Saveable):
     _freq_dependent: bool = True
     _coord_dependent: bool = False
-
+    skip_fields = ("_func","_fmax")
     def __init__(self, 
                  scalar: Callable | None = None,
                  vector: Callable | None = None,
@@ -115,8 +121,6 @@ class FreqDependent(MatProperty):
             vector (Callable | None, optional): The diagonal rank-2 tensor function returning a (3,) array. Defaults to None.
             matrix (Callable | None, optional): The rank-2 tensor function returning a (3,3) array. Defaults to None.
 
-        Returns:
-            _type_: _description_
         """
         if scalar is not None:
             def _func(f: float) -> np.ndarray:
@@ -144,12 +148,17 @@ class FreqDependent(MatProperty):
         data[:,:,self._apply_to] = np.repeat(self._func(f)[:,:,np.newaxis], self._apply_to.shape[0], axis=2)
         return data
     
+    @property
+    def value(self) -> float:
+        raise ValueError('Frequency dependent material properties have no fixed value. Use the scalar(f) method to get the value at a specific frequency.')
+    
     def scalar(self, f: float):
         return self._func(f)[0,0]
     
-class CoordDependent(MatProperty):
+class CoordDependent(MatProperty,Saveable):
     _freq_dependent: bool = False
     _coord_dependent: bool = True
+    skip_fields = ("_func","_fmax")
     def __init__(self, 
                  max_value: float,
                  scalar: Callable | None = None,
@@ -173,9 +182,7 @@ class CoordDependent(MatProperty):
             scalar (Callable | None, optional): The scalar value function returning a float/complex. Defaults to None.
             vector (Callable | None, optional): The diagonal rank-2 tensor function returning a (3,) array. Defaults to None.
             matrix (Callable | None, optional): The rank-2 tensor function returning a (3,3) array. Defaults to None.
-
-        Returns:
-            _type_: _description_
+            
         """
         
         if scalar is not None:
@@ -211,13 +218,17 @@ class CoordDependent(MatProperty):
         data[:,:,self._apply_to] = self._func(self._x, self._y, self._z)
         return data
 
+    @property
+    def value(self) -> float:
+        return self._func(0,0,0)[0,0]
+    
     def scalar(self, f: float):
         return self._func(0,0,0)[0,0]
     
-class FreqCoordDependent(MatProperty):
+class FreqCoordDependent(MatProperty, Saveable):
     _freq_dependent: bool = True
     _coord_dependent: bool = True
-    
+    skip_fields = ("_func","_fmax")
     def __init__(self, 
                  max_value: float,
                  scalar: Callable | None = None,
@@ -240,8 +251,6 @@ class FreqCoordDependent(MatProperty):
             vector (Callable | None, optional): The diagonal rank-2 tensor function returning a (3,) array. Defaults to None.
             matrix (Callable | None, optional): The rank-2 tensor function returning a (3,3) array. Defaults to None.
 
-        Returns:
-            _type_: _description_
         """
         if scalar is not None:
             def _func(f, x, y, z) -> np.ndarray:
@@ -277,10 +286,14 @@ class FreqCoordDependent(MatProperty):
         data[:,:,self._apply_to] = self._func(f,self._x, self._y,self._z)
         return data
     
+    @property
+    def value(self) -> float:
+        raise ValueError('Frequency and coordinate dependent material properties have no fixed value. Use the scalar(f) method to get the value at a specific frequency and coordinate.')
+    
     def scalar(self, f: float):
         return self._func(f, 0,0,0)[0,0]
 
-class Material:
+class Material(Saveable):
     """The Material class generalizes a material in the EMerge FEM environment.
 
     If a scalar value is provided for the relative permittivity or the relative permeability
@@ -295,7 +308,7 @@ class Material:
     
     """
     _pickle_exclude = {"_neff"}
-    
+    skip_fields = ("_neff",)
     def __init__(self,
                  er: float | complex | np.ndarray | MatProperty = 1.0,
                  ur: float | complex | np.ndarray | MatProperty = 1.0,
@@ -330,10 +343,12 @@ class Material:
             self._neff: Callable = lambda f: np.sqrt(self.ur._fmax(f)*self.er._fmax(f))
         else:
             self._neff: Callable = lambda f: _neff
-        hex_str = self.color.lstrip('#')
-        self._color_rgb = tuple(int(hex_str[i:i+2], 16)/255.0 for i in (0, 2, 4))
         self._metal: bool = _metal
     
+    @property
+    def _color_rgb(self) -> tuple[float,float,float]:
+        return tuple(int(self.color.lstrip('#')[i:i+2], 16)/255.0 for i in (0, 2, 4))
+    
     def __getstate__(self):
         state = self.__dict__.copy()
         for k in self._pickle_exclude:
@@ -382,18 +397,14 @@ class Material:
     @property
     def frequency_dependent(self) -> bool:
         """If The material property are at all frequency dependent.
-
-        Returns:
-            bool: _description_
+        
         """
         return self.er._freq_dependent or self.ur._freq_dependent or self.tand._freq_dependent or self.cond._freq_dependent
     
     @property
     def coordinate_dependent(self) -> bool:
         """If the material properties are at all coordinate dependent
-
-        Returns:
-            bool: _description_
+        
         """
         return self.er._coord_dependent or self.ur._coord_dependent or self.tand._coord_dependent or self.cond._coord_dependent
     

{emsutil-0.2.2 → emsutil-0.4.0}/src/emsutil/themes.py

@@ -184,8 +184,66 @@ class _Document(EMergeTheme):
             "#4DBEEEFF",
         ]
 
+class Stylish(EMergeTheme):
+    """ A custom EMerge theme. """
+    def define(self):
+        self.backgroung_grad_1 = "#FFFFFF"
+        self.backgroung_grad_2 = "#FFFFFF"
+        self.grid_color = "#676767FF"
+        self.brightness = 1.0
+        
+        self.label_color = "#FFFFFF"
+        self.text_color = "#000000FF"
+        self.render_metal = True
+        self.line_width = 3.0
+        
+        self.geo_edge_width = 3.0
+        self.geo_edge_color = "#000000ff"
+        
+
+        self.draw_xax = False
+        self.draw_yax = False
+        self.draw_zax = False
+        self.draw_xplane = False
+        self.draw_yplane = False
+        self.draw_zplane = False
+        self.draw_xgrid = False
+        self.draw_ygrid = False
+        self.draw_zgrid = False
+        
+        self.axis_x_color = "#FF0000FF"
+        self.axis_y_color = "#00FF00FF"
+        self.axis_z_color = "#0000FFFF"
 
+        self.aa_active = True
+        self.aa_samples = 5
+        self.cmap_npts = 64
+        
+        self.draw_pvgrid = False
+        self.render_shadows = True
+        # Basic clear academic scales
+        # Amplitude is Jet
+        # Wave is blue to transparent to red
+        
+        self.colormaps = {
+            'amplitude': (("#0000FF", "#00FFFF", "#00FF00", "#FFFF00", "#FF0000"),(0.0, 0.25, 0.5, 0.75, 1.0)),
+            'wave': (("#FF0000", "#FFAAAA00","#0000FF00", "#0000FF"), (0.0, 0.49, 0.51, 1.0)),
+        }
+        
+        self.render_style = 'surface'
+        
+        # Clear high contrast matlab colors for paper
+        self.line_color_cycle = [
+            "#003E67FF",
+            "#7F2600FF",
+            "#5B09A7FF",
+            "#0C5D14FF",
+            "#A35100FF",
+            "#6C006AFF",
+        ]
+        
 VaporWave = _VaporWave()
 Vintage = _Vintage()
 Tron = _Tron()
 Document = _Document()
+Stylish = Stylish()

{emsutil-0.2.2 → emsutil-0.4.0}/uv.lock

@@ -279,11 +279,13 @@ wheels = [
 
 [[package]]
 name = "emsutil"
-version = "0.2.2"
+version = "0.4.0"
 source = { editable = "." }
 dependencies = [
     { name = "loguru" },
     { name = "matplotlib" },
+    { name = "msgpack" },
+    { name = "msgpack-numpy" },
     { name = "numpy" },
     { name = "pyvista" },
     { name = "scipy", version = "1.15.3", source = { registry = "https://pypi.org/simple" }, marker = "python_full_version < '3.11'" },
@@ -294,6 +296,8 @@ dependencies = [
 requires-dist = [
     { name = "loguru", specifier = ">=0.7.3" },
     { name = "matplotlib", specifier = ">=3.8.0" },
+    { name = "msgpack", specifier = ">=1.1.2" },
+    { name = "msgpack-numpy", specifier = ">=0.4.8" },
     { name = "numpy", specifier = ">=1.24,<2.3" },
     { name = "pyvista", specifier = ">=0.46.4" },
     { name = "scipy", specifier = ">=1.14.0" },
@@ -560,6 +564,80 @@ wheels = [
     { url = "https://files.pythonhosted.org/packages/73/e4/6d6f14b2a759c622f191b2d67e9075a3f56aaccb3be4bb9bb6890030d0a0/matplotlib-3.10.8-pp311-pypy311_pp73-manylinux2014_x86_64.manylinux_2_17_x86_64.whl", hash = "sha256:1ae029229a57cd1e8fe542485f27e7ca7b23aa9e8944ddb4985d0bc444f1eca2", size = 8713867, upload-time = "2025-12-10T22:56:48.954Z" },
 ]
 
+[[package]]
+name = "msgpack"
+version = "1.1.2"
+source = { registry = "https://pypi.org/simple" }
+sdist = { url = "https://files.pythonhosted.org/packages/4d/f2/bfb55a6236ed8725a96b0aa3acbd0ec17588e6a2c3b62a93eb513ed8783f/msgpack-1.1.2.tar.gz", hash = "sha256:3b60763c1373dd60f398488069bcdc703cd08a711477b5d480eecc9f9626f47e", size = 173581, upload-time = "2025-10-08T09:15:56.596Z" }
+wheels = [
+    { url = "https://files.pythonhosted.org/packages/f5/a2/3b68a9e769db68668b25c6108444a35f9bd163bb848c0650d516761a59c0/msgpack-1.1.2-cp310-cp310-macosx_10_9_x86_64.whl", hash = "sha256:0051fffef5a37ca2cd16978ae4f0aef92f164df86823871b5162812bebecd8e2", size = 81318, upload-time = "2025-10-08T09:14:38.722Z" },
+    { url = "https://files.pythonhosted.org/packages/5b/e1/2b720cc341325c00be44e1ed59e7cfeae2678329fbf5aa68f5bda57fe728/msgpack-1.1.2-cp310-cp310-macosx_11_0_arm64.whl", hash = "sha256:a605409040f2da88676e9c9e5853b3449ba8011973616189ea5ee55ddbc5bc87", size = 83786, upload-time = "2025-10-08T09:14:40.082Z" },
+    { url = "https://files.pythonhosted.org/packages/71/e5/c2241de64bfceac456b140737812a2ab310b10538a7b34a1d393b748e095/msgpack-1.1.2-cp310-cp310-manylinux2014_aarch64.manylinux_2_17_aarch64.manylinux_2_28_aarch64.whl", hash = "sha256:8b696e83c9f1532b4af884045ba7f3aa741a63b2bc22617293a2c6a7c645f251", size = 398240, upload-time = "2025-10-08T09:14:41.151Z" },
+    { url = "https://files.pythonhosted.org/packages/b7/09/2a06956383c0fdebaef5aa9246e2356776f12ea6f2a44bd1368abf0e46c4/msgpack-1.1.2-cp310-cp310-manylinux2014_x86_64.manylinux_2_17_x86_64.manylinux_2_28_x86_64.whl", hash = "sha256:365c0bbe981a27d8932da71af63ef86acc59ed5c01ad929e09a0b88c6294e28a", size = 406070, upload-time = "2025-10-08T09:14:42.821Z" },
+    { url = "https://files.pythonhosted.org/packages/0e/74/2957703f0e1ef20637d6aead4fbb314330c26f39aa046b348c7edcf6ca6b/msgpack-1.1.2-cp310-cp310-musllinux_1_2_aarch64.whl", hash = "sha256:41d1a5d875680166d3ac5c38573896453bbbea7092936d2e107214daf43b1d4f", size = 393403, upload-time = "2025-10-08T09:14:44.38Z" },
+    { url = "https://files.pythonhosted.org/packages/a5/09/3bfc12aa90f77b37322fc33e7a8a7c29ba7c8edeadfa27664451801b9860/msgpack-1.1.2-cp310-cp310-musllinux_1_2_x86_64.whl", hash = "sha256:354e81bcdebaab427c3df4281187edc765d5d76bfb3a7c125af9da7a27e8458f", size = 398947, upload-time = "2025-10-08T09:14:45.56Z" },
+    { url = "https://files.pythonhosted.org/packages/4b/4f/05fcebd3b4977cb3d840f7ef6b77c51f8582086de5e642f3fefee35c86fc/msgpack-1.1.2-cp310-cp310-win32.whl", hash = "sha256:e64c8d2f5e5d5fda7b842f55dec6133260ea8f53c4257d64494c534f306bf7a9", size = 64769, upload-time = "2025-10-08T09:14:47.334Z" },
+    { url = "https://files.pythonhosted.org/packages/d0/3e/b4547e3a34210956382eed1c85935fff7e0f9b98be3106b3745d7dec9c5e/msgpack-1.1.2-cp310-cp310-win_amd64.whl", hash = "sha256:db6192777d943bdaaafb6ba66d44bf65aa0e9c5616fa1d2da9bb08828c6b39aa", size = 71293, upload-time = "2025-10-08T09:14:48.665Z" },
+    { url = "https://files.pythonhosted.org/packages/2c/97/560d11202bcd537abca693fd85d81cebe2107ba17301de42b01ac1677b69/msgpack-1.1.2-cp311-cp311-macosx_10_9_x86_64.whl", hash = "sha256:2e86a607e558d22985d856948c12a3fa7b42efad264dca8a3ebbcfa2735d786c", size = 82271, upload-time = "2025-10-08T09:14:49.967Z" },
+    { url = "https://files.pythonhosted.org/packages/83/04/28a41024ccbd67467380b6fb440ae916c1e4f25e2cd4c63abe6835ac566e/msgpack-1.1.2-cp311-cp311-macosx_11_0_arm64.whl", hash = "sha256:283ae72fc89da59aa004ba147e8fc2f766647b1251500182fac0350d8af299c0", size = 84914, upload-time = "2025-10-08T09:14:50.958Z" },
+    { url = "https://files.pythonhosted.org/packages/71/46/b817349db6886d79e57a966346cf0902a426375aadc1e8e7a86a75e22f19/msgpack-1.1.2-cp311-cp311-manylinux2014_aarch64.manylinux_2_17_aarch64.manylinux_2_28_aarch64.whl", hash = "sha256:61c8aa3bd513d87c72ed0b37b53dd5c5a0f58f2ff9f26e1555d3bd7948fb7296", size = 416962, upload-time = "2025-10-08T09:14:51.997Z" },
+    { url = "https://files.pythonhosted.org/packages/da/e0/6cc2e852837cd6086fe7d8406af4294e66827a60a4cf60b86575a4a65ca8/msgpack-1.1.2-cp311-cp311-manylinux2014_x86_64.manylinux_2_17_x86_64.manylinux_2_28_x86_64.whl", hash = "sha256:454e29e186285d2ebe65be34629fa0e8605202c60fbc7c4c650ccd41870896ef", size = 426183, upload-time = "2025-10-08T09:14:53.477Z" },
+    { url = "https://files.pythonhosted.org/packages/25/98/6a19f030b3d2ea906696cedd1eb251708e50a5891d0978b012cb6107234c/msgpack-1.1.2-cp311-cp311-musllinux_1_2_aarch64.whl", hash = "sha256:7bc8813f88417599564fafa59fd6f95be417179f76b40325b500b3c98409757c", size = 411454, upload-time = "2025-10-08T09:14:54.648Z" },
+    { url = "https://files.pythonhosted.org/packages/b7/cd/9098fcb6adb32187a70b7ecaabf6339da50553351558f37600e53a4a2a23/msgpack-1.1.2-cp311-cp311-musllinux_1_2_x86_64.whl", hash = "sha256:bafca952dc13907bdfdedfc6a5f579bf4f292bdd506fadb38389afa3ac5b208e", size = 422341, upload-time = "2025-10-08T09:14:56.328Z" },
+    { url = "https://files.pythonhosted.org/packages/e6/ae/270cecbcf36c1dc85ec086b33a51a4d7d08fc4f404bdbc15b582255d05ff/msgpack-1.1.2-cp311-cp311-win32.whl", hash = "sha256:602b6740e95ffc55bfb078172d279de3773d7b7db1f703b2f1323566b878b90e", size = 64747, upload-time = "2025-10-08T09:14:57.882Z" },
+    { url = "https://files.pythonhosted.org/packages/2a/79/309d0e637f6f37e83c711f547308b91af02b72d2326ddd860b966080ef29/msgpack-1.1.2-cp311-cp311-win_amd64.whl", hash = "sha256:d198d275222dc54244bf3327eb8cbe00307d220241d9cec4d306d49a44e85f68", size = 71633, upload-time = "2025-10-08T09:14:59.177Z" },
+    { url = "https://files.pythonhosted.org/packages/73/4d/7c4e2b3d9b1106cd0aa6cb56cc57c6267f59fa8bfab7d91df5adc802c847/msgpack-1.1.2-cp311-cp311-win_arm64.whl", hash = "sha256:86f8136dfa5c116365a8a651a7d7484b65b13339731dd6faebb9a0242151c406", size = 64755, upload-time = "2025-10-08T09:15:00.48Z" },
+    { url = "https://files.pythonhosted.org/packages/ad/bd/8b0d01c756203fbab65d265859749860682ccd2a59594609aeec3a144efa/msgpack-1.1.2-cp312-cp312-macosx_10_13_x86_64.whl", hash = "sha256:70a0dff9d1f8da25179ffcf880e10cf1aad55fdb63cd59c9a49a1b82290062aa", size = 81939, upload-time = "2025-10-08T09:15:01.472Z" },
+    { url = "https://files.pythonhosted.org/packages/34/68/ba4f155f793a74c1483d4bdef136e1023f7bcba557f0db4ef3db3c665cf1/msgpack-1.1.2-cp312-cp312-macosx_11_0_arm64.whl", hash = "sha256:446abdd8b94b55c800ac34b102dffd2f6aa0ce643c55dfc017ad89347db3dbdb", size = 85064, upload-time = "2025-10-08T09:15:03.764Z" },
+    { url = "https://files.pythonhosted.org/packages/f2/60/a064b0345fc36c4c3d2c743c82d9100c40388d77f0b48b2f04d6041dbec1/msgpack-1.1.2-cp312-cp312-manylinux2014_aarch64.manylinux_2_17_aarch64.manylinux_2_28_aarch64.whl", hash = "sha256:c63eea553c69ab05b6747901b97d620bb2a690633c77f23feb0c6a947a8a7b8f", size = 417131, upload-time = "2025-10-08T09:15:05.136Z" },
+    { url = "https://files.pythonhosted.org/packages/65/92/a5100f7185a800a5d29f8d14041f61475b9de465ffcc0f3b9fba606e4505/msgpack-1.1.2-cp312-cp312-manylinux2014_x86_64.manylinux_2_17_x86_64.manylinux_2_28_x86_64.whl", hash = "sha256:372839311ccf6bdaf39b00b61288e0557916c3729529b301c52c2d88842add42", size = 427556, upload-time = "2025-10-08T09:15:06.837Z" },
+    { url = "https://files.pythonhosted.org/packages/f5/87/ffe21d1bf7d9991354ad93949286f643b2bb6ddbeab66373922b44c3b8cc/msgpack-1.1.2-cp312-cp312-musllinux_1_2_aarch64.whl", hash = "sha256:2929af52106ca73fcb28576218476ffbb531a036c2adbcf54a3664de124303e9", size = 404920, upload-time = "2025-10-08T09:15:08.179Z" },
+    { url = "https://files.pythonhosted.org/packages/ff/41/8543ed2b8604f7c0d89ce066f42007faac1eaa7d79a81555f206a5cdb889/msgpack-1.1.2-cp312-cp312-musllinux_1_2_x86_64.whl", hash = "sha256:be52a8fc79e45b0364210eef5234a7cf8d330836d0a64dfbb878efa903d84620", size = 415013, upload-time = "2025-10-08T09:15:09.83Z" },
+    { url = "https://files.pythonhosted.org/packages/41/0d/2ddfaa8b7e1cee6c490d46cb0a39742b19e2481600a7a0e96537e9c22f43/msgpack-1.1.2-cp312-cp312-win32.whl", hash = "sha256:1fff3d825d7859ac888b0fbda39a42d59193543920eda9d9bea44d958a878029", size = 65096, upload-time = "2025-10-08T09:15:11.11Z" },
+    { url = "https://files.pythonhosted.org/packages/8c/ec/d431eb7941fb55a31dd6ca3404d41fbb52d99172df2e7707754488390910/msgpack-1.1.2-cp312-cp312-win_amd64.whl", hash = "sha256:1de460f0403172cff81169a30b9a92b260cb809c4cb7e2fc79ae8d0510c78b6b", size = 72708, upload-time = "2025-10-08T09:15:12.554Z" },
+    { url = "https://files.pythonhosted.org/packages/c5/31/5b1a1f70eb0e87d1678e9624908f86317787b536060641d6798e3cf70ace/msgpack-1.1.2-cp312-cp312-win_arm64.whl", hash = "sha256:be5980f3ee0e6bd44f3a9e9dea01054f175b50c3e6cdb692bc9424c0bbb8bf69", size = 64119, upload-time = "2025-10-08T09:15:13.589Z" },
+    { url = "https://files.pythonhosted.org/packages/6b/31/b46518ecc604d7edf3a4f94cb3bf021fc62aa301f0cb849936968164ef23/msgpack-1.1.2-cp313-cp313-macosx_10_13_x86_64.whl", hash = "sha256:4efd7b5979ccb539c221a4c4e16aac1a533efc97f3b759bb5a5ac9f6d10383bf", size = 81212, upload-time = "2025-10-08T09:15:14.552Z" },
+    { url = "https://files.pythonhosted.org/packages/92/dc/c385f38f2c2433333345a82926c6bfa5ecfff3ef787201614317b58dd8be/msgpack-1.1.2-cp313-cp313-macosx_11_0_arm64.whl", hash = "sha256:42eefe2c3e2af97ed470eec850facbe1b5ad1d6eacdbadc42ec98e7dcf68b4b7", size = 84315, upload-time = "2025-10-08T09:15:15.543Z" },
+    { url = "https://files.pythonhosted.org/packages/d3/68/93180dce57f684a61a88a45ed13047558ded2be46f03acb8dec6d7c513af/msgpack-1.1.2-cp313-cp313-manylinux2014_aarch64.manylinux_2_17_aarch64.manylinux_2_28_aarch64.whl", hash = "sha256:1fdf7d83102bf09e7ce3357de96c59b627395352a4024f6e2458501f158bf999", size = 412721, upload-time = "2025-10-08T09:15:16.567Z" },
+    { url = "https://files.pythonhosted.org/packages/5d/ba/459f18c16f2b3fc1a1ca871f72f07d70c07bf768ad0a507a698b8052ac58/msgpack-1.1.2-cp313-cp313-manylinux2014_x86_64.manylinux_2_17_x86_64.manylinux_2_28_x86_64.whl", hash = "sha256:fac4be746328f90caa3cd4bc67e6fe36ca2bf61d5c6eb6d895b6527e3f05071e", size = 424657, upload-time = "2025-10-08T09:15:17.825Z" },
+    { url = "https://files.pythonhosted.org/packages/38/f8/4398c46863b093252fe67368b44edc6c13b17f4e6b0e4929dbf0bdb13f23/msgpack-1.1.2-cp313-cp313-musllinux_1_2_aarch64.whl", hash = "sha256:fffee09044073e69f2bad787071aeec727183e7580443dfeb8556cbf1978d162", size = 402668, upload-time = "2025-10-08T09:15:19.003Z" },
+    { url = "https://files.pythonhosted.org/packages/28/ce/698c1eff75626e4124b4d78e21cca0b4cc90043afb80a507626ea354ab52/msgpack-1.1.2-cp313-cp313-musllinux_1_2_x86_64.whl", hash = "sha256:5928604de9b032bc17f5099496417f113c45bc6bc21b5c6920caf34b3c428794", size = 419040, upload-time = "2025-10-08T09:15:20.183Z" },
+    { url = "https://files.pythonhosted.org/packages/67/32/f3cd1667028424fa7001d82e10ee35386eea1408b93d399b09fb0aa7875f/msgpack-1.1.2-cp313-cp313-win32.whl", hash = "sha256:a7787d353595c7c7e145e2331abf8b7ff1e6673a6b974ded96e6d4ec09f00c8c", size = 65037, upload-time = "2025-10-08T09:15:21.416Z" },
+    { url = "https://files.pythonhosted.org/packages/74/07/1ed8277f8653c40ebc65985180b007879f6a836c525b3885dcc6448ae6cb/msgpack-1.1.2-cp313-cp313-win_amd64.whl", hash = "sha256:a465f0dceb8e13a487e54c07d04ae3ba131c7c5b95e2612596eafde1dccf64a9", size = 72631, upload-time = "2025-10-08T09:15:22.431Z" },
+    { url = "https://files.pythonhosted.org/packages/e5/db/0314e4e2db56ebcf450f277904ffd84a7988b9e5da8d0d61ab2d057df2b6/msgpack-1.1.2-cp313-cp313-win_arm64.whl", hash = "sha256:e69b39f8c0aa5ec24b57737ebee40be647035158f14ed4b40e6f150077e21a84", size = 64118, upload-time = "2025-10-08T09:15:23.402Z" },
+    { url = "https://files.pythonhosted.org/packages/22/71/201105712d0a2ff07b7873ed3c220292fb2ea5120603c00c4b634bcdafb3/msgpack-1.1.2-cp314-cp314-macosx_10_13_x86_64.whl", hash = "sha256:e23ce8d5f7aa6ea6d2a2b326b4ba46c985dbb204523759984430db7114f8aa00", size = 81127, upload-time = "2025-10-08T09:15:24.408Z" },
+    { url = "https://files.pythonhosted.org/packages/1b/9f/38ff9e57a2eade7bf9dfee5eae17f39fc0e998658050279cbb14d97d36d9/msgpack-1.1.2-cp314-cp314-macosx_11_0_arm64.whl", hash = "sha256:6c15b7d74c939ebe620dd8e559384be806204d73b4f9356320632d783d1f7939", size = 84981, upload-time = "2025-10-08T09:15:25.812Z" },
+    { url = "https://files.pythonhosted.org/packages/8e/a9/3536e385167b88c2cc8f4424c49e28d49a6fc35206d4a8060f136e71f94c/msgpack-1.1.2-cp314-cp314-manylinux2014_aarch64.manylinux_2_17_aarch64.manylinux_2_28_aarch64.whl", hash = "sha256:99e2cb7b9031568a2a5c73aa077180f93dd2e95b4f8d3b8e14a73ae94a9e667e", size = 411885, upload-time = "2025-10-08T09:15:27.22Z" },
+    { url = "https://files.pythonhosted.org/packages/2f/40/dc34d1a8d5f1e51fc64640b62b191684da52ca469da9cd74e84936ffa4a6/msgpack-1.1.2-cp314-cp314-manylinux2014_x86_64.manylinux_2_17_x86_64.manylinux_2_28_x86_64.whl", hash = "sha256:180759d89a057eab503cf62eeec0aa61c4ea1200dee709f3a8e9397dbb3b6931", size = 419658, upload-time = "2025-10-08T09:15:28.4Z" },
+    { url = "https://files.pythonhosted.org/packages/3b/ef/2b92e286366500a09a67e03496ee8b8ba00562797a52f3c117aa2b29514b/msgpack-1.1.2-cp314-cp314-musllinux_1_2_aarch64.whl", hash = "sha256:04fb995247a6e83830b62f0b07bf36540c213f6eac8e851166d8d86d83cbd014", size = 403290, upload-time = "2025-10-08T09:15:29.764Z" },
+    { url = "https://files.pythonhosted.org/packages/78/90/e0ea7990abea5764e4655b8177aa7c63cdfa89945b6e7641055800f6c16b/msgpack-1.1.2-cp314-cp314-musllinux_1_2_x86_64.whl", hash = "sha256:8e22ab046fa7ede9e36eeb4cfad44d46450f37bb05d5ec482b02868f451c95e2", size = 415234, upload-time = "2025-10-08T09:15:31.022Z" },
+    { url = "https://files.pythonhosted.org/packages/72/4e/9390aed5db983a2310818cd7d3ec0aecad45e1f7007e0cda79c79507bb0d/msgpack-1.1.2-cp314-cp314-win32.whl", hash = "sha256:80a0ff7d4abf5fecb995fcf235d4064b9a9a8a40a3ab80999e6ac1e30b702717", size = 66391, upload-time = "2025-10-08T09:15:32.265Z" },
+    { url = "https://files.pythonhosted.org/packages/6e/f1/abd09c2ae91228c5f3998dbd7f41353def9eac64253de3c8105efa2082f7/msgpack-1.1.2-cp314-cp314-win_amd64.whl", hash = "sha256:9ade919fac6a3e7260b7f64cea89df6bec59104987cbea34d34a2fa15d74310b", size = 73787, upload-time = "2025-10-08T09:15:33.219Z" },
+    { url = "https://files.pythonhosted.org/packages/6a/b0/9d9f667ab48b16ad4115c1935d94023b82b3198064cb84a123e97f7466c1/msgpack-1.1.2-cp314-cp314-win_arm64.whl", hash = "sha256:59415c6076b1e30e563eb732e23b994a61c159cec44deaf584e5cc1dd662f2af", size = 66453, upload-time = "2025-10-08T09:15:34.225Z" },
+    { url = "https://files.pythonhosted.org/packages/16/67/93f80545eb1792b61a217fa7f06d5e5cb9e0055bed867f43e2b8e012e137/msgpack-1.1.2-cp314-cp314t-macosx_10_13_x86_64.whl", hash = "sha256:897c478140877e5307760b0ea66e0932738879e7aa68144d9b78ea4c8302a84a", size = 85264, upload-time = "2025-10-08T09:15:35.61Z" },
+    { url = "https://files.pythonhosted.org/packages/87/1c/33c8a24959cf193966ef11a6f6a2995a65eb066bd681fd085afd519a57ce/msgpack-1.1.2-cp314-cp314t-macosx_11_0_arm64.whl", hash = "sha256:a668204fa43e6d02f89dbe79a30b0d67238d9ec4c5bd8a940fc3a004a47b721b", size = 89076, upload-time = "2025-10-08T09:15:36.619Z" },
+    { url = "https://files.pythonhosted.org/packages/fc/6b/62e85ff7193663fbea5c0254ef32f0c77134b4059f8da89b958beb7696f3/msgpack-1.1.2-cp314-cp314t-manylinux2014_aarch64.manylinux_2_17_aarch64.manylinux_2_28_aarch64.whl", hash = "sha256:5559d03930d3aa0f3aacb4c42c776af1a2ace2611871c84a75afe436695e6245", size = 435242, upload-time = "2025-10-08T09:15:37.647Z" },
+    { url = "https://files.pythonhosted.org/packages/c1/47/5c74ecb4cc277cf09f64e913947871682ffa82b3b93c8dad68083112f412/msgpack-1.1.2-cp314-cp314t-manylinux2014_x86_64.manylinux_2_17_x86_64.manylinux_2_28_x86_64.whl", hash = "sha256:70c5a7a9fea7f036b716191c29047374c10721c389c21e9ffafad04df8c52c90", size = 432509, upload-time = "2025-10-08T09:15:38.794Z" },
+    { url = "https://files.pythonhosted.org/packages/24/a4/e98ccdb56dc4e98c929a3f150de1799831c0a800583cde9fa022fa90602d/msgpack-1.1.2-cp314-cp314t-musllinux_1_2_aarch64.whl", hash = "sha256:f2cb069d8b981abc72b41aea1c580ce92d57c673ec61af4c500153a626cb9e20", size = 415957, upload-time = "2025-10-08T09:15:40.238Z" },
+    { url = "https://files.pythonhosted.org/packages/da/28/6951f7fb67bc0a4e184a6b38ab71a92d9ba58080b27a77d3e2fb0be5998f/msgpack-1.1.2-cp314-cp314t-musllinux_1_2_x86_64.whl", hash = "sha256:d62ce1f483f355f61adb5433ebfd8868c5f078d1a52d042b0a998682b4fa8c27", size = 422910, upload-time = "2025-10-08T09:15:41.505Z" },
+    { url = "https://files.pythonhosted.org/packages/f0/03/42106dcded51f0a0b5284d3ce30a671e7bd3f7318d122b2ead66ad289fed/msgpack-1.1.2-cp314-cp314t-win32.whl", hash = "sha256:1d1418482b1ee984625d88aa9585db570180c286d942da463533b238b98b812b", size = 75197, upload-time = "2025-10-08T09:15:42.954Z" },
+    { url = "https://files.pythonhosted.org/packages/15/86/d0071e94987f8db59d4eeb386ddc64d0bb9b10820a8d82bcd3e53eeb2da6/msgpack-1.1.2-cp314-cp314t-win_amd64.whl", hash = "sha256:5a46bf7e831d09470ad92dff02b8b1ac92175ca36b087f904a0519857c6be3ff", size = 85772, upload-time = "2025-10-08T09:15:43.954Z" },
+    { url = "https://files.pythonhosted.org/packages/81/f2/08ace4142eb281c12701fc3b93a10795e4d4dc7f753911d836675050f886/msgpack-1.1.2-cp314-cp314t-win_arm64.whl", hash = "sha256:d99ef64f349d5ec3293688e91486c5fdb925ed03807f64d98d205d2713c60b46", size = 70868, upload-time = "2025-10-08T09:15:44.959Z" },
+]
+
+[[package]]
+name = "msgpack-numpy"
+version = "0.4.8"
+source = { registry = "https://pypi.org/simple" }
+dependencies = [
+    { name = "msgpack" },
+    { name = "numpy" },
+]
+sdist = { url = "https://files.pythonhosted.org/packages/08/94/61e8aee142733ebfdc400a05bdac6e1763c4514bba3b42743d223f388450/msgpack-numpy-0.4.8.tar.gz", hash = "sha256:c667d3180513422f9c7545be5eec5d296dcbb357e06f72ed39cc683797556e69", size = 10923, upload-time = "2022-06-09T03:43:08.739Z" }
+wheels = [
+    { url = "https://files.pythonhosted.org/packages/9b/5d/f25ac7d4fb77cbd53ddc6d05d833c6bf52b12770a44fa9a447eed470ca9a/msgpack_numpy-0.4.8-py2.py3-none-any.whl", hash = "sha256:773c19d4dfbae1b3c7b791083e2caf66983bb19b40901646f61d8731554ae3da", size = 6919, upload-time = "2022-06-09T03:43:06.82Z" },
+]
+
 [[package]]
 name = "numpy"
 version = "2.2.6"

emsutil-0.2.2/src/emsutil/inexport/ffdata.py

@@ -1,58 +0,0 @@
-import numpy as np
-from .strutil import arry_to_line, matrix_to_lines, arry_to_fwl
-from ..emdata import EHFieldFF
-
-def export_ffdata(filename: str, 
-                  thetas: np.ndarray,
-                  phis: np.ndarray,
-                  frequencies: np.ndarray,
-                  fields: list[EHFieldFF],
-                  precision: int = 4) -> None:
-    
-    lines = []
-    lines.append('% Theta (deg)')
-    lines.append(arry_to_line(thetas, precision=precision))
-    lines.append('% Phi (deg)')
-    lines.append(arry_to_line(phis, precision=precision))
-    lines.append('')
-    nF = frequencies.shape[0]
-    
-    actual_frequencies = []
-    blocks = []
-    
-    T, P = np.meshgrid(thetas, phis, indexing='ij')
-    
-    thetal = T.flatten()
-    phil = P.flatten()
-    
-    for iF in range(nF):
-        freq = frequencies[iF]
-        actual_frequencies.append(freq)
-        
-        farfield = fields[iF]._E
-        Fx = farfield[0,:,:].squeeze().flatten()
-        Fy = farfield[1,:,:].squeeze().flatten()
-        Fz = farfield[2,:,:].squeeze().flatten()
-        
-        block_lines = []
-        
-        block_lines.append(f'# {freq} (Hz)')
-        block_lines.append('$ Theta(deg); Phi(deg); E_x[re](V/m); E_x[im](V/m); E_y[re](V/m); E_y[im](V/m); E_z[re](V/m); E_z[im](V/m)')
-        positions = (0, 14, 24, 38, 52, 66, 80, 94)
-        for th, ph, ex, ey, ez in zip(thetal, phil, Fx, Fy, Fz):
-            re_ex = np.real(ex)
-            im_ex = np.imag(ex)
-            re_ey = np.real(ey)
-            im_ey = np.imag(ey)
-            re_ez = np.real(ez)
-            im_ez = np.imag(ez)
-            line_values = np.array([th, ph, re_ex, im_ex, re_ey, im_ey, re_ez, im_ez])
-            line = arry_to_fwl(line_values, positions, precision=precision)
-            block_lines.append(line)
-
-        lines.extend(block_lines)
-        lines.extend('')
-    
-    text = '\n'.join(lines)
-    with open(filename, 'w') as f:
-        f.write(text)