3x Speed-up

* Use 1D Gaussian filters, apply twice, along each axis. This is what scipy's gaussian_filter does, and is much faster, yielding the same result
* correctly set mu in relation to gaussian_kernel_width; not always 5
* use -= where appropriate

These changes result in a ~3x performance improvement (on my machine) without any changes to the API or the results.