Is your feature request related to a problem? Please describe.

The VQT function (which powers cqt and the like) recursively downsamples the input signal by factors of 2 to compute one octave worth of frequencies at a time at the lowest feasible sampling rates (subject to integral hop length constraints, etc):

There's an assumption in CQT/VQT that all octaves have a fixed number of bins (except maybe the first/last), but this really isn't necessary.

Describe the solution you'd like
We could instead compute a set of feasible sampling rates, and then divide up the basis frequencies into whatever sampling rate range is most reasonable. This change should be invisible to the user for the currently supported use cases.

However, restructuring the code in this way would relax the constraint that each octave has an equal number of frequencies, which in turn should allow for basically arbitrary frequency sets. This, I think, would fully resolve #18 .