Jørgen Schartum Dokken

Note to self, test on DOLFINx output

For this to work on codim-1 meshes, I believe that the way forward is to: 1. Expose `dolfinx::fem::CoordinateElement::permute_subentity_closure` to be able to map the input geometry dofmap correctly (https://github.com/FEniCS/dolfinx/blob/main/cpp/dolfinx/fem/CoordinateElement.h#L103-L123) 2....

I would probably refer to: http://jsdokken.com/FEniCS23-tutorial/src/ufl_formulation.html#functionals-and-derivatives which goes through the functional version of linear elasticity.

I get that you want that included in this demo. However, the functional principle is the same for both linear and the non-linear problem: 1. Formulate an energy functional 2....

Do you mean in equation 10? Is this what you would like: $a(u, v) = \int_{\Omega} \nabla u \cdot \nabla v~\mathrm{d}x + \int_{\partial\Omega}\left\[-(n \cdot\nabla u) v - (n \cdot \nabla...

In the case above, we are constraining a single constant value to an interface. Then we only supply a single space in the `locate_dofs_topological` function. When we want to use...

@tyagi-ankit I've now added a proposed text to correct this. It refers to the more detailed discussion I've made in https://jsdokken.com/FEniCS-workshop/src/deep_dive/mixed_problems.html#dirichlet-conditions-in-mixed-spaces

Following up on this, this means that no package can really use `adios2`'s Python bindings and be a functioning spack package (without the workaround of Jack).

Fixed with https://github.com/spack/spack-packages/pull/880

Further reduced to: ```python import numpy as np # p0 = np.array([[44.2919, 39.7508, 18.8208]]) p0 = np.array([[44.2919, 39.7508, 18.8208]]) p1 = np.array( [ [46.0659, 38.6251, 16.3637], [44.2764, 42.3867, 18.239], [50.0397,...