Connect 1D domains

The way to do this is through the usage of submeshes (using dolfinx.fem.create_submesh) and create a mixed variational form for these.

See for instance:

or

(The latter link doesn’t cover the coupling from the first post, but it illustrates the idea with more documentation).