Heat solver works in parallel with MPI