GitHub - tak0kada/cnthd: Halfedge data structure for triangular mesh ハーフエッジデータ構造のC++ライブラリ。面の向き付けがバラバラなファイルに対応しようと思っていたがバグを出してしまって(しかもそれほど使わない機能)書き直す時間が惜しいのでまともなファイルが来るのが前提になっている。テストもいい加減で目視で確かめただけである。以下は使用例。
- example (test_mesh.cpp)
#include <iostream> #include <string> #include "mesh.hpp" int main(int argc, char* argv[]) { using namespace cnthd; if (argc == 1) { std::cout << "./a.out input.obj output.obj" << std::endl; return 0; } std::string infile{argv[1]}; std::string outfile{argv[2]}; Mesh mesh = read_obj(infile); std::cout << mesh << std::endl; write_obj(mesh, outfile); return 0; }
cube_small.obj
v 0 0 0 v 0 0 1 v 0 1 0 v 0 1 1 v 1 0 0 v 1 0 1 v 1 1 0 v 1 1 1 f 1 5 2 f 2 5 6 f 2 6 8 f 2 8 4 f 5 7 8 f 5 8 6 f 1 3 5 f 5 3 7 f 1 2 3 f 2 4 3 f 3 4 8 f 3 8 7
- stdout
nV: 8, nE: 18, nF: 12, g = 1 - 0.5 * (nV - nE + nF) = 0 vertex: v1: {0.000000, 0.000000, 0.000000} v2: {0.000000, 0.000000, 1.000000} v3: {0.000000, 1.000000, 0.000000} v4: {0.000000, 1.000000, 1.000000} v5: {1.000000, 0.000000, 0.000000} v6: {1.000000, 0.000000, 1.000000} v7: {1.000000, 1.000000, 0.000000} v8: {1.000000, 1.000000, 1.000000} edge: e1: 1-5 e2: 5-2 e3: 2-1 e4: 5-6 e5: 6-2 e6: 6-8 e7: 8-2 e8: 8-4 e9: 4-2 e10: 5-7 e11: 7-8 e12: 8-5 e13: 1-3 e14: 3-5 e15: 3-7 e16: 2-3 e17: 4-3 e18: 8-3 halfedge: he: 1->5 he: 5->2 he: 2->1 he: 2->5 he: 5->6 he: 6->2 he: 2->6 he: 6->8 he: 8->2 he: 2->8 he: 8->4 he: 4->2 he: 5->7 he: 7->8 he: 8->5 he: 5->8 he: 8->6 he: 6->5 he: 1->3 he: 3->5 he: 5->1 he: 5->3 he: 3->7 he: 7->5 he: 1->2 he: 2->3 he: 3->1 he: 2->4 he: 4->3 he: 3->2 he: 3->4 he: 4->8 he: 8->3 he: 3->8 he: 8->7 he: 7->3 face: f1: 1-5-2 f2: 2-5-6 f3: 2-6-8 f4: 2-8-4 f5: 5-7-8 f6: 5-8-6 f7: 1-3-5 f8: 5-3-7 f9: 1-2-3 f10: 2-4-3 f11: 3-4-8 f12: 3-8-7 adjacent matrix: 0 1 1 0 1 0 0 0 1 0 1 1 1 1 0 1 1 1 0 1 1 0 1 1 0 1 1 0 0 0 0 1 1 1 1 0 0 1 1 1 0 1 0 0 1 0 0 1 0 0 1 0 1 0 0 1 0 1 1 1 1 1 1 0 adjacent list: {{4, 1, 2}, {4, 0, 5, 7, 3, 2}, {0, 4, 6, 1, 3, 7}, {7, 1, 2}, {0, 1, 5, 6, 7, 2}, {4, 1, 7}, {4, 7, 2}, {5, 1, 3, 6, 4, 2}}
null.obj
nV: 0, nE: 0, nF: 0, g = 1 - 0.5 * (nV - nE + nF) = 1 vertex: edge: halfedge: face: adjacent matrix: adjacent list: {}