-
Notifications
You must be signed in to change notification settings - Fork 2
/
GraysKlein.cpp
63 lines (52 loc) · 1.37 KB
/
GraysKlein.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
#include "GraysKlein.h"
#include "Point.h"
#include <glm/gtc/constants.hpp>
#include <glm/vec3.hpp>
#include <glm/glm.hpp>
GraysKlein::GraysKlein(QOpenGLShaderProgram* prog, float radius, unsigned int nSlices, unsigned int nStacks, unsigned int sMax, unsigned int tMax) :
ParametricSurface(prog, nSlices, nStacks, sMax, tMax),
_A(2),
_M(1),
_N(2),
_radius(radius)
{
setAutoIncrName("Gray's Klein Bottle");
buildMesh();
}
GraysKlein::~GraysKlein()
{
}
TriangleMesh* GraysKlein::clone()
{
return new GraysKlein(_prog, _radius, _slices, _stacks, _sMax, _tMax);
}
float GraysKlein::firstUParameter() const
{
return 0.0;
}
float GraysKlein::lastUParameter() const
{
return 4 * glm::pi<float>();
}
float GraysKlein::firstVParameter() const
{
return 0.0;
}
float GraysKlein::lastVParameter() const
{
return glm::two_pi<float>();
}
Point GraysKlein::pointAtParameter(const float& u, const float& v)
{
Point P;
float x, y, z;
// Gray's Klein bottle
//0 <= u <= 4 PI
//0 <= v <= 2 PI
//A = 2, N = 2, M = 1
x = _radius * (_A + cos(_N * u / 2.0) * sin(v) - sin(_N * u / 2.0) * sin(2 * v)) * cos(_M * u / 2.0);
y = _radius * (_A + cos(_N * u / 2.0) * sin(v) - sin(_N * u / 2.0) * sin(2 * v)) * sin(_M * u / 2.0);
z = _radius * sin(_N * u / 2.0) * sin(v) + cos(_N * u / 2.0) * sin(2 * v);
P.setParam(x, y, z);
return P;
}