Subject: CGAL users discussion list
List archive
- From: Sebastien Loriot <>
- To:
- Subject: Re: [cgal-discuss] 3D triangulation: absolute minimum of facets
- Date: Tue, 5 Jul 2022 17:21:46 +0200
- Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
- Ironport-data: A9a23:up9zy6lXJsM+RqxK3XPKI47o5gyzIERdPkR7XQ2eYbSJt1+Wr1Gzt xIdUD+Eb/2MNjDyeY13O4TgphsO6pXcytdgTgc4ri1jFFtH+JHPbTi7BhepbnnKdqUvb2o+s p5AMoGYRCwQZiWBzvt4GuG59RGQ7YnRGvykTrSs1hlZHWeIcg944f5Ys7N/09QAbeSRWVvX4 4ut+5CHYjdJ5hYtWo4qw/LbwP9QlK+q0N8olgRWiSdj4TcyP1FMZH4uDfnZw0nQGuG4LcbmL wr394xVy0uCl/sb5nxJpZ6gGqECaua60QFjERO6UYD66vRJjnRaPqrWqJPwZG8P4whlkeydx /0KkM2JagoEMZfXndQEfBdiTitVGLZZreqvzXiX6aR/zmXDenrohuprVQQ4YdFe9eFwDmVDs /cfLVjhbDjZ37PwkO/9E7M8wJh/RCXoFNt3VnVI1jXeF/chXdbGR43F4NZZ2HE7gcUm8fP2P ppHMWszMESojxtnKA4xOIM8tsqSh1LxUG1/q1evvLYH7D2GpOB2+OG1bIC9lsaxbc5al0Ldq mPd9HniGTkBJdmHwHyE9Gitj6nBh0vGtJk6EbS58rt1hQTWyDBPThIRUlS/rL+yjUvWt89jx 1I8vQ9/lIs/2UORVdzNRjK2miOOlTsHRI8FewEl0z2lxq3R6gefI2ELSD9dddAr3PPaoxR6h jdlePu5VVRSXK2ppWG1rejL8GvjUcQBBSpTOn9eFFptD8zL+dlr1nryosBf/LlZZ+AZ9Bn1y jGO6Tc73vAd0ZJN2KK88lTKxTmro/AlrzLZBC2HDgpJDSsjPOZJgrBED3CFtJ6sy67HFzG8U IAswZT20Qz3JcjleNaxaOsMBqq1wP2OLSfRh1Vid7F4qWnzqiPyIdsLu24vTKuMDiriUW+5C KM0kVMBjKK/wFP3BUOKS9nsVJhylvSI+SrND6mNMoAmjmdNmP+vpXkyPyZ8Lkjil08jlaxXB HtoWZfEMJruMow+lGDeb75Fj9cDn3lirUuOG82T50n4idK2OS/NIZ9YYQDmRr1ovMus/l+Jm /4BbJDi40sEAIXDjtz/q9F7waYidihlW/gbaqV/Koa+H+aRMDh+VKCLme1+IdwNcmY8vr6gw 0xRk3RwkDLX7UAr4y3QApy6QL+wD5t5s1whOikgYQSh13Q5MNSg6a4ec908erx+rL5vyvt9T v8kfcScA6QXGm6XpWhFNZSt/pZ/cBmLhB6VO3X3bTU6ealmTVOb99LheDzp6yRTXDG8stEzo uH72w6CGcgDSg1uAdz4cvWqy1/t73ERlPgjDUTNK9hXPk7r9dEyeSD2i/Y2JeAKKAnClmPKj VbIXU9AqLCU8YEv8dTPiaSVlKuTErNzThhAAm3WzbeqLi2FrGeuxIl3VuzXLz3QUWXD/rr7O bdYwvT6B/0wnFhQtr16Hbs2n7k14MHipuMDwwlpQCfLYlCsBu8yK3WKx5MU5KhEx7scpg7vH 0zTpp9VPrKGPM6jG1kUfVJ3YuOG3PASuz/T8fVlfxmgtXEvpOKKARdIIh2BqC1BN78pYokr9 uEs5ZwN4Aulhxt2b9uL0nJO+2KXIiBSWqkrrMtBUoriiw5u1VQbJJKAVmn555aAb9gKOU4ve 2fGiK3HjrVa50zDb3tjSiSXjLQF3cwD6EJQ0VsPB1WVgd6Z1PU56xtcrGYsRQNPwxQbju9+N wCH7aGuyXliItupuCRCY4xoMwRIBRnc5UmojlVVyCvWSE6nUmGLJ2o4UQpIEIb17EoEFgW3P pnBoIombdouVM701yo2H0VirpQPiPRvoxbalpnP89utRvEHjPmMvkNqTWUNohrjR8g2gSUrY AWsEPlYMcXGCMLbn0H350R2G1jdpNBo6VGumc1cwZ4=
- Ironport-hdrordr: A9a23:KKBHpq5vuLjl9wdgIwPXwOvXdLJyesId70hD6qkDc20yTiX+rb HToB17726ItN91YhsdcL+7Scu9qB/nhOVICMwqTM2ftWrdyQ6VxeNZnOjfKlTbckWUltK1l5 0QC5SWYOeQMbEQt7ec3ODXKadY/DDKytHNuQ4c9RtQpMNRBp2IIz0XNu9TKCNLeDU=
- Ironport-phdr: A9a23:nOyU5BDBIvqkuA4qYgBWUyQU9kkY04WdBeb1wqQuh78GSKm/5ZOqZ BWZua82ygaXBM6Lu7ptsKn/i+jYQ2sO4JKM4jgpUadncFs7s/gQhBEqG8WfCEf2f7bAZi0+G 9leBhc+pynoeUdaF9zjaFLMv3a88SAdGgnlNQpyO+/5BpPeg9642uys9ZDfeRhEiTugbb99L xi7rhjau9ULj4dlNqs/0AbCrGFSe+RRy2NoJFaTkAj568yt4pNt8Dletuw4+cJYXqr0Y6o3T bpDDDQ7KG81/9HktQPCTQSU+HQRVHgdnwdSDAjE6BH6WYrxsjf/u+Fg1iSWIdH6QLYpUjmk8 qxlSgLniD0fOjA6/m/YhcJ/gr9FrhKvpxJwwZLbbo6OOfpifa7QZ88WSXZPU8tTUSFKH4Oyb 5EID+oEJetWspfyqEAUohulGQmsBf3gyjlVjXHw2q06yeUhEQba3Ac9GN8Oq3DUrNLvNKcTV uC10LPEzTDHb/5N1jf97ZLHchElof2WQb1wds/RxFApGgjYgVqetZbrMCmJ1uQRrWeb9exgW PqthmM6pQx8oCWiy8YwhoTUiIwYy1/J+CZ7zYsoO9C1TEF2b96kHpZSsyyXK4R7T8EsTWxnp is217MLt5CmcSUOzpks2hDRa/uCc4eS4xLjUv6cITZ+hHJ/ebKwnQy+8Ua6yuD6S8K6005Ko zJbntXQsn0BzR/e58idRvdg4Eus2iyD2g/c5+1cIk05lLDXJpwgz7IqiJYev1rPEjXrlEj2g 6KbcFgv9PKy5OT9eLrmo4eRN49qhQH6NaQjgsm/Dv45MggKRmSa9+G82KD68U33T7hHgec6k qbesJDdKsQborC2DxVJ3YYk7hazFzam0NIGknkbNF9JZg6LgozzN1zNIP30F+qzj0qvnTtxy P3LMaXtApDXIXjClLfhc6x960lZyAcr1tBf45VUCqsAIPL0Rk/9rtnYDhsiPgy7xuboEtR91 ocEVW2TBa+ZNbvesVmT6e03OemDeJcVuCrhK/gi//PikGc1mUUBcqmxwZsXdHe4E+x6LEWWe 3Xshs4NHnoLvgolUOPqlUaCUSVIanapX6M84yk7B5i8AYfCQICtmr2B0z2hEp1YfGAVQm2KC mriIoWYR+8XOmXVOd5kijVCVL67SoZn2wvprx7/07MgL+za/WoTupvnkdR0/ObOjgphyDpvE s6h3nGRGmFognsTFXhxx7F6uUU7y1GZ0KE+jeYfDs1W//oOUwE0MtnXwOV+TtzzQQncZczad FHzSdqvBXQ9T8k63sQVS0d7AdSryB7ZjASwBLpArLGBDYco86/al1z2Pcd61z6Sz60mlVghX o1KMUWpg6d+807YAIufwBbRrLqjaalJhH2Fz2yE12fb5CmwMSZ1WKTBBjUEY1fO6M/+/gXER qOvDrIuNk1Az9SDI+1EcI6hlk1IEdHkPtmWeGetgyGoHx/d3beLdovtZyMY2A3SDUEFl0YY+ nPVfRMmCHKZqnnFRCdrCUqpZkrt9edkr3buVUsz1QCNcwtk05K6/xcUgbqXTPZAlqkctnIHr DN5VE24w8qQC9eEoF95e75AZNom/Fpd/WfQtgg4IZn5aq4+2hgRdAN4u06o3BJyYmlZue4tq n5ijA97KKbDlUhEayvdxpfofLveNmj1+hmrLa/QwFDXltiMqO8J77wjplPvsRvMdAJq+mh70 9RTz3qX54nbRAsUX5XrV08r9h9877jEayg57onQ2DViK66x+jPF3tsoAqMixHPCN59ENKSeF QjuVcgeL8erIe0u3VOua1NMPexf8rI1I9Lzb+GPi+agOOdtmi7jjHwSutgslBLRsXMmGqiUg 8VgobnQxAaMWjbigU30t8n2ndsBfjQOBi+kzjCiAodNZ6p0dIJNCGG0IsTxyM8t4vylE3Ne6 lOnAEsLncGzfh/HJUf52hdR0lhRpHiPli6xzjgymDYs5Pn6vmSG06H5eRwLN3QeDnJmik3tJ pTyitQyU02hbgxvnxygrxWfpeATtOF0KG/dRl1Ndi79IjR5U6e+gbGFZtZG9JIitSg/vP2UW VmBUfa9phIb13imBG5C3HUgcCnsvJzlnht8gWbbLXBprXOfd9sijRvY4dXdQ7ZW0F9kDGFjj T7NB1+gedyt1dqRnpbH9Ou5UiqtW4ZSfi/i0Y6b/HHjtCs6XFvmxaD1xoGvGBNyySLh0th2S SjEyXS0KpLm0ai3K6Mvf0VlAkP999svH4h/loUqg5RDkXMeh5iT4T8GiTKpaYQdif+4NSBdA 2JXkLu3qED/1UZuL2yE3df8X3SZmY56YsWiJ3gRwmQ75txLD6Gd6PpFmzF0qxy2t1G0A7A1k zEDxP8p8HNfjfsOvV9n1SGaGLEVAQ9dOQTjkh2J65a1q6AdNwPNOfCgkVFzm9ysFuTIuQVbQ nf+ZtEnGQd/68x+NBTH13i5ueSGMJHAKNkUsBOTiRLJieNYfYkwmvQ9ji1iIWvhvHchxr1zn Vl00Jq9poTCN3R18ffzHEtDLjOsLZB2mHmlneNEk82RxYzqAphxBmBBQs7zVfzxWDMK6aa8a kDXQWV68CvEX+KYR1PX6V86/SyTVcrwbDfOeiFflZI7FXz/bARemFxGAmt8x8ZjUFjsnIu7K A94/mxDuAC+8EcdjLIwcUG4CD+XpR/0OGhuDsHDakMHtEcaoB6FVK7WpuNrQ3MHotv49lHLc irDIF0WRWARBh7dXwClZ+bxo4mGq6/CX6K/N6ecOOrV77UPC7HQg8roi9UDnX7EN93TbCM6X rtrhwwaByA/Q4OAxH0OU3BFzXuTKZPL4k7toGsv6ZnuuPXzBFC1vNXJUeADd441vUjx2PbmV abYki99LXwwOook43jOxfBf2VcTj3orbDyxCfEbsjaLSqvMm6hRBhpdaiVpNcIO4bhulg9Kc dXWjN/4zNsaxrY8FktFWFr9m8qoed1CImezM0nCDVqKM7LOLCPCwsX+a6exAbNKi+Acuxq1s DedW0jtW1bL3yHuTAyqOPpQgTuzORVfvMSid08oBzS/Ct3hbRK/PZl8ijh3ibw4i3XWNHINZ Dhxd0Qey9/YpShcg/h5BylA9i8/dbjCy3vftrOIbMpK4p4JSmxumulX4Wo30e5Q5SBAH7lun TfK68Vpuxegm/WOzTxuVFxPrCxKjcSFpxYHW+2R+59eVHLD5B9I43+XDkFAvN9iENzop+ZVz vDAkavyLHFJ9Neerq5+T4DEbdmKNnYsK0+jADnPEA4MViKmL0najk1Z1e6Xrzib8sR8pZ/rl 54DDLRcURZmc5FSQlQgF9sELpBtWzoimrPOl88E6021qxzJTdlbtJTKPhp9KfrqITedy7JDY klRqVsdBYsWP4z/nUdlbwsi9GwrM0/ZXNQIvSg4KwFo/ANC939xSmB10EXgOFvF3Q==
- Ironport-sdr: 46fZkuPZrI9dZ4CW5YzRjjY64EOzgpCZUX0b1dlNV0Gwr/2JuMXdUOt8H7B+veBOQmAMFn4/+8 s/HGoPOAA0gEXeppUXQqJScDFlMH6tPfukpxIKm7QoXcEoI4GOV1llhvZ0Fm8N2r2Rr33BLXKE E7HddR/2XGhVjUA5s6Ftvl6689DLLSLCid7s2iMCL+TB2aluPj5Fb4o8L1fcwrHlWoJIaEaJM1 3ilio1MrtBqbnu2pQyiK6M1H28+U3Et69W5L35ksQ3bI9AbOc+BdpUrBbem5M/jbL8SduKDkoE h7FBelqhpSXVmihKXFGmWHNI
I'm not sure to understand what is your problem. A minimal self contained example that we can run and that shows your problem would
be perfect.
What I can tell from the code I see:
You are computing a 3 triangulations that consists of tetrahedra.
`finite_facets()` provides the range of face that are not connected
to the infinite vertex. A facet is given a std::pair of a Cell_handle
(tetrahedron) and the id of the vertex in the cell opposite to the
facet. About the access to the vertices of a facet, the %4 is here so
that the index stays in [0,3] while ignoring pair.second (that might or
might not be infinite).
HTH,
Sebastien.
On 7/5/22 14:46, KL ( via cgal-discuss Mailing List) wrote:
Hi,
The issue: while the code below works, it returns too many triangles for a simple shape.
So, for a square, instead of two triangles made by one of the diagonal of the square, I get four tris, because the two diagonals are connecting.
I get my vertices using:
using K = CGAL::Exact_predicates_inexact_constructions_kernel;
using Triangulation = CGAL::Triangulation_3 <K>;
using Cell_handle = Triangulation::Cell_handle;
using Vertex_handle = Triangulation::Vertex_handle;
using Locate_type = Triangulation::Locate_type;
using Point_CGAL = Triangulation::Point;
...
Triangulation tri;
...
// load data
...
//loop
for ( auto &e : tri.finite_facets() ){
auto &vert = *( e.first );
auto starter_id = e.second;
for ( u32 x = 1; x < 4; x += 1 ){
auto &pt = *( vert.vertex( ( starter_id + x ) % 4 ) ); // saw that line on the stackoverflow -> 3 finite vertices and 1 infinite vertice, hence the 4.
// store my result
}
}
I need the triangulation to be the minimum number of tris for any given shape. A square should have 2 tris.
I read the doc, but it is not among the easiest to wrap one's head around.
What should I modify, please?
Thanks
--
MKL
--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss
- [cgal-discuss] 3D triangulation: absolute minimum of facets, KL, 07/05/2022
- Re: [cgal-discuss] 3D triangulation: absolute minimum of facets, Sebastien Loriot, 07/05/2022
- Re: [cgal-discuss] 3D triangulation: absolute minimum of facets, KL, 07/05/2022
- Re: [cgal-discuss] 3D triangulation: absolute minimum of facets, Sebastien Loriot, 07/05/2022
- Re: [cgal-discuss] 3D triangulation: absolute minimum of facets, KL, 07/05/2022
- Re: [cgal-discuss] 3D triangulation: absolute minimum of facets, Sebastien Loriot, 07/05/2022
Archive powered by MHonArc 2.6.19+.