Allocating different variables to different subdomains

[dokken]

Using @jpdean’s implementation of mixed assembly this is very doable.

I have an example of this kind of setup at:

Wells-Group/TEAM30/blob/main/team30_A_phi.py#L150

Although I do admit that it is not the easiest example to start with. There will be more demos of this in the future.

I have also covered a similar example at:

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…