Using dolfinx_mpc in generalised plain strain formulation

dolfinx
4

Depending on your DOLFINx version, you need to:

Stable/0.9.0

  1. Load your geometry as usual (you have omitted this part of the code)
  2. Create the two function spaces V, R
  3. Create MPC on V
  4. Use the assembly procedure from dolfinx_mpc/python/demos/demo_stokes_nest.py at v0.9.2 · jorgensd/dolfinx_mpc · GitHub
  5. Create your ksp object as in: dolfinx_mpc/python/demos/demo_stokes_nest.py at v0.9.2 · jorgensd/dolfinx_mpc · GitHub

Main branch
Step 1-3 are the same.
4-5. Use dolfinx_mpc.LinearProblem or dolfinx_mpc.NonLinearProblem as described in: dolfinx_mpc/python/demos/demo_stokes_nonlinear_nest.py at main · jorgensd/dolfinx_mpc · GitHub