Changed the ADI2D to BE2D in diffusion.cpp, the implementation is slow, but should be correct
??