How to import xdmf mesh generated by gmsh into fenics

[dokken]

As i do not know what meshing tool you use to go from an image to a mesh, it is hard to give you any specific tips.

If you use Gmsh for this, there might be a tutorial on their webpage:
https://gmsh.info/

@taytay33

2nd option is quite complex.

The second option doesn’t take many lines of Python to program. See for instance:

How to create lists of cell neighbours?

Instead of using the filter function, you can simply unwrap it in multiple loops: from dolfin import * mesh = UnitSquareMesh(2, 2) # Init facet-cell connectivity tdim = mesh.topology().dim() mesh.init(tdim - 1, tdim) # For every cell, build a list of cells that are connected to its facets # but are not the iterated cell cell_neighbours = {} for cell in cells(mesh): index = cell.index() cell_neighbours[index] = [] for facet in facets(cell): facet_cells = facet.entities(tdi…