random number generator seeding should be better
Added one step of diffusion equation
Merged branch from intensity-diffusion