I'd like to solve two pdes on two different meshes that are coupled using boundary conditions. Right now I have an iterative process set-up to alternate between equations , but I think solving both equations simultaneously would be a superior method (and is something that I know works in other FEM software).
Does anyone have an idea how to implement it in FEniCS?
Thanks