I am using DG for a two-domain problem.

The domain is a small square inclosed in a large square.

The question is when I do symbolic integration

a += symbolicBFI(u*v, VOL, skeleton=True)
does this skeleton integrator include the integration on the interface of the two domains?
It seems I can integrate on the interface boundary as follows:
a += symbolicBFI(u*v, skeleton=True, definedon=mesh.Boundaries(“interface_id”))

I am quite confused…