simpleGMap3.cpp 5.62 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 33

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

37 38
	CellMarker mE(myMap,EDGE);

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 62 63
//	for(Dart d = myMap.begin() ; d != myMap.end() ; myMap.next(d))
//		mE.mark(d);

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

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

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

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

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

85 86 87 88 89 90
	myMap.check();
//
	myMap.splitFace(d,myMap.phi1(myMap.phi1(myMap.phi1(d))));

	myMap.check();

91 92
	for(unsigned int i = position.begin() ; i != position.end() ; position.next(i))
		position[i] += VEC3(0,2,0);
93

94 95 96 97 98 99 100 101
	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();
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
}

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);
128
//	Algo::Render::GL1::renderTopoMD3<PFP>(myMap, position, true, true, true, 0.9f, 0.9f, 0.9f);
129 130 131 132 133 134

	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);
135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154
}


/**********************************************************************************************
 *                                      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() ;
}