Odd bugs with Constructive Solid Geometry 2D

I've noticed some odd behavior when generating 2D meshes using Constructive Solid Geometry. Consider the following code, which should generate the intersection of two overlapping circles:

We instead get a single circle centered at at the origin. The same thing happens with a union. An intersection yields an empty mesh; it's like the two circles are the same. If we shift the center a small amount, the code works as expected.

I've done this using a conda install on my Macbook Pro as well as with an apt-get install on Ubuntu 20.04 LTS, both with the same result.

Hello,

I cannot reproduce the issue with Netgen 6.2.2009-69-gea7f6c1e.
Which version are you running? If it's older than 6.2.2009, I suggest an update.

Best,
Matthias

matthiash wrote: I cannot reproduce the issue with Netgen 6.2.2009-69-gea7f6c1e.
Which version are you running? If it's older than 6.2.2009, I suggest an update.

Thanks for the response. I am running 6.2.2009 on both my Mac and on my Ubuntu server. I've converted a Jupyter notebook to a webpage so you can see my version and full code here:

marksmath.org/ProblematicMeshIllustration.html

Here's a similar problem. To produce this, I set up a virtual machine running Ubuntu 20.04 LTS on Digital Ocean. Python version is 3.8. I then logged in as root and followed almost exactly the Ubuntu instructions on this page. - thus:

(Since I was root, I didn't need sudo.)

I then used my favorite, console based text editor to create a file named "ng_bad_script.py" with the following content:

If I run this script, meshing then fails:

The situation is worse if I delete the point (0.5,0.5) from the Solid2d object and store the program in "ng_really_bad_script.py", for then I enter an infinite loop:

Note that loclines.Size=1 printout repeats until I terminate.

Again, this is a totally fresh install and the exact same thing happens on my Mac with an Anaconda based installation.

To be clear, though, I'm generally quite impressed with the software and I'm looking forward to using it to generate examples for the PDE class that I'll be teaching this semester. Examples where the points are in generic position tend to work; thus, problems like the above tend to go away if the points are perturbed a bit.

Thanks for the detailed information. I discovered the reason for this bug in the source code. I will let you know, when it's fixed.

Best,
Matthias

The issues are fixed on the master branch. Nightly builds will be available tomorrow.

Best,
Matthias