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
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
|
// -----------------------------------------------------------------------------
//
// Gmsh C++ tutorial 9
//
// Plugins
//
// -----------------------------------------------------------------------------
// Plugins can be added to Gmsh in order to extend its capabilities. For
// example, post-processing plugins can modify views, or create new views based
// on previously loaded views. Several default plugins are statically linked
// with Gmsh, e.g. Isosurface, CutPlane, CutSphere, Skin, Transform or Smooth.
//
// Plugins can be controlled through the API functions in the `gmsh::plugin'
// namespace, or from the graphical interface (right click on the view button,
// then `Plugins').
#include <set>
#include <gmsh.h>
int main(int argc, char **argv)
{
gmsh::initialize();
gmsh::model::add("t9");
// Let us for example include a three-dimensional scalar view:
try {
gmsh::merge("../view3.pos");
} catch(...) {
gmsh::logger::write("Could not load post-processing views: bye!");
gmsh::finalize();
return 0;
}
std::vector<int> v;
gmsh::view::getTags(v);
if(v.size() != 1) {
gmsh::logger::write("Wrong number of views!", "error");
gmsh::finalize();
return 1;
}
// We then set some options for the `Isosurface' plugin (which extracts an
// isosurface from a 3D scalar view), and run it:
gmsh::plugin::setNumber("Isosurface", "Value", 0.67); // Iso-value level
gmsh::plugin::setNumber("Isosurface", "View", 0); // Source view is View[0]
int v1 = gmsh::plugin::run("Isosurface"); // Run the plugin!
// We also set some options for the `CutPlane' plugin (which computes a
// section of a 3D view using the plane A*x+B*y+C*z+D=0), and then run it:
gmsh::plugin::setNumber("CutPlane", "A", 0);
gmsh::plugin::setNumber("CutPlane", "B", 0.2);
gmsh::plugin::setNumber("CutPlane", "C", 1);
gmsh::plugin::setNumber("CutPlane", "D", 0);
gmsh::plugin::setNumber("CutPlane", "View", 0);
int v2 = gmsh::plugin::run("CutPlane");
// Add a title (By convention, for window coordinates a value greater than
// 99999 represents the center. We could also use `General.GraphicsWidth / 2',
// but that would only center the string for the current window size.):
gmsh::plugin::setString("Annotate", "Text", "A nice title");
gmsh::plugin::setNumber("Annotate", "X", 1.e5);
gmsh::plugin::setNumber("Annotate", "Y", 50);
gmsh::plugin::setString("Annotate", "Font", "Times-BoldItalic");
gmsh::plugin::setNumber("Annotate", "FontSize", 28);
gmsh::plugin::setString("Annotate", "Align", "Center");
gmsh::plugin::setNumber("Annotate", "View", 0);
gmsh::plugin::run("Annotate");
gmsh::plugin::setString("Annotate", "Text", "(and a small subtitle)");
gmsh::plugin::setNumber("Annotate", "Y", 70);
gmsh::plugin::setString("Annotate", "Font", "Times-Roman");
gmsh::plugin::setNumber("Annotate", "FontSize", 12);
gmsh::plugin::run("Annotate");
// We finish by setting some options:
gmsh::view::option::setNumber(v[0], "Light", 1);
gmsh::view::option::setNumber(v[0], "IntervalsType", 1);
gmsh::view::option::setNumber(v[0], "NbIso", 6);
gmsh::view::option::setNumber(v[0], "SmoothNormals", 1);
gmsh::view::option::setNumber(v1, "IntervalsType", 2);
gmsh::view::option::setNumber(v2, "IntervalsType", 2);
// Launch the GUI to see the results:
std::set<std::string> args(argv, argv + argc);
if(!args.count("-nopopup")) gmsh::fltk::run();
gmsh::finalize();
return 0;
}
|
