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
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
|
// Gmsh - Copyright (C) 1997-2021 C. Geuzaine, J.-F. Remacle
//
// See the LICENSE.txt file for license information. Please report all
// issues on https://gitlab.onelab.info/gmsh/gmsh/issues.
#include <string.h>
#include "GmshConfig.h"
#include "CutSphere.h"
#include "Context.h"
#if defined(HAVE_OPENGL)
#include "drawContext.h"
#endif
StringXNumber CutSphereOptions_Number[] = {
{GMSH_FULLRC, "Xc", GMSH_CutSpherePlugin::callbackX, 0.},
{GMSH_FULLRC, "Yc", GMSH_CutSpherePlugin::callbackY, 0.},
{GMSH_FULLRC, "Zc", GMSH_CutSpherePlugin::callbackZ, 0.},
{GMSH_FULLRC, "R", GMSH_CutSpherePlugin::callbackR, 0.25},
{GMSH_FULLRC, "ExtractVolume", GMSH_CutSpherePlugin::callbackVol, 0.},
{GMSH_FULLRC, "RecurLevel", GMSH_CutSpherePlugin::callbackRecur, 4},
{GMSH_FULLRC, "TargetError", GMSH_CutSpherePlugin::callbackTarget, 0.},
{GMSH_FULLRC, "View", nullptr, -1.}};
extern "C" {
GMSH_Plugin *GMSH_RegisterCutSpherePlugin()
{
return new GMSH_CutSpherePlugin();
}
}
void GMSH_CutSpherePlugin::draw(void *context)
{
#if defined(HAVE_OPENGL)
GLint mode[2];
glGetIntegerv(GL_POLYGON_MODE, mode);
glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
glColor4ubv((GLubyte *)&CTX::instance()->color.fg);
glLineWidth((float)CTX::instance()->lineWidth);
drawContext *ctx = (drawContext *)context;
ctx->drawSphere(
CutSphereOptions_Number[3].def, CutSphereOptions_Number[0].def,
CutSphereOptions_Number[1].def, CutSphereOptions_Number[2].def, 40, 40, 1);
glPolygonMode(GL_FRONT_AND_BACK, mode[1]);
#endif
}
double GMSH_CutSpherePlugin::callback(int num, int action, double value,
double *opt, double step, double min,
double max)
{
switch(action) { // configure the input field
case 1: return step;
case 2: return min;
case 3: return max;
default: break;
}
*opt = value;
GMSH_Plugin::setDrawFunction(draw);
return 0.;
}
double GMSH_CutSpherePlugin::callbackX(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[0].def,
CTX::instance()->lc / 100., -2 * CTX::instance()->lc,
2 * CTX::instance()->lc);
}
double GMSH_CutSpherePlugin::callbackY(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[1].def,
CTX::instance()->lc / 100., -2 * CTX::instance()->lc,
2 * CTX::instance()->lc);
}
double GMSH_CutSpherePlugin::callbackZ(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[2].def,
CTX::instance()->lc / 100., -2 * CTX::instance()->lc,
2 * CTX::instance()->lc);
}
double GMSH_CutSpherePlugin::callbackR(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[3].def,
CTX::instance()->lc / 100., 0., 2 * CTX::instance()->lc);
}
double GMSH_CutSpherePlugin::callbackVol(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[4].def, 1., -1.,
1.);
}
double GMSH_CutSpherePlugin::callbackRecur(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[5].def, 1, 0,
10);
}
double GMSH_CutSpherePlugin::callbackTarget(int num, int action, double value)
{
return callback(num, action, value, &CutSphereOptions_Number[6].def, 0.01, 0.,
1.);
}
std::string GMSH_CutSpherePlugin::getHelp() const
{
return "Plugin(CutSphere) cuts the view `View' with the "
"sphere (X-`Xc')^2 + (Y-`Yc')^2 + (Z-`Zc')^2 = `R'^2.\n\n"
"If `ExtractVolume' is nonzero, the plugin extracts "
"the elements inside (if `ExtractVolume' < 0) or "
"outside (if `ExtractVolume' > 0) the sphere.\n\n"
"If `View' < 0, the plugin is run on the current view.\n\n"
"Plugin(CutSphere) creates one new list-based view.";
}
int GMSH_CutSpherePlugin::getNbOptions() const
{
return sizeof(CutSphereOptions_Number) / sizeof(StringXNumber);
}
StringXNumber *GMSH_CutSpherePlugin::getOption(int iopt)
{
return &CutSphereOptions_Number[iopt];
}
double GMSH_CutSpherePlugin::levelset(double x, double y, double z,
double val) const
{
double a = CutSphereOptions_Number[0].def;
double b = CutSphereOptions_Number[1].def;
double c = CutSphereOptions_Number[2].def;
double r = CutSphereOptions_Number[3].def;
return (x - a) * (x - a) + (y - b) * (y - b) + (z - c) * (z - c) - r * r;
}
PView *GMSH_CutSpherePlugin::execute(PView *v)
{
int iView = (int)CutSphereOptions_Number[7].def;
_ref[0] = CutSphereOptions_Number[0].def;
_ref[1] = CutSphereOptions_Number[1].def;
_ref[2] = CutSphereOptions_Number[2].def;
_extractVolume = (int)CutSphereOptions_Number[4].def;
_recurLevel = (int)CutSphereOptions_Number[5].def;
_targetError = CutSphereOptions_Number[6].def;
_valueIndependent = 1;
_valueView = -1;
_valueTimeStep = -1;
_orientation = GMSH_LevelsetPlugin::SPHERE;
PView *v1 = getView(iView, v);
if(!v1) return v;
return GMSH_LevelsetPlugin::execute(v1);
}
|
