First working iteration, hardcoded Gaussian kernel

Need to refactor convolution code; allow customizing Gaussian kernel
radius