Creating a matrix of Gaussian wavelets at dyadic scales
In this case the problem was that I was working in the frequency domain, but implementing an algorithm in the time domain. 

I've removed the pre-multiplications of the Fourier matrix in the time domain and everything works fine.