I just went through the Fenics book to get the nodal displacements and tried this code ....u_nodal_values = U.vector()u_array = u_nodal_values.array()coor = mesh.coordinates()center = (0.5,0.5,0.5)print U(center)if mesh.num_vertices()==len(u_array):   for i in range(mesh.num_vertices()):        print U(%8g,%8g,%8g) = %g% (coor[i][0], coor[i][1], u_array[i])but i don't get the expected result..can anyone suggest me the modification ??
Fenics : looking forward to compute displacements at nodes
python;fenics
I think the problem is that you are relying on a specific DOF ordering, which was changed in release 1.1. However, the vertex to DOF map which you need to solve your problem can be obtained by the following codemap = V.dofmap().vertex_to_dof_map(mesh)