c++ solver for 2d diffusion equation