mshr with overlaping subdomains missing edges

Question

Hi,

In the following code the vertical subdomain is missing some edges. It gives the expected mesh if the size of r2 in the x direction is greater than 0.1. For smaller values the edges of r2 inside r1 disappear. Any help is appreciated.

Thanks,
Chaffra

from dolfin.cpp.mesh import Point
from dolfin import parameters, plot

from mshr import generate_mesh, Rectangle
import numpy as np


r1 = Rectangle(Point(0.0,0.0), Point(2.0,1.0))
r2 = Rectangle(Point(0.0,0.0), Point(0.1,2.0))

geometry = r1+r2
geometry.set_subdomain(1,r1)
geometry.set_subdomain(2,r2)
mesh = generate_mesh(geometry, 1.0, 'cgal')

plot(mesh)