Parallel numerical simulations of anisotropic and heterogeneous diffusion equations with GPGPU