1 #ifndef INCLUDE_UNITTESTS_TRIMESH_COLLAPSE_HH
2 #define INCLUDE_UNITTESTS_TRIMESH_COLLAPSE_HH
4 #include <gtest/gtest.h>
5 #include <Unittests/unittests_common.hh>
14 virtual void SetUp() {
18 virtual void TearDown() {
41 Mesh::VertexHandle vhandle[4];
43 vhandle[0] = mesh_.add_vertex(Mesh::Point(0, 0, 0));
44 vhandle[1] = mesh_.add_vertex(Mesh::Point(0, 1, 0));
45 vhandle[2] = mesh_.add_vertex(Mesh::Point(1, 1, 0));
46 vhandle[3] = mesh_.add_vertex(Mesh::Point(1, 0, 0));
49 std::vector<Mesh::VertexHandle> face_vhandles;
51 face_vhandles.push_back(vhandle[0]);
52 face_vhandles.push_back(vhandle[1]);
53 face_vhandles.push_back(vhandle[2]);
54 mesh_.add_face(face_vhandles);
56 face_vhandles.clear();
58 face_vhandles.push_back(vhandle[0]);
59 face_vhandles.push_back(vhandle[2]);
60 face_vhandles.push_back(vhandle[3]);
61 mesh_.add_face(face_vhandles);
63 face_vhandles.clear();
65 face_vhandles.push_back(vhandle[2]);
66 face_vhandles.push_back(vhandle[1]);
67 face_vhandles.push_back(vhandle[3]);
68 mesh_.add_face(face_vhandles);
70 face_vhandles.clear();
72 face_vhandles.push_back(vhandle[3]);
73 face_vhandles.push_back(vhandle[1]);
74 face_vhandles.push_back(vhandle[0]);
75 mesh_.add_face(face_vhandles);
77 mesh_.request_vertex_status();
78 mesh_.request_edge_status();
79 mesh_.request_face_status();
81 Mesh::HalfedgeHandle v0v1 = mesh_.halfedge_handle(0);
82 Mesh::HalfedgeHandle v1v0 = mesh_.opposite_halfedge_handle(v0v1);
84 Mesh::HalfedgeHandle v1vL = mesh_.next_halfedge_handle(v0v1);
85 Mesh::HalfedgeHandle vLv1 = mesh_.opposite_halfedge_handle(v1vL);
86 Mesh::HalfedgeHandle vLv0 = mesh_.next_halfedge_handle(v1vL);
87 Mesh::HalfedgeHandle v0vL = mesh_.opposite_halfedge_handle(vLv0);
89 Mesh::HalfedgeHandle vLvR = mesh_.next_halfedge_handle(v0vL);
90 Mesh::HalfedgeHandle vRvL = mesh_.opposite_halfedge_handle(vLvR);
92 Mesh::HalfedgeHandle v0vR = mesh_.next_halfedge_handle(v1v0);
93 Mesh::HalfedgeHandle vRv0 = mesh_.opposite_halfedge_handle(v0vR);
94 Mesh::HalfedgeHandle vRv1 = mesh_.next_halfedge_handle(v0vR);
95 Mesh::HalfedgeHandle v1vR = mesh_.opposite_halfedge_handle(vRv1);
99 Mesh::VertexHandle v0 = mesh_.from_vertex_handle(v0v1);
100 Mesh::VertexHandle v1 = mesh_.to_vertex_handle(v0v1);
101 Mesh::VertexHandle vL = mesh_.to_vertex_handle(mesh_.next_halfedge_handle(v0v1));
102 Mesh::VertexHandle vR = mesh_.to_vertex_handle(mesh_.next_halfedge_handle(v1v0));
108 EXPECT_TRUE( mesh_.is_collapse_ok(v0v1) ) <<
"Collapse not ok for halfedge 0";
109 EXPECT_TRUE( mesh_.is_collapse_ok(v1v0) ) <<
"Collapse not ok for opposite of halfedge 0";
112 EXPECT_EQ(0, v0.idx() ) <<
"Index wrong for from vertex of collapse halfedge";
113 EXPECT_EQ(1, v1.idx() ) <<
"Index wrong for to vertex of collapse halfedge";
114 EXPECT_EQ(2, vL.idx() ) <<
"Index wrong for left vertex of collapse halfedge";
115 EXPECT_EQ(3, vR.idx() ) <<
"Index wrong for right vertex of collapse halfedge";
118 EXPECT_EQ(0, v0v1.idx() ) <<
"Index wrong for collapse halfedge";
119 EXPECT_EQ(1, v1v0.idx() ) <<
"Index wrong for opposite collapse halfedge";
121 EXPECT_EQ(2 , v1vL.idx() ) <<
"Index wrong for v1vL halfedge";
122 EXPECT_EQ(3 , vLv1.idx() ) <<
"Index wrong for vLv1 halfedge";
123 EXPECT_EQ(4 , vLv0.idx() ) <<
"Index wrong for vLv0 halfedge";
124 EXPECT_EQ(5 , v0vL.idx() ) <<
"Index wrong for v0vL halfedge";
126 EXPECT_EQ(6 , vLvR.idx() ) <<
"Index wrong for vLvR halfedge";
127 EXPECT_EQ(7 , vRvL.idx() ) <<
"Index wrong for vRvL halfedge";
129 EXPECT_EQ(8 , vRv0.idx() ) <<
"Index wrong for vRv0 halfedge";
130 EXPECT_EQ(9 , v0vR.idx() ) <<
"Index wrong for v0vR halfedge";
132 EXPECT_EQ(10 , v1vR.idx() ) <<
"Index wrong for v1vR halfedge";
133 EXPECT_EQ(11 , vRv1.idx() ) <<
"Index wrong for vRv1 halfedge";
139 mesh_.collapse(v0v1);
160 EXPECT_EQ(4u , mesh_.n_faces() ) <<
"Wrong number of faces (garbage collection not executed!)";
163 EXPECT_TRUE( mesh_.status(mesh_.face_handle(0)).deleted() ) <<
"Face 0 not deleted";
164 EXPECT_FALSE( mesh_.status(mesh_.face_handle(1)).deleted() ) <<
"Face 1 deleted";
165 EXPECT_FALSE( mesh_.status(mesh_.face_handle(2)).deleted() ) <<
"Face 2 deleted";
166 EXPECT_TRUE( mesh_.status(mesh_.face_handle(3)).deleted() ) <<
"Face 3 not deleted";
169 Mesh::FaceHandle fh_1 = mesh_.face_handle(1);
170 Mesh::FaceHandle fh_2 = mesh_.face_handle(2);
174 EXPECT_EQ(1 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 1 of face 1";
176 EXPECT_EQ(2 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 2 of face 1";
178 EXPECT_EQ(3 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 3 of face 1";
180 fv_it = mesh_.fv_begin(fh_2);
181 EXPECT_EQ(2 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 1 of face 2";
183 EXPECT_EQ(1 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 2 of face 2";
185 EXPECT_EQ(3 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 3 of face 2";
188 Mesh::HalfedgeHandle fh_1_he = mesh_.halfedge_handle(fh_1);
190 EXPECT_EQ(11 , fh_1_he.idx() ) <<
"Index wrong for first halfedge of face 1";
191 EXPECT_EQ(1 , mesh_.to_vertex_handle(fh_1_he).idx() ) <<
"First halfedge inside face 1 pointing to wrong vertex";
193 Mesh::HalfedgeHandle next = mesh_.next_halfedge_handle(fh_1_he);
194 EXPECT_EQ(2 , next.idx() ) <<
"Index wrong for second halfedge inside face 1 ";
195 EXPECT_EQ(2 , mesh_.to_vertex_handle(next).idx() ) <<
"second halfedge inside face 1 pointing to wrong vertex ";
197 next = mesh_.next_halfedge_handle(next);
198 EXPECT_EQ(6 , next.idx() ) <<
"Index wrong for third halfedge inside face 1 ";
199 EXPECT_EQ(3 , mesh_.to_vertex_handle(next).idx() ) <<
"Third halfedge inside face 1 pointing to wrong vertex ";
202 Mesh::HalfedgeHandle fh_2_he = mesh_.halfedge_handle(fh_2);
204 EXPECT_EQ(7 , fh_2_he.idx() ) <<
"Index wrong for first halfedge of face 2";
205 EXPECT_EQ(2 , mesh_.to_vertex_handle(fh_2_he).idx() ) <<
"First halfedge inside face 2 pointing to wrong vertex";
207 next = mesh_.next_halfedge_handle(fh_2_he);
208 EXPECT_EQ(3 , next.idx() ) <<
"Index wrong for second halfedge inside face 2";
209 EXPECT_EQ(1 , mesh_.to_vertex_handle(next).idx() ) <<
"second halfedge inside face 2 pointing to wrong vertex ";
211 next = mesh_.next_halfedge_handle(next);
212 EXPECT_EQ(10 , next.idx() ) <<
"Index wrong for third halfedge inside face 2";
213 EXPECT_EQ(3 , mesh_.to_vertex_handle(next).idx() ) <<
"Third halfedge inside face 2 pointing to wrong vertex ";
217 EXPECT_EQ(10 , voh_it.handle().idx() ) <<
"Index wrong for first outgoing halfedge of vertex 1";
219 EXPECT_EQ(2 , voh_it.handle().idx() ) <<
"Index wrong for second outgoing halfedge of vertex 1";
221 EXPECT_EQ(10 , voh_it.handle().idx() ) <<
"Index wrong for third(one lap) outgoing halfedge of vertex 1";
224 voh_it = mesh_.voh_begin(mesh_.vertex_handle(2));
225 EXPECT_EQ(3 , voh_it.handle().idx() ) <<
"Index wrong for first outgoing halfedge of vertex 2";
227 EXPECT_EQ(6 , voh_it.handle().idx() ) <<
"Index wrong for second outgoing halfedge of vertex 2";
229 EXPECT_EQ(3 , voh_it.handle().idx() ) <<
"Index wrong for third(one lap) outgoing halfedge of vertex 2";
232 voh_it = mesh_.voh_begin(mesh_.vertex_handle(3));
233 EXPECT_EQ(11 , voh_it.handle().idx() ) <<
"Index wrong for first outgoing halfedge of vertex 3";
235 EXPECT_EQ(7 , voh_it.handle().idx() ) <<
"Index wrong for second outgoing halfedge of vertex 3";
237 EXPECT_EQ(11 , voh_it.handle().idx() ) <<
"Index wrong for third(one lap) outgoing halfedge of vertex 3";
242 mesh_.garbage_collection();
262 EXPECT_EQ(2u , mesh_.n_faces() ) <<
"Wrong number of faces (garbage collection executed!)";
265 Mesh::FaceHandle fh_0 = mesh_.face_handle(0);
266 fh_1 = mesh_.face_handle(1);
268 fv_it = mesh_.fv_begin(fh_0);
270 EXPECT_EQ(2 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 1 of face 0 after garbage collection";
272 EXPECT_EQ(1 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 2 of face 0 after garbage collection";
274 EXPECT_EQ(0 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 3 of face 0 after garbage collection";
276 fv_it = mesh_.fv_begin(fh_1);
277 EXPECT_EQ(1 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 1 of face 1 after garbage collection";
279 EXPECT_EQ(2 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 2 of face 1 after garbage collection";
281 EXPECT_EQ(0 , fv_it.handle().idx() ) <<
"Index wrong for Vertex 3 of face 1 after garbage collection";
284 Mesh::HalfedgeHandle fh_0_he = mesh_.halfedge_handle(fh_0);
286 EXPECT_EQ(5 , fh_0_he.idx() ) <<
"Index wrong for first halfedge of face 0";
287 EXPECT_EQ(2 , mesh_.to_vertex_handle(fh_0_he).idx() ) <<
"First halfedge inside face 0 pointing to wrong vertex";
289 next = mesh_.next_halfedge_handle(fh_0_he);
290 EXPECT_EQ(3 , next.idx() ) <<
"Index wrong for second halfedge inside face 0 ";
291 EXPECT_EQ(1 , mesh_.to_vertex_handle(next).idx() ) <<
"second halfedge inside face 0 pointing to wrong vertex ";
293 next = mesh_.next_halfedge_handle(next);
294 EXPECT_EQ(0 , next.idx() ) <<
"Index wrong for third halfedge inside face 0 ";
295 EXPECT_EQ(0 , mesh_.to_vertex_handle(next).idx() ) <<
"Third halfedge inside face 0 pointing to wrong vertex ";
298 fh_1_he = mesh_.halfedge_handle(fh_1);
300 EXPECT_EQ(1 , fh_1_he.idx() ) <<
"Index wrong for first halfedge of face 1";
301 EXPECT_EQ(1 , mesh_.to_vertex_handle(fh_1_he).idx() ) <<
"First halfedge inside face 1 pointing to wrong vertex";
303 next = mesh_.next_halfedge_handle(fh_1_he);
304 EXPECT_EQ(2 , next.idx() ) <<
"Index wrong for second halfedge inside face 1 ";
305 EXPECT_EQ(2 , mesh_.to_vertex_handle(next).idx() ) <<
"second halfedge inside face 1 pointing to wrong vertex ";
307 next = mesh_.next_halfedge_handle(next);
308 EXPECT_EQ(4 , next.idx() ) <<
"Index wrong for third halfedge inside face 1 ";
309 EXPECT_EQ(0 , mesh_.to_vertex_handle(next).idx() ) <<
"Third halfedge inside face 1 pointing to wrong vertex ";
313 voh_it = mesh_.voh_begin(mesh_.vertex_handle(0));
314 EXPECT_EQ(1 , voh_it.handle().idx() ) <<
"Index wrong for first outgoing halfedge of vertex 0";
316 EXPECT_EQ(5 , voh_it.handle().idx() ) <<
"Index wrong for second outgoing halfedge of vertex 0";
318 EXPECT_EQ(1 , voh_it.handle().idx() ) <<
"Index wrong for third(one lap) outgoing halfedge of vertex 0";
321 voh_it = mesh_.voh_begin(mesh_.vertex_handle(1));
322 EXPECT_EQ(0 , voh_it.handle().idx() ) <<
"Index wrong for first outgoing halfedge of vertex 1";
324 EXPECT_EQ(2 , voh_it.handle().idx() ) <<
"Index wrong for second outgoing halfedge of vertex 1";
326 EXPECT_EQ(0 , voh_it.handle().idx() ) <<
"Index wrong for third(one lap) outgoing halfedge of vertex 1";
329 voh_it = mesh_.voh_begin(mesh_.vertex_handle(2));
330 EXPECT_EQ(3 , voh_it.handle().idx() ) <<
"Index wrong for first outgoing halfedge of vertex 2";
332 EXPECT_EQ(4 , voh_it.handle().idx() ) <<
"Index wrong for second outgoing halfedge of vertex 2";
334 EXPECT_EQ(3 , voh_it.handle().idx() ) <<
"Index wrong for third(one lap) outgoing halfedge of vertex 2";
336 EXPECT_FALSE( mesh_.is_collapse_ok(mesh_.halfedge_handle(0)) ) <<
"Collapse should be not ok for halfedge 0";
337 EXPECT_FALSE( mesh_.is_collapse_ok(mesh_.halfedge_handle(1)) ) <<
"Collapse should be not ok for halfedge 1";
338 EXPECT_FALSE( mesh_.is_collapse_ok(mesh_.halfedge_handle(2)) ) <<
"Collapse should be not ok for halfedge 2";
339 EXPECT_FALSE( mesh_.is_collapse_ok(mesh_.halfedge_handle(3)) ) <<
"Collapse should be not ok for halfedge 3";
340 EXPECT_FALSE( mesh_.is_collapse_ok(mesh_.halfedge_handle(4)) ) <<
"Collapse should be not ok for halfedge 4";
341 EXPECT_FALSE( mesh_.is_collapse_ok(mesh_.halfedge_handle(5)) ) <<
"Collapse should be not ok for halfedge 5";
352 Mesh::VertexHandle vhandle[5];
355 vhandle[0] = mesh_.add_vertex(Mesh::Point(0 , 0, 0));
356 vhandle[1] = mesh_.add_vertex(Mesh::Point(1 , 0, 0));
357 vhandle[2] = mesh_.add_vertex(Mesh::Point(0 ,-1, 0));
358 vhandle[3] = mesh_.add_vertex(Mesh::Point(0 , 1, 0));
359 vhandle[4] = mesh_.add_vertex(Mesh::Point(-1, 0, 0));
362 std::vector<Mesh::VertexHandle> face_vhandles;
364 face_vhandles.push_back(vhandle[0]);
365 face_vhandles.push_back(vhandle[4]);
366 face_vhandles.push_back(vhandle[2]);
367 mesh_.add_face(face_vhandles);
369 face_vhandles.clear();
371 face_vhandles.push_back(vhandle[3]);
372 face_vhandles.push_back(vhandle[4]);
373 face_vhandles.push_back(vhandle[0]);
374 mesh_.add_face(face_vhandles);
376 face_vhandles.clear();
378 face_vhandles.push_back(vhandle[2]);
379 face_vhandles.push_back(vhandle[4]);
380 face_vhandles.push_back(vhandle[3]);
381 mesh_.add_face(face_vhandles);
383 face_vhandles.clear();
385 face_vhandles.push_back(vhandle[2]);
386 face_vhandles.push_back(vhandle[1]);
387 face_vhandles.push_back(vhandle[0]);
388 mesh_.add_face(face_vhandles);
390 face_vhandles.clear();
392 face_vhandles.push_back(vhandle[0]);
393 face_vhandles.push_back(vhandle[1]);
394 face_vhandles.push_back(vhandle[3]);
395 mesh_.add_face(face_vhandles);
397 mesh_.request_vertex_status();
398 mesh_.request_edge_status();
399 mesh_.request_face_status();
406 Mesh::HalfedgeHandle heh_collapse1 = mesh_.halfedge_handle(0);
408 EXPECT_EQ(4, mesh_.to_vertex_handle(heh_collapse1).idx() ) <<
"To vertex of collapse halfedge 1 is wrong";
409 EXPECT_EQ(0, mesh_.from_vertex_handle(heh_collapse1).idx() ) <<
"from vertex of collapse halfedge 1 is wrong";
411 EXPECT_TRUE( mesh_.is_collapse_ok(heh_collapse1) ) <<
"Collapse not ok for collapse first halfedge (0)";
412 mesh_.collapse(heh_collapse1);
414 Mesh::HalfedgeHandle heh_collapse2 = mesh_.halfedge_handle(2);
416 EXPECT_EQ(2, mesh_.to_vertex_handle(heh_collapse2).idx() ) <<
"To vertex of collapse halfedge 2 is wrong";
417 EXPECT_EQ(4, mesh_.from_vertex_handle(heh_collapse2).idx() ) <<
"from vertex of collapse halfedge 2 is wrong";
419 EXPECT_TRUE( mesh_.is_collapse_ok(heh_collapse2) ) <<
"Collapse not ok for collapse second halfedge (2)";
420 mesh_.collapse(heh_collapse2);
422 Mesh::HalfedgeHandle heh_collapse3 = mesh_.halfedge_handle(6);
424 EXPECT_EQ(2, mesh_.to_vertex_handle(heh_collapse3).idx() ) <<
"To vertex of collapse halfedge 3 is wrong";
425 EXPECT_EQ(3, mesh_.from_vertex_handle(heh_collapse3).idx() ) <<
"from vertex of collapse halfedge 3 is wrong";
427 EXPECT_FALSE( mesh_.is_collapse_ok(heh_collapse3) ) <<
"Collapse not ok for collapse third halfedge (6)";
441 Mesh::VertexHandle vhandle[7];
443 vhandle[0] = mesh_.add_vertex(Mesh::Point( 0, 1, 0));
444 vhandle[1] = mesh_.add_vertex(Mesh::Point( 1, 0, 0));
445 vhandle[2] = mesh_.add_vertex(Mesh::Point( 2, 1, 0));
446 vhandle[3] = mesh_.add_vertex(Mesh::Point( 0,-1, 0));
447 vhandle[4] = mesh_.add_vertex(Mesh::Point( 2,-1, 0));
448 vhandle[5] = mesh_.add_vertex(Mesh::Point(-1, 0, 0));
449 vhandle[6] = mesh_.add_vertex(Mesh::Point( 3, 0, 0));
452 std::vector<Mesh::VertexHandle> face_vhandles;
454 face_vhandles.push_back(vhandle[0]);
455 face_vhandles.push_back(vhandle[5]);
456 face_vhandles.push_back(vhandle[1]);
457 mesh_.add_face(face_vhandles);
459 face_vhandles.clear();
461 face_vhandles.push_back(vhandle[1]);
462 face_vhandles.push_back(vhandle[5]);
463 face_vhandles.push_back(vhandle[3]);
464 mesh_.add_face(face_vhandles);
466 face_vhandles.clear();
468 face_vhandles.push_back(vhandle[0]);
469 face_vhandles.push_back(vhandle[1]);
470 face_vhandles.push_back(vhandle[2]);
471 mesh_.add_face(face_vhandles);
473 face_vhandles.clear();
475 face_vhandles.push_back(vhandle[1]);
476 face_vhandles.push_back(vhandle[3]);
477 face_vhandles.push_back(vhandle[4]);
478 mesh_.add_face(face_vhandles);
480 face_vhandles.clear();
482 face_vhandles.push_back(vhandle[2]);
483 face_vhandles.push_back(vhandle[1]);
484 face_vhandles.push_back(vhandle[4]);
485 mesh_.add_face(face_vhandles);
487 face_vhandles.clear();
489 face_vhandles.push_back(vhandle[2]);
490 face_vhandles.push_back(vhandle[4]);
491 face_vhandles.push_back(vhandle[6]);
492 mesh_.add_face(face_vhandles);
504 mesh_.request_vertex_status();
505 mesh_.request_edge_status();
506 mesh_.request_face_status();
512 Mesh::HalfedgeHandle heh_collapse;
515 for ( Mesh::HalfedgeIter he_it = mesh_.halfedges_begin() ; he_it != mesh_.halfedges_end() ; ++he_it ) {
516 if ( mesh_.from_vertex_handle(he_it).idx() == 1 && mesh_.to_vertex_handle(he_it).idx() == 4 )
517 heh_collapse = he_it;
521 EXPECT_EQ(4, mesh_.to_vertex_handle(heh_collapse).idx() ) <<
"To vertex of collapse halfedge is wrong";
522 EXPECT_EQ(1, mesh_.from_vertex_handle(heh_collapse).idx() ) <<
"from vertex of collapse halfedge is wrong";
523 EXPECT_TRUE( mesh_.is_collapse_ok(heh_collapse) ) <<
"Collapse not ok for collapse first halfedge (0)";
526 Mesh::VertexHandle vh_from = mesh_.from_vertex_handle(heh_collapse);
527 Mesh::VertexHandle vh_to = mesh_.to_vertex_handle(heh_collapse);
530 mesh_.collapse(heh_collapse);
532 EXPECT_TRUE( mesh_.status(vh_from).deleted() ) <<
"From vertex not deleted";
533 EXPECT_FALSE( mesh_.status(vh_to).deleted() ) <<
"To Vertex deleted";
538 #endif // INCLUDE GUARD