I have been trying to work with H1Curl elements for a FEEC simulation and ran a simple code to plot the mesh coordinates of the Finite Elements to confirm on the literature and found the following results.
Is the method for mesh coordinates wrong?
The H1Curl and H1Div Elements are giving results for degree=1 but are not giving any results for degree>=2
Expected…
from fenics import *
import matplotlib.pyplot as plt
mesh = UnitSquareMesh(1,1)
El_0form = FiniteElement('P-', cell='triangle', degree=2, form_degree=1)#, form_degree=1)
Q = FunctionSpace(mesh, El_0form)
plt.figure(), plt.grid()
plt.scatter(Q.tabulate_dof_coordinates()[:,0], Q.tabulate_dof_coordinates()[:,1])
plot(mesh)
I am getting a result,
RuntimeError: tabulate_dof_coordinates is not defined for this element
Any help is appreciated.