How to plot VectorFunctionSpace with ndim=2 over IntervalMesh

[dokken]

The built in plotting in dolfin is quite limited, and I would suggest you write the data to file and visualize with paraview, or create your own custom visualization extracting the data as explained in:

How to get the solution obtained by fenics as a vector?

Please make a minimal example, as your original question is not clear on what kind of vector you would like to use. What you can do is to get each of the components of a VectorFunctionSpace by collapsing each of its components. from dolfin import * import numpy as np mesh = UnitSquareMesh(10, 10) V = VectorFunctionSpace(mesh, "CG", 1) u = Function(V) u.interpolate(Expression(("3*x[0]", "2*x[1]"), degree=1)) num_dofs_per_component = int(V.dim()/V.num_sub_spaces()) num_sub_spaces = V.num_sub_spac…