Add 2D conv

[ericphanson]

It’s tricky to do efficiently so we should have an implementation here (even though it doesn’t need a new atom). We should also speed up the sparse matrix creation in the 1D convolution.

I implemented both here: https://discourse.julialang.org/t/convex-jl-slow-in-creating-constraints/34775/31

We could also let users reuse the matrix representation as is done there.

[RoyiAvital]

The paper link you used is no longer available.
This is its archive: https://web.archive.org/web/20170809080755/https://sites.ualberta.ca/~mostafan/Files/Papers/md_convolution_TLE2009.pdf.