Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] Alpha shape 3D : Are Exterior , interior, singular, regular exclusive

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] Alpha shape 3D : Are Exterior , interior, singular, regular exclusive


Chronological Thread 
  • From: <>
  • To:
  • Subject: Re: [cgal-discuss] Alpha shape 3D : Are Exterior , interior, singular, regular exclusive
  • Date: Sun, 09 Jul 2023 20:41:36 +0200
  • Authentication-results: mail2-smtp-roc.national.inria.fr; spf=None ; spf=None ; spf=Pass
  • Ironport-data: A9a23:iiTzSa6kyq1OCvLA8x0n3AxRtBXDchMFZxGqfqrLsTDasY5as4F+v mZMUDzUM/jZamH1fdF/b9uw8ksHsZKGzdYyQQFqrS5kZn8b8sCt6faxfh6hZXvKRiHgZBs6t JtGMoGowOQcFCK0SsKFa+C5xZVE/fjUAOC6UoYoAwgpLSd8UiAtlBl/rOAwh49skLCRDhiE0 T/Ii5S31GSNhXgsawr414rZ8Ek05Kuq6WtC1rADTakjUGH2xyF94K03fvnZw0vQGuF8AuO8T uDf+7C1lkuxE8AFV7tJOp6iGqE7aua60Tqm0hK6aID+6vR2nRHe545gXBYqhei7vB3S9zx54 I0lWZVd0m7FNIWU8AgWe0Ew/y2TocSqUVIISJSymZX78qHIT5fj68lAB1oaH64yw/RMD31e6 PgCciAPPw/W0opawJrjIgVtrt4uKMD6YsYT/HRpzDWfAv8gTZGFRajWjTNa9G1p2YYRRbCHN 5FfMGs/BPjDS0Un1lM/DZ87luqyrnLyaCZRuRSSv6VfD237lVUvj+e9b4eNEjCMbfh8o0S7m UDKxjnoAQsnMtDO+zTYzW3504cjmguhCdpNSODgnhJwu3WYyWUXTREXTlCmuuKRkV+7Q9sZK koO+yNoo7JayaCwZsPwWxSp/juI+BsVWt4WHOQ85ADLzKfIi+qEOoQaZgJvY8cKluwcfC4R2 EOZgtfzXyxt4aLAHBpx6YyohT+1PCEUK0oLaikFURYJ7rHfTGcb002nojFLTf/dszHlJd3j6 2za/HVu3d3/meZRjP7rojgrlhrx/sChc+Ij2unAdkyfhj6Viaa/aoih+Aae4bBFJYedCFaIu nQF3c6T8Iji7K1hdgTSEI3h/5nzvZ5p1QEwZ3YzQ/HNEBzxpBaekXh4um0WGauQGp9slcXVS EHSoxhNw5RYIWGna6R6C6roVZR6nfS5TY+4CqqEBjarXnSXXFLXlM2JTRDJt10BbGB2+U3CE cvCIJbwZZrkIf08lWHeqxghPU8DnHFmmj+OHfgXPjyo2LySZTaOQKxtDbd9Rr5R0U9wmy2Mq 4w3H5LTk313CbSiCgGJq9J7BQ5RdhATW8umw/G7g8bfeWKK7kl6U6CJqV7gEqQ595loehDgp yDiBhYIlwOl1BUq62yiMxheVV8mZr4nxVpTAMDmFQ/AN6ELMdbxvpQMPYA6Z6cm/+FFxPt5B atNMcaZD/gFDnyN9z0BZNOv5MZvZTa6tzKoZiCFWTkYe4I/Zgrr/tS/QBDj2hNTBQWKtOw/g Yaa6CXlfbQ5ST9PMuPqedO07lbovXEiiON4BETJBd9IeXTTyothKg2vr/kVJ8sjB0jlx2KZ8 wPPAR4nuvT/+aEowvbKmKqBkd+IEvR/LGVeDWL0/ba7DgiE32uBkKtrcveEQiDZb0zwoJ6dX ORyy+ruFdE+h3NYmtNYP5cy6oxm/PrphbtR7jo8LUXxd17xV49ReCiX7/dApohm5+F/uzLvf mmt59MDG7GCGP29IW4rPAB/M9iyj6AFqALzs8YwDl7xvhJs3by9Vk5XARmApQpdIJZxM6In2 e0Rg9EX2SPulisVNsu6sQ4M+1SuNnAgV4AVhqMeCqLviSspzQhMWoyDKynU5JrUVc5AHHN3K RCphY3DpY9m+Gz8T1QJG0Pg59Ft3aY1hEgSzXspBUi4pd7et/pmgDxT6WsWSypW/DVm0sVyG GpZB2tnFJWN4j0zps5iXn+tKS5FFha26kz88HpXtWz7Hm2DdH3BE305AsmJpHsmymN7eiNK2 rOyxELOcyfYTOuo0gQcAUda+uHeF/pv/Qj8qeWbNsWiHahiRwH6g6WrNFE6mzG+Dewf3ET49 PRXpsBuYqjGNAkVka0xK6+e8Z8yEBmkBmhzcctNzZMzP1P3WW+NgGCVCkWLZMlyCeTA8ha4B +xQN8t/bUmC+xjUnA8LJ5wnAuFSrKYy6csga4HbAzcMk4GiowpDtLPS8SnDh1EXfeh+rPZlF KTvc2OtL2/Bo1pVhG7Hk+dcMEWaf9QvRVPxzcK1wsozBrMBt+BdK38g2eGa4lGQYRpV7jOPn Qb5f6SN5fdT+YdtuIrNE6t4GAS/L+3oZtmI6Ayet9dvb8vFFMX/qDMutVjsOjpJMYsrW9hYk aqHtPj11hjnuIkafn/4mZ7bMYV0/uS3AfRqN/zoIElgnSetXNHm5z0B8TuaLb1LiNZs2dm1d TCnacefdc8nZPkF/Sd7MxNhKhc6D7j7SozCpinn9vSFNUU74Dz9ddii8SfkUHFfeio2IKbBM w7Tucj/1uAA+c4ITFUBCup9CpB1HE77VOF0P5ftvD2fFS+zjknEprLmkgE65CrWDmWfVvz3+ o/BWgO0YSHaVHskFz2Fm9caUtwr4HdBbS0YZE8Z8s8ozjz8CWcHKaISOJMKC9dSnzCaOFTQe mTWdGV7YcnidW0sTPk+yI2LssSj6igmMdriPTU4uUSOZE9awW9G7KRJrk9dDrQfRtcn5OqqM soX6zv3JBfZLlSFgwoMzqTTvNqLDc82Cp7FFY4RXiAy79sj7W026UFc
  • Ironport-hdrordr: A9a23:lyqNo6NZW0ErvMBcTsOjsMiBIKoaSvp037BL7S9MoHNuE/Bw+/ rCoB1f73PJYVUqKRIdcLK7Scu9qBrnnPYe3WB7B8bGYOCJgguVxc1ZnO/fK0yKIUPD398Y+4 BMN5J3EtiYNzNHsfo=
  • Ironport-phdr: A9a23:K/9E+BBZL26z0bngtVj+UyQUQUkY04WdBeb1wqQuh78GSKm/5ZOqZ BWZua8zygWUFtyBu7oE07OQ7/q9HzRYoN6oizMrTt9lb1c9k8IYnggtUoauKHbQC7rUVRE8B 9lIT1R//nu2YgB/Ecf6YEDO8DXptWZBUhrwOhBoKevrB4Xck9q41/yo+53Ufg5EmCexbal9I RmorAjdrNcajIttJqos1BfEoWZDd+BKyW91P16ekAjw68eq8JJ/7yhcvu8q+tJdX6n9Y6k3Q rtUASg8PWso/sPrrx7DTQWO5nsYTGoblwdDDhbG4h/nQJr/qzP2ueVh1iaUO832Vq00Vi+57 6h3Uh/oiTwIOCA//WrKl8F/lqNboBampxxi347ZZZyeOfRicq/Be94RWGxMVdtTWSNcGIOxd 4sBAOgdMuhWsYfyqFgAowawCgS3Huzj1iVFi2Xq0aAg0eksFxzN0gw6H9IJtXTZtND1O7kXU e+u0qbI1y7OYO9I1jfm84jHaBAgreuLXbJqfsra1E4iFwbbgVWVsoPqISmV2/8RvGeF9OpgT vigim49pwF2uDivyd4hh4/UjY0a1l7K7z92wJopJdKmUk57Z8apHpteuSyVM4Z6XN0uTn9nt ig0ybALtpq2cigUxZg72xLTd/+Kf5aH7B/+W+icPyt0iWxqdr+jiBi/8EmtxO3iWsSyzV1Eo C1FktzWuXAM0Rze8tSISvRn8UelwzqP0BrT6u5cLUA1k6rbMZ8hzaQqmpodq0TOGDL9lkbuj KKOaEko5+ql5/7lb7n6vJOQKox5hw/kPqgznsGyAuI1ORUUUWeB4+Szzrjj8FX5QLpUiv02l bHUsIrEKsQfvKG4Aw5Y3p8+5BqlCjepytUYnX0dIFJKYh2HiZXmO1DULPD+E/i/n0ygkDZtx //YIr3sGovBImXBnbv7f7tx91RQxBcwwN1R/Z5YF7UMLfHrVk/0rtPYDxs5MwKuw+bgDdVwz oYeVnyKAqCCLqzSq0SF6f4xLOmXf48aojb8J+Ii5/70gn81gVkdcrKx0psNdH+4BuhmI1meY Xf0n9sNCX0KsRYmTOz2lF2CViZeaGq9X60m4jE3EZ+pDYbYRo+xnbyBxzy7E4ZNZmFGD1CMC W3ne5+FW/cKciKSI9VuniYKVbi7GMcc0wqzvlr616Z/NbiTvTYJsIrqktlz/eza0x8osidlC tyUlGCLQWYzlWwBQ3o63bt0vFdmmWqFyrVyo+BdEYlT++9RSVV9coXNyvRzTdH0QAPIONmTD 023R82vRjA3QNV2yNAHZwNxGs6pkwvYjBesGKIfq7GbGMk07r7ExCq2YN1szm7PkqgnlVgvB MVVcna3g7Z2sAnVCYmOmEqQk+OmdL8XwTXWp1uEmGGBtUUdXA9rWrjeRlgeYFHXpJL3/BDsV birXLEqNwFM0+aJJ7BRY8yvik9JDKPmNdHUbniZkGqrHh+ZgLSWY9y5KC0mwCzBBR1cwEgo9 nGcOF1mbs/Ai2fXDTg0UEnqf1up6+517nWyUk4zyQiOKUxnzbu8vBAP1rSHU/1G+LUCtW86r ilsWk6n1ofTBt+KqhBJd6hHetosplBd2jGRrBRza7qnKa0qnVsCa0JytkLq2Q9wD9BFnM0go WkCxwdpMquJllRbeGDQxoj+b4XeMXK65xWzc+jW11XZhc6R4bsK4e8kpk/LpgytH1t8tnkh1 tBU1z2T75PGDUwUXI6Zvl8f0R98qvmaZyA849mRznhwKeyutTSE3ds1BewjwxLmftFFMarCG hWgW8sdT9OjLuAngT3LJloNIfxS+agoPsina+rO2ailO/xllS6nimIP6Z5000aF/S5xAuDS2 JNNz/ad1wqBHzDy6TXp+sX6nIpJfxkZE3GkwDmiA5ReJ+VzcYsNFWayMpiv3Nwtz5XpWnNe6 BuiHwZcgZDvJkLUNgymm1QBhRdywzTvgyazwj1qnit8q6Oe2HaL2OH+bF8dPWUNQmB+jFDqK IzyjtYAXUHuYRJ68XntrUv82aVfo7xyamfJRkIdNS37KmtvSIO7u6GZat8J74kn+3YyMqz0c RWBR7jxrgFPmSbuGmFT3hgweiqxs4S/mAZ1wjHVPDN4q3zXftt1zBHU6YnHRPJf6TEBQTFxl TjdAlXU08CB29KSmt+Dt+m/UzjkTZhPaWzxyojGsiKn5GpsCBn5nvapm9ShHxJomSn8n8JnU ynFtnOeKsHiyri6POR7f0JpGE606sx0HZt7m5cxg5dY0GYTh5Gc934K2WnpNtAT1aX7ZXsLD TkFprydqA3o3UduNFqCypniX2TbxdFuJpG7bm4QxiMh/pVPAaaQ4qZDmHgQwBLwpgbQbP5h2 zYFnKF3tzhF26dQ6FFrl3TNUdVwVQFCMCfhlgqF9YW7paRTPiO0dKSokVB5lpanBa2DpQdVX DD4fI0jFGl+9JYaUhqE3Xvt54XjYNSVY8gUs0jelhrFheVPAJk8jeEDlGxgI2237hhHg6Yry Adj2523pt3NLmho9q+jKhVVLSb8doUU4D6n3u5O282R2Y6oBJBoHD4GCYDpQfyfGzUXrf37N gyKHVXQs1+jEKHEVU+a4UZi9TfUFoyzcmuQLz8fxMljQx+UIApehhoVVXM0hMxxGgevzc3nO EB3g1JZrkb/sQdJw/l0OgPXSGrbrRb2LztyTZGeKFxZ5wdO5gHTPNDW4u9oHi5e94Gst0TXc zPdOFwOVDFSHBXVVzWBdvGn/pHY/vKdB/aiIveGeriIpeFEFr+JyZ+pzop67mOMO8GIbTFpC /w23FYGXGgsQZ2G3WpXDXVHzGSUPp7IwXX0sjd6pc2+7vnxDQfm5I/VTqBXLc0q4Re9x6GKK ++Xgi99bzdezJIFg3HSm919lBYfjT9jcz61HPEOrynIGejVl6pdAgQzYC5oL8JZqaUm0UMeX KyTwsOwzbN+gvMvXh1dUkf9n8iyecEQC3uwMFrWWgOHcrGPJDmNzMjxbaL6T7BMxrYx1VX4q XOQFEntOS6Gnj/iWkW0MO1CuyqcOQRXpIC3dhs+QXimVt/tbQe3dcNmlTBji6NhnWvEbCRPV Fo0O1MItLCb6jlUx+lyC3AUpGQwNvGKwm6Qp6zRLpJc2ROOKi95jPlT8TI90bQHtEmsq9R+l TDJp88opEupwLHnIttPShdKoycVwY7NuExjPePW/59MWDDC8Q5ftA2t
  • Ironport-sdr: 64aaff60_oUOEQABc5tYAORYHSEYFKbj7O263r9h/UZERPuGH1wEmgK1 k4kfimEa95+mb4ofeMvElb6MtJpHndn9lWfLXGQ==

Sebastien,
Thanks for your response!
There is no problem on CGAL side. 
I found the problem. The number of singular faces in this case was signficanly small than reported. 

On Friday, July 07, 2023 00:05 PDT, "Sebastien Loriot" ( via cgal-discuss Mailing List) <> wrote:
 
">It should indeed be partition as the classify() function for each facet
gives a unique result.

Maybe if you want a complete example I can compile and run, I can have a
look and tell you if it is a bug in CGAL or in your code.

Best,

Sebastien.

On 7/7/23 02:43, wrote:
> It is not clear to me of a facet belongs exclusively to one of the four
> types.
> I have
> read https://doc.cgal.org/latest/Alpha_shapes_3/index.html#Alpha_shapes_3Definitions and section 2.3 of it. From that it appears it should be exclusive.
> However the sum of number of diffrent types is not eqaul to the number
> when no condtion is used. I modified the code below and executed it 5
> times, four times for each type and then once without the "if" condition
> so that all facets are counted.
> ​​​​​​ Do I need to use infinite precesion for exact match?
> Component 1 , general mode
> Exterior 7208 ALL 67754
> Interior 56668
> 63876 diff 3878
> Regular 3872
> Singular 5450
> 9322
>
>
>
>   double alpha = 0;
>     Alpha_shape_3 as(points.begin(), points.end(), alpha, mode);
>     Alpha_iterator opt = as.find_optimal_alpha(num_of_components);
>     as.set_alpha(*opt);
> *    Classification_type in_class_type = Classification_type::SINGULAR;*
>     for (auto it = as.finite_facets_begin();
>         it != as.finite_facets_end();
>         it++)
>     {
>         auto class_type  = as.classify(*it);
> *        if (in_class_type == class_type)*
>         {
>
> --
> 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

 



 


Archive powered by MHonArc 2.6.19+.

Top of Page