def METHOD_NAME():
    intervals = librosa.plimit_intervals(primes=[3, 5], bins_per_octave=24, sort=False)
    intervals_s = librosa.plimit_intervals(primes=[3, 5], bins_per_octave=24, sort=True)
    assert np.allclose(sorted(intervals), intervals_s)
    assert np.allclose(
        intervals,
        [
            1,          # Unison
            3 / 2,      # P5
            4 / 3,      # P4
            9 / 8,      # major wholetone
            5 / 4,      # Pt M3
            15 / 8,     # Pt M7
            5 / 3,      # Pt M6
            45 / 32,    # Pt TT
            8 / 5,      # Pt m6
            6 / 5,      # Pt m3
            16 / 15,    # major diatonic semitone
            9 / 5,      # Pt m7
            16 / 9,     # Py m7
            27 / 16,    # Py M6
            10 / 9,     # minor wholetone
            64 / 45,    # Pt dim5
            32 / 27,    # Py m3
            25 / 16,    # Pt aug5
            75 / 64,    # Pt aug2 
            25 / 24,    # minor chromatic semitone
            135 / 128,  # major limma
            225 / 128,  # Pt aug6
            40 / 27,    # Pt narrow 5
            25 / 18,    # Rameau's tritone
        ],
    )