Things to note
Look at Linear solver
Maybe use mumps
oomph-lib has multi-domain and multi-field approaches for mixed physics
Multifield mix the elements together one mesh
Multidomain each equation on there own mesh and they can (do not have to) overlap.
Foeppl-von Karman equations can now to the stress free case but not merged into the trunk.