Jørgen Schartum Dokken

The `subspace_slave` and `subspace_master` is just for convenience for cases where you can isolate the components. If you do not set these, it will work with dofs on the full...

You would need to code each coefficient for each degree of freedom separately, i.e. for each (x1, y1) we have a slave `dof_0(x1, y1)_x`, masters dof_1(map(x1, y1))_x, dof_1(map(x1,y1))_y, we would...

Maybe the `create_general_constraint` is not the right abstraction for your problem (I need to think about it) and revisit the code. The most general way of adding a constraint is...

[result.bp.zip](https://github.com/user-attachments/files/15573779/result.bp.zip) Here is a zip file of an even smaller mesh that exhibits the same behavior. Here the mesh is only 1x1 unit square with two triangles, with ```bash bpls...

The 69 cell tag is the arbitrary order Lagrange cells. We use them to write our meshes to file, as we support arbitrary order geometry. Ref: ``` // Arbitrary order...

> I dove into the code for the reader. It looks like there is indeed an issue with reading cell data that changes over time. For some reason, it is...

I think the issue is this around these lines ( https://github.com/FEniCS/dolfinx/blame/f1c5f307c70f414872a090e349f1d88211f54b58/python/dolfinx/fem/forms.py#L369), as it seems like the order of the integrals here are different that the one in the generated code....

So I think we can simplify the Python code to: ```python # Extract subdomain ids from ufcx_form integral_offsets = [ufcx_form.form_integral_offsets[i] for i in range(4)] for i in range(3): integral_type =...

Pinned down to being an ADIOS2 issue: https://github.com/ornladios/ADIOS2/issues/4287

@minrk, do you have any idea? (Being the ipyparallel wizard!)