Can I plot the obtained H(curl) and H(div) basis function on the edges of the cells?
For example,
gdim = 2 # Geometric dimension
shape = “triangle”
degree = 1 # Degree of the polynomial for the Lagrange element
cell = ufl.Cell(shape, geometric_dimension=gdim)
domain = ufl.Mesh(ufl.VectorElement(“Lagrange”, cell, degree))
x = np.array([[0.0, 0.0], # Vertex 0
[0.0, 1.0], # Vertex 1
[1.0, 0.0]]) # Vertex 2
cells = np.array([[0, 1, 2]], dtype=np.int32) # One triangle composed of the above three vertices
mesh = dolfinx.mesh.create_mesh(MPI.COMM_WORLD, cells, x, domain)
V = FunctionSpace(mesh, (“Nedelec 1st kind H(curl)”, 1))
Can I plot and obtain the value of trial function as a numpy array at any point?
Thank you