Hi,
I am trying to use FEniCS to create the nodal interpolation operator $\Pi_h$ used to in a multigrid method for $H(curl)$ problems. The operator takes vector nodal functions to edge functions (Nedelec functions). The operator is of the form
$$
\Pi_h \textbf{v} = \sum_{e \in E} \int_e \textbf{v}\cdot d\overrightarrow{s}\cdot \textbf{b}_e,
$$
where $\textbf{v}$ is a vector nodal function, $\textbf{b}_e$ are edge element basis functions and $E$ is the set of edges.
I would like to be able to create a matrix $P$ that does this operation for an order elements (in the vector nodal and edge space), how would I go about this?
