Subject: CGAL users discussion list
List archive
- From: Sebastien Loriot <>
- To:
- Subject: Re: [cgal-discuss] Create a 3D mesh as coarse as possible
- Date: Tue, 5 Jul 2022 11:58:58 +0200
- Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
- Ironport-data: A9a23:o5Uku6Alp91uQBVW/xvlw5YqxClBgxIJ4kV8jS/XYbTApDIg12YPn 2tLW2mCPKyKYWb0L91zatjnp0oDuZLUmodnOVdlrnsFo1Bi+ZOUX4zBRqvTF3rPdZObFBoPA +E2MISowBUcFyeEzvuVGuG96yE6j8lkf5KkYAL+EnkZqTRMFWFw03qPp8Zj2tQy2YfhXlvU0 T/Pi5S31GGNi2Yc3l08sPrrRCNH5JwebxtF1rCWTakjUG72zxH5PrpHTU2CByeQrr1vIwKPb 72rIIdVUY/u10xF5tuNyt4Xe6CRK1LYFVDmZnF+A8BOjvXez8A/+v5TCRYSVatYoxSJwvN4w dRyia2TF1w3YJ+VxcsQTzANRkmSPYUekFPGCX22sMjW0EifNnWxmLNhC0Y5OYBe8eFyaY1M3 aZAeXZdM1ba3bLwnenTpupE3qzPKOHwOIQFu3Z8izTdJfkjSJHHBa7N4Le02R9p1p0RRqyON 6L1bxIyXRjtXw8XZG0NN7Yfhdjxqyfkdxtx/Qf9Sa0fujCPlmSdyoPFO9XcfpmGRN5eg12Dj nnX+nzwRBAcLt2WjzSfmk9AncfKlCL/HZ0XTfi2q6Usj1qUyWgeThYRUDNXvMVVlGafW498I kUVxRYfjvZp6hKZF4TSUgGR9SvsUgEnZ/JcFOgz6Qeow6XS4hqECmVsctKnQIx23CPRbWx6v mJlj+8FFhQ07+LIESP1GqO86GLtaXJMfAfucAddFVNdi+QPtr3fmf4mczqOOKu8j9mwBjSph j7X924xgLIcicNN3KK+lbwmv95OjsiZJuLWzl+PNo5A0u+fTND+D2BPwQaAhcus1K7DEjG8U IEswqByFtwmA5CXjzCqS+4QBryv7PvtGGSC3AA+Qsh9rG/2pC/LkWVsDNdWdBcB3iEsKW+BX aMvkV45CGJ7Zyb1MfcsP+pd9exzlPG8T7wJqcw4nvIXOsQrHON21C5pYkGU0gjQfLsEwMkC1 WOgWZ/0Vx4yUPw5pBLvHrt1+eJ1m0gWmDyLLbimnkzP+efPPxa9FOZVWHPQNLxRxP3f+239r Y0PX/ZmPj0FD4USlAGModBNRb3LRFBnba3LRzt/L7LcfFM7QDt+U5c8A9oJIuRYokicrc+Ql lnVZ6OS4ACXaaTvJVrYZ3Z9RqnoWJoj/3s3MTZ9b1mt0nknJ42o6f5HJZcweLAm8s1lzOJ1H 6FVIZXeXqwXR2SV4SkZYLn8sJdmK0amiAeICCyvP2oycptmcArW94K2ZQDo7iQPUnG6uJJm8 b2t3w/WW7QZQAFmAJqEYf6j1QLjsn0UmeY0VEzNe4EBdELp+YlsCirwkv5ne5FWeUufnmOXj l/EDw0ZqO/Bp54O3OPI3a3U/Z20F+ZeH1ZBGzaJ4LuzMx7c9DXxzIJFVtGOYmmBBm75/aOVZ dJVwev5B/sJkQsYqIF7Cbtqkfsz6taz9b9XygNoQCfCY1ixUOgyJ3CH2YxQs/QIyOIC/wSxX U2L959RPrDQYJHpF1sYJQwEaOWf1KFLxmODs6xteEiqtjVq+LenUFlJO0XegiJqKrYoYpgux v0suZJL5gHj2AAmNM2K0nJd+2iWdCdSVqwmstQFBdauhFdwjF5FZpPYB2n955TWM4dANUwjI zm1gqvehuQDmhCTLSJrTXWdj/BAgZkuuQxRyANQLVq+nNeY1OQ82wdc8GhqQwlYpvmdPzmf5 oS224xJyaSyE/NAgcFCWyWzHlgECkDCvEP2zFQNmSvSSEzAuqkh6oEiEb7lwazb2zs0kvtnE HWwx2PsUDKsd8b0tsf3cVAwsOTtFLSd6SWb8P1K3K25810Sbj/sg6vobm0Nw/ci7QXdm2Wfz dRXEC1MhWEX+MLeT2DXy2VX6FjIdC25GQ==
- Ironport-hdrordr: A9a23:rbUlhqg9kuJb1WJOKAXYFUX6vHBQXs0ji2hC6mlwRA09TyX4ra CTdZsgvyMc5Ax9ZJhCo7G90cu7Lk80nKQdieIs1NyZMjUOzVHYSb2KmrGSpwEIdReTygck79 YDT0EoMqySMWRH
- Ironport-phdr: A9a23:ANkdkhNQQglSqw4fA6sl6nZXBBdPi9zP1u491JMrhvp0f7i5+Ny6Z QqDv64r1QGQFt+Bo9t/yMPu+5j6XmIB5ZvT+FsjS7drEyE/tMMNggY7C9SEA0CoZNTjbig9A dgQHAQ9pyLzPkdaAtvxaEPPqXOu8zESBg//NQ1oLejpB4Lelcu62/689pHJfglFhSexbbxsI BmosA7cqtQYjYx+J6gr1xDHuGFIe+NYxWNpIVKcgRPx7dqu8ZBg7ipdpesv+9ZPXqvmcas4S 6dYDCk9PGAu+MLrrxjDQhCR6XYaT24bjwBHAwnB7BH9Q5fxri73vfdz1SWGIcH7S60/VDK/5 KlpVRDokj8KOT03/m7YhMN+kbxVoByhqRN934Hab5qYNOZ9c67HYd8WWWRMU8RXWidcAo28d YwPD+8ZMOhFsYbyvEEBrQCiBQmtGuzvzCJDiHnx3a0g0uQhFh3G1xEnEtIVtHTUq8/1NKYVU eCp1qXIyDHDYuhZ2Tf48ofIcxQhreuQUrJ3dMrc0E8iHB7KgVuMs4LqJS+V1vgTvGiB6eptT f6jhnAmpg1trTaixMMhh4vJi48IxV7I6Tl0zYUoKNC6SEB2f96pHYdTuiyaOYZ6X8MvTW9ot Som1LEIuZ62cicMxZ86yRDfbPmHfJKJ4hLlTOuRLjZ4hG5leLKinBm+61Svyur5VsWszFlKs zFKncfWunAK2RzT5NCLRedy/ki73TaDzQfT6vtLIUwslKrbLYQhwr8qmZYJv0TDBCn7k1j1g q+Obkgo5PSk5uD9brjlppKQLZJ4hh/9P6g0msGyD/w0PwsPUmSB++mx1afv8VP3TblUlPE6j 6nUvZ/cKM8GvKC2GRVV3Zwm6xunDzepztAYnX4fIVJAYh2HjozpN0jAIfD9EPuzmlqsnTdly vzcMb3hBZLNLnfHkLj/Z7py90lcyA8rwdBe4ZJbFK0BLertVkPtsNHUFBw0Pgyuz+r6Ftlw1 ZkSVGKBD6OBNaPdq16I5uYhI+mWY48VvS7wK+Mh5/7pl385mFgdfbWz0ZsSbHC1BftmI0CDb nrthtcNC3sFvg07TODyjl2NSiZcZ2yuUKIk+jE7FIWmAJ/eSoy1mryOwD+7HoFKZmBBEl2DD XjoeJ+AW/sVdS2SI9RhnSAZVbi6UI8szgqutQ//y7p/NOXY4CwYtZT51Nh0/eLfjx8y9SYnR /mbyHyHGmFognsTFXhxx7F6uUU7y1GZ0KE+jeYfDs1W//oOUwE0MtnXwOV+TtzzQQncZczad VC9X9+aDCEtG9It38cVMQE6AMSnlhmF3iywArZTmabMH40x6qub3n7/IIF2xH/CkaUglFI7W dAcCWrzjaF2803fBpXCjl6Cv6esb6UVmiDXp0mZym/bh0xSWRVsUKjDFVQYfErRsZyt/UfFV b6pFfIiNiNOzMeDLu1Bbdi/3gYOf+vqJNmLOzH5oGy3HxvdnttkDaLvcmQZh2DGDVQc1hsU5 TCAPBQ/ASGopyTfCiZvHBTheRCk6vFw/VW8SEJ81ASWdwt5zbPg4RoSn/2bV7UW2poLvS4gr 3N/G1Pul8nOBY+4rhF6NL5Zfct75V5G0WzDsAkoJZinNaFlmhgbdyx4ukrv01N8DYAT2dMyo iYMywx/YbmdzEsHdz6c2sXoPabLL2Dp4B21Q6vf21Wbz9PPv6lSs7I3rFLsuAzvHU0nm5l++ /9S1XbUppDDDQ5JFIn0Tl5y7B9i4bfTfig64YrQk3xqK6i99DHYiZovA6M+xxCscs06UuvMH RLuE8AcG8mlKfA701mvYBUeOelO9akyd8q4fvqC0aSvMa5uhjWjxWhA5Yl81AqL+U8eAqbT2 5EfwvaEmA6Dfzj5hVal9MvwnMEMZD0fGHa+1TmxHJRYNeV5eYcGD3vrItXinI0vwc6wHScCp Bj+XABjuofhYxeZYl3j0BcF0E0WpSfigi6k13lvlClvqKOD3SvIyuCkdRwdO2cNSnMx6DWka YWykd0eW1ClKgYzkx7wr13+wLJao7g5Km37TkJBfiywJGZnGPjV1PLKc4tU5ZUkvD8CGv+4Z kqbTaK7ph8y3CbqHm8Yzzc+PWLP2N2xj1lxj2STK2x2pXzSdJRrxBvR09fbQOZYwjsMQCQQZ SD/PlGnJJHp+NyVk8yGqeWiTye6UYUVdyD3zISGvS/95GtwABT5kerh0tHgFAE712f82bwIH W3TqBHmY473kaG+Ge1id0hsQlT77oJ2F5p/nY05mJwLkSJC19PFoDxeyze1bI4T0Lm2dHcXQ D8X39PZhWqtkFZuKH6E3cOxV3mQxNdge8jvZ2oX3iwn6MUZQKyQ7bFCgW50ug/i9VOXMaU7x 21Nj6Z1uxt4y6kTtQEgzzuQGOUXFEhcZmn3kgiQqsu5tONRbXqudr650Ax/m8qgBfeMuFI5O j6xd5E8ECt39sg6PkjL1Si58YXjYtjXcZQWsjWblh7Bi65eL5d7xZ9ozWJ3fHnwu3Eo0btxl hhpx5C9oM6CL01i+au4BlhTMTi/NKZxsnn9yK1ZmMiRxYWmGJ5sTy4KUJXfRvWtCDsOtP7jO lXGAHgmp3ycA7aaAR6H5RIssSfUC57yfSLyRjFR3ZB4SRKaPkAanA0EQGBwgMsiDg7zjM35L BUiu3ZItwa+8EcTjLovbUW3U3+D9ln0LG1vE97GckIQtkYbtiK3eYSf9r4hQX8ep8X76lTLc ivBP0xJFT1bBBLCXQyyeOn2o4GHqbDQB/LifaSUJ+zS77UPDbHQgsv/t+kutzeUapfQYj87V aB9ggwbGiknU8XBx2dWE3xRzn2SKZ7d/FDmo2V2tpztqauwHlu+uc3XTeMVaIsKmVj+gL/fZ bTI1WApdHABjMNKnTiRl/AexAJA0XgwMWT9V+1R72iVC/uB0q5PU0xBMn01bpAZqflmmFEKY J++6Ju9wLd8irRd505tc1vnl4npYMULJzr4L1bbHAOQM6zAIzTXwsbxaKf6SLtKjewSuQfi8 TCcW1TuODiOjVyLH1imLP1MgSeHPRdfpJD1cxBjDnLmRc7nbRvzOcF+jDk/y7k5zn3QMmtUP T95ekJL5rqeiEEQyu14AHBE52F5IPOsni+Y66zHL89Tv6YwRCtzkO1e7TIxzL4UpCBISfppm TfD+95jp1b194vHgjFjURdItnNKnNfR5RQkafifr8MQHyufr3dvpS2KBh8HpsVoEIjqsqFUk Z3Ukb7rbSxF65TS9NcdAM7dLISGNmAgOFznAm2xbkNNQDi1OGXYn0EYnuuV8yjftZw9sJnrh NwLTpdUUVU0ErURDUEvT7lgaN9nGygpl7KWlptC/X2lsBzYX9lXpLjCX/OWROrqcXOX0OICa BwPzrf1a48UM8eovi4qIkk/l4PMFU3KWNlLqSA0dQ44rnJG931mR3Ey0UboAutIyHAWHP+w2 BUxj1kmCQzM3DLl4lNyOV+T4SVtyw8+ntLqhT3Xezn0fv/YtW5+BC/9tkx3OZT+EV8dUA==
- Ironport-sdr: DFtxEB1Z/ADHEnlolDwtVNY3wMPjA88m6PhfOhHZI24qBn6R8A2MzwPyGf4c4QsbRCVNtffLQY I0G9ZtoU17GaVHq7/cETS8JLszX+5D1PeJDz0+0huCbPbzFN/S04MXD1Vik7c+FgdX9ZDtvzhU ItBvRjGFz9cvMQzVNAyaa44IZfbSvXWiD+RzLsj64HScGfIZXkdirXK3AChvvbfIM5DLKSF/IX QRdOjFLH0tlWZ9lUH6jMae8IO/Ux4GTip4nDmJeHxJOC75wG+n0UNCXSSOYiI0ZyAQoyUA6P8k HRVl4SuqGaffgc+f9jC4C1l9
Sorry for the late reply.
Unfortunately, you won't be able to get much less that 52 cells.
In this particular case, you have some vertices on edges added so
that they are in the output and then initialization (some points
thrown in the domain). For this particular case, we could probably
not add any extra initialization points thanks to protected edges
but in general this step is essential to get a valid initial
triangulation to be refined.
Maybe what you need is more something like a 3D constrained Delaunay
triangulation that CGAL does not provide yet.
Best,
Sebastien.
On 5/10/22 15:55, Bob Bill ( via cgal-discuss Mailing List) wrote:
Here's the link to the minimal working example on Gist: https://gist.github.com/boblukas/603c16f4aceca0d7a5b46fa8150fae29 <https://gist.github.com/boblukas/603c16f4aceca0d7a5b46fa8150fae29> (The cube.off is taken from CGAL repo: https://github.com/CGAL/cgal/blob/master/Data/data/meshes/cube.off <https://github.com/CGAL/cgal/blob/master/Data/data/meshes/cube.off> )
The output shows 52 cells, which are way too many for my application.
Thanks a lot for your time,
Bob
Il martedì 10 maggio 2022, 15:27:41 CEST, Sebastien Loriot" ( via cgal-discuss Mailing List) <> ha scritto:
The number of input triangles does no affect the number of output cells.
Note that a simple working example that we can compile and run on
gist.github.com would most probably understand your problem.
Right now I can only guess that the number of output cells is the one
needed to preserve the topological criteria.
Best,
Sebastien.
On 5/10/22 15:23, Bob Bill ( <> via cgal-discuss Mailing
List) wrote:
> Unfortunately the result is more or less the same, but I think I
> identified the issue: my original surface mesh is not triangulated, so
> maybe CGAL::triangulate_faces() is adding too many triangles, right?
> The following snippet highlights what I am referring to. Is there any
> way to have a "minimal" triangulation for the faces when calling
> triangulate_faces()?
>
> Best,
> Bob
>
> CGAL::Polygon_mesh_processing::triangulate_faces(surface_mesh);
> Mesh_domain domain(surface_mesh);
> domain.detect_features();
> Mesh_criteria criteria(CGAL::parameters::facet_size = 0,
> CGAL::parameters::facet_distance = 0,
> CGAL::parameters::cell_radius_edge_ratio = 2,
> CGAL::parameters::cell_size = 0);
> C3t3 triangulation = CGAL::make_mesh_3<C3t3>(domain,
> criteria,
> CGAL::parameters::no_perturb(),
> CGAL::parameters::no_exude());
>
> Il martedì 10 maggio 2022, 09:29:46 CEST, Sebastien Loriot"
> ( <> via cgal-discuss Mailing List)
> < <>> ha scritto:
>
>
> If you put facet_size=0, facet_distance=0, cell_radius_edge_ratio=2 and
> cell_size=0 then there is no geometric criteria applied to refine the
> mesh (only the topological one). If you mesh is too coarse then you
> can change those values to match the level of detail you are expecting.
>
> Best,
>
> Sebastien.
>
> On 5/9/22 14:51, Bob Bill ( <>
> < <>> via cgal-discuss Mailing
> List) wrote:
> > Dear all,
> >
> > I successfully applied the function make_mesh_3() to fill a 3D mesh. So
> > far the result is too refined, and I'd like to let my C3t3 to be as
> > coarse as possible, i.e. I want only few tets inside.
> >
> > If I use the default constructor
> >
> > Mesh_criteria criteria;
> >
> > then the result is non-sense. Is there any way to have something like
> > the bottom left figure here
> > (https://doc.cgal.org/latest/Mesh_3/index.html#fig__figureparameters <https://doc.cgal.org/latest/Mesh_3/index.html#fig__figureparameters >
> <https://doc.cgal.org/latest/Mesh_3/index.html#fig__figureparameters <https://doc.cgal.org/latest/Mesh_3/index.html#fig__figureparameters>>) ,
> > but without specifying the parameters for every situation?
> >
> > Bob
> >
> > --
> > You are currently subscribed to cgal-discuss.
> > To unsubscribe or access the archives, go to
> > https://sympa.inria.fr/sympa/info/cgal-discuss <https://sympa.inria.fr/sympa/info/cgal-discuss >
> <https://sympa.inria.fr/sympa/info/cgal-discuss <https://sympa.inria.fr/sympa/info/cgal-discuss>>
>
> >
>
> --
> You are currently subscribed to cgal-discuss.
> To unsubscribe or access the archives, go to
> https://sympa.inria.fr/sympa/info/cgal-discuss <https://sympa.inria.fr/sympa/info/cgal-discuss >
> <https://sympa.inria.fr/sympa/info/cgal-discuss <https://sympa.inria.fr/sympa/info/cgal-discuss>>
>
>
>
> --
> You are currently subscribed to cgal-discuss.
> To unsubscribe or access the archives, go to
> https://sympa.inria.fr/sympa/info/cgal-discuss <https://sympa.inria.fr/sympa/info/cgal-discuss>
>
--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss <https://sympa.inria.fr/sympa/info/cgal-discuss>
--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss
- Re: [cgal-discuss] Create a 3D mesh as coarse as possible, Sebastien Loriot, 07/05/2022
- <Possible follow-up(s)>
- Re: [cgal-discuss] Create a 3D mesh as coarse as possible, Sebastien Loriot, 07/05/2022
Archive powered by MHonArc 2.6.19+.