c++ - CGAL:添加功能失败
问题描述
任务
很简单:我想创建一个四面体网格,我想向它添加几个特征。
例子
作为输入,获取cube.off
可以在CGAL-4.11/examples/Mesh_3/data
. 我想添加的功能(只是立方体的十二条边)保存在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 的澄清。我还尝试先将域分成几个补丁。但是,补丁的边界会发生变化。编辑:保留它们的一种方法是将补丁边界添加为特征。
反转表面的方向(并交换
0
和1
)incident_subdomains
:没有可见的变化。更改要素线的顺序
cube.edges
:无可见变化。与旧版本的 Boost 链接:没有明显的变化。
解决方案
该课程Polyhedral_complex_mesh_domain_3
是一年前在 CGAL-4.11 中添加的新课程。你遇到的问题是,它从来没有在没有调用的情况下测试过domain.detect_features()
,并且代码有一个错误:数据成员patch_id_to_polyhedron_id
只填写了detect_features()
. 这就解释了为什么当您调用detect_features()
而不是add_features(..)
.
编辑:我今天已经修复了这个错误,请参阅CGAL 的 PR #3393。该修复将在未来的错误修复版本 CGAL-4.12.2 和 CGAL-4.13.1 中进行。
推荐阅读
- python - 不知道为什么这个函数不能执行正确的数学运算
- docker - 如何在 Kubernetes 中为容器设置 VPN 和反向代理?
- docker - Wireguard 无法解析 Docker 上的本地 DNS
- python - 通过在 2 个字段上应用聚合在 Django 中创建评估字段
- html - 输入字段中的占位符只显示给我,不显示给其他人
- reactjs - 在所需位置显示 recharts xaxis 刻度的问题
- flutter - WidgetSpan 访问父 TextSpan fontSize
- python - 无法创建文件错误,文件不存在
- mysql - SQL - 如何从过去 5 年中选择日期以包括最早年份的 1 月 1 日
- javascript - 带有自定义字体系列的 Syncfusion JavaScript (ES5) 导出图表