dynamic_range()

std::vector< double > MayaFlux::Kinesis::Discrete::dynamic_range	(	std::span< const double >	data,
		size_t	n_windows,
		uint32_t	hop_size,
		uint32_t	window_size
	)

Dynamic range in dB per window.

Computes 20*log10(max_abs / max(min_abs, 1e-10)) per window.

Parameters

data	Input span
n_windows	Pre-computed window count from num_windows()
hop_size	Samples between window starts
window_size	Samples per window

Returns

Per-window dynamic range values in dB

Definition at line 89 of file Analysis.cpp.

{
    std::vector<double> out(n_windows);
    std::vector<size_t> idx(n_windows);
    std::iota(idx.begin(), idx.end(), 0);
 
    Parallel::for_each(Parallel::par_unseq, idx.begin(), idx.end(),
        [&](size_t i) {
            const size_t start = i * hop_size;
            auto w = data.subspan(start, std::min<size_t>(window_size, data.size() - start));
            double mn = std::numeric_limits<double>::max();
            double mx = std::numeric_limits<double>::lowest();
            for (double s : w) {
                double a = std::abs(s);
                mn = std::min(mn, a);
                mx = std::max(mx, a);
            }
            out[i] = 20.0 * std::log10(mx / std::max(mn, k_epsilon));
        });
 
    return out;
}