Subject: CGAL users discussion list
List archive
- From: Sebastien Loriot <>
- To:
- Subject: Re: [cgal-discuss] Issue with computing volume of a surface mesh
- Date: Mon, 11 Jul 2022 08:45:34 +0200
- Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=Pass ; spf=None
- Ironport-data: A9a23:GG7SK6NtzkmXlR7vrR2XkcFynXyQoLVcMsEvi/4bfWQNrUpx0DMEy WoXXG6BPavZM2rxeogjPNywpE1Q6pLSnIJkTnM5pCpnJ55ogZqcVI7Bdi8cHAvLc5adFBo/h yk6QoOdRCzhZiaE/n9BCpC48T8kk/vgqoPUUIYoAAgoLeNfYHpn2EgLd9IR2NYy24DmWlrV4 LsenuWGULOb824sWo4rw/nbwP9flKyaVOQw4zTSzdgS1LPvvyF94KA3fcldHFOkKmVgJdNWc s6YpF2PEsw1yD92Yj+tuu6TnkTn2dc+NyDW4pZdc/DKbhSvOkXe345jXMfwZ3u7hB2MuPpV2 MREu6ePRFkUZvWQteUjTBxhRnQW0a1uoNcrIFC6uM2XilLDKj7in6goA0YxMokVvO1wBAmi9 9RCcGFLPk3F3rjmhu7mIgVvrpxLwM3DJ4cYoHBn0XfcCd4pRJnCR+PB4towMDIY2p0WQK+GN 5dxhTxHQgbZMzx3AVUtLI8SmryvgHv4Syd8twfAzUYwyzGLkFYZPKLWGNHac9jPScROlVuDv Urd7mHhC1cbMsaewHyL6BqRavTnmCr6XMcNEeT9+KI3xlKUwWMXBVsdUl7TTeSFZlCWXvV1D 1wF8yQVkqVsyRTzEvT8YQK9mSvR1vIDYOZ4H+o/4QCL76Pb5QeFG2QJJgKtjvR25KfaohR6h je0c8PV6S9H6+LKFCrMnluAhXbjZnhPdD5qiTosFFNdu7HeTJcPYgUjp+uP/Yawh9zxXC71m nWE8HV4iLIUgsoGka68+DgrYg5ARLCZH2bZBS2NBgpJCz+Vgqb7PuREDnCFtJ59wH6xFAXpg ZT9s5H2ABoyJZ+MjjeRZ+4GAauk4f2IWBWF3wMwQ8Z+p27wpyD8FWy13N2YDBc5WirjUW+5C HI/RSsMjHOuFCD3MvMtMtrZ5zoClPCwRIiNug/ogipmO8AtLmdrDQlhYkmf222FraTfufBXB HtvSu71VSxyIf0/klKeHr5BuZd2mH1W7T6MHfjTkkX/uZLDNSX9YepUYDOmMLplhIva+1692 4gEZ6O3J+B3CrKWjt//qt5Nczjn7BETWfjLliCgXrTdeFU+RDpwVqS5LHFIU9UNopm5X9zgp hmVMnK0AnKl7ZEeAQnVOH1ldp31WpNz8SAyMSA2bASn3nEiZcCk66JGL8k7erwu9epCy/9oT qldK5/QXKgXEjmXqS4AaZTdrZB5cEv5iA+LOR2jamdtcpNlQTvP5dK5LBDk8zMDD3bsuMZn+ ++g2wrXTIAtXQNnCMqKOvujw0ng73cYkeN2GUDPJ4ALKknr9YFrLQ33j+M2c5lcc0WdmmPC2 l/PUxkCpOTLr4sky/XzhPiJ/9WzDu9zPktGBG2Ev7u7MC/t+GD8k4JNVeC/ey+ECDH59aCkU uViz//mNcoBklsX4ZF3FKxmzP5n6tbi++1awwBjECmZZlinEOk7cHyP3M0KqagUg7EF5E25X UWA/tQcMrKMYZu3HFkULQsjT+KCyfBExWWIvKpteB33tH1t4b6KcUROJB3Q2iZTG714bdE+y uA7tc9KtgGy1kgwPtCdgnwG/miANCZbAaAut5VfH4qyzwR3lxdNZpvTDiKw65aKMo0ePk4vK z6SpazDm7UMmRaYIiRrTSDAjbhHmJADmBFW11tedV6HrdzI260s1xpL/DVrEwlYw32rCQ6o1 rSH6qG0GUmPw9usrM1KXmTpBAQYQRPFpwr+zFwGkGCfRE6tPoAIwKvRJs7VlH31MUoFFtSYw F1c4GngWDfuOsr220PenGZ7/uf7Q4UZGhLqwaia8gfsI3X+STXgi66qI2EPrnMLxC/3aFLv/ YFXwQq7VUE32eP8bUH250l2GIn8kCy5GVE=
- Ironport-hdrordr: A9a23:kH9uRapEAGmdYwpFpN2s7r4aV5oIeYIsimQD101hICG9Afbo8f xG+85rqSMc6QxhJU3I/OrrBEDuewK/yXcY2/hrAV7mZnidhILKFvAG0WKB+UyCJ8SWzIc0vs ddmupFeaXN5DBB7PoSizPIdOrIruP3iZxAyd2uq0uFjjsHV0ij1WpE48qgfXGejTMpOaYE
- Ironport-phdr: A9a23:uZ9nKxTp8W8sSaVF62L4L/p3LtpsoiaVAWYlg6HPa5pwe6iut67vI FbYra00ygOTBsOCtKgP0rOP+4nbGkU4qa6bt34DdJEeHzQksu4x2zIaPcieFEfgJ+TrZSFpV O5LVVti4m3peRMNQJW2aFLduGC94iAPERvjKwV1Ov71GonPhMiryuy+4ZLebxtGiTanfL9/L gu6oQrMusULnIdpN7o8xAbOrnZUYepd2HlmJUiUnxby58ew+IBs/iFNsP8/9MBOTLv3cb0gQ bNXEDopPWY15Nb2tRbYVguA+mEcUmQNnRVWBQXO8Qz3UY3wsiv+sep9xTWaMMjrRr06RTiu8 6FmQwLuhSwaNTA27XvXh9RzgqxVrx2uqQFxzZDaYI+VNvVxYqzTcMgGRWdCRMtdSzBND42+Y oYJEuEPPfxYr474p1YWsBWxGxKjC/n0yj9Uh3/5w6s62PknHwHBxgwgHtUOsGjJp9jyL6cSV ee1zLPJzTXEc/NawzL96IzTfxAupPGDR7Nwcc7LxUYzEAPFi0ydpIr4NDyayuoDqXKU7/Z8V e2xkW4nrRl8rzqyy8oiiITEm4MYxkzZ+Ct3wYs4Kt61RUxmbNO6FJZduCGUOYt0T84mX2xlt yQ3xLMJtJOlYSQHxpopygPbZvGBboOG7BXjVOOLLjd5gnJoYLO/hxCo8Uih0OLwTMe00ExMo yFYkdfMrmgA2wLP5sWDUPdw/Ues1SyS2wzN6exIO085mKrDJ5Miw7M8jJgevEXZEiPolkj6k rWaelgk9+Wn7ensf6/oqYWGN4BujwHzKqQuldK7AeQ/KgUOWnKU+eW41LH680z5RahGguQ4k qTWsJ3WP8sbpqm+Aw9a1oYs9QyzACuh0NQdhXUHLVRFdwybj4XxJV3CPPT1Ae28jlmsijtn2 u7KM777DpjCLnXPiLLhcqx8605Yxgoz19df55dMB7EaJPLzXVXxtN3fDh8lKQC0xuLnBc5y1 oMbQ22PA6uZPLnOvl+P4+IjO++Ma5QNtzbnN/cl/+LujWM+mVIFYKWlxYEXZ2ygHvR6P0WZZ mLhjcsOEWgQugoyVfHliFyZUTFPenayRLk86yogCIOmCIfDXpqijKaA3Ce9BJ1WZ3pJBkqCE XfyJM25XewRYneSPtN5iW5DEqOwTpcokxCorg7zjbR9afHF/zUR8pPl2t8y7OLakVQ+9Cd/E t+GgFyLVHx+om4YW2o2wLxnuh46jUyS1LBxxf1eD91aof1TFRwrMIbViO18Bde1UQ3IepKFS U2tX869Ug02G9k+yttLb0dmEMi5lTjC2TCrCvkbje+lHpsxp5nR1XHqO8dwzT7i07MgiEVuF tBLMnengbI58gz7CIvAkkHfnKGvI/dPlBXR/XuOmDLd9HpTVxR9BP2ttRE3Y0LXqY68/UbeV /q0Dr9hNAJdyMmEI68Ma9vzjFwASu2wcM/GbTeXnGG9TQ2N2qvKdJDjLn4Z2z/cD1RClgQ78 nOPNAx4DSCk8CrFFDI7LVv0eAv39PVm7nayT0s61QaPOlZl0KC0/QJTgPi0RPYa37ZCsyAk+ H1vBFjo+dXQBpKbohZ5OqVRZdRo+FBcyWfQrBBwJLSlJqFmw0EbKkF54x6o2BJwBYFN18Mtq RvG1SJULqSVmBNEfjKch9XrP6HPb3P191apYrLX3VfX1JCX/L0O4bI2sQerugbhDUck/3h9t rsdm3KB+pXHChYTWpPtQw428RZ9vbTTfig64cvdy3RtNaC+tjKK1cguAaMpzROpftEXN63hd ke6CMkdHcmpNKovn3CmaxsFOKZZ86t1d8Kqev2a2bK6af56lWHuhmBG7YZhl0OUonAkG6iYg tBcmaHeglbUMlW0xE2suc32h41eMDQbH27ljDPhGJYUfKp5O4ACFWapJcSzgNR4nZ/kHXBCp zvBTxsL3tGkfR2KYhnzxwpVgA4MpXu9mCykiTlwuz4sp6ubmifJxq6xEXhPcn4OX2RkgVr2d MKvi9cAXU+0KQ0tvBSg7Ef+gaNcoe4sSgubCVcNdC/wIWZ4V6K2vbfXeM9D5qQjtiBPWfi9a 1SXIlLkiyMTyDirX25XxTRhMiqvpo28hRtxzmSUMHd0qnPdP8B23xbWotLGF7Zd2T8PRS8wj jeyZBD0JN2u59STi9HGtsixUmugUttYdiyjwY6btSS97HFnGlXlx6H1yoChSFJql3Omn9BxH T3FthP9fpXm28HYeap8c09kCUW9o8t2F4dik5ch0ZQZ2HwUnJKQrjIMlWb+N8ke2LqrNiJcA 25WhYePsE65hR4GTDrB3Y/yW3SDz9E0YtC7ZjlTwSch94VQD6zS6rVYnCxzq17+rATLYPE7k C1OrJlmoHMcnewNvxIgiyuHBbVHV1JcOjbtkAjO6tSWo6BeZWLper+1nhkb/5jpHPSZrwdQV WysMIwmGjVx6dk5NVbk33j664Wic97VJ4F21FXcg1LLiO5bL4g0n/wBiH98OG7zinYizvYyk R1k2Zzp9JjCMWhm+7i1RwJJLjCgLd1G4Snj1OwN+6Tel5DqBJhqHS8HGYflXe79WixHruzpb k6PCGFu8SrdQOuHW1XDtwE+6CiTW5GzayPJeD9Dlo4kHUfFYhQY2VFxPn1yn4ZlRF70gpW5K gEhoGhWvAawqwMQmLw2cUOjAyGP/EHwLW1sAJmHcEgJtEcbuwGMYJbYtqUqT0Q6ttWgtFDfd TDdPl4VSzlPAgvdWRjiJuX8vIGQtbHHWazuaaOJO+zGqPQCBa7XntT2g9cgp3DUcZzRWxsqR /wjhhgZBSE/S5mfwm9fDXRQzn2FbtbH9k3lpGsq/obmoa6tAEW2tMOOE+cAa4wxvUrt0OHYb ajIw38oTFQQnoUFwXuCoFQG9HgVjSwmNzykELBa8DXIULqVgKhcSRgSdyJ0MsJMqaM6xAhEf 8DB2Jvz0fZjg/g5Bk0gNxSpk9y1ZcEMP2C2NU/WTEeNOrOcIDTXwsbxKaqiQLxUhe9QulW+o zGeW0PkOz2CkXHuWXXNealUizqHORVFpIynWhNkCGymU9e/LxPnbZl4ijo5xbByjXTPdCYdP TV6b0JRv+iQ4Cdf0ZAdUyRK6ntoK/XBmj7MtbGJbMZL96IxUmIoyrE/gjxy0bZe4SBaSeYgn SLTqoQruFS6iqyVzTEhVhNSqzFNjYbNvEN4OKyf+IMTPBSMtB8L82iUDAwH4tV/Ddi68btUz cLOk761LTNq/Nfd/M9aDM/RYpHiUjJpIV/yFTjYARFQByasLn3ajldBneu683SUqt0lr8Gpl sZTDLBcU1MxG7URDUEvT7lgaN9nGzgjl7CclssB43Gz+QLQSMtttZfCTvuOAP/rJV5xaJFLY hIJxfXzKoFBbuUTOmRnY1h+2ZvPQg/eAY0LrSpmYQs55k5K9SomJoXc803gYwKppnQUEKzt9 iM=
- Ironport-sdr: MwHJxF8uYFmrPPCIJg9Ljd3U73dtGpEfRf0Q7H7CaGFcghcuxuEWMSOSEqL7xEckS+xhVk92Kf BVTM/sMm4hqBzNnd131v5WfUr4WUFGGlWWfl7QPxd8oVbh4Kdv16h6zi7IDURs2sqyFl8N9F9o CIQAMNnf2voSlfR9qNz8GKu2if1bOKdU5JHW44CBp0u3UTKJ9b60U+rqtf3xT9i5IQBkKX9Qqv 2SSIclf7lH7pDdB5WbEfmMeFtzyYv2mgmlYRbEj3vXsWJZCTJsmrfBFJUBjVfiBBnm/nDXQe8M RNdchRrnPEOCxLs5cTBgvsMh
The 3D triangulation of the points of a surface is a partition of the
convex-hull of the points into tetrahedra. It is not a 3D constrained
triangulation restricted to the inside of the volume bounded by your
input mesh. Such a triangulation is not yet available in CGAL.
Best,
Sebastien.
On 7/9/22 10:26, Bob Bill ( via cgal-discuss Mailing List) wrote:
Okay, I discovered what is the problem.
One of the faces of my original surface_mesh is *not* straight. So, when I construct the triangulation giving only the set of points, this feature is not present in the resulting triangulation. See the attached image: in blue I have the original (triangulated) surface mesh, while white line is coming from the Triangulation3 I created using the vertices of the surface only. So, the question now is:
How can I get a triangulation that is truly representing also the boundary of the surface mesh? I think this should be possible, but I've not been able to find an example about how to do this. This should be possible. Also, my surface_mesh is a triangulated surface mesh, as it's the result of a BooleanOperation.
Best,
Bob
Il venerdì 8 luglio 2022, 13:37:52 CEST, Bob Bill" ( via cgal-discuss Mailing List) <> ha scritto:
Dear all,
I notice that if I compute the volume bounded by a Surface_mesh with CGAL::Polygon_mesh_processing::volume(surface_mesh)
I have a result that is different w.r.t what I get if I first mesh the inside and then sum the volume of each tetrahedron, i.e.
using K = CGAL::Exact_predicates_inexact_constructions_kernel;
using Triangulation_3 = CGAL::Triangulation_3<K>;
double sum =0.;
Triangulation_3 tria;
tria.insert(sm.points().begin(), sm.points().end());
for (const auto &c : tr_q.finite_cell_handles())
{
const auto &tet = tr_q.tetrahedron(c);
sum += std::abs(tet.volume());
}
What is the reason for that? I spent a whole day but I really cannot figure it out.
Best, BoB
--
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
- [cgal-discuss] Issue with computing volume of a surface mesh, Bob Bill, 07/08/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Bob Bill, 07/09/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Sebastien Loriot, 07/11/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Bob Bill, 07/11/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Sebastien Loriot, 07/11/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Bob Bill, 07/11/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Sebastien Loriot, 07/11/2022
- Re: [cgal-discuss] Issue with computing volume of a surface mesh, Bob Bill, 07/09/2022
Archive powered by MHonArc 2.6.19+.