Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Create a coarse surface mesh produces areas of very small elements

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Create a coarse surface mesh produces areas of very small elements


Chronological Thread 
  • From: Andreas Fabri <>
  • To: "Sebastian Bachmann ( via cgal-discuss Mailing List)" <>
  • Subject: Re: [cgal-discuss] Create a coarse surface mesh produces areas of very small elements
  • Date: Thu, 25 May 2023 12:51:46 +0200
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-data: A9a23:0SDlBqwy1SUDR04PUnZ6t+fXwirEfRIJ4+MujC+fZmUNrF6WrkUCz DNNXm6CPfqCazDwfYsjO4rgp0JSu5fTz4BlQAY4qlhgHilAwSbnLYTAfx2oZ0t+DeWaERk5t 51GAjXkBJppJpMJjk71atANlVEliefSAOCU5NfsYkhZXRVjRDoqlSVtkus4hp8AqdWiCmthg /uryyHkEAHjg2Yc3l48sfrZ80sy5qiq5Fv0g3RnDRx1lA+G/5UqJMlHTU2BByOQapVZGOe8W 9HCwNmRlo8O10pF5nuNy94XQ2VSKlLgFVDmZkl+B8BOtiN/Shkaic7XAhazhXB/0F1ll/gpo DlEWAfZpQ0BZsUgk8xFO/VU/r0X0QSrN9YrLFDm2fF/wXEqfFPM56RtNlgsFrQC8+tMLU5j7 /UKMxYSO0Xra+KemNpXS8E2wMEkcZm3eoYWu3UlyizFS/E4QdbFTrmiCd1wgGl235gTW6mEI ZFHAdZsREyojxlnM1EbEp83h6Gmj3PlchVXpVWQqLYt8mbaxxB2yqmrO93QEjCPbZgLzxnH/ jqdl4j/KjMaGuTA1SWryUq12cWRxz3HZIItP5Tto5aGh3XJnjROWXX6T2CTqvawjguyWslUN lcP0jE/qLA7sk2tVNj0GROiyENopTZGHd9XTr1mrgSEy66R5BuFQG8aTnhHZcBOWNIKqSIC9 nvYnOnwHxdWqrSYcl6y262dvyiUJn1ARYMdXhMsQQwA6tjlhYg8iBPTU9pueJJZaPWoRVkcJ BjX90ADa6UvYd0ji/XroACd6964jsKSFV5qjunCdj/9tmtEiJiZi5uAx2Kz0BqtBNzfS17c7 CRCnsGf6KUJEI3LkzKNBuMAANlFBspp0hWC3TaD/LF7r1xBHkJPm6gMvVmSw281bq45lcfBO hO7hO+ozMY70IGWRaF2eZmtLM8h0LLtE9/oPtiNMIoSM8UqLFTWp34zDaJ144wLuBh3+U3YE cfAGftA8V5LVPg5pNZLb7xMieF7mXtWKZ37HMqmlE//uVZhWJJlYe1daAbfBgzIxKSaoE3O+ s03Cid540U3bQEKWQGOqdR7BQlSfRATXMmmw+QJLLLrClQ8RwkJVaSLqZt/INMNokigvr2Tl p1LchMAlwWXaLyuAVniV02Pn5u2A8gv/SxibHR0VbtqslB6CbuSAG4kX8NfVdEaGCZLlJaYl tFUI5vSMecFUTnd5TUWYL/0qYEoJlzhhhuDM2DhKHIzdoJpDV6BsNL1XBrdxA9XBAqOtOw6v +KB0CHfSsE9XAhMNpvdR8+u6FKThkIjvtxOcXHGGfRpQ3W0wrNWc3TwqtQVP/AzLQ7yw2rG9 gSOXjYdi+r/g64019jrg5K7/pubI9FvE3VoH1iA16i9MyrHzzDy3a5BT+e6UjTPX0zk+KiZR Lt0ztOtFNYljVp1o45HPLIz9p0H5vzrvK58zC5oOF7pfmabIOpsDVfe1PYeq5AX4KFSvDWHf 36m+/5YCO2vA9zkGltAHzgVRL2P+t9MkwaD8MlvBlvx4RJ227+1UU9yGR2oow4FJZtXNLIV+ 8sQiPQ02See1CVzas2nixpK/VuiNnYDCqUrlq8LCb/R1zYE9AtwXoz+OATXvre/M85BI2s7E A+y3aDivYlR9mDGUng0FEXO49ZjuIQziEhK4mMvd1WtsfjZt8AzxyxUoGgWTBwK7xBp0NBTG 2lMNmxkAaDX4RFUi5BnWjmsEj5eGBfCxBTVzgpTnUaEUnuYb23pB08+MNar40o23T99fD9a3 bfA00fjc2/gU//Q1xsIe3xOiqLcX/lu0DbdieabH8ihNLsrUwrP26OBSzIBlEr6PJkXmkbCm 9hPwM9xTq/KbQgrvKwxDtih54Q6ERyrCjRLfqB8wfkvA2rZRTCV3AqOIWCXfudmBaTD0W28O vxUCvN/bTaM/wfQkWlDHo8JGaF+o9Ax7tlber/LG38Pg4HCkhVX6qDv5grMr051ZeU2isssC JLjRxTbGEyquHZksWvsrs5FB2mGXec5dDDMhN6SzuFYOK8A4cdNcF4z2ISap3+6EhVq1DPKs RLhZ5359f1Dy4NtrdHSDK98OV2/L9f2ZduuqArokYlqaPXRO/zvrCIQkEHsZC5NDIsSWvN2t LWDi8H210X7p4QLU3jVtp2CNqtR7+CwbfVnCd32J3xkgieyYs/gzB8d8WSeK5YSstdiyuS4Z gm/Mu2cSMU0XopD+XhrdCRuKRYRJKDpZKPGpynmjfCtCAAY4DPXPuGc6n7lQmFKRBAmY6SkJ FfPhM+vwdREoKBnJhwOXahmCqAlBm7TY/ItctmpuASICmWtvEi5hYLjshgeuAH7UizOVI6w5 J/eXRHxeSijoKyCnpkTr4V2uQZRF3pnx/U5ekUG4dNtljSmFyg8IP8ANYkdQIRh+sAoOEoUu BmWBIfjNcn8YdiAWU+65dO+DlvZA+UPPpL+OyBv+F2ULSG7GOtsxVenGjhIux9LlvnLlYlL6 u3yPlX/Mx+0z41zVOga7eC8m/YhzfTfrp7N0V6oiNT8Wn7yHp1TvEGM32NxuejvHMzKkUjXP 3k7TGtYR1uqD0X2FK6MvpKT9A4x5FvS8tnjUctDLBszdWlWIC2sBcAT493O74A=
  • Ironport-hdrordr: A9a23:TWyirKAkYaEhMGjlHelK55DYdb4zR+YMi2TDt3ocdfUnSL39qy nOppomPHDP6Ar5NEtQ4OxofZPwOU80lqQFgrX5X43CYOCOggLBR+sO0WKh+UyHJ8SUzJ856U 4PSdkYNDS9NykZse/KpC+GEpILxdmK972JnI7loEtFfEVNdKAl0hxwAAqBVm16LTM2fqYRJd 6i3Y58nhKFU1wqVOjTPAh4Y8Hz4+Lx0KjeSTo9LTgdxGC1/E6VwY+/LCPd5C43dwlz6YoPmF K16TDE2g==
  • Ironport-phdr: A9a23:TDRljhbHcG/v1E58NwJup1T/LTEr2YqcDmcuAnoPtbtCf+yZ8oj4O wSHvLMx1gKPBtqEoKkcw8Pt8IneGkU4qa6bt34DdJEeHzQksu4x2zIaPcieFEfgJ+TrZSFpV O5LVVti4m3peRMNQJW2aFLduGC94iAPERvjKwV1Ov71GonPhMiryuy+4ZLebxtIiTanb75/L he7oQrfu8UKnIBvNrs/xhzVr3RHfOhb2XlmLk+JkRbm4cew8p9j8yBOtP8k6sVNT6b0cbkmQ LJBFDgpPHw768PttRnYUAuA/WAcXXkMkhpJGAfK8hf3VYrsvyTgt+p93C6aPdDqTb0xRD+v4 btnRAPuhSwaLDMy7n3ZhdJsg6JauBKhpgJww4jIYIGOKfFyerrRcc4GSWZdW8pcUTFKDIGhY IsVF+cOIeVWoYfyqFQAsBWwGwesCuPtxTJTnHD6wbc33v49HQzI3gEtGc8FvnTOrNXyMacfS fy4zK/VwjXZcfxZwyz945XIfB87uv6MR7NwetfWxEY1EQPKk0+QppL/MDOOzOgNvXKU4PF6V eKukGMnpB9+oj+qxsgylonJnIIUyl7e+ipl3IY1Ise4SU5hbdG4F5tQsjiXOo1rScwtX29mo jw1yqEauZGlZigKzowqywLCZvKJcYWF4BzuWumMLDp2mH9pZa+yiwux/EWjy+DxUs253UpXo idYjtTAqmwB2h7O58WZSPZw+kiv1DKR2w3V9+pKL0c0la/BJJ4gxL49joYcsUvZES/ygkn2i 7WWdko89uin6uTnfrXmpoKHOINulg7+NaEul8qiCuo7KggDR2eW9OWm2LH+4UH0QK9Gg/8on qXDvp3XJNwXq6qlDwNPz4ov9hWyAy273NkYk3QLNkxJdRyDgoTxNVzDIur0APe7jli2jTtmy O7KMqflD5jOKHXIja3vcqxn60FGzQo+1dBf6IxQCrEGOP//RFX9tNnfDhI5MAG5zPzrCNBn2 YMfXWKDG6iZML7OsVOS4eIvPemMaJUTuDnjN/gp///ujXknll8BZaSk0poaZGqlEvh6I0iVe 3jhjskbHWoEvwcyVOnqh0eDUT5XaXayRaU85jQjBY28C4fMW5qtj6Kb0yegBJFWaHpJBUuWE XfpbIiEVOkDaCWXIsB/jjwET7igS4E72h6yqA/6zKBnL+rO9iIEspLj0cB55+vImh4q7zx7E 9yd032RT2Fzhm4EWic63Lpjrkxl1leDza94juREGtxc/fxJVh42OoPdz+xhF9/yRxnBfsySR Va9QtSmBCkxQcgrz98PZUZ9AdSigQrZ0yqkGb9G34CMUdYP/6jVxGTwLsA14X/N1KwoiV9uA stCMWmhrrRy8BiVDpKBjUSQ362lfKBZ0C/W/3qY1iKrtUBcUQh/GYHERmxeLAOCptvw4gbOT qSlFK88GgpH08+LbKVQPI7Hl1JDEd7lNM7Tann5tW66HxfAkrqKYJDncn5b0izXEkksnA0U+ HuaLxkwDyy9pHjPSjdpEAS8MAvX7eBipSbjHQcPxAaQYhgkjuLtkvZ0rfmVSvdIm6kBpD9ks DJsWlC0w9PRDdOE4QtnZqRVJ90nsx9czWyMkQt7M9S7Krx6wEYEelF+uU/02hxsTItJm9Irh Hwnyw9/NbiJ3lpKazSCzNb7PbiEYnLq8kWXYrXNkkrbzM7Q/64O7PoirFC2tQelDEcr6DNp2 tNP0lOT65LPARYIQJz4WVow7Qk8rLbfMWEm/42B83pqPOGvtyPandIkAOxw0hG7Y9JWK7+JD ieoVcgdWpjybukjmlzsaQ8YeudM9OgyMt/Om+Ku/qmtMa4gmTuniT8C+4VhygeX8DI6TOfU3 pEDyvXe3w2dVj66gk3z+sbw0ZtJYz0fBA/dgWDtGZJRa6tufI0KFXbmIsu5wc97joLsXHgQ/ UCqBlcP0sukMRSIaFm10QpV3EUR6Xup/EnwhzVynyskp7HZ0iXE2eXKexcAP2NXXnhshFzwJ pKlydsdWQngbgQkkgek+Vev37JS98EdZyHYRUZFeTSzLnk3DvX27+LfJZEVrsh24kA1GKymb FuXS6DwuU4f2iLnRC5FwSwjMiqtotP/lgB7j2SUKDByqmDYcId+30S6hpSUSPhP0z4BXCQ9h yPQAw32Nt+l59idjNHNu+qkVkquW5pWfDXx3I2Jvze8/3wsChq61aPW+JWvAU0h3Cn32sM/H yzOoAzxa5Kt2a2wK+NPcURvAVLg8dt0E4pilZEhwpoX3DJJ4/fdtWpCmmD1P9JB3Kv4Z3dYX j8Hzene5w390VFiJHaEr27gfk2U2dApJ9yzY2dNnzk489gPE6CMqrpNgSpypFO86wPXe/l02 DkHm7Mi73sTgudBvwRIrG3VC7YeB0RRJmrinh6S7vixoaJSaXq1YLa531Z5h8HnB7aH6g1RQ 3f2fJ4+EDQ4sps5aQqRliKprNq6MNDLCLBb/gWZiRLBk/RYJNoqm/wGiDAmcWPxsHs5yvIq2 Bln3JW0po+CeC1m+KO0BAIdNyWgO51Vo2mwy/kE2J3IjOXNVt17Fz4GXYXlV6etGTMW7rH8M hqWVSY7sjGdEKbeGgmW7AFnqWjOGtalLSLyRjFRwNN8SR2aPEEajhoTWWBwl5g1DA2j34rve U1j5xge61n9pwdW2+xhPAX4SHaZrwCtIGRRKtDXPF9N4wdO6l2Ae8WY4vh+Fjoe8JSrtg2lJ WGcYglUF3AHU0eYAEryeLKp4JOTlorQTvr7JPzIb7KUrOVYXPrd3pOj3Lxt+DOUP9mONH1vX LUrn1BOVndjF4HFii0CHmYJwjnVYZfR93LesmVn69qy+/PxVEfz6JuTXvFMZM539Un+mbuZZ euZhSI8LitEkJQSxTnOxaRXxEII2SRndj3rFqoA8yjAUMey0udWF0JJO2Uob40RtP16h1EQf pWBwtLtiuwh36RzWgYDCgOnypjuPpZCdGi5MBmv6F+jErOAKHWLxsj2Zfj5UrhMlKBPsAX2v z+HEkjlNzDFlj/zVhnpP/sexCecdAdTvo2waHMPQSDqUc7mZxunMdR2kSx+wLs6gWnPPHIdN j40el1Er7mZ5ydVyvtlHGkJ4n1gJOiC0yGXiouQYo4RquduCz9omvhyujI/zOYFsGdBTf1x3 SzPsphpvVHglOSTi3JmXBdItjdXldeLsEFlasC7vtFLXXfJ+g5I7H3FV0pM/oM4TIK1/fsOm b2t3OrpJTxP8szZ554ZDsnQc4ecNWY5dAHuAHjSBRcESjiiMSfegVZcmbec7C7wzNByp57yl Z4JUrIeWkYyE6ZQA0JpBtELO9F5Vzk+kJaUgckN6GakvRfYT9ldpIGBXfWXS6aKSn7Rnfxfa h0EzKmtZ5wULZH+0ld+Z0NSxsLPHBeAAZZIqyxlKwgpvANK7ng4SGAvkRGACEvl8DoYEvi6m QQzgw11bLE29Tvi1FwwI0LDuCo6lERZcTTN23aUd2OhcOG1VIBSTi3pqw42LJO9RQtpP1Xad alMOzrDQr9NlapufGt3jxXN/5BIHKwEJUWrSBAdwvSae+85311XtiK91AlM4u6XUPNf
  • Ironport-sdr: 646f3dc5_jdntdG01sk/7t2g5UlxMn7qFHTIqNJLSQojuPrawJwBcl5W gT4JOaFj3733WmQAFHLfj8oE4nxLr6/JaXpmZcw==


Hi Sebastian,

Even 5.5.2 needs a call to custom initialization.  

https://doc.cgal.org/latest/Mesh_3/Mesh_3_2mesh_3D_image_with_custom_initialization_8cpp-example.html

So you must add

#include <CGAL/Mesh_3/initialize_triangulation_from_labeled_image.h>

and then call

 C3t3 c3t3;
  initialize_triangulation_from_labeled_image(c3t3,
                                              domain,
                                              image,
                                              criteria,
                                              static_cast<unsigned char>(0));
  CGAL::refine_mesh_3(c3t3, domain, criteria);

Best,
Andreas

On 5/25/2023 9:15 AM, Sebastian Bachmann ( via cgal-discuss Mailing List) wrote:
I tested this now in 5.5.2 (according to github it should be in there as well?) but I cannot see any difference to 5.4 - except when using the custom init (see attached image). Code is still this one: https://gist.github.com/reox/0df1f446f0ff9d46e0688b60e8e858b8

Best,
Sebastian

On 25.05.2023 08:32, Sebastian Bachmann ( via cgal-discuss Mailing List) wrote:
Hi!
Thanks for the hint as well!
I'm using  5.4 / 1050401000.
I'll also test with 5.4.4 (looks like that is the version where it was merged?)

However, also the custom initialization works really good so far.
I fixed the seed for the random number generator to get deterministic results, which works very well.

Best,
Sebastian


On 25.05.2023 08:15, Jane Tournois wrote:
Hello,

this local over-refinement issue has been fixed in the master branch by the following pull-request :
https://github.com/CGAL/cgal/pull/6978

Which version of CGAL are you using?

Best,

Jane.



-- 
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project



Archive powered by MHonArc 2.6.19+.

Top of Page