simpleGMap3.cpp 5.53 KB
Newer Older
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
/*******************************************************************************
* CGoGN: Combinatorial and Geometric modeling with Generic N-dimensional Maps  *
* version 0.1                                                                  *
* Copyright (C) 2009, IGG Team, LSIIT, University of Strasbourg                *
*                                                                              *
* This library is free software; you can redistribute it and/or modify it      *
* under the terms of the GNU Lesser General Public License as published by the *
* Free Software Foundation; either version 2.1 of the License, or (at your     *
* option) any later version.                                                   *
*                                                                              *
* This library is distributed in the hope that it will be useful, but WITHOUT  *
* ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or        *
* FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License  *
* for more details.                                                            *
*                                                                              *
* You should have received a copy of the GNU Lesser General Public License     *
* along with this library; if not, write to the Free Software Foundation,      *
* Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301 USA.           *
*                                                                              *
20
* Web site: http://cgogn.unistra.fr/                                  *
21
22
23
24
25
26
27
28
* Contact information: cgogn@unistra.fr                                        *
*                                                                              *
*******************************************************************************/

#include "simpleGMap3.h"
#include "Utils/GLSLShader.h"
#include "Algo/Geometry/boundingbox.h"
#include "Algo/Modelisation/primitives3d.h"
untereiner's avatar
untereiner committed
29
#include "Algo/Modelisation/subdivision3.h"
30
31
32

SimpleGMap3::SimpleGMap3()
{
33
34
35
	position = myMap.addAttribute<VEC3, VERTEX>("position");
	normal = myMap.addAttribute<VEC3, VERTEX>("normal");
	volume = myMap.addAttribute<VEC3, VOLUME>("volume");
36

Thomas's avatar
Thomas committed
37
	CellMarker<EDGE> mE(myMap);
38

Thomas's avatar
Thomas committed
39
	Algo::Modelisation::Primitive3D<PFP> primCat(myMap,position);
40
41
	Dart d = primCat.hexaGrid_topo(3,1,1);
	primCat.embedHexaGrid(2,1,1);
Thomas's avatar
Thomas committed
42
	myMap.check();
43
44
45
//	DartMarker markOrient(myMap);
//	std::vector<Dart> orient;
//	FunctorStore fs(orient);
Pierre Kraemer's avatar
Pierre Kraemer committed
46
	
47
48
//	d = 49;
//	myMap.foreach_dart_of_oriented_volume(d, fs);
49

50
51
//	for(std::vector<Dart>::iterator it = orient.begin() ; it != orient.end() ; ++it)
//		markOrient.mark(*it);
52

53
54
55
//	SelectorMarked sm(markOrient);
//	std::cout << "AAA"<< std::endl;
//	Algo::Modelisation::catmullClarkVol<PFP,PFP::TVEC3,PFP::VEC3>(myMap, position, sm);
56

57
58
59
//	Geom::Plane3D<PFP::REAL> pl(VEC3(0.5,0.14,0.5),VEC3(1.5,0.45,0.5),VEC3(0.5,0.15,1.5));
	Geom::Plane3D<PFP::REAL> pl(VEC3(-1,-0.5,-0.5),VEC3(-1,-0.5,0.5),VEC3(1,0.5,0.5));
	Algo::Modelisation::sliceConvexVolume<PFP>(myMap, position, d, pl);
60

61
	myMap.check();
62

63
64
	for(unsigned int i = position.begin() ; i != position.end() ; position.next(i))
		position[i] += VEC3(2,0,0);
65

66
	Algo::Modelisation::Primitive3D<PFP> prim(myMap, position);
67
68
	d = prim.hexaGrid_topo(2,2,1);
	prim.embedHexaGrid(1,1,1);
Thomas's avatar
Thomas committed
69

70
	Dart d1 = myMap.phi1(myMap.phi1(myMap.phi2(myMap.phi1(myMap.phi1(d)))));
71
	VEC3 mid0 = (position[d1] + position[myMap.phi1(d1)]) / 2.0f;
72
73
74
	myMap.cutEdge(d1);
	position[myMap.phi1(d1)] = mid0;

75
	VEC3 mid1 = (position[d] + position[myMap.phi1(d)]) / 2.0f;
76
77
78
79
	myMap.cutEdge(d);
	position[myMap.phi1(d)] = mid1;

	d = myMap.phi1(myMap.phi1(myMap.phi2(myMap.phi1(myMap.phi1(d)))));
80
	VEC3 mid = (position[d] + position[myMap.phi1(d)]) / 2.0f;
81
82
83
	myMap.cutEdge(d);
	position[myMap.phi1(d)] = mid;

84
85
86
87
	myMap.splitFace(d,myMap.phi1(myMap.phi1(myMap.phi1(d))));

	myMap.check();

88
89
	for(unsigned int i = position.begin() ; i != position.end() ; position.next(i))
		position[i] += VEC3(0,2,0);
90

91
92
93
94
95
96
97
98
	Algo::Modelisation::Primitive3D<PFP> prim2(myMap,position);
	d = prim2.hexaGrid_topo(2,1,1);
	prim2.embedHexaGrid(1,1,1);

	d = myMap.phi2(myMap.phi1(myMap.phi1(myMap.phi2(d))));
	myMap.unsewVolumes(d);

	myMap.check();
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
}

void SimpleGMap3::initGUI()
{
}

void SimpleGMap3::cb_initGL()
{
	Utils::GLSLShader::setCurrentOGLVersion(1) ;

	Geom::BoundingBox<PFP::VEC3> bb = Algo::Geometry::computeBoundingBox<PFP>(myMap, position) ;
	VEC3 gPosObj = bb.center() ;
	float tailleX = bb.size(0) ;
	float tailleY = bb.size(1) ;
	float tailleZ = bb.size(2) ;
	float gWidthObj = std::max<float>(std::max<float>(tailleX, tailleY), tailleZ) ;
	setParamObject(gWidthObj, gPosObj.data());
}

void SimpleGMap3::cb_redraw()
{
	glDisable(GL_LIGHTING);
	glLineWidth(1.0f);
	Algo::Render::GL1::renderTopoGMD3<PFP>(myMap, position, true, true, true, true, 0.9f, 0.9f, 0.9f, 0.9f);
123
//	Algo::Render::GL1::renderTopoMD3<PFP>(myMap, position, true, true, true, 0.9f, 0.9f, 0.9f);
124

125
126
127
128
129
//	glDisable(GL_LIGHTING);
//	glColor3f(1.0f, 1.0f, 1.0f);
//	glLineWidth(1.0f);
//	glPolygonMode(GL_FRONT_AND_BACK, GL_LINE);
//	Algo::Render::GL1::renderTriQuadPoly<PFP>(myMap, Algo::Render::GL1::LINE, 1.0,position, normal);
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
}

/**********************************************************************************************
 *                                      MAIN FUNCTION                                         *
 **********************************************************************************************/

int main(int argc, char **argv)
{
	QApplication app(argc, argv) ;

	SimpleGMap3 sqt ;
	sqt.setGeometry(0, 0, 1000, 800) ;
 	sqt.show() ;

	sqt.initGUI() ;

	return app.exec() ;
}