MatlEditor/BreatherSurface.cpp at master · sharjith/MatlEditor · GitHub

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
#include "BreatherSurface.h"
#include "Point.h"

#include <glm/gtc/constants.hpp>
#include <glm/vec3.hpp>
#include <glm/glm.hpp>


BreatherSurface::BreatherSurface(QOpenGLShaderProgram* prog, GLfloat radius, GLuint nSlices, GLuint nStacks) :
	ParametricSurface(prog, nSlices, nStacks),
	_radius(radius)
{
	_name = "Breather Surface";
	buildMesh(nSlices, nStacks);
}


BreatherSurface::~BreatherSurface()
{
}

float BreatherSurface::firstUParameter() const
{
	return -13.2f;
}

float BreatherSurface::lastUParameter() const
{
	return 13.2f;
}

float BreatherSurface::firstVParameter() const
{
	return -37.4f;
}

float BreatherSurface::lastVParameter() const
{
	return 37.4f;
}

Point BreatherSurface::pointAtParameter(const float& u, const float& v)
{
	Point P;
	float x, y, z;

	//http://xahlee.info/surface/breather_p/breather_p.html

	float b = 0.4f;
	float r = 1 - b * b;
	float w = sqrt(r);
	float denom = b * pow((w*cosh(b*u)), 2) + pow((b*sin(w*v)), 2);
	x = _radius * (-u + (2 * r*cosh(b*u) * sinh(b*u)) /denom);
	y = _radius * ((2 * w*cosh(b*u) * (-(w*cos(v) * cos(w*v)) - sin(v) * sin(w*v))) / denom);
	z = _radius * ((2 * w*cosh(b*u) * (-(w*sin(v) * cos(w*v)) + cos(v) * sin(w*v))) / denom);


	P.setParam(x, y, z);
	return P;
}