Difference in `fft` implementation between stan and numpy/scipy

Thank you for the link. They mention scaling which should not affect the results up to a global rescaling. The other difference is how they handle real FFTs where only half the spectrum is required. Neither should lead to the differences observed above, I think.

1 Like