On the implementation of linear finite element method for nonlocal diffusion model over 2D domain

We propose an implementation of linear finite element method for nonlocal diffusion problem in 2D space. In the implementation, we reduce the integral from 4D to 2D which would simplify the computation significantly.