Surface flux calculation doesn't compile

I wanted to compute the RT1 field that has the same fluxes as the centred flux for a P0 finite volume scheme (i.e., a flux reconstruction). This requires to assemble the matrix from the form

\int_{\Gamma} w.n u.n dS, for all w in RT1, and with u in RT1

for the edge space \Gamma. Assembling this form for a sphere mesh fails to compile due to various out of scope variables. I think this is a use case that we didn't check when we did the assembly on manifolds stuff.

Python file and input mesh attached.

Comments (2)