How to define spatially varying Young’s modulus from a NumPy array in FEniCSx?

[dokken]

Please reduce the amount of code to a minimal reproducible example, as this has way too much information in it.

It seems like what you want to do is to read in some vertex data from file, and attach it to each of the vertices of your mesh?

This is for instance covered in:

How to apply a boundary condition to a range of DOFs?

Do I really need to re-read your code and compare it to the one I created above? Why not simply use the code I supplied, i.e. the function: There is a reason for me adding all the steps above, hich is that you need to go through: Of course, if you have a higher order mesh, and want to map nodes from say a midpoint of a facet to the dof at a midpoint of a facet, something slightly different should be done. """ Compute local dof indices of global input node indices Author: Jørgen S. Dokken …