E = field("hcurl", [2,3])
E_re = E.harmonic(2)
E_im = E.harmonic(3);
(I think) the primal and dual parts are E and J. So I defined two ports
Ep = port()
Jp = port()
But the fact of the matter is I want to set them only in the z direction. And this is illegal in Sparselizard:
# Include excitation part of the integral
maxwell += integral(feed,-(1/C0)*sigma*dt(dof(E))*tf(E))