Can I do interpolate expression for Nedelec element?

Hello I want to ask a question about the Nedelec element

Let us consider the Gauss’s law in 2D:

\pmb{E} = \nabla\rho

If I define the term \rho in a H1 space, then it naturally holds true that \pmb{E} should be in a H curl space. So I define my function spaces as follows:

H1:  "Lagrange", "quadrilateral", 1, LagrangeVariant.equispaced)
Hcurl: "Nedelec 1st kind H(curl)", "quadrilateral", 1, LagrangeVariant.equispaced)

My question is: If I have the function \rho, can I interpolate \nabla\rho and get \pmb{E} directly. For normal Lagrange element, I know that I can do it by defining expressions but this requires tabulating the dofs in reference element; whille for Nedelec element, is it still possible? I tried to use tabulate_dofs_coordinates, but get the error:

 Cannot evaluate dof coordinates - this element does not have pointwise evaluation.

Does anybody have any suggestions? Thanks in advance.

Expressions require the interpolation_points, not the tabulated dof coordinates.
This is for instance illustrated in: dolfinx/python/test/unit/fem/test_interpolation.py at fb0b99447b96434a48b44c4a256c29330ffac7c6 · FEniCS/dolfinx · GitHub