Heat PDE solver, iterative methods