A class of evolution equations with nonlocal diffusion is considered in this work. These are integro-differential arising as models propagation phenomena continuum media interactions including neural tissue, porous media, peridynamics, and fractional diffusion, well limits interacting dynamical systems. The principal challenge numerical integration systems stems from the lack spatial regularity...