-
Notifications
You must be signed in to change notification settings - Fork 2
/
SaddleTorus.cpp
61 lines (50 loc) · 1.49 KB
/
SaddleTorus.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
#include "SaddleTorus.h"
#include "Point.h"
#include <glm/gtc/constants.hpp>
#include <glm/vec3.hpp>
#include <glm/glm.hpp>
SaddleTorus::SaddleTorus(QOpenGLShaderProgram* prog, float radius, unsigned int nSlices, unsigned int nStacks, unsigned int sMax, unsigned int tMax) :
ParametricSurface(prog, nSlices, nStacks, sMax, tMax),
_radius(radius)
{
setAutoIncrName("Saddle Torus");
buildMesh();
}
SaddleTorus::~SaddleTorus()
{
}
TriangleMesh* SaddleTorus::clone()
{
return new SaddleTorus(_prog, _radius, _slices, _stacks, _sMax, _tMax);
}
float SaddleTorus::firstUParameter() const
{
return 0.0;
}
float SaddleTorus::lastUParameter() const
{
return glm::two_pi<float>();
}
float SaddleTorus::firstVParameter() const
{
return 0.0;
}
float SaddleTorus::lastVParameter() const
{
return glm::two_pi<float>();
}
Point SaddleTorus::pointAtParameter(const float& u, const float& v)
{
Point P;
float x, y, z;
//Saddle torus
// Where F(s) = 1 - cos2(s) - cos2(s + 2 PI / 3)
// 0 <= u <= 2 PI, 0 <= v <= 2 PI
auto sign = [](auto a) { if (a == 0) return 0; else if (a < 0) return -1; else return 1; };
auto F = [](auto s) {return 1 - pow(cos(s), 2) - pow(cos(s + 2 * glm::pi<float>() / 3), 2); };
x = _radius * (2 + cos(u)) * cos(v);
y = _radius * (2 + cos(u + 2 * glm::pi<float>() / 3)) * cos(v + 2 * glm::pi<float>() / 3);
z = _radius * (2 + sign(F(u)) * sqrt(abs(F(u)))) * sign(F(v)) * sqrt(abs(F(v)));
P.setParam(x, y, z);
return P;
}