My problem involves ufl.MixedFunctionSpace with individual function spaces defined on multiple submeshes. I need to use create_contact_inelastic_condition to impose continuity on a non-conformal interface found in one of the submeshes. I noticed that dolfinx_mpc.problem.LinearProblem does not admit entity maps in the arguments. Adding it by hand and updating the relevant _fem.form operations does not fix the problem. I still get an error related to MPI issues. I would appreciate any advice in case the functionality is there but I am not using it right or if it is presently impossible to make it work for my problem.
It will soon be compatible, ref: Couple constraints across nest sub matrices/functionspaces(meshes) by jorgensd · Pull Request #263 · jorgensd/dolfinx_mpc · GitHub
@sbhasan it is now supported on the main branch ; Accept submeshes and entity maps coupling by jorgensd · Pull Request #263 · jorgensd/dolfinx_mpc · GitHub
Thank you much! Is there a possibility to use it with the current 0.11 or it relies on features exclusive to dolfinx on the main branch?
I’ve not had the bandwidth to port it back to 0.11.
I’ll hopefully post a docker image very soon with DOLFINx and DOLFNx_MPC with this support that you could use.