I am having some trouble figuring out how the member `holes` (and the
method `set_holes()`) of `meshpy.triangle.MeshInfo` objects works.
Referring to the example file "examples_test_triangle.py" that computes a
mesh in a domain comprised between a rectangle and an ellipse, the steps to
compute such mesh are:
a. define points and facets for the inner boundary (rectangle);
b. define points and facets for the outer boundary (ellipse). Extend the
arrays containing points and facets at a. with the points and facets
computed at b.;
c. create a `MeshInfo` object, set its points, facets and holes and then
build the mesh.
So my questions are:
1. how does MeshPy distinguish the points (and the facets) computed at a.
and b. as belonging to different boundaries if they are contained in the
2. what is the argument that is given to the `set_holes` method?
3. is it possible to add multiple holes in the outer domain? If so, how?
Thanks in advance for any kind of input.
P.S. Thanks to Andreas Klöckner for this great tool: it is very powerful
and practical to be used from Python, too.
Show replies by date