ezmsg-sigproc 2.13.1__py3-none-any.whl → 2.15.0__py3-none-any.whl

ezmsg/sigproc/__version__.py

@@ -28,7 +28,7 @@ version_tuple: VERSION_TUPLE
 commit_id: COMMIT_ID
 __commit_id__: COMMIT_ID
 
-__version__ = version = '2.13.1'
-__version_tuple__ = version_tuple = (2, 13, 1)
+__version__ = version = '2.15.0'
+__version_tuple__ = version_tuple = (2, 15, 0)
 
 __commit_id__ = commit_id = None

ezmsg/sigproc/affinetransform.py

@@ -23,6 +23,8 @@ from ezmsg.baseproc import (
 from ezmsg.util.messages.axisarray import AxisArray, AxisBase
 from ezmsg.util.messages.util import replace
 
+from ezmsg.sigproc.util.array import array_device, is_float_dtype, xp_asarray, xp_create
+
 
 def _find_block_diagonal_clusters(weights: np.ndarray) -> list[tuple[np.ndarray, np.ndarray]] | None:
     """Detect block-diagonal structure in a weight matrix.
@@ -245,13 +247,25 @@ class AffineTransformTransformer(
 
             self._state.new_axis = replace(message.axes[axis], data=np.array(new_labels))
 
-        # Convert to match message.data namespace for efficient operations in _process
+        # Convert to match message.data namespace and device for _process.
+        # Weights are stored as numpy float64 after cluster detection; some
+        # devices (e.g. MPS) don't support float64, so we downcast weight
+        # arrays to the message's dtype when the message is floating-point.
         xp = get_namespace(message.data)
+        dev = array_device(message.data)
+        msg_dt = message.data.dtype
+        # Downcast weights dtype only for float message data (avoids casting
+        # float weights to integer when message data happens to be int).
+        w_dt = msg_dt if is_float_dtype(xp, msg_dt) else None
         if self._state.weights is not None:
-            self._state.weights = xp.asarray(self._state.weights)
+            self._state.weights = xp_asarray(xp, self._state.weights, dtype=w_dt, device=dev)
         if self._state.clusters is not None:
             self._state.clusters = [
-                (xp.asarray(in_idx), xp.asarray(out_idx), xp.asarray(sub_w))
+                (
+                    xp_asarray(xp, in_idx, device=dev),
+                    xp_asarray(xp, out_idx, device=dev),
+                    xp_asarray(xp, sub_w, dtype=w_dt, device=dev),
+                )
                 for in_idx, out_idx, sub_w in self._state.clusters
             ]
 
@@ -345,7 +359,7 @@ class AffineTransformTransformer(
 
         # Pre-allocate output (omitted channels stay zero)
         out_shape = data.shape[:-1] + (self._state.n_out,)
-        result = xp.zeros(out_shape, dtype=data.dtype)
+        result = xp_create(xp.zeros, out_shape, dtype=data.dtype, device=array_device(data))
 
         for in_idx, out_idx, sub_weights in self._state.clusters:
             chunk = xp.take(data, in_idx, axis=data.ndim - 1)
@@ -371,7 +385,10 @@ class AffineTransformTransformer(
                 # The weights are stacked A|B where A is the transform and B is a single row
                 #  in the equation y = Ax + B. This supports NeuroKey's weights matrices.
                 sample_shape = data.shape[:axis_idx] + (1,) + data.shape[axis_idx + 1 :]
-                data = xp.concat((data, xp.ones(sample_shape, dtype=data.dtype)), axis=axis_idx)
+                data = xp.concat(
+                    (data, xp_create(xp.ones, sample_shape, dtype=data.dtype, device=array_device(data))),
+                    axis=axis_idx,
+                )
 
             if axis_idx in [-1, len(message.dims) - 1]:
                 data = xp.matmul(data, self._state.weights)

ezmsg/sigproc/rollingscaler.py

@@ -12,8 +12,6 @@ from ezmsg.baseproc import (
 from ezmsg.util.messages.axisarray import AxisArray
 from ezmsg.util.messages.util import replace
 
-from ezmsg.sigproc.sampler import SampleMessage
-
 
 class RollingScalerSettings(ez.Settings):
     axis: str = "time"
@@ -168,8 +166,8 @@ class RollingScalerProcessor(BaseAdaptiveTransformer[RollingScalerSettings, Axis
 
         self._state.samples.append((n_b, mean_b, M2_b))
 
-    def partial_fit(self, message: SampleMessage) -> None:
-        x = message.sample.data
+    def partial_fit(self, message: AxisArray) -> None:
+        x = message.data
         self._add_batch_stats(x)
 
     def _process(self, message: AxisArray) -> AxisArray:

ezmsg/sigproc/sampler.py

@@ -21,7 +21,7 @@ from ezmsg.util.messages.util import replace
 
 from .util.axisarray_buffer import HybridAxisArrayBuffer
 from .util.buffer import UpdateStrategy
-from .util.message import SampleMessage, SampleTriggerMessage
+from .util.message import SampleTriggerMessage
 from .util.profile import profile_subpub
 
 
@@ -75,17 +75,7 @@ class SamplerState:
 
 
 class SamplerTransformer(BaseStatefulTransformer[SamplerSettings, AxisArray, AxisArray, SamplerState]):
-    def __call__(self, message: AxisArray | SampleTriggerMessage) -> list[SampleMessage]:
-        # TODO: Currently we have a single entry point that accepts both
-        #  data and trigger messages and we choose a code path based on
-        #  the message type. However, in the future we will likely replace
-        #  SampleTriggerMessage with an agumented form of AxisArray,
-        #  leveraging its attrs field, which makes this a bit harder.
-        #  We should probably force callers of this object to explicitly
-        #  call `push_trigger` for trigger messages. This will also
-        #  simplify typing somewhat because `push_trigger` should not
-        #  return anything yet we currently have it returning an empty
-        #  list just to be compatible with __call__.
+    def __call__(self, message: AxisArray | SampleTriggerMessage) -> list[AxisArray]:
         if isinstance(message, AxisArray):
             return super().__call__(message)
         else:
@@ -109,7 +99,7 @@ class SamplerTransformer(BaseStatefulTransformer[SamplerSettings, AxisArray, Axi
             self._state.triggers = deque()
         self._state.triggers.clear()
 
-    def _process(self, message: AxisArray) -> list[SampleMessage]:
+    def _process(self, message: AxisArray) -> list[AxisArray]:
         self._state.buffer.write(message)
 
         # How much data in the buffer?
@@ -119,7 +109,7 @@ class SamplerTransformer(BaseStatefulTransformer[SamplerSettings, AxisArray, Axi
         )
 
         # Process in reverse order so that we can remove triggers safely as we iterate.
-        msgs_out: list[SampleMessage] = []
+        msgs_out: list[AxisArray] = []
         for trig_ix in range(len(self._state.triggers) - 1, -1, -1):
             trig = self._state.triggers[trig_ix]
             if trig.period is None:
@@ -152,13 +142,13 @@ class SamplerTransformer(BaseStatefulTransformer[SamplerSettings, AxisArray, Axi
             # Note: buffer will trim itself as needed based on buffer_dur.
 
             # Prepare output and drop trigger
-            msgs_out.append(SampleMessage(trigger=copy.copy(trig), sample=buff_axarr))
+            msgs_out.append(replace(buff_axarr, attrs={**buff_axarr.attrs, "trigger": copy.copy(trig)}))
             del self._state.triggers[trig_ix]
 
         msgs_out.reverse()  # in-place
         return msgs_out
 
-    def push_trigger(self, message: SampleTriggerMessage) -> list[SampleMessage]:
+    def push_trigger(self, message: SampleTriggerMessage) -> list[AxisArray]:
         # Input is a trigger message that we will use to sample the buffer.
 
         if self._state.buffer is None:
@@ -198,7 +188,7 @@ class Sampler(BaseTransformerUnit[SamplerSettings, AxisArray, AxisArray, Sampler
     SETTINGS = SamplerSettings
 
     INPUT_TRIGGER = ez.InputStream(SampleTriggerMessage)
-    OUTPUT_SIGNAL = ez.OutputStream(SampleMessage)
+    OUTPUT_SIGNAL = ez.OutputStream(AxisArray)
 
     @ez.subscriber(INPUT_TRIGGER)
     async def on_trigger(self, msg: SampleTriggerMessage) -> None:
@@ -228,7 +218,8 @@ def sampler(
 
     Returns:
         A generator that expects `.send` either an :obj:`AxisArray` containing streaming data messages,
-        or a :obj:`SampleTriggerMessage` containing a trigger, and yields the list of :obj:`SampleMessage` s.
+        or a :obj:`SampleTriggerMessage` containing a trigger, and yields the list of :obj:`AxisArray` s
+        with trigger info stored in ``attrs["trigger"]``.
     """
     return SamplerTransformer(
         settings=SamplerSettings(

ezmsg/sigproc/util/array.py

@@ -0,0 +1,59 @@
+"""Portable helpers for Array API interoperability.
+
+These utilities smooth over differences between Array API libraries
+(NumPy, PyTorch, MLX, CuPy, etc.) — in particular around ``device``
+placement and ``dtype`` introspection, which are not uniformly supported.
+"""
+
+import numpy as np
+
+
+def array_device(x):
+    """Return the device of an array, or ``None`` for device-less libraries."""
+    try:
+        from array_api_compat import device
+
+        return device(x)
+    except (AttributeError, TypeError):
+        return None
+
+
+def xp_asarray(xp, obj, *, dtype=None, device=None):
+    """Portable ``xp.asarray`` that omits unsupported kwargs.
+
+    Some Array API libraries (e.g. MLX) don't accept a ``device`` keyword.
+    This helper builds the kwargs dict dynamically so that only supported
+    arguments are forwarded.
+    """
+    kwargs = {}
+    if dtype is not None:
+        kwargs["dtype"] = dtype
+    if device is not None:
+        kwargs["device"] = device
+    return xp.asarray(obj, **kwargs)
+
+
+def xp_create(fn, *args, dtype=None, device=None, **extra):
+    """Call a creation function (``zeros``, ``ones``, ``eye``) portably.
+
+    Omits ``device`` if it is ``None`` (for libraries that don't support it).
+    """
+    kwargs = dict(extra)
+    if dtype is not None:
+        kwargs["dtype"] = dtype
+    if device is not None:
+        kwargs["device"] = device
+    return fn(*args, **kwargs)
+
+
+def is_float_dtype(xp, dtype) -> bool:
+    """Check whether *dtype* is a real floating-point type, portably."""
+    try:
+        return xp.isdtype(dtype, "real floating")
+    except AttributeError:
+        pass
+    # Fallback for libraries without isdtype (e.g. MLX).
+    try:
+        return xp.issubdtype(dtype, xp.floating)
+    except (AttributeError, TypeError):
+        return np.issubdtype(np.dtype(dtype), np.floating)

{ezmsg_sigproc-2.13.1.dist-info → ezmsg_sigproc-2.15.0.dist-info}/METADATA

@@ -1,13 +1,13 @@
 Metadata-Version: 2.4
 Name: ezmsg-sigproc
-Version: 2.13.1
+Version: 2.15.0
 Summary: Timeseries signal processing implementations in ezmsg
 Author-email: Griffin Milsap <griffin.milsap@gmail.com>, Preston Peranich <pperanich@gmail.com>, Chadwick Boulay <chadwick.boulay@gmail.com>, Kyle McGraw <kmcgraw@blackrockneuro.com>
 License-Expression: MIT
 License-File: LICENSE
 Requires-Python: >=3.10.15
 Requires-Dist: array-api-compat>=1.11.1
-Requires-Dist: ezmsg-baseproc>=1.1.0
+Requires-Dist: ezmsg-baseproc>=1.3.0
 Requires-Dist: ezmsg>=3.6.0
 Requires-Dist: numba>=0.61.0
 Requires-Dist: numpy>=1.26.0

{ezmsg_sigproc-2.13.1.dist-info → ezmsg_sigproc-2.15.0.dist-info}/RECORD

@@ -1,8 +1,8 @@
 ezmsg/sigproc/__init__.py,sha256=8K4IcOA3-pfzadoM6s2Sfg5460KlJUocGgyTJTJl96U,52
-ezmsg/sigproc/__version__.py,sha256=LrOgVsOxSmuj7RbXPLo3yIvC77lH9VlW4tk7Ihs46rY,706
+ezmsg/sigproc/__version__.py,sha256=a2xWIti3jyJFUZUTWeA586-E0pruFwlKCcTQl-q5UQQ,706
 ezmsg/sigproc/activation.py,sha256=83vnTa3ZcC4Q3VSWcGfaqhCEqYRNySUOyVpMHZXfz-c,2755
 ezmsg/sigproc/adaptive_lattice_notch.py,sha256=ThUR48mbSHuThkimtD0j4IXNMrOVcpZgGhE7PCYfXhU,8818
-ezmsg/sigproc/affinetransform.py,sha256=mjA21DRVYm0kS2tK7dNR_mU5XAxDbJGuuhnuzz0gtw4,21679
+ezmsg/sigproc/affinetransform.py,sha256=rnbQ2QtLq15eoFrqL-ij12yXSoLLq5ncOwRRzQgB1lA,22574
 ezmsg/sigproc/aggregate.py,sha256=7Hdz1m-S6Cl9h0oRQHeS_UTGBemhOB4XdFyX6cGcdHo,9362
 ezmsg/sigproc/bandpower.py,sha256=dAhH56sUrXNhcRFymTTwjdM_KcU5OxFzrR_sxIPAxyw,2264
 ezmsg/sigproc/base.py,sha256=SJvKEb8gw6mUMwlV5sH0iPG0bXrgS8tvkPwhI-j89MQ,3672
@@ -33,8 +33,8 @@ ezmsg/sigproc/merge.py,sha256=LmuN3LDIZF7DynMcjLp7eGc2G3Yxks9Zd8-luSqqXuA,15436
 ezmsg/sigproc/messages.py,sha256=KQczHTeifn4BZycChN8ZcpfZoQW3lC_xuFmN72QT97s,925
 ezmsg/sigproc/quantize.py,sha256=y7T4_67BHZluX3gyl2anp8iL6EEI6JvsK7Pmp1vapsk,2268
 ezmsg/sigproc/resample.py,sha256=3mm9pvxryNVhQuTCIMW3ToUkUfbVOCsIgvXUiurit1Y,11389
-ezmsg/sigproc/rollingscaler.py,sha256=GcLctVAWTmx9J39r0-dt3e7C_hs25s7M0dDnKiGhkC4,8955
-ezmsg/sigproc/sampler.py,sha256=iOk2YoUX22u9iTjFKimzP5V074RDBVcmswgfyxvZRZo,10761
+ezmsg/sigproc/rollingscaler.py,sha256=1cJf-wESetWJihrIWKJUycSkH_5OG3K5JncxjM-1TTI,8895
+ezmsg/sigproc/sampler.py,sha256=IgQ4d1VE-aaf0U2vkGDXtIwsKfpoW9ZR_s8Fi9VW3Kc,10103
 ezmsg/sigproc/scaler.py,sha256=kVXjRbqoxJ5yJICGsGagRXYIDW3-oihSnbBj-n3s55o,6816
 ezmsg/sigproc/signalinjector.py,sha256=mB62H2b-ScgPtH1jajEpxgDHqdb-RKekQfgyNncsE8Y,2874
 ezmsg/sigproc/singlebandpow.py,sha256=BVlWhFI6zU3ME3EVdZbwf-FMz1d2sfuNFDKXs1hn5HM,4353
@@ -55,6 +55,7 @@ ezmsg/sigproc/math/log.py,sha256=JhjSqLnQnvx_3F4txRYHuUPSJ12Yj2HvRTsCMNvlxpo,202
 ezmsg/sigproc/math/pow.py,sha256=0sdlXFUEBXmpEV_i75oshGRjMguv8L13nLt7hlvdX3E,1284
 ezmsg/sigproc/math/scale.py,sha256=4_xHcHNuf13E1fxIF5vbkPfkN4En6zkfPIKID7lCERk,1133
 ezmsg/sigproc/util/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+ezmsg/sigproc/util/array.py,sha256=S6wTGRbtJpU5I9u2t0MTUAMSKkm1PJref145s87xHLg,1831
 ezmsg/sigproc/util/asio.py,sha256=aAj0e7OoBvkRy28k05HL2s9YPCTxOddc05xMN-qd4lQ,577
 ezmsg/sigproc/util/axisarray_buffer.py,sha256=TGDeC6CXmmp7OUuiGd6xYQijRGYDE4QGdWxjK5Vs3nE,14057
 ezmsg/sigproc/util/buffer.py,sha256=83Gm0IuowmcMlXgLFB_rz8_ZPhkwG4DNNejyWJDKJl8,19658
@@ -62,7 +63,7 @@ ezmsg/sigproc/util/message.py,sha256=ppN3IYtIAwrxWG9JOvgWFn1wDdIumkEzYFfqpH9VQkY
 ezmsg/sigproc/util/profile.py,sha256=eVOo9pXgusrnH1yfRdd2RsM7Dbe2UpyC0LJ9MfGpB08,416
 ezmsg/sigproc/util/sparse.py,sha256=NjbJitCtO0B6CENTlyd9c-lHEJwoCan-T3DIgPyeShw,4834
 ezmsg/sigproc/util/typeresolution.py,sha256=fMFzLi63dqCIclGFLcMdM870OYxJnkeWw6aWKNMk718,362
-ezmsg_sigproc-2.13.1.dist-info/METADATA,sha256=mjyCiq1zy3JbgayLzq3BQLB-i1LOouUeGHe0tn3v7MY,1909
-ezmsg_sigproc-2.13.1.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
-ezmsg_sigproc-2.13.1.dist-info/licenses/LICENSE,sha256=seu0tKhhAMPCUgc1XpXGGaCxY1YaYvFJwqFuQZAl2go,1100
-ezmsg_sigproc-2.13.1.dist-info/RECORD,,
+ezmsg_sigproc-2.15.0.dist-info/METADATA,sha256=WHVPI25XYKw4BAzF8749TZ5k-rg0Ik1iiH-p5ze8byc,1909
+ezmsg_sigproc-2.15.0.dist-info/WHEEL,sha256=WLgqFyCfm_KASv4WHyYy0P3pM_m7J5L9k2skdKLirC8,87
+ezmsg_sigproc-2.15.0.dist-info/licenses/LICENSE,sha256=seu0tKhhAMPCUgc1XpXGGaCxY1YaYvFJwqFuQZAl2go,1100
+ezmsg_sigproc-2.15.0.dist-info/RECORD,,