[Gmsh] Bug? Points in Volume only meshed for first volume

[Christophe Geuzaine]

Dear Gregor,

Indeed, this was only supported with our new Delaunay algorithm. I've fixed the old algorithm so that it works there, too.

Thanks for the report,

Christophe


> On 14 Nov 2016, at 13:17, Gregor Mitscha-Baude <gregor.mitscha-baude at gmx.at> wrote:
> 
> Dear all,
> 
> I was happy to see that the new gmsh version (2.14.1) includes the long-missing "Point In Volume" feature.
> However, if multiple volumes are created in the .geo file, with points added to each volume by "Point In Volume", then only the points lying in the *first* volume are meshed in the output .msh file.
> 
> I attached a simple example file input.geo with two cubes stacked on top of each other, divided by the z=0 plane. vol1 is the lower cube (z<0) and vol2 the upper one (z>0). A series of points of the form (0, 0, z) with z negative and positive is added with Point in Volume. The output file was obtained by running
> 
> >> gmsh -3 input.geo -o out.msh
> 
> and is also attached. As can be seen, all the embedded points with positive z values are missing from the .msh file. Hence, only the points embedded in vol1 were meshed. I observed this also on more complicated examples with many different volumes: only the points embedded in volume 1 are meshed.
> Is this a bug? Is there any workaround to mesh embedded points in multiple volumes?
> 
> Cheers, and thanks a lot,
> Gregor Mitscha-Baude
> <out.msh><input.geo>_______________________________________________
> gmsh mailing list
> gmsh at onelab.info
> http://onelab.info/mailman/listinfo/gmsh

-- 
Prof. Christophe Geuzaine
University of Liege, Electrical Engineering and Computer Science 
http://www.montefiore.ulg.ac.be/~geuzaine

Free software: http://gmsh.info | http://getdp.info | http://onelab.info