Быстрая свертка дискретных сигналов, как? Опции

[phantom]

Надо написать алгоритм КИХ фильтрации. Особенность процесса такова, что сигнал выдается фрагментами (начиная от 64 отсчетов).Светрку реализовили через БПФ в нескольких вариантах (overlap&save, overlap&add), но необходимой производительности достичь не удалось (получается в 2-8раз медленнее чем у существующих программ). Причем у программ-аналогов загрузка процессора практически не зависит от длины импульсной характеристики, а у нас - чем больше длина ИХ, тем больше времени требуется на обработку. Думали может БПФ медленное, заменили на Intel FFT которое в 10 раз быстрее общепринятых алгоритмов, но это существенно не помогло. Осталось попаробовать метод блоков переменной длины, но не совсем понятно как это работает. Что можно было бы сделать?

[DRUID3]

Надо написать алгоритм КИХ фильтрации. Особенность процесса такова, что сигнал выдается фрагментами (начиная от 64 отсчетов).

...ну и что? Вы что, после каждого приращения буфера на 64 сэмпла делаете новое FFT на 4096(допустим) отсчета?

Светрку реализовили через БПФ в нескольких вариантах (overlap&save, overlap&add), но необходимой производительности достичь не удалось (получается в 2-8раз медленнее чем у существующих программ).

А откуда Вам знать, что там именно так сделано?

Причем у программ-аналогов загрузка процессора практически не зависит от длины импульсной характеристики, а у нас - чем больше длина ИХ, тем больше времени требуется на обработку.

Даже не знаю, что сказать ...

Думали может БПФ медленное, заменили на Intel FFT которое в 10 раз быстрее общепринятых алгоритмов...

Из цикла "аФФтАр пЕши ИСчО" ...

Забыл добавить - ИХ длинная >130тис. Выиграш есть конечно по сравнению с простыми алгоритмами - но все равно не то.

На 130 000 выигрыш FFT очень ощутим. Вот получили свои 130 000 и сделали один раз FFT...

[phantom]

[quote name='DRUID3' post='711474' date='Jan 31 2010, 17:42']...ну и что? Вы что, после каждого приращения буфера на 64 сэмпла делаете новое FFT на 4096(допустим) отсчета?

А откуда Вам знать, что там именно так сделано?
>>>>>>>>>А я и не знаю. Я только проверяю время обработки для своего алгоритма и программ-аналогов.


Из цикла "аФФтАр пЕши ИСчО" ...
>>>>>>>>> Ну это я не понял, но так для общей инфы: FFT которые описаны в общедоступной литературе, проигрывают в производительности некоторым оптимизированным алгоритмам, например FFTW, либо Intel MKL от 8 до 15 раз.

а Вы CUDA попробуйте.
>>>>> да я в курсе, но надо на СPU а не на GPU