Add diffusion as test case

* Pythran performs well on the pure python code but not on the numpy version.