How to refine a mesh (through the script) imported from GMSH?

[dokken]

You can use the refine and adapt commands, as used in:

Mesh refinement gives solution with zero values

I strongly doubt that, as you can clearly see from just printing the markers before and after refinement: from dolfin import * from mshr import * from math import sin, cos, pi a = 1.0 b = 1.2 c_1 = 0.8 c_2 = 1.4 r = 0.1 R = 5.0 n = 10 domain = Circle(Point(0, 0), R) cylinder = Circle(Point(0, 0), b) - Circle(Point(0, 0), a) angles_N = [i*2*pi/n for i in range(n)] angles_S = [(i + 0.5)*2*pi/n for i in range(n)] wires_N = [Circle(Point(c_1*cos(v), c_1*sin(v)), r) for v in ang…