circulation.cpp 1.3 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162
  1. #include <test_common.h>
  2. #include <igl/circulation.h>
  3. #include <igl/edge_flaps.h>
  4. #include <igl/unique_edge_map.h>
  5. #include <igl/matlab_format.h>
  6. TEST(circulation,single_edge)
  7. {
  8. // 7
  9. // /₆|₇\
  10. // 4 - 5 - 6
  11. // |₂/₃|₄\₅|
  12. // 1 - 2 - 3
  13. // \₀|₁/
  14. // 0
  15. const Eigen::MatrixXi F = (Eigen::MatrixXi(8,3)<<
  16. 0,2,1,
  17. 0,3,2,
  18. 1,5,4,
  19. 1,2,5,
  20. 2,3,5,
  21. 3,6,5,
  22. 4,5,7,
  23. 5,6,7).finished();
  24. Eigen::MatrixXi E,uE;
  25. Eigen::VectorXi EMAP;
  26. std::vector<std::vector<int> > uE2E;
  27. igl::unique_edge_map(F, E, uE, EMAP, uE2E);
  28. Eigen::MatrixXi EI,EF;
  29. {
  30. const auto & cuE = uE;
  31. const auto & cEMAP = EMAP;
  32. igl::edge_flaps(F,cuE,cEMAP,EF,EI);
  33. }
  34. // Find (2,5) in uE
  35. int ei = 0;
  36. for(;ei<E.rows();ei++)
  37. {
  38. if(uE(ei,0) == 2 && uE(ei,1) == 5) break;
  39. if(uE(ei,1) == 2 && uE(ei,0) == 5) break;
  40. }
  41. Eigen::VectorXi Nccw;
  42. igl::circulation(ei,true,EMAP,EF,EI,Nccw);
  43. const Eigen::VectorXi Nccwgt =
  44. (Eigen::VectorXi(6)<<
  45. 4,
  46. 5,
  47. 7,
  48. 6,
  49. 2,
  50. 3).finished();
  51. test_common::assert_eq(Nccw,Nccwgt);
  52. Eigen::VectorXi Ncw;
  53. igl::circulation(ei,false,EMAP,EF,EI,Ncw);
  54. const Eigen::VectorXi Ncwgt =
  55. (Eigen::VectorXi(4)<<
  56. 4,
  57. 1,
  58. 0,
  59. 3).finished();
  60. test_common::assert_eq(Ncw,Ncwgt);
  61. }