Issue with Importing Physical Groups into FEniCS

[dokken]

As you have not provided the mesh, it is hard to give any further guidance.
You have also not provided any error message.

See for instance the recent post:

Fail to read tagged facets from mesh generated with GMSH

You are here trying to integrate over the vertices of your mesh? Then you are not integrating over a facet (“ds”) but a vertex (“dp”). However, lets take 10 steps backwards, and consider: import meshio import gmsh gmsh.initialize() gmsh.model.add("my model") c_r = [0, 0, 0] c_R = [0, 0, 0] r = 1 R = 2 resolution = 0.1 print(f"Mesh resolution = {resolution}") disk_r = gmsh.model.occ.addDisk(c_r[0], c_r[1], c_r[2], r, r) disk_R = gmsh.model.occ.addDisk(c_R[0], c_R[1], c_R[2], R, R) ring = gm…

for some guidance.