FFT - faq

From

Evgeny Sharandin (2:5020/755.12)

To

Nick Poroshin

Date

2002-12-13T02:06Z

Area

RU.ALGORITHMS

Reply-To: shar@nep.cplire.ru

Привет Nick!

10 декабря 2002 года (а было тогда 20:40)
Nick Poroshin в своем письме к All писал:

 NP> Пpедлагаю подпpавленный ваpиант-ляпы подпpавлены и сделано, чтобы
 NP> массивы начинались с 0(а не 1)-так имо удобнее.

 NP> Кстати, он медленнее ооуpовского cdft в 1.5-1.6 pаза, а его
 NP> машинно-оптимизиpованная мной веpсия - в 1.1-1.2 pаза. На это можно
 NP> посмотpеть двояко. С одной стоpоны - медленнее значит медленнее. С
 NP> дpугой - 1.2 не так уж много

В качестве теста использовалась задача дифракции монохроматического пучка света
на диафрагме (fft, домножение спектра на параболу, rfft). Размер массива
8388608 точек. Athlon-XP1500+. Компилятор gcc 3.20 с ключами

-O3 -fno-force-mem -mcpu=athlon-xp -march=athlon-xp -fssa
-momit-leaf-frame-pointer

Ооуровский справился за 9.8с с максимальным отклонением полученного решения
относительно аналитического 2.6e-5 %.

Предлагаемый - 24.7с и 7.3e18 %, соответственно. Вторая цифра совершенно
неприемлима ;).

 NP> пpотив соотношения pазмеpа исходников ~2kb/~40-80kb (т.к. иногда
 NP> большой объём исходников нежелателен)

Так в эти 80К засунуто 6 вариантов ft. Хотя действительно, объем исходников
заметно больше.

С уважением, Evgeny                           13 декабря 2002 года

---
 * Origin: LID (2:5020/755.12)