For my problem I have two meshes, finer mesh by splitting the coarse one by half.
I want to do an intalledges over the coarse mesh. As one field is defined on the fine mesh, I need the quadrature to be applied in a piecewise way on each half piece of the edge.
mesh Th1 = trunc(Th,1,split=2);
fespace Vh(Th1,P1);
fespace Qh(Th,P1);
Vh fh = x;
varf RHS(p,q) = intalledges(Th,qforder=2)( fh*q );
real[int] rhs = RHS(0,Qh);
The above won’t work properly as fh*q has two pieces on the edge and we should apply the quadrature on each piece. Is there any straightforward way to make this work?