secs-dsp 0.7.0

f64 DSP primitives and optimizers shared by the engine, the mdat builder, the CamillaDSP writer and the measurement path.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28

//! Adapted from Impulcifer-py313 `crates/impulcifer-dsp/src/peaks.rs` (MIT, Copyright (c) 2024- 115dkk).
//! See docs/THIRD_PARTY.md for attribution and licenses.
//! Height-only peak selection (`scipy.signal.find_peaks` as SECS 2's `mimo` and
//! `mimo_iir` call it).

/// Mirrors `scipy.signal.find_peaks(x, height=height)` of SciPy 1.18.1,
/// pinned by the `dsp` golden family's peak cases.
/// Returns ascending indices, inclusive height, floor midpoint for plateaus,
/// and excludes endpoints. NaN comparisons are false, so NaNs break candidates;
/// an NaN height excludes all peaks. Infinities retain ordinary comparison rules.
pub fn find_peaks(x: &[f64], height: Option<f64>) -> Vec<usize> {
    let mut peaks = Vec::new();
    let mut i = 1;
    while i + 1 < x.len() {
        if x[i - 1] < x[i] {
            let mut right = i;
            while right + 1 < x.len() && x[right + 1] == x[i] {
                right += 1;
            }
            if right + 1 < x.len() && x[right + 1] < x[i] && height.is_none_or(|h| x[i] >= h) {
                peaks.push(i + (right - i) / 2);
            }
            i = right;
        }
        i += 1;
    }
    peaks
}