О функции FFT

Question

Кто-нибудь знает, какой алгоритм используется в Юлии для выполнения быстрого преобразования Фурье?Документация гласит только:

...
A one-dimensional FFT computes the one-dimensional discrete Fourier transform (DFT) as defined by

\operatorname{DFT}(A)[k] =
  \sum_{n=1}^{\operatorname{length}(A)}
  \exp\left(-i\frac{2\pi
  (n-1)(k-1)}{\operatorname{length}(A)} \right) A[n].
...

В частности, у меня есть расхождение в моих преобразованных данных, то есть эти преобразованные данные «сдвинуты», как мне кажется, pi.Существует ли соглашение об исправлении этой глобальной фазы?

РЕДАКТИРОВАТЬ: Возможно, стоит сказать, что если я выполняю обратное FFT, то несоответствие в фазе исправляется.

Answer 1

Я полагаю, что Джулия использует библиотеку FFTW, которая использует несколько вариантов алгоритма Кули-Тьюки, как описано в приведенной ниже ссылке.

http://www.fftw.org/fftw-paper-ieee.pdf