I am personally not very familiar with solving PDEs defined on different domains. The functionality for this is implemented in the development version of fenics, with the MeshView feature.
Maybe take a look here: https://arxiv.org/pdf/1911.01166.pdf
and here MeshView and MixedFunctionSpace
Also, I would suggest that you supply a minimal working example (e.g. using your iterative solver)