Hi
Consider this:
from dolfin import *
from mshr import *
#Making a cylindrical geometry
cylinder = Cylinder(Point(0, 0, -1), Point(0, 0, 1), 0.5, 0.5)
geometry = cylinder
# Making Mesh (40 corresponds to the mesh density)
mesh = generate_mesh(geometry, 40)
# Save the mesh
File("cylinder.pvd") << mesh
Here is the mesh: