Search papers, labs, and topics across Lattice.
This paper reformulates cascaded second-order recursive filtering into a block-tridiagonal linear system and introduces two parallel solution algorithms, achieving over 600 Megasamples per second on a single SIMD core. It addresses the challenge of serializing outputs from multi-core CPUs and GPUs by employing superposition and a divide-and-conquer cyclic reduction method, allowing for efficient asynchronous state propagation. The resulting implementations demonstrate significant performance improvements, with a wavefront pipeline achieving 3.95x scaling on six cores and GPU kernels reaching up to 38.2 Gigasamples per second, outperforming existing parallel recurrence baselines.
Achieving 85% of the memory-bandwidth limit, the proposed methods redefine the performance landscape for real-time and batched recursive filtering on multi-core and GPU architectures.
The companion of this paper reformulated cascaded second-order (biquad) recursive filtering as a block-tridiagonal linear system and developed two parallel solution algorithms, PH factorization and cyclic reduction, reaching over 600 Megasamples per second on a single SIMD core. This paper scales that framework to multi-core CPUs and GPUs, where a new obstacle appears: the terminal outputs of each signal block group are the initial conditions of the next, so naively distributed groups serialize. The dependency is resolved by superposition -- each group's output splits into a zero-state response, computable immediately, and a homogeneous correction applied when the state arrives -- and by a divide-and-conquer form of cyclic reduction that exposes both terminal blocks before back substitution, as asynchronous state propagation requires. Two implementations pair the two dominant deployment scenarios with opposite treatments of the dependency. For real-time streaming, a wavefront pipeline realized with TBB flow graphs parallelizes across cascade sections, preserves first-in-first-out order, and achieves 3.95x scaling on six performance cores, about 2.4 Gigasamples per second for a 16th-order filter. For batched processing, a single-kernel GPU implementation carries each group through the entire cascade in registers and parallelizes across groups with a decoupled look-back protocol; a communication-based cost model, comprising a memory roof, a barrier price, and a latency-hiding floor, reduces tuning to two parameters and predicts the measured behavior across two GPU generations. The best kernels reach 38.2 Gigasamples per second for a single second-order section on an RTX~3060, 85% of the memory-bandwidth roof, exceed the strongest published parallel recurrence baseline at every filter order, and remain numerically valid at order 16, where the direct-form baseline fails.
