ng_vol.cpp


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73

/*
  original from  Joachim Schoeberl
  ..../netgen-4.9.13/nglib/ng_vol.cpp
  License: LGPL
*/
#include <fstream>
#include <iostream>

using namespace std;

namespace nglib {
#include <nglib.h>
}

int main(int argc, char **argv) {
  using namespace nglib;

  cout << "Netgen tet-mesher" << endl;

  if (argc < 2) {
    cerr << "use: ng_vol filename" << endl;
    return 1;
  }

  Ng_Mesh *mesh;

  Ng_Init();

  // creates mesh structure
  mesh = Ng_NewMesh();

  int i, np, nse, ne;
  double point[3];
  int trig[3], tet[4];

  // reads surface mesh from file
  ifstream in(argv[1]);

  in >> np;
  cout << "Reading " << np << " points...";
  cout.flush();
  for (i = 1; i <= np; i++) {
    in >> point[0] >> point[1] >> point[2];
    Ng_AddPoint(mesh, point);
  }
  cout << "done" << endl;

  in >> nse;
  cout << "Reading " << nse << " faces...";
  cout.flush();
  for (i = 1; i <= nse; i++) {
    in >> trig[0] >> trig[1] >> trig[2];
    Ng_AddSurfaceElement(mesh, NG_TRIG, trig);
  }
  cout << "done" << endl;

  // generate volume mesh
  Ng_Meshing_Parameters mp;
  mp.maxh = 1e6;
  mp.fineness = 1;
  mp.second_order = 0;

  if (argc == 3)
    mp.maxh = atof(argv[2]);

  cout << "start meshing" << endl;
  Ng_GenerateVolumeMesh(mesh, &mp);
  cout << "meshing done" << endl;

  Ng_SaveMesh(mesh, "test.vol");

  return 0;
}