исходник свёртки

google fftw.c

Мне не fft, мне convolution. гугление по convolution/convolver source/src ничег оне дало.

Свертку можно делать через FFT (Fast convolution). А самый простой пример вероятно так:

def conv(d, h):
  res = []
  for i in range(len(d)-len(h)):
    s = 0
    for j in range(len(h)):
      s += d[i+j]*h[j]
    res.append(s)

  return res

print conv(range(100), [1,2,3,4])

Не проверял, могут быть ошибки.

Что-то я не понимаю, что значит "перемножение двух массивов". Массивы обычно содержат какие-то определённые данные, например коэффициенты полиномов или координаты векторов (это если рассматривать наипростейшие случаи), и для разных типов данных будут разные алгоритмы перемножения.

Сигнал и импульсный отклик системы. Запись звука и реверберация помещения, если угодно.

> s += d[i+j]*h[j]

А не i-j? Там же второй массив в обратном порядке перемножать надо...

> А не i-j? Там же второй массив в обратном порядке перемножать надо...

а там от оператора зависит. А вообще там прям по теории простейший цикл записывается.