I have encountered an error when requesting multiple harmonic terms from PyNAFF while setting getFullSpectrum=False. The issue seems straightforward, I think the orthanormalization step gramsc returns a complex array which then conflicts with the np.fft.rfft call in the main naff loop without the full spectrum flag. The attached script reproduces the behavior.
pynaffError.txt
I am afraid I do not fully understand the orthanormalization step in naff, and so this issue may be a fundamental mathematical misunderstanding. If that is the case, it may be useful to throw a descriptive warning for this configuration, since it is not immediately obvious that the different spectral lines of a strictly real signal cannot be extracted.
I have encountered an error when requesting multiple harmonic terms from PyNAFF while setting
getFullSpectrum=False. The issue seems straightforward, I think the orthanormalization stepgramscreturns a complex array which then conflicts with thenp.fft.rfftcall in the main naff loop without the full spectrum flag. The attached script reproduces the behavior.pynaffError.txt
I am afraid I do not fully understand the orthanormalization step in naff, and so this issue may be a fundamental mathematical misunderstanding. If that is the case, it may be useful to throw a descriptive warning for this configuration, since it is not immediately obvious that the different spectral lines of a strictly real signal cannot be extracted.