Реализация КИХ-фильтра с разрядностью 32 бит, на Blackfin и С Опции

[K A A]

Пишу программу вычисления действующего значения сигнала с помощью КИХ-фильтрации.
Алгоритм следующий: 16-бит выборки возвожу в квадрат, получаю 32 бит, затем умножаю на весовой коэффициент (32 бит) и прибавляю к аккумулятору. При такой разрядности нужно использовать все 40 бит аккумулятора MAC, и то этого мало. Но никак не пойму, как это сделать на С. Функции из fract_math.h реализуют умножение 32x32 -> 32 и сложение 32+32 -> 32, как обрабатывать дополнительные 8 бит? Или писать на ассемблере?

ЗЫ: почему такая разрядность - просто погрешность результата должна быть на уровне 0.001%. Почти получается, когда алгоритм был реализован с использованием вещественной арифметики. Для улучшения стабильности нужно увеличить размер окна в несколько раз (сейчас 4096, хорошо бы раз в 16 больше, путем интерполяции на лету), соответственно увеличить частоту - тут уже проц не поспевает. Еще возможные проблемы с float при большом окне - ограничение разрядности, когда сумма накопится большая, а слагаемые очень маленькие, при этом они будут игнорироваться.

Так как это сделать на С?

[fontp]

Вообще-то их и делать особенно не нужно. Так подправить...
Взять исходники библиотечных функций fir_fr16.asm или аналогичной, поменять атрибуты инструкции делающей mac, добавить выходную нормировку. Там почти везде два мака на такт