Stiff matrix for interior/exterior nodes

Hello, I am trying to implement a research paper (https://arxiv.org/pdf/2310.19571) where I need to separately compute the stiffness/mass matrix for nodes which are inside the domain and for nodes which are on the boundary, respectively.

Any help is more than appreciated!

https://multiphenics.github.io/ can help in doing that, see for instance tutorial 3.

dolfinx 0.9.0 has intregated support for mixed domain assembly FEniCS release notes | FEniCS Project , so that would be another route.

1 Like