Subject: CGAL users discussion list
List archive
- From: Andreas Fabri <>
- To:
- Subject: Re: [cgal-discuss] Is the following CDT buggy?
- Date: Tue, 15 Mar 2022 17:20:04 +0100
- Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=Pass ; spf=None
- Ironport-data: A9a23:kQtEc6Ky1W8vE/p/FE+RbZMlxSXFcZb7ZxGr2PjKsXjdYENS12EDx 2UaUT+Ca/7bYGekctF/a4vjoRsDvMTSx9QwQQMd+CA2RRqmi+KVXIXDdh+Y0wC6d5CYEho/t 63yTvGacajYm1eF/k/F3oDJ9CU6jefRLlbFILas1hpZHGeIcw98z0M78wIFqtQw24LhWFrQ4 YiaT/D3YTdJ5RYkagr41IrY8HuDjNyq0N/PlgFWiVhj5TcyplFNZH4tDfnZw0jQHuG4KtWHq 9Prl9lVyI92EyAFUbtJmp6jGqEDryW70QKm0hK6UID66vROS7BbPqsTbJIhhUlrZzqhmNB9k P9ouZWKFgYZbr/um7pEDghdKnQrVUFG0OevzXmX6IqWyBCYKz3pyvRqSUYrIcsf5OYxB2xSn RAaAGxTKEnb2qTom/TiELAEascLdKEHOKsWvH16wDjCS/gvS4rCa6rH4tpVwC0hiMlFAfHEd owSbj8HgBHoOUMVaw1IVMhWcOGAqH7aTB9DkH2sqqMa33HX1ikryKrhP4+AEjCNbZ8NxhvJz o7cxEzyDRgecdCe0jGY6WmEnf7KhSq9WYQIFbT++OQCvbGI7jVVDRpLCwP9pPC4jgu5Rs4ZL FIUvCwjscDe6XBHUPH6XxHgiluItyRDWuoAAc85tljQ26/ttlPx6nc/chZNb9kvtckTTDMs1 0OUk96BOdCJmODPIZ563urNxQ5eKRT5PkdeNX9bHVFtD83L+thj0EmnosNLTfbt5uAZDw0c1 BiklkDSbZ0pkMMK0+Ce9FTCijO3oZ6houUdt1+KBjrNAu9Rz2ihasmk5B319/9GJ4vxc7Vsl CFCns/DsrhICJiMkGmKXflLG6+pofCILFUwYGKD/bF/p1xBGFb5Jui8BQ2Swm81bK7onhe1P SfuVft5vsM7AZdTRfYfj3iNI8or17P8Mt/uS+rZaNFDCrAoKlPYpnA2OR7Ngj+1+KTJrU3ZE cfAGSpLJSZAYZmLMBLnGb5GuVPW7nthmD+JHcqTI+qPi+LDPiLIIVv6DLd+RrphtP7VyOkk2 9FBN4ORwg43bQENSnm/zGLnFnhTdSJTLcmu9aR/L7ffSiI7Rj1JI6KPkNsJJtw094wIx7yg1 izsCidwlgCg7VWZcl/iQi44N9vSsWNX9i9T0doEZgrzhBDOoO+HsM8iSnfAVeJ+pbI5kqIsE aVtlgfpKq0ndwkrMg81NfHVxLGOvjzy7e5XFyb6MjU5YbB6QAnFpo3tcgf1rXhcASewsMYiu ayu3wjHRoASAQ9lCZ+OOv6oylqwu1kbmf5zDhqUfIYIIhW8/dg4MTH1g982P9oIdUfJyAyF2 lvEGhwfv+TM/dM4qYGbmaCeoo61OOJiBU4GTXLD5LO7OHCCrGquyINNSsiSejXZWD+m8amuf 7wEnf7xN/gKgExbvYN3D7Fx3OQ14N62/+1WyQFtHXPqaVW3C+IxfifZhpMf5qAUn+1Xowq7X E6L6+J2A7TRNZO3CkMVKSokcv+HiaMelA7N4KlnO079/iJ2ouaKCB0AIxmWhSVBB7JpK4d5k /w5scsb5gHj2BonNtGK0nJd+2iWdCBSCuMiscxEWMnugwsvj1ZffdrbFC+w5pyTMo0ePk4vK z6SpazDm7UMlxWYLCpvTSDAjbhHmJADmBFW11tcdV6HrdrygKNl1hNm9zlqHB9eyQ9K0r4oN 2U3ZVd5I76Co2VhiMRZBTz+QkdEA0DHogr0wloN0WrEUw+vSGyLKmAhYL7f8Ecc+mNaXz5a4 LDImTm/Dma2JJn8jnkoREpoi/3/VtgvpAfMr8CQGZjXFZcNZze40LSlYnAFqke8DM487KEdS TKGIAqkhWzH2S8sT2kTUs+f0ulOEFaBLW1GBPZ84OUOA2GafjyusdRLx4ZdZesVT8EmM2fhY yCtGi6Lfxuz0yODsioKC6cHP7hujbgi499qlnbDOzsdq7XGxtZ2mMu4y8U97VPHh/1hnME4J 5/LZjyLGXCXn2oSkGjIxCWB1qxUfvFcDDDBMCuJHCnl2n7NXCyAsa3/71dsg0ioDQ==
- Ironport-hdrordr: A9a23:j4tXDqrYlZus5OOdPl92dgIaV5rPeYIsimQD101hICG9Evb0qy ncpoV/6faUslkssR4b+OxoVJPwJ080sKQFmLX5Xo3DYOCFggWVxehZhOOIqQEIWReOkdK1vp 0QEZSWZueAaWSSwfyKgzWQIpIH/9jC2qWpieDC724Fd21XQpAlyRh5TjuBGkd7X01nCPMCZe Khz/sCmAHlQ00+Su6SOkMoNtKz2uHjpdbYWlo9PjYC0iOpth/A0s+BLzGomigkFwl36ZMZzE D5v2XCl96ej80=
- Ironport-phdr: A9a23:bhj6+hPRY//XePtBjMYl6nb1BBdPi9zP1u491JMrhvp0f7i5+Ny6Z QqDv68r1QSRFtWGo9t/yMPu+5j6XmIB5ZvT+FsjS7drEyE/tMMNggY7C9SEA0CoZNTjbig9A dgQHAQ9pyLzPkdaAtvxaEPPqXOu8zESBg//NQ1oLejpB4Lelcu62/6s95HJfglFiiaxbbxxI Ri2sA7cqtQYjYx+J6gr1xDHuGFIe+NYxWNpIVKcgRPx7dqu8ZBg7ipdpesv+9ZPXqvmcas4S 6dYDCk9PGAu+MLrrxjDQhCR6XYaT24bjwBHAwnB7BH9Q5fxri73vfdz1SWGIcH7S60/VDK/5 KlpVRDokj8KOCI6/m/ZhcN/gr9VrgyvpxJ/wYHZe5uaOOZkc67HYd8XS2hMU8BMXCJBGIO8a I4PAvIdMOZesob9vUUBogGlBQKxBO3g0CRHhmX33aYn1OkuCwfG0xE6H90QqnvUt8/5NKkIX u+u1qnIzC7Ob/xI1jf67YjHbAwhoeuMXLJ+a8Xe1VUvFwTfjlSQs4DqIzSV1uEUvmWd8uFvW v6hhXQ9pAFtvjig2N0sio/Ri48Vy13K+zt0zoYpKNCmVEJ2ZdCqHZpeuSybOYZ6XMcvT391t Cs11rAKpIK3ciYWxJkkyRPSaeKLfoaG7x7/SuqcJypzinxieLK6nRmy8E6gx/XzVsm1zFZKr jdFncLWun8R0BzT79CLRedy/ki73TaDzQfT6vtLIUwslKrbLYQhwr8qmZYJv0TDBCn7k1j1g q+Obkgp9eml5/76brn6ppKQLY55hw/kPqgzhMCzHPw0PhYUU2SG+emwzqDv8E77TblQkPE7k qbUvZbHLsoBvKG5GRVa0oM75ha/ETim1NMYkGEaIF5ffRKHiIjkNlDSLPziCPewmVWskDNyy /DBI73hHo/BLn/dkLfhZ7p97VBTyBYrwdBe4ZJbFK0BLer8V0LwrtDUEwM1Pgyuz+vmBthxz J4SVXySDqOBNaPdq16I5uYhI+mWY48VvS7wJOIj5/Hzl382h0USfbKz0ZsZcn+3Be5mI0SCb HrjnNgBFmIKsRQiTOzwiFyCXyRcaGquX6I7/DE6CJiqApvERoComLCB2jy0EYVMZm9aElCMD Wvod4KcVvsRcC6SONNukiQYVbi9TI8szQ2htADgxLpjN+bb5ywYtYn/29ho/O3Tjgo/9Sd0D sSYy2GCVXt4nmIORz8s3aBwu1ZxylmZ0fswvvpDCNYG5+9VShxoctnH3uliApbzXBjAd5GHU hG9U9C+CHYwSNw2hNQBakI4F9S5hQ3YxHmXBaQInZyXAZhh8r7Ax2OjYIFm2nPe3e8giUMnS 41BLyq9l6tn/k/SAYDO1E6WnqLveaUH1zPW7zS+y3GTtn1VQBIlUbnZRWtNIQzNvNHh7wXDS aWvAPIpKExa2MuaI+xLbNPuyl5JTfOmNNXFaH+qgDSMAkODybqIKYbrYG4Axz71CU4ekglV8 2zVGxI5A3KIpW7EATV1XXzmaVntub13rn+hQ0YviQ+DZVdg/7W4/RsYmeaNRfoYwrUeqWEqr DAiTwX15M7fF9fV/1kpR65be95ouD+vtErcvg15ZNm7Krx6w0QZeEJxtl/v0BN+DsNBl9Irp TUk1lk6MrqWhXVGcT7QxpXsIvvPMGCn9RajcaPfwRfQ2dyM+48A5fM9pkn5rQ+gHVYl6WQh2 N5QgDOH/pufNAMJStrqV1ovsR1zprXUeC44sorS2WdhOLLyvDbIwdMBC+Yizxu8ZcZROaiYE xXjVcYdApvmM/Qkzn6uaB9MJ+VO7Og0MsehIuOBw7KuNf18kSiO1CJC590jigSJ/it4D+nVw 9AC3fHe2AabP9vlpHGmtM2/2YVNZDVIW3G61TChHolaIKt7YYcMD26qZcyx3NR3wZD3CTZe8 xa4ClUK1dXMG1LaZkHh3QBWyUUcoGC20Sq+wTtuljg1r62ZlCXQyuXmfRADNyZFXm5nxVvrJ IG1iZgdUi3KJ0Agkhe/6E/hgaZSrr5+B2beRkJFYzLnIWhpTqyqp/yJZMsOoJIkvCNLUfitN EiAQ+2YwVNS2CfiEm1CgTEjIm30/Myhwlok0yTEdi8gyRiRMdt9zhre+tHGEPtY3z5cATJ9l SGSHF+3ed+g4dSTkZ7H9OG4TWOoEJNJIkyJhcuNsjW24WpyDFixhfe2z5fuFwQg3CbgkddjX z/JhBnxZYzmy76rP+tsYk5yFRn37M8wSeQc2sMgwYod33QXnMDf9HwLi2r6KpNV0KjkbVICS DkOzsLP8Qbs01FkNGPPzIX8HCb4oIMpd5yxZWUY3Tg45sZBBfKP7bBKqiByp0KxsQPbZfUVc i418fI18zZahugIvFBo1SCBGvUIGkIeOyXwlhOO5tT4raNNZW/pf6LinEZ5mNmgCvmFrGQ+E D7wfJs4ECZrqMt2OkjN+HL+7YTpZMPBY9sYqhqOglHLiO0dJJ8ql/UMjDZqIiqh4zt8k7N90 0cohMvj9IGcTgcltLq0GBtZKiH4a4sI9zfhgLwf1sea0oazH4lwTzACXZ/mV/WtQ3oZsfXqM RrLESVp9indQOKCW1HErh489iGqcdjjLXycKXgHwM83QRCcIBYamwUIRHAgmZV/EAm2xcvne UM/5zYL51e+pAEfr4AgfxT5TGrbox+lLzkuT53KZhNQ4htP7lyTP8WU9ONbECxf+5C9th2DI 2eHYB5ZS2oOXwbXYjKrdqnr/tTG/+WCU6C3KffUbLyV7+JXXe2J7Z2i1Y5r4yydOMyEIn54H rsw3U8JDhUbU4zJ3j4ITSIQjSfEacWW8Qy99iNApce66P33WQjr6NjHG/5IPN5o4RzznbabO rvamnNiMTgBnMBppzeA2P0F0VUVkS0razS9DeFKq3vWVKyJ0rdLFUweYiJ3ccxV8+Q7wAkLP 8PHwsvpzOx9h/8xTVBbVRrnl93MB4RCInnhZgOBXx3Nbe/ZY2STnoasO+u9UeEC1bQO8Uzh/ 2nEVR+6eW/Z327gUxTlWQ1VpAecOhEW+IS0cxI3THPmUMqjcBqwdtl+kTwxx7QwwHLML28Vd zZmIQtLqfWL4CVUj+8aeSQJ52d5LeSChyeS7vXJYpcQv/xxBy1oluVcqH0kwrpR5StASbR7g izX5tJpplinlKGIxF8FGFJWrS1XgYuQoUh4Efyf854dBy6C+RsM6SCXFghMoMVlTNvip+EYy 9TCkr7yNCYX89/Q+ph5ZYCcI8aGPXw9dBvxTWKKSlJdEHj1bieG3RUO9ZPavmeYpZU7tJX2z Z8HS7sBEUcwCutfEENuWtoLPJZwWDog17+dls8Bo3Sk/3yzDI1XuI7KUvWKDLDhMjGc2PNNY RcSzLrja4oaPJf680NvbVx3gJ7bFUPbQddXs2tqaQp+8yAvuDBuC3Y+3U7ocFbn+HgICfu9h QI7kCMnJ+EqqWm0pVI+J16PozYs1k4vmZPjjCzbI1uTZO+gGIpRDSTzrU04NJj2Fh10YQOFl kthLD7YRrhVgtOIlEgyzgrb5MUXX/tVTKkBZwIMg/aJZ7Mu3EgO8k1PKmdI6OzABId4hQUje oKrtWMG0AVmPoddzUn4K6dOy1VMnLOAtyS00foghgQZIhRVmF4=
- Ironport-sdr: qCM/7OfA8xmtD0n7393q1VaRnTSrF7A4Rr4Ihk0CQXizSdzsvcW/BSI7QOWNoAkdXyNGbPc6mh GTfIBffcuRWKm2GLJN3+S/cybtvSIGbzIA4ZFnZwuee/clhQefOQdl90q4wIHIKfzLkNWe+npr bFzDFxfwN9CNS51id4HlYJCTFwg3ZJg9nft0IipEWNCZaphzBB5KkBvUL+h0n5ezxAJ0+fTbmq gwb7zfFnEudD4XJrYF0WyahWi+1451rZ9pfVj2xKF52hh/ItVBe7lSOWLoHi5cjr7/3CcLsuYp g0uz1ju9C8lgrmI1kP0luhhv
Hi Bob
one PR later I now understand where your initial problem comes
from.
You commented refine_Delaunay_mesh_2().
This function internally does a traversal
as the mark_domain() function.
When not called there is no domain marked.
Best,
Andreas
ps.: The above mentioned PR is here https://github.com/CGAL/cgal/pull/6412
On 3/14/2022 12:46 PM, Andreas Fabri
wrote:
This should do it, although not yet elegant
https://gist.github.com/afabri/82e9680181eea814fb488f96ceca8c12
On 3/14/2022 12:09 PM, Bob Bill (
via cgal-discuss Mailing List) wrote:
Thanks Sebastien and Andreas,
I see what Andreas means, indeed if I modify my example accordingly and then call refine_Delaunay_mesh_2() I have the expected error " no member named 'set_in_domain'...".
I tried to adapt it by adding typedef CGAL::Delaunay_mesh_face_base_2<K> Fb; in place of typedef CGAL::Constrained_triangulation_face_base_2<K,Fbb> Fb, but I get an error "no member named 'info' in..." coming from the start->info(), so I don't really know how to adapt it in order to have both functions.
Should I only change the typedefs?
Best, Bob
Il lunedì 14 marzo 2022, 08:59:26 CET, Andreas Fabri <>
ha scritto:
Hello,
Note that write_VTU is documented in the CGAL package 2D Conforming
Triangulations and Meshes.
When you call a function like refine_Delaunay_mesh() you define the
domain (for example the interior
of your polygon, as a domain, otherwise the convex hull of the points is
the domain.
When you follow the example Sebastien just gave you must adapt it, as
the Delaunay_mesh_face_base_2
class has a function to set the domain.
https://doc.cgal.org/latest/Mesh_2/classDelaunayMeshFaceBase__2.html#ac4a80317799961e7830b9690121d6bdc
Best,
Andreas
On 3/14/2022 8:40 AM, Sebastien Loriot ( via
cgal-discuss Mailing List) wrote:
> A triangulation of a point set is a triangulation of the convex hull
> of the points. If some constrain edges are added, you then have the
> guarantee that those edges will be present in the triangulation of the
> convex hull of your points.
>
> If your constraints define a closed shape then you need a flooding like
> algorithm to only pick the triangles that are inside the shape, like in
> this example:
>
> https://doc.cgal.org/latest/Triangulation_2/Triangulation_2_2polygon_triangulation_8cpp-example.html
>
>
> Best,
>
> Sebastien.
>
> On 3/13/22 02:46, Bob Bill ( via cgal-discuss
> Mailing List) wrote:
>> I am given the vertices of a polygon, and I need to mesh it. I am
>> following one of the first examples in the CDT package documentation.
>>
>> The code is the following:
>>
>> #include <CGAL/Exact_predicates_inexact_constructions_kernel.h> >> #include <CGAL/Constrained_Delaunay_triangulation_2.h> >> #include <CGAL/Delaunay_mesher_2.h> >> #include <CGAL/Delaunay_mesh_face_base_2.h>
>> #include <CGAL/Delaunay_mesh_size_criteria_2.h> >> #include <CGAL/IO/write_VTU.h> >>
>> #include <iostream>
>> typedef CGAL::Exact_predicates_inexact_constructions_kernel K; >> typedef CGAL::Triangulation_vertex_base_2<K> Vb; >> typedef CGAL::Delaunay_mesh_face_base_2<K> Fb; >> typedef CGAL::Triangulation_data_structure_2<Vb, Fb> Tds; >> typedef CGAL::Constrained_Delaunay_triangulation_2<K, Tds> CDT;
>> typedef CGAL::Delaunay_mesh_size_criteria_2<CDT> Criteria; >> typedef CDT::Vertex_handle Vertex_handle; >> typedef CDT::Point Point;
>>
>> int main()
>> {
>> CDT cdt;
>> Vertex_handle va = cdt.insert(Point(-1,-1)); >> Vertex_handle vb = cdt.insert(Point(-1,-0.25714));
>> Vertex_handle vc = cdt.insert(Point(-1.5,-1.4));
>> Vertex_handle vd = cdt.insert(Point(1.3,-1.2));
>> Vertex_handle ve = cdt.insert(Point(1.3,0.5));
>> Vertex_handle vf = cdt.insert(Point(1, 0.45714));
>> Vertex_handle vg = cdt.insert(Point(1, -1)); >>
>> cdt.insert_constraint(va, vb); >> cdt.insert_constraint(vb, vc); >> cdt.insert_constraint(vc, vd); >> cdt.insert_constraint(vd, ve); >> cdt.insert_constraint(ve,vf); >> cdt.insert_constraint(vf,vg); >> cdt.insert_constraint(vg,va); >> std::cout << "Number of vertices: " << cdt.number_of_vertices() << >> std::endl;
>> std::cout << "Meshing the triangulation..." << std::endl;
>> // CGAL::refine_Delaunay_mesh_2(cdt, Criteria(0.125, 0.5));
>> std::cout << "Number of vertices: " << cdt.number_of_vertices() << >> std::endl;
>>
>> std::ofstream output("intersected_mesh.vtu");
>> CGAL::write_vtu(output, cdt); >>
>> //Show the vertices of the triangulation >> for(auto it =
>> cdt.finite_faces_begin();it!=cdt.finite_faces_end();++it){ >> std::cout << cdt.triangle(it) <<'\n';
>> }
>> }
>>
>> If I look at the .vtu file, it seems to have the correct shape (see
>> attached image). However, after I printed the vertices for each
>> triangle, it appears that one element of the triangulation is also
>> the one with vertices
>> (-1, -0.25714)
>> (1, -1)
>> (1, 0.45714)
>>
>> which should not be there since is not inside that polygon. Isn't
>> this a bug? If not, how can I create a mesh for which this is
>> avoided? I mean, I'd like to have all the triangles inside the domain.
>>
>> Best regards,
>> Bob
>>
>> --
>> You are currently subscribed to cgal-discuss. >> To unsubscribe or access the archives, go to >> https://sympa.inria.fr/sympa/info/cgal-discuss >>
>
--
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project
--
You are currently subscribed to cgal-discuss. To unsubscribe or access the archives, go to https://sympa.inria.fr/sympa/info/cgal-discuss
Note that write_VTU is documented in the CGAL package 2D Conforming
Triangulations and Meshes.
When you call a function like refine_Delaunay_mesh() you define the
domain (for example the interior
of your polygon, as a domain, otherwise the convex hull of the points is
the domain.
When you follow the example Sebastien just gave you must adapt it, as
the Delaunay_mesh_face_base_2
class has a function to set the domain.
https://doc.cgal.org/latest/Mesh_2/classDelaunayMeshFaceBase__2.html#ac4a80317799961e7830b9690121d6bdc
Best,
Andreas
On 3/14/2022 8:40 AM, Sebastien Loriot ( via
cgal-discuss Mailing List) wrote:
> A triangulation of a point set is a triangulation of the convex hull
> of the points. If some constrain edges are added, you then have the
> guarantee that those edges will be present in the triangulation of the
> convex hull of your points.
>
> If your constraints define a closed shape then you need a flooding like
> algorithm to only pick the triangles that are inside the shape, like in
> this example:
>
> https://doc.cgal.org/latest/Triangulation_2/Triangulation_2_2polygon_triangulation_8cpp-example.html
>
>
> Best,
>
> Sebastien.
>
> On 3/13/22 02:46, Bob Bill ( via cgal-discuss
> Mailing List) wrote:
>> I am given the vertices of a polygon, and I need to mesh it. I am
>> following one of the first examples in the CDT package documentation.
>>
>> The code is the following:
>>
>> #include <CGAL/Exact_predicates_inexact_constructions_kernel.h> >> #include <CGAL/Constrained_Delaunay_triangulation_2.h> >> #include <CGAL/Delaunay_mesher_2.h> >> #include <CGAL/Delaunay_mesh_face_base_2.h>
>> #include <CGAL/Delaunay_mesh_size_criteria_2.h> >> #include <CGAL/IO/write_VTU.h> >>
>> #include <iostream>
>> typedef CGAL::Exact_predicates_inexact_constructions_kernel K; >> typedef CGAL::Triangulation_vertex_base_2<K> Vb; >> typedef CGAL::Delaunay_mesh_face_base_2<K> Fb; >> typedef CGAL::Triangulation_data_structure_2<Vb, Fb> Tds; >> typedef CGAL::Constrained_Delaunay_triangulation_2<K, Tds> CDT;
>> typedef CGAL::Delaunay_mesh_size_criteria_2<CDT> Criteria; >> typedef CDT::Vertex_handle Vertex_handle; >> typedef CDT::Point Point;
>>
>> int main()
>> {
>> CDT cdt;
>> Vertex_handle va = cdt.insert(Point(-1,-1)); >> Vertex_handle vb = cdt.insert(Point(-1,-0.25714));
>> Vertex_handle vc = cdt.insert(Point(-1.5,-1.4));
>> Vertex_handle vd = cdt.insert(Point(1.3,-1.2));
>> Vertex_handle ve = cdt.insert(Point(1.3,0.5));
>> Vertex_handle vf = cdt.insert(Point(1, 0.45714));
>> Vertex_handle vg = cdt.insert(Point(1, -1)); >>
>> cdt.insert_constraint(va, vb); >> cdt.insert_constraint(vb, vc); >> cdt.insert_constraint(vc, vd); >> cdt.insert_constraint(vd, ve); >> cdt.insert_constraint(ve,vf); >> cdt.insert_constraint(vf,vg); >> cdt.insert_constraint(vg,va); >> std::cout << "Number of vertices: " << cdt.number_of_vertices() << >> std::endl;
>> std::cout << "Meshing the triangulation..." << std::endl;
>> // CGAL::refine_Delaunay_mesh_2(cdt, Criteria(0.125, 0.5));
>> std::cout << "Number of vertices: " << cdt.number_of_vertices() << >> std::endl;
>>
>> std::ofstream output("intersected_mesh.vtu");
>> CGAL::write_vtu(output, cdt); >>
>> //Show the vertices of the triangulation >> for(auto it =
>> cdt.finite_faces_begin();it!=cdt.finite_faces_end();++it){ >> std::cout << cdt.triangle(it) <<'\n';
>> }
>> }
>>
>> If I look at the .vtu file, it seems to have the correct shape (see
>> attached image). However, after I printed the vertices for each
>> triangle, it appears that one element of the triangulation is also
>> the one with vertices
>> (-1, -0.25714)
>> (1, -1)
>> (1, 0.45714)
>>
>> which should not be there since is not inside that polygon. Isn't
>> this a bug? If not, how can I create a mesh for which this is
>> avoided? I mean, I'd like to have all the triangles inside the domain.
>>
>> Best regards,
>> Bob
>>
>> --
>> You are currently subscribed to cgal-discuss. >> To unsubscribe or access the archives, go to >> https://sympa.inria.fr/sympa/info/cgal-discuss >>
>
--
Andreas Fabri, PhD
Chief Officer, GeometryFactory
Editor, The CGAL Project
--
You are currently subscribed to cgal-discuss. To unsubscribe or access the archives, go to 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
-- Andreas Fabri, PhD Chief Officer, GeometryFactory Editor, The CGAL Project
--
You are currently subscribed to cgal-discuss.
To unsubscribe or access the archives, go to
https://sympa.inria.fr/sympa/info/cgal-discuss
-- Andreas Fabri, PhD Chief Officer, GeometryFactory Editor, The CGAL Project
- [cgal-discuss] Is the following CDT buggy?, Bob Bill, 03/13/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Sebastien Loriot, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Andreas Fabri, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Bob Bill, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Andreas Fabri, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Andreas Fabri, 03/15/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Andreas Fabri, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Bob Bill, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Andreas Fabri, 03/14/2022
- Re: [cgal-discuss] Is the following CDT buggy?, Sebastien Loriot, 03/14/2022
Archive powered by MHonArc 2.6.19+.