Home Explore Help
Sign In
cpp
/
libigl
mirror of https://github.com/libigl/libigl.git
2
0
Fork 0
Files Issues 0 Wiki
Branch: alecjacobson/parallel_for-no-nesting
Branches Tags
libigl
/
include
/
igl
/
copyleft
/
cgal
/
fast_winding_number.cpp

fast_winding_number.cpp 2.2 KB
Permalink History Raw

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263 
  1. #include "fast_winding_number.h"
    2. 

  2. #include "../../fast_winding_number.h"
    3. 

  3. #include "../../octree.h"
    4. 

  4. #include "../../knn.h"
    5. 

  5. #include "../../parallel_for.h"
    6. 

  6. #include "point_areas.h"
    7. 

  7. #include <vector>
    8. 

  8. 

  9. template <
    10. 

  10.   typename DerivedP, 
    11. 

  11.   typename DerivedN, 
    12. 

  12.   typename DerivedQ,
    13. 

  13.   typename BetaType, 
    14. 

  14.   typename DerivedWN>
    15. 

  15. IGL_INLINE void igl::copyleft::cgal::fast_winding_number(
    16. 

  16.   const Eigen::MatrixBase<DerivedP>& P,
    17. 

  17.   const Eigen::MatrixBase<DerivedN>& N,
    18. 

  18.   const Eigen::MatrixBase<DerivedQ>& Q,
    19. 

  19.   const int expansion_order,
    20. 

  20.   const BetaType beta,
    21. 

  21.   Eigen::PlainObjectBase<DerivedWN>& WN)
    22. 

  22. {
    23. 

  23.   typedef typename DerivedWN::Scalar real;
    24. 

  24.         
    25. 

  25.   std::vector<std::vector<int> > point_indices;
    26. 

  26.   Eigen::Matrix<int,Eigen::Dynamic,8> CH;
    27. 

  27.   Eigen::Matrix<real,Eigen::Dynamic,3> CN;
    28. 

  28.   Eigen::Matrix<real,Eigen::Dynamic,1> W;
    29. 

  29.   Eigen::MatrixXi I;
    30. 

  30.   Eigen::Matrix<real,Eigen::Dynamic,1> A;
    31. 

  31.   
    32. 

  32.   octree(P,point_indices,CH,CN,W);
    33. 

  33.   knn(P,21,point_indices,CH,CN,W,I);
    34. 

  34.   point_areas(P,I,N,A);
    35. 

  35.   
    36. 

  36.   Eigen::Matrix<real,Eigen::Dynamic,Eigen::Dynamic> EC;
    37. 

  37.   Eigen::Matrix<real,Eigen::Dynamic,3> CM;
    38. 

  38.   Eigen::Matrix<real,Eigen::Dynamic,1> R;
    39. 

  39.   
    40. 

  40.   igl::fast_winding_number(
    41. 

  41.     P,N,A,point_indices,CH,expansion_order,CM,R,EC);
    42. 

  42.   igl::fast_winding_number(
    43. 

  43.     P,N,A,point_indices,CH,CM,R,EC,Q,beta,WN);
    44. 

  44. }
    45. 

  45.       
    46. 

  46. template <
    47. 

  47.   typename DerivedP, 
    48. 

  48.   typename DerivedN, 
    49. 

  49.   typename DerivedQ, 
    50. 

  50.   typename DerivedWN>
    51. 

  51. IGL_INLINE void igl::copyleft::cgal::fast_winding_number(
    52. 

  52.   const Eigen::MatrixBase<DerivedP>& P,
    53. 

  53.   const Eigen::MatrixBase<DerivedN>& N,
    54. 

  54.   const Eigen::MatrixBase<DerivedQ>& Q,
    55. 

  55.   Eigen::PlainObjectBase<DerivedWN>& WN)
    56. 

  56. {
    57. 

  57.   fast_winding_number(P,N,Q,2,2.0,WN);
    58. 

  58. }
    59. 

  59. 

  60. #ifdef IGL_STATIC_LIBRARY
    61. 

  61. // Explicit template instantiation
    62. 

  62. template void igl::copyleft::cgal::fast_winding_number<Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, Eigen::Matrix<double, -1, -1, 0, -1, -1>, double, Eigen::Matrix<double, -1, -1, 0, -1, -1> >(Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, Eigen::MatrixBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> > const&, int, double, Eigen::PlainObjectBase<Eigen::Matrix<double, -1, -1, 0, -1, -1> >&);
    63. 

  63. #endif
    64.