Skip to Content.
Sympa Menu

cgal-discuss - Re: [cgal-discuss] On cell deletion, do neighboring cells need neighbors to be specified?

Subject: CGAL users discussion list

List archive

Re: [cgal-discuss] On cell deletion, do neighboring cells need neighbors to be specified?


Chronological Thread 
  • From: Sebastien Loriot <>
  • To:
  • Subject: Re: [cgal-discuss] On cell deletion, do neighboring cells need neighbors to be specified?
  • Date: Tue, 5 Jul 2022 10:37:35 +0200
  • Authentication-results: mail3-smtp-sop.national.inria.fr; spf=None ; spf=Pass ; spf=None
  • Ironport-data: A9a23:3Q8CTaK/AMlMf9N+FE+R4pMlxSXFcZb7ZxGr2PjKsXjdYENShmAGn WFOXmuPM/aMa2DxLY9/Oo2/9hsCucfTn9EyGQYd+CA2RRqmi+KVXIXDdh+Y0wC6d5CYEho/t 63yTvGacajYm1eF/k/F3oDJ9CU6jefSLlbFILas1hpZHGeIcw98z0M58wIFqtQw24LhXFvQ4 YqaT/D3YTdJ5RYkagr41IrY8HuDjNyq0N/PlgFWiVhj5TcyplFNZH4tDfnZw0jQHuG4KtWHq 9Prl9lVyI92EyAFUbtJmp6jGqEDryW70QKm0hK6UID66vROS7BbPqsTbJIhhUlrZzqhwf9M6 YpmqKeLGV14GIj8mqcBUQtFKnQrVUFG0OevzXmXtMWSywjZdiKpzag2Sk4xOoIc96B8BmQmG f4wcmhcKEDewbvskPTmEoGAhex7RCXvFJgbvWthyi2fCPIOTpXKQqGM7thdtNs1rpkeQKiHO JdJAdZpRCTdcQceIFVHM74ns+ysul/HIzJEjl3A8MLb5ECKlFAruFT3C/LeddWOAMlUhU2Fv Xnu5HX8GhhcNdqFyDPD/GjEuwPUtSbyWYZXC7jhs/A32BucwWscDBBQXly+yRWktqKgc4xwF m1L0TcVlvhopWqOftu+Qx+qvUfR63bwROFsO+E97QiMzI/d7ACYGnUIQ1Z9hDoO5J9eqdsCh g/hoj/5OdB8mObKFi/Fp994uRv3aHdFdzZTDcMRZVJdu4GLnW0lsv7Yoj9e/EOdi9T0HXTvx mnPonVmwboUisEP2uOw+lWvb9OQSnrhHl5dCub/BDrNAuZFiGiNOdXABb/zs68oEWphZgPd1 EXoYuDHhAz0MbmDlTaWXMIGF6yz6vCOPVX02AAyQch8pmj2oCX8LOi8BQ2Swm85Yq7onhe5M CfuVf95ufe/wVPxMPArO9ntYyjU5fG5TIi/PhwrUja+SsEpKFXvENBGakmX0GTg+HXAYoluU ap3hf2EVC5AYYw+lGTeb75EjdcDm35jrUuOGsiT50n2idK2OS/OIZ9YYQvmRr5jsMus/l6Om /4BbJvi9vmqeLehCsUh2dVDcw5iwLlSLcyelvG7gcbZfVA4Rzl/VaC5LHFIU9UNopm5X9zgp hmVMnK0AnKm7ZEeAQnVOH1ldp31WpNz8SAyMSA2bASn3nEiZcCk66JGL8k7erwu9epCy/9oT qldK5/QXKgXEjmXqS4AaZTdrZB5cEv5iA+LOR2jamdtcpNlQTvP5dK5LBDk8zMDD3bsuMZn+ ++g2wrXTIAtXQNnCMqKOvujw0ng73cYkeN2GUDPJ4ALKknr9YFrLQ33j+M2c5lcc0WdmmPC2 l/PUxkCpOTLr4sky/XzhPiJ/9WzDu9zPktGBG2Ev7u7MC/t+GD8k4JNVeC/ey+ECDH59aCkU uViz//mNcoBklsX4ZF3FKxmzP5n6tbi++1awwBjECmZZlinEOk7cHyP3M0KqagUg7EE4U25X UWA/tQcMrKMYZu3HFkULQsjT+KCyfBExWWIvKpteB33tH1t4b6KcUROJB3Q2iZTG714bdE+y uA7tc9KtgGy1kgwPtCdgnwG/miANCZbAaAut5VfH4qyzwR2khdNZpvTDiKw65aKMo0ePk4vK z6SpazDm7UMmRaYIiRrTSDAjbhHmJADmBFW11tedV6HrdzI260s1xpL/DVrEwlYw32rCQ6o1 rSH6qG0GUmPw9usrM1KXmTpBAQYQRPEogr+zFwGkGCfRE6tPoAIwKvRJs7VlH31MUoFFtSYw F1c4GngWDfuOsr220PenGZ7/uf7Q4UZGhLqwaia8gfsI3X+STXgi66qI2EPrnMLxC/3aFLv/ YFXwQq7VUE32eP8bUH250l2GIn8kCy5GVE=
  • Ironport-hdrordr: A9a23:oi9lTaygKP8O3d/ZxA2QKrPwD71zdoMgy1knxilNoG9uA6qlfq eV7YgmPH7P+UsssRQb8+xoV5PwI080maQFmrX5eI3SJjUO21HYSb2Kj7GSoAEIcheWnoU86U 4jSdkHNDSZNzlHZK3BkW6F+rgbsaC6GeyT9IPjJrRWIT2CqZsM0+60MGmm+4RNKjV7OQ==
  • Ironport-phdr: A9a23:o7dgbh2FYpCNBX7fsmDO+Q0yDhhOgF0UFjAc5pdvsb9SaKPrp82kY BaEo6wz1RSRA83y0LFts6LuqafuWGgNs96qkUspV9hybSIDktgchAc6AcSIWgXRJf/uaDEmT owZDAc2t360PlJIF8ngelbcvmO97SIIGhX4KAF5Ovn5FpTdgsip2e2+4YDfbgRViDayfb9/L Qi9oBnMuMURnYZsMLs6xAHTontPdeRWxGdoKkyWkh3h+Mq+/4Nt/jpJtf45+MFOTav1f6IjT bxFFzsmKHw65NfqtRbYUwSC4GYXX3gMnRpJBwjF6wz6Xov0vyDnuOdxxDWWMMvrRr0yRD+s7 bpkSAXwhSkEOTA2/2/YhNFugqJVoByupQBwzYHPbYGJN/dzZL/Rcc8USGdDWMtaSixPApm7b 4sKF+cPIPhXoJLnqFsVsRu+BBejBOXywTFOgX/2x6w60+M9EQrb2wEgBNYOsHPOrNX0M6cST fq5zKzMzTjYYPNW3C3y6InMchw7vf6MWrdwfNPXxEIyGAzLkk+eppb5PzOJyOsNqW6b4vJ9W e6zhWMppQF/riWty8oojoTEm4EYx1DY+SlnwIg5ON61RkFnbNOkH5ZdtT2XOYR4T84jTGxmu Ck0x7MFtJO1eiUB1Zopxxnaa/OdcoiI5AruVOmQITdkhHJlZamwiwyu/kinz+3wTtO730tQo SpAndnDqGwN2ADT6seZTvt9+l2t2TGV1wDc8u1EIEY0lbDaK5E72LIwmIATvELeFSH1gEX7l LGaelkg9+Sy6OnqYq/qqoKCO4J3kA3zMrgiltS8DOk5KAQDXWmW9f6i2LDt40H1WqhGg/05n 6TfrZvUP94UprSjDA9Qyosj6wiwDzOh0NkAmHkINlNFeBaeg4nnPlHCPen0Dfm/jliwijtry PfGPrruApXJMHfPiqvufbF460JEyQozy85Q545MB70fPP7+XlX9ud/YAxMjLgC43eXqBM992 48CQW6PB7WWMKLWsV+G/OIvJOyMaZcOuDb5LPgl/f/ujXglmVMHeKmp2YAaaH+9Hvt8IkWZZ WDgjcsGEWcPpgY+Svflh0WYUTFPf3ayQ7485jYjBY67FYvDXJyigKSd3CenGZ1bfnxJCl+WH nfsbomLRvYMaDmOLc98iTwETqOsS5Qh1BGrrA/10aBrLuvS+i0CtJLszsJ55+PJlUJ6yDpvE s7I032RV3om2SQTVjou1eZ+p1Z8wxGNy+9jkvlAHJtS4f1OFQw1PJqZw+1hAM3pQVH8eMyUQ mqrUsnzAS0tVsljhJgVckNlEpOjiArC1mykGfgOhrmTDds19KzbmHP+LsI4x3fd37Q6lAoaR Z5EOmSiw6J+7APOHJXhkkODlq/se75P8jTK8TK4wGCHp11ZXQg4daLfXHcDLh/NqdPj50TeC bqqIbsiOwpFj8WFL/0ZOZXSkVxaSaK7a5zlaGWrljLobf7p7raFbY6xPn4YwD2YE08c1QYa4 XeBMwE6QCanuWPXSjJ0Rhr0e0251+54pTugS1MsiRmQZhh6072v+xkJw/mYY/wW17MA/iwmr mY8B06ziurfEMHIvA99ZONZaNI57k1A0DfCswtnP5u8aaVmrlEbegVz+Ujp0kY/EZ1OxO4tq n5i1w9uMeSY3VdGIiufxoz1M6bLJ3Pa+Rmub+vJ3wib3ovOvKgI7/s8phPouwTB+lMK1XJh3 pEV1nKd4s6PFw8OSdfqVV5x8RFmprbcay175oXO1HQqP7Pm+jnFk8kkAuco0HPCN59WLb+EG QnuEsYbG9nmKeoknEKsZw4FO+Ya/bA9Psevff+LkKCxO+MokDWjhGVBqIdztyDEvzF4TfTJ2 IpDxvWw0Q6OVjO6h1Ck84j2lY1CeTAODz+n0yG3TIVVZ6B0YcMKETL0e5zxlog43sS9HSMAp zvBTxsc1cSkeASfdQn41AxUjgEMpGC/3DC/13pymi0oqayW2GrPxf7jfVwJID0uJiEqgFHyL Iyzl91fUlKval1jjx+i/0f92+5eoIxwKmDSRQFDeC28fAQAGuOg86GPZcJC8sZiqiFQSuW7f RafTpbyphIb12XoGG4Ul1VZP3m6/57+mRJ9kmeUKn1++WHYdc9HzhDa/NXARPRV01LqXQFAg CLMThi5NtitppCPkovb9/u5TySnX4FSdi/iycWBsjG67CtkG0/3k/e2k9zhWQ80tE2zn8JuU j/Jqwq6Z43D2KGzMOYhdU5tTFPx8Mt1HIhin5B43slBnyhHwMzPrTxayS/6Kp1D1Lj7bWYRS DJuoZad+wXj1EB5bzqIy4//SnSB055kbti+bHkR33F15MRLBaGIqb1cyHEt8xzo8ESIO6g7w 29OrJlmoGQXiOwIpgc3myCUA7RJWFJdITSpjBOQqda3sKRQYm+rN7m2zkt32967X9Tg6klRX mj0fpA6EGp+9MJ6ZRjX1Hrp64b4PtzURd0WvxyQ1RzHiqIGTfB53upPnidhNW/n6Dc+z+khj Btym5S+lIeCImRpuqm+B1QLU1+9L9NW8Tbrg6FEm8+Q1I36BZRtFAIAW57wRO6pGjYf5rz3c hyDGzompjKHCKLSSEWBvVx+oSuFQPXJfzmHYWMUxtJ4SFyBKVxD1UoKCS4ikMdxFxj2lpe8N h4ovnZLugG+8lwWlqppL0WtDDuZ/lzzLG5qEN7Ha0MHi2MKr0bNbZ7AsKQqR3septv561bVY m2DO1YWUydTBh3CVwilZv70vZHB67TKWbD4dqeIOOTU77QZDqftp9rn05M6rWnQcJzVYz86S aV8gxQLXGglSZ2By3NWFHNRx2SVKJTC7Bakpn8u8Zv5qaW3HlqpvczWVd4weZ1u40zk2//Sc b7Nwn8jeXAAkcpTjX7QlOpFhQBU1nEoLmj3V+xH7H+FTbqMyPUOUVhBMHI1b5EOt+Vlj2wvc YbNg9fxnNaUl9YTDFFIHRzkk8CtPokRJn2lcUjAHACNPaiHIjvCx4f2Z7m9QPtelrccsRr4o judH0L5W1bL3zD0SxCiN/1NhyCHLVRfvo+6aBNkFWnkSprvdBS6NNZ9iTB+z6czgzvGMmsVM D40dE0ozPXY9SRDnvB2AHBM9FJgJOiA3jmctqzWd85QvvxsDSB50eld5TVyyrdY6j1FWO0gm CbWqY0LwRnume2Oxzx7FRtW/2wT1cTb4AM7Y/Wfq8ETPBSMtAgA5miRFRkQ8t5sC9m0/rtV1 sCKj6XrbjFL79PT+8IYQcnSMsOOdnQ7Yn+LUHbZChUISTmzOCTRnUtYxbuJ8nqPr58m7J3os JULQ75fElcyE7lJbyYtVMxHO5pxUj4+xPSDi9UU4HOlsBTLbMBTv5SCSffLRPuycHCWir5LY xZOyrT9Z9d2VMWzywlpbV90m57PEkzbUIVWoyFvWQQzpV1E7Hl0SmBbM6PNbwqs5DoMEKfxk ENuzAR5Zusp+XHn5FJlfjIiSwM/lUAwnZPuhjXDKVYZyY+/WIhXD2z/sE1javvG
  • Ironport-sdr: PQT5aaGK1jN0CLVMVbsvf1qqL50cpGn/J7nvYHFkiRqPiWFmRJMGRdhHEOikJ9JeYipOvMh20c SVJfeV6yBwjm8R6J8GNpMUiFOm50f/d8KjjP3LP9ufAQx2eDQcvra6W3Hu/W6I6oM2sn8iYbK6 dEvukgXK+P1jjStQbT+e0dTYRRRt17XxlMSPVBgHmBIos9T8HxekAfA9dcU+4O+XeGezd0p7gf mdNog5ALuh+FVUxBIzX0jmW6d0ASlEIpV2PnjvgDM/RFmG21Dc81liKAsxgJEvFUsI6n4u8E0v kULdhN2FDy8yGX+Hll3s+C2M

Hello Adam,

did you already have a look at CGAL/Tetrahedral_remeshing/internal/flip_edges.h

flip_n_to_m() seems to be what you are looking for.

Jane is not around for a few weeks. I'm sure she will get back to
you when she's back.

Best,

Sebastien.

On 6/29/22 22:52, Adam Getchell ( via cgal-discuss Mailing List) wrote:
Hello all,

I’m working on a bistellar flip for 3D triangulations, briefly as follows:

Step 1: Find a candidate edge for which a circulator around that edge obtains 4 incident cells.

Step 2: Get the 2 vertices comprising the candidate edge.

Step 3: Get the remaining 4 vertices in the 4-cell complex (there are 6 vertices in the cell complex).

Step 4: Obtain a new pivot edge between 2 different vertices for which a circulator around that pivot edge would obtain 4 new incident cells comprised of the 6 vertices.

Step 5: Find the 8 neighboring cells of the original 4-cell complex.

Step 6: Delete the old cells.

Step 7: Create the new cells using the same 6 vertices.

Step 8: Assign the 8 neighboring cells to the cells in the new 4-cell complex.

Details here: https://github.com/acgetchell/bistellar-flip <https://github.com/acgetchell/bistellar-flip>

The issue is at the conclusion of Step 8, all cells have been verified as valid, cell orientations have been verified/reoriented, but the triangulation itself is not valid.

The error is: “neighbor of c has not c as neighbor”

My question is, do I also need to track the 32 neighboring cells of the 8 neighbors of the 4-cell complex being flipped? (And in my test case, the 8 neighbors are all infinite.)

Is there some other mistake?

Thanks for your advice,
Adam

P.S. I have looked at the Tetrahedral Remeshing package at Sebastien Loriot’s suggestion. I’m presently unsure if any of those functions (e.g. construct_opposite_vector_3_object()) are relevant.
--
Adam Getchell
https://adamgetchell.org <https://adamgetchell.org>




--
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