dokken
June 27, 2023, 8:05pm
2
First of all, please format all code with
```python
#add code here
```
Secondly, you would have to either supply the xdmf-file and the binary h5 file, a pure ascii xdmf file, or the original mesh file, and the functions that have been used to convert it to xdmf.
Also see:
Oh that was my bad. I had not read the previous answers. As @dokken mentioned in a previous reply, since you don’t have any facet regions with a physical tag:
mvc = MeshValueCollection("size_t", mesh, 1)
with XDMFFile("mf.xdmf") as infile:
infile.read(mvc, "line")
mf = cpp.mesh.MeshFunctionSizet(mesh, mvc)
This part of the code would not be required as it is basically reading the tagged boundaries into the object mf using the MeshValueCollection and MeshFunction functions.
Please remove …