It seems like some versions of Fenics do not support interactive() command anymore, also there are some issues with vtkfile creation, e.g.
u = Function(V) t = 0 for n in range(num_steps): t += dt u_y.t = t solve(a == L, u, bc) #plot(u) vtkfile = File("test.pvd") vtkfile << (u, t) u_n.assign(u) #interactive()
This creates two files pvd & vtu files, when I have opened pvd file in ParaView, then some errors come up, as well as no results can be visualized; when opening pvd file, it looks like xml file
<?xml version="1.0"?> <VTKFile type="Collection" version="0.1"> <Collection> <DataSet timestep="1.9999999999999998" part="0" file="test000000.vtu" /> </Collection> </VTKFile>
Has anyone faced such issue before ?