ezmsg-sigproc 1.7.1__tar.gz → 1.8.2__tar.gz

{ezmsg_sigproc-1.7.1 → ezmsg_sigproc-1.8.2}/PKG-INFO

@@ -1,12 +1,13 @@
 Metadata-Version: 2.4
 Name: ezmsg-sigproc
-Version: 1.7.1
+Version: 1.8.2
 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>
 License-Expression: MIT
 License-File: LICENSE.txt
 Requires-Python: >=3.10.15
 Requires-Dist: ezmsg>=3.6.0
+Requires-Dist: numba>=0.61.0
 Requires-Dist: numpy>=1.26.0
 Requires-Dist: pywavelets>=1.6.0
 Requires-Dist: scipy>=1.13.1

{ezmsg_sigproc-1.7.1 → ezmsg_sigproc-1.8.2}/pyproject.toml

@@ -12,6 +12,7 @@ requires-python = ">=3.10.15"
 dynamic = ["version"]
 dependencies = [
     "ezmsg>=3.6.0",
+    "numba>=0.61.0",
     "numpy>=1.26.0",
     "pywavelets>=1.6.0",
     "scipy>=1.13.1",

{ezmsg_sigproc-1.7.1 → ezmsg_sigproc-1.8.2}/src/ezmsg/sigproc/__version__.py

@@ -17,5 +17,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '1.7.1'
-__version_tuple__ = version_tuple = (1, 7, 1)
+__version__ = version = '1.8.2'
+__version_tuple__ = version_tuple = (1, 8, 2)

{ezmsg_sigproc-1.7.1 → ezmsg_sigproc-1.8.2}/src/ezmsg/sigproc/downsample.py

@@ -14,7 +14,7 @@ from .base import GenAxisArray
 
 @consumer
 def downsample(
-    axis: str | None = None, target_rate: float | None = None
+    axis: str | None = None, target_rate: float | None = None, factor: int | None = None
 ) -> typing.Generator[AxisArray, AxisArray, None]:
     """
     Construct a generator that yields a downsampled version of the data .send() to it.
@@ -28,6 +28,7 @@ def downsample(
             Note: The axis must exist in the message .axes and be of type AxisArray.LinearAxis.
         target_rate: Desired rate after downsampling. The actual rate will be the nearest integer factor of the
             input rate that is the same or higher than the target rate.
+        factor: Explicitly specify downsample factor.  If specified, target_rate is ignored.
 
     Returns:
         A primed generator object ready to receive an :obj:`AxisArray` via `.send(axis_array)`
@@ -39,7 +40,7 @@ def downsample(
     msg_out = AxisArray(np.array([]), dims=[""])
 
     # state variables
-    factor: int = 0  # The integer downsampling factor. It will be determined based on the target rate.
+    q: int = 0  # The integer downsampling factor. It will be determined based on the target rate.
     s_idx: int = 0  # Index of the next msg's first sample into the virtual rotating ds_factor counter.
 
     check_input = {"gain": None, "key": None}
@@ -61,19 +62,21 @@ def downsample(
             check_input["key"] = msg_in.key
             # Reset state variables
             s_idx = 0
-            if target_rate is None:
-                factor = 1
+            if factor is not None:
+                q = factor
+            elif target_rate is None:
+                q = 1
             else:
-                factor = int(1 / (axis_info.gain * target_rate))
-            if factor < 1:
+                q = int(1 / (axis_info.gain * target_rate))
+            if q < 1:
                 ez.logger.warning(
                     f"Target rate {target_rate} cannot be achieved with input rate of {1/axis_info.gain}."
                     "Setting factor to 1."
                 )
-                factor = 1
+                q = 1
 
         n_samples = msg_in.data.shape[axis_idx]
-        samples = np.arange(s_idx, s_idx + n_samples) % factor
+        samples = np.arange(s_idx, s_idx + n_samples) % q
         if n_samples > 0:
             # Update state for next iteration.
             s_idx = samples[-1] + 1
@@ -92,7 +95,7 @@ def downsample(
                 **msg_in.axes,
                 axis: replace(
                     axis_info,
-                    gain=axis_info.gain * factor,
+                    gain=axis_info.gain * q,
                     offset=axis_info.offset + axis_info.gain * n_step,
                 ),
             },
@@ -107,6 +110,7 @@ class DownsampleSettings(ez.Settings):
 
     axis: str | None = None
     target_rate: float | None = None
+    factor: int | None = None
 
 
 class Downsample(GenAxisArray):
@@ -116,5 +120,7 @@ class Downsample(GenAxisArray):
 
     def construct_generator(self):
         self.STATE.gen = downsample(
-            axis=self.SETTINGS.axis, target_rate=self.SETTINGS.target_rate
+            axis=self.SETTINGS.axis,
+            target_rate=self.SETTINGS.target_rate,
+            factor=self.SETTINGS.factor,
         )

{ezmsg_sigproc-1.7.1 → ezmsg_sigproc-1.8.2}/tests/test_downsample.py

@@ -21,7 +21,8 @@ from ezmsg.util.debuglog import DebugLog
 
 @pytest.mark.parametrize("block_size", [1, 5, 10, 20])
 @pytest.mark.parametrize("target_rate", [19.0, 9.5, 6.3])
-def test_downsample_core(block_size: int, target_rate: float):
+@pytest.mark.parametrize("factor", [None, 1, 2])
+def test_downsample_core(block_size: int, target_rate: float, factor: int | None):
     in_fs = 19.0
     test_dur = 4.0
     n_channels = 2
@@ -60,7 +61,7 @@ def test_downsample_core(block_size: int, target_rate: float):
     in_msgs = list(msg_generator())
     backup = [copy.deepcopy(msg) for msg in in_msgs]
 
-    proc = downsample(axis="time", target_rate=target_rate)
+    proc = downsample(axis="time", target_rate=target_rate, factor=factor)
     out_msgs = []
     for msg in in_msgs:
         res = proc.send(msg)
@@ -70,7 +71,7 @@ def test_downsample_core(block_size: int, target_rate: float):
     assert_messages_equal(in_msgs, backup)
 
     # Assert correctness of gain
-    expected_factor: int = int(in_fs // target_rate)
+    expected_factor: int = int(in_fs // target_rate) if factor is None else factor
     assert all(msg.axes["time"].gain == expected_factor / in_fs for msg in out_msgs)
 
     # Assert messages have the correct timestamps
@@ -132,7 +133,13 @@ class DownsampleSystem(ez.Collection):
 
 @pytest.mark.parametrize("block_size", [10])
 @pytest.mark.parametrize("target_rate", [6.3])
-def test_downsample_system(block_size: int, target_rate: float, test_name: str | None = None):
+@pytest.mark.parametrize("factor", [None, 2])
+def test_downsample_system(
+    block_size: int,
+    target_rate: float,
+    factor: int | None,
+    test_name: str | None = None,
+):
     in_fs = 19.0
     num_msgs = int(4.0 / (block_size / in_fs))  # Ensure 4 seconds of data
 
@@ -146,7 +153,7 @@ def test_downsample_system(block_size: int, target_rate: float, test_name: str |
             fs=in_fs,
             dispatch_rate=20.0,
         ),
-        down_settings=DownsampleSettings(target_rate=target_rate),
+        down_settings=DownsampleSettings(target_rate=target_rate, factor=factor),
         log_settings=MessageLoggerSettings(output=test_filename),
         term_settings=TerminateTestSettings(time=1.0),
     )
@@ -160,7 +167,7 @@ def test_downsample_system(block_size: int, target_rate: float, test_name: str |
     ez.logger.info(f"Analyzing recording of { len( messages ) } messages...")
 
     # Check fs
-    expected_factor: int = int(in_fs // target_rate)
+    expected_factor: int = int(in_fs // target_rate) if factor is None else factor
     out_fs = in_fs / expected_factor
     assert np.allclose(
         np.array([1 / msg.axes["time"].gain for msg in messages]),