三角メッシュ用のライブラリ書いた

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:
{}