Basic 1D convolution function

convolve.cpp

#include <vector>

template<typename T>
auto convolve(std::vector<T> v1, std::vector<T> v2){
    std::vector<T> res(v1.size() + v2.size() - 1);

    for(int i = 0; i < v1.size(); i++){
        int j = i;
        for(const auto &ev2 : v2){
            res[j++] += v1[i]*ev2;
        }
    }

    return res;
}