Automatic order selection #98
Description
We could select the harmonic order automatically with an information criterion: fit models with different orders, then compare the model fit and choose one that is simplest. This would slow down runtime but might help against overfitting.