Hang (?) in constructive solid geometry support

[Patrick Farrell]

Consider the following code:

from dolfin import *

r = 1.0/3.0
info_blue("Creating domain")
domain = Rectangle(0, 0, 1, 1) - Circle(0, 0, r) - Circle(0, 1, r) - Circle(1, 0, r) - Circle(1, 1, r)

info_blue("Meshing ...")
mesh = Mesh(domain, 1)
plot(mesh, interactive=True)

I expected this to work near-instantaneously, but it appears to hang; it's been running for several hours without completing.

This is with

 *** 1.2.0+git20131015~ppa1~precise1-1 0
        500 http://ppa.launchpad.net/fenics-packages/fenics-dev/ubuntu/ precise/main amd64 Packages

[Benjamin Dam Kehlet]

This is fixed in this branch: https://bitbucket.org/benjamik/dolfin/branch/benjamik/csg-2d

(This branch also adds support for subdomains in 2D CSG which will be marked the resulting mesh). However, apparently Dolfin is not capable of distributing the mesh data in parallel. It fails with this message: ------------------------------------------------------------------------- Error: Unable to extract local mesh data. Reason: Marked subdomains are not yet supported. Where: This error was encountered inside LocalMeshData.cpp. Process: 0 -------------------------------------------------------------------------

Any chance this can be fixed soon? (Anders?)

[Benjamin Dam Kehlet]

This should be fixed in master now.