Geometry with a subdomain not meshed correctly
A Rectangle geometry with a Circle subtracted is not meshed correctly when a subdomain is set.
In one case the mesh is as expected (see "good.png"), but changing the size and location of the Rectangle and the Circle generates an erroneous mesh (see "bad.png").
The subdomain extends outside the domain in both cases, is that a problem? No warning is given.
Without the subdomain, the mesh is correct.
See attached files.
Comments (5)
-
-
I'm having trouble reproducing this. What version of mshr are you using? I get this with the current dev version.
-
reporter I was using the Docker image jjan/fenics-mooc:test (https://hub.docker.com/r/jjan/fenics-mooc/tags/) and I see that it uses dolfin 1.6.0 and mshr version 1.4.0.
I tried with the Docker image quay.io/fenicsproject/stable:1.6.0 instead, which uses dolfin 1.6.0 and mshr 2016.2.0.dev0. Then it works.
Maybe the problem is solved for mshr version > 1.4.0.
Sorry for this.
-
No problem. Yes, 1.4.0 is pretty old. A lot have been fixed since then.
-
- changed status to invalid
- Log in to comment
Hi Frida. Thanks for reporting. I am looking into it now.
No, it is not a problem that the subdomain geometry is not completely within the domain. The domain still defines the region that is meshed (so the area marked is the intersection of the domain and the subdomain)