[Octave-bug-tracker] [bug #39566] Suspect interpft is wrong regarding treatment of Nyquist frequency

[Mike Miller]

Update of bug #39566 (project octave):

                  Status:               Confirmed => In Progress            
             Assigned to:                    None => mtmiller               

    _______________________________________________________

Follow-up Comment #2:

I've got a patch for both problems you described, please test the attached if
you are able to and let me know if this does what you expect.

I'd like to add some test cases to this change as well but haven't done those
yet.

Also if you have access to Matlab, which it seems you do, it would be helpful
to post some examples of interpft on complex input, since the Nyquist bin will
not be purely real and the FFT is not conjugate symmetric. For example, what
do you get with


> fft (interpft ((1+j)*[1:4]', 6))


(file #28731)
    _______________________________________________________

Additional Item Attachment:

File name: interpft.diff                  Size:0 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?39566>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/