Could different Functionspaces be defined on different subdomains in one mesh?

[dokken]

For legacy dolfin. See MeshView and the corresponding paper and examples:

Abstractions and Automated Algorithms for Mixed Domain Finite Element Methods...

GitHub - cdaversin/mixed-dimensional-examples: Code to reproduce numerical examples presented in...

For similar abstractions in DOLFINx, see for instance

Discontinuity at interface using mixed domains

I think this is the variational form that you would like: # SPDX-License-Identifier: MIT from mpi4py import MPI import dolfinx import dolfinx.fem.petsc import ufl import numpy as np from petsc4py import PETSc class NewtonSolver: max_iterations: int bcs: list[dolfinx.fem.DirichletBC] A: PETSc.Mat b: PETSc.Vec J: dolfinx.fem.Form b: dolfinx.fem.Form dx: PETSc.Vec def __init__( self, F: list[dolfinx.fem.form], J: list[list[dolfinx.fem.for…