任务

很简单：我想创建一个四面体网格，我想为其添加几个特征。

例

作为输入，请获取可在CGAL-4.11/examples/Mesh_3/data中找到的文件cube.off。我想添加的功能(仅立方体的十二条边(保存在cube.edges：

2 -1 -1 -1 -1 1 -1
2 -1 -1 -1 1 -1 -1
2 -1 -1 -1 -1 -1 1
2 -1 1 -1 1 1 -1
2 -1 1 -1 -1 1 1
2 1 1 -1 1 -1 -1
2 1 1 -1 1 1 1
2 1 -1 -1 1 -1 1
2 -1 -1 1 -1 1 1
2 -1 -1 1 1 -1 1
2 -1 1 1 1 1 1
2 1 1 1 1 -1 1

代码的 MWE：

#include <CGAL/Exact_predicates_inexact_constructions_kernel.h>
#include <CGAL/Mesh_triangulation_3.h>
#include <CGAL/Mesh_complex_3_in_triangulation_3.h>
#include <CGAL/Mesh_criteria_3.h>
#include <CGAL/Polyhedral_complex_mesh_domain_3.h>
#include <CGAL/make_mesh_3.h>
// From CGAL-4.11/examples/Mesh_3 but for simplicity copied to the current folder.
#include "read_polylines.h"
typedef CGAL::Exact_predicates_inexact_constructions_kernel K;
typedef CGAL::Mesh_polyhedron_3<K>::type Polyhedron;
typedef CGAL::Polyhedral_complex_mesh_domain_3<K> Mesh_domain;
typedef CGAL::Sequential_tag Concurrency_tag;
typedef CGAL::Mesh_triangulation_3<Mesh_domain,CGAL::Default,Concurrency_tag>::type Tr;
typedef CGAL::Mesh_complex_3_in_triangulation_3<
Tr,Mesh_domain::Corner_index,Mesh_domain::Curve_segment_index> C3t3;
typedef CGAL::Mesh_criteria_3<Tr> Mesh_criteria;
typedef K::Point_3 Point;
typedef std::vector<std::vector<Point> > Polylines;
int main()
{
// Read the (one) patch.
std::vector<Polyhedron> patches(1);
std::ifstream input("cube.off");
input >> patches[0];
const std::pair<int, int> incident_subdomains[] = { std::make_pair(1,0) };
Mesh_domain domain(patches.begin(), patches.end(), incident_subdomains, incident_subdomains+1);
// Read the features.
std::string feature_edges="cube.edges";
Polylines polylines;
read_polylines<Point>(feature_edges.c_str(), polylines);
domain.add_features(polylines.begin(), polylines.end());
// Create the mesh.
Mesh_criteria criteria(CGAL::parameters::edge_size = 0.25);
C3t3 c3t3 = CGAL::make_mesh_3<C3t3>(domain, criteria);
// Write it (if it hasn't crashed before).
std::ofstream medit_file("out.mesh");
c3t3.output_to_medit(medit_file, false, true);
}

这(使用选项-DCGAL_MESH_3_VERBOSE编译时(崩溃并显示以下错误消息：

Start volume scan...Scanning triangulation for bad cells (sequential)... terminate called after throwing an instance of 'CGAL::Assertion_exception'
what():  CGAL ERROR: assertion violation!
Expr: patch_id > 0 && std::size_t(patch_id) < patch_id_to_polyhedron_id.size()
File: /yatmp/scr1/ya10813/cgal-install/include/CGAL/Polyhedral_complex_mesh_domain_3.h
Line: 457
Aborted

问题

我错过了什么？它一定是非常基本的东西。CGAL 是一个非常可靠的库，可以很好地处理复杂数据;我不相信一个有 12 条边的立方体足以阻止它。

我也尝试过的事情

• 当我将行domain.add_features(polylines.begin(), polylines.end());替换为domain.detect_features();时，程序正确终止。令人困惑的是，检测到的特征正是我尝试添加的特征(我知道，因为我在Mesh_domain_with_polyline_features_3中制作了一个为我打印边缘的功能;如有必要，我可以在这里分享(。

• 当我使用Polyhedral_mesh_domain_with_features_3而不是Polyhedral_complex_mesh_domain_3时，程序正确终止。这些特征似乎保留在生成的几何中。但是，它只是一个补丁(out.mesh中的所有三角形都具有相同的最后一个数字，而在这种情况下我希望它们具有数字 0 到 11(。编辑：为了实现这一点，必须使用detect_features();而不是add_features(...);，或者将域拆分为补丁，让CGAL从中创建复合体。感谢@lrineau的澄清。

• 我还尝试先将域拆分为几个补丁。但是，随后补丁的边界会发生变化。编辑：保留它们的一种方法是将面片边界添加为要素。

• 反转表面的方向(并在incident_subdomains中交换0和1(：没有明显的变化。

• 更改特征线的顺序cube.edges：无可见变化。

• 与旧版本的 Boost 链接：没有明显的变化。