Y(k) ← H(k) · Y(k), H(k) derived from σ / cutoff / order).
ODE coupling then iterates in the same FFT space: for each bin k, push
P(k) toward jω·Y(k), and pull Y(k) back along
jω. Strength α blends these into one operator.