#555

Request to merge #555 into the master branch.
It's generally a bad idea to compare floats with `==` or `!=` due to the likelihood of slight rounding errors causing things to not be equal.  Probably this should be changed to something like `if abs(IPC_kernel.sum() - 1.0) > 1.e-8:` to allow for numerical rounding errors.