Using fitted FEM to solve the interface problem, how to calculate the curve integral？

[dokken]

I don’t understand where you are trying to enforce the jump discontinuity of u in your code.

I would strongly suggest looking at

Imposing a discontinuity at interface using DG method

Hi all, This is a simpler version of the problem stated in Imposing a discontinuity in the solution, on an internal interface I’d like (in dolfinx) to impose a discontinuity at the interface between two domains. I understood that MeshView would be the best option since we would have CG elements everywhere except at the interface. However it’s not yet implemented in FEniCSx I’d like to give the DG method a go! This is a simple poisson problem and I’d like to enforce the following equation at…

as your formulation will surely be of a similar flavor, either with a pure DG method, or with the mixed approach I presented in DOLFINx