1D линейная свертка в коде ANSI C?

Вместо того, чтобы изобретать колесо, мне интересно, может ли кто-нибудь направить меня к фрагменту кода 1D линейной свертки в ANSI C? Я искал в google и в stack overflow, но не смог найти ничего на C, что я мог бы использовать.

Например, для массивов A, B и C, все с двойной точностью, где A и B - входы, а C - выход, имеющих длины len_A, len_B и len_C = len_A + len_B - 1, соответственно.

Размеры моих массивов невелики, поэтому никакого увеличения скорости при реализации быстрой свертки с помощью БПФ не требуется. Ищу прямое вычисление.