def METHOD_NAME(self):
    f0 = 110
    freqs = [2.9, 3.1]
    mags = [1.0, 0.5]
    freqs = [freq*f0 for freq in freqs]
    f, m = HarmonicPeaks(maxHarmonics=3)(freqs, mags, f0)
    self.assertAlmostEqualVector(f, [110, 220, 319])
    self.assertEqualVector(m, [0, 0, 1])