('heat_gaussian/solution.pvd') Inside the time loop, we may then append the solution values to this file: vtkfile << (u, t) This line is called in each time step, resulting in the creation of a new file with suffix .vtu containing all data for the time step (the mesh and the ...