解决cgal中的surface_mesh_deformation,代码copy到别处报错

我只想说我遇到过多少奇葩的bug啊。。。这么奇葩的错误我竟然能debug出来真是醉醉的


cgal4.5之后,有实现surface_mesh_deformation的example。

1.如何正常编译example?

最好是在编译cgal的时候,cmake中的entry就有,勾选上with_eigen,然后配置include目录即可。

这样以后省事,所以编译cgal的时候,要好好检查外置的库,会用到的都先加上。

如果不想重新编译cgal,在重新用cmake编译surface_mesh_deformation这个example和demo时,选择with_eigen即可。


2.当你正常运行example之后,把代码原样copy到自己的project中,编译会报错TTTTTTT

1>C:\Program Files\CGAL-4.6.1\include\CGAL/Surface_mesh_deformation.h(212): error C2039: 'Matrix' : is not a member of 'CGAL::Default' 
1>          C:\Program Files\CGAL-4.6.1\include\CGAL/Default.h(31) : see declaration of 'CGAL::Default' 
1>          ModelDecompoent.cpp(245) : see reference to class template instantiation 'CGAL::Surface_mesh_deformation<HG>' being compiled 
1>          with 
1>          [ 
1>              HG=Polyhedron 
1>          ] 

......
一堆错

翻墙说加上#define CGAL_EIGEN3_ENABLED
因为看CGAL/Surface_mesh_deformation.h代码中:
#if defined(CGAL_EIGEN3_ENABLED)
#include <CGAL/Eigen_solver_traits.h>  // for sparse linear system solver
#include <CGAL/Deformation_Eigen_polar_closest_rotation_traits_3.h>  // for 3x3 closest rotation computer
#endif
就是说CGAL_EIGEN3_ENABLED被define了才可以。

然而我加上并没有用!!!

为什么呢?
因为顺序反了= =

正确的顺序:
#define CGAL_EIGEN3_ENABLED
#include <CGAL/Surface_mesh_deformation.h>


错误的顺序:
#include <CGAL/Surface_mesh_deformation.h>
#define CGAL_EIGEN3_ENABLED //这样子都include完了,define还有什么用。。。

真是奇葩的bug

要在 `CGAL::Surface_mesh` 上挖洞,您需要执行以下步骤: 1. 确定要挖去的区域,并创建表示该区域的多边形。 2. 使用 `CGAL::Polygon_mesh_processing::triangulate_hole()` 函数将多边形三角化,生成一个新的三角网格。 3. 使用 `CGAL::Polygon_mesh_processing::stitch_borders()` 函数将新的三角网格与原始网格连接起来,生成一个没有挖空区域的新网格。 下面是一个简单的示例代码,用于在 `CGAL::Surface_mesh` 上挖一个圆形洞: ```cpp #include <CGAL/Simple_cartesian.h> #include <CGAL/Polygon_mesh_processing/triangulate_hole.h> #include <CGAL/Polygon_mesh_processing/stitch_borders.h> #include <CGAL/Surface_mesh.h> typedef CGAL::Simple_cartesian<double> Kernel; typedef Kernel::Point_3 Point_3; int main() { CGAL::Surface_mesh<Point_3> mesh; // 添加一些三角形到网格... // 创建表示要挖去的区域的多边形 std::vector<Point_3> hole; double cx = 0.0, cy = 0.0, r = 1.0; int n = 20; for (int i = 0; i < n; ++i) { double theta = 2 * M_PI * i / n; double x = cx + r * std::cos(theta); double y = cy + r * std::sin(theta); hole.push_back(Point_3(x, y, 0.0)); } // 三角化洞 CGAL::Polygon_mesh_processing::triangulate_hole(mesh, hole.begin(), hole.end()); // 连接边界 CGAL::Polygon_mesh_processing::stitch_borders(mesh); // 现在,mesh 已经成功挖去了一个圆形洞 return 0; } ``` 请注意,上述代码仅是一个示例,并且可能需要根据您的具体应用进行修改。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值