16.09.2011, Изучение быстродействия и оптимизация алгоритма БПФ
Korogodin (обсуждение | вклад) |
Admin (обсуждение | вклад) |
||
| (не показана 1 промежуточная версия 1 участника) | |||
| Строка 1: | Строка 1: | ||
| − | |||
<summary [ hidden ]>Попытки ускорить БПФ на ARM</summary> | <summary [ hidden ]>Попытки ускорить БПФ на ARM</summary> | ||
| Строка 69: | Строка 68: | ||
== См. также == | == См. также == | ||
* [http://fpga.parallel.ru/fft/ Основная схема быстрого преобразования Фурье] | * [http://fpga.parallel.ru/fft/ Основная схема быстрого преобразования Фурье] | ||
| + | |||
| + | [[Категория:НИИ КП]] | ||
{{wl-publish: 2011-09-16 14:37:37 +0400 | Korogodin }} | {{wl-publish: 2011-09-16 14:37:37 +0400 | Korogodin }} | ||
Текущая версия на 02:44, 18 марта 2013
Под исходные коды заведен проект fft-for-arm-search.
Исследуется БПФ от 2048 точек.
Согласно изученной литературе, минимальное количество операций достижимо для БПФ размером 
("четной степени двойки"). Для него:
- число комплексных умножений;
- число комплексных сложений.
Остаются два рычага власти: варьировать N и менять разрядность и тип переменных в функции (уменьшать время умножения и сложения).
По произведенным оценкам ARM'у на одно комплексное умножение (float+jfloat)*(float+jfloat) надо около 2.5-2.8 мкс (при расчете пришлось принять гипотезу исчезающе малого влияния суммирования на время вычисления). Отсюда, в изначальном варианте оценка свертки мс-ого сигнала для 10 частот - 300 мс.
| Ревизия | ARM, ms | Pentium, ms | Расположение в приемнике | Примечание |
|---|---|---|---|---|
| 2 | 23.3 | 0.3 | Внешнее ОЗУ | Исходный алгоритм с коэф, вх. и вых. данными во float |
| 0.3 | Внешнее ОЗУ | Входные данные int: изменений не замечено | ||
| 3 | 84.5 | 0.3 | Внешнее ОЗУ | Все float переведены в int. Грубая оценка быстродействия целочисленного алгоритма. |
| 4 | 66 | - | Внешнее ОЗУ | Все int переведены в int16. Грубая оценка быстродействия целочисленного алгоритма с int16. |
| 5 | 83 | - | Внешнее ОЗУ | Все int переведены в int16, но значения ограничены 8 битами. Грубая оценка быстродействия целочисленного алгоритма с int16, но маленькими числами. |
Перенос программы и данных в TCM выигрыша в быстродействии не дал.
[править] Бонус
In this figure we compare the performance in MFLOPS of a complex-complex 5*n*log2(n)/t, where n is the length of the sequence and t is the time for one transform. The computations are done on one processor of the IBM Winterhawk II with 375-MHz Power3-II processors. We compare FFTPACK 5 by P. Swarztrauber, NCAR (FFTPACK), the new Spectral Toolkit by Rodney James, NCAR (STK), FFTW by M. Frigio and S. Johnson, MIT (FFTW), and the IBM ESSL Library (ESSL).
[ Хронологический вид ]Комментарии
Войдите, чтобы комментировать.