Discrete convolution (source code)

= Discrete convolution

For sequences $a,b$, their discrete convolution is
$$
(a*b)_n=\sum_m a_{n-m}b_m.
$$