tilings.cpp 12.2 KB
Newer Older
Lionel Untereiner's avatar
Lionel Untereiner committed
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
/*******************************************************************************
* CGoGN: Combinatorial and Geometric modeling with Generic N-dimensional Maps  *
* version 0.1                                                                  *
* Copyright (C) 2009-2012, 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.           *
*                                                                              *
* Web site: http://cgogn.unistra.fr/                                           *
* Contact information: cgogn@unistra.fr                                        *
*                                                                              *
*******************************************************************************/

untereiner's avatar
untereiner committed
25
#include "tilings.h"
Lionel Untereiner's avatar
Lionel Untereiner committed
26 27 28

#include "Geometry/vector_gen.h"

untereiner's avatar
untereiner committed
29 30 31 32 33
#include "Algo/Import/import.h"
#include "Algo/Geometry/boundingbox.h"
#include "Utils/GLSLShader.h"
//#include "Algo/Geometry/area.h"
#include "Algo/Geometry/normal.h"
untereiner's avatar
untereiner committed
34 35
#include "Geometry/matrix.h"
#include "Geometry/transfo.h"
Lionel Untereiner's avatar
Lionel Untereiner committed
36

untereiner's avatar
untereiner committed
37 38
#include "Utils/cgognStream.h"
#include "Utils/chrono.h"
Lionel Untereiner's avatar
Lionel Untereiner committed
39 40 41 42 43 44 45 46 47

using namespace CGoGN ;

/**
 * Struct that contains some informations about the types of the manipulated objects
 * Mainly here to be used by the algorithms that are parameterized by it
 */
struct PFP: public PFP_STANDARD
{
48 49
	// definition of the map
	typedef EmbeddedMap2 MAP;
Lionel Untereiner's avatar
Lionel Untereiner committed
50 51
};

Pierre Kraemer's avatar
Pierre Kraemer committed
52 53 54
typedef PFP::MAP MAP;
typedef PFP::MAP::IMPL MAP_IMPL;
typedef PFP::VEC3 VEC3;
untereiner's avatar
untereiner committed
55 56

PFP::MAP myMap;
Pierre Kraemer's avatar
Pierre Kraemer committed
57 58 59 60

VertexAttribute<VEC3, MAP_IMPL> position;
VertexAttribute<VEC3, MAP_IMPL> position2;
VertexAttribute<VEC3, MAP_IMPL> normal;
untereiner's avatar
untereiner committed
61 62 63

void MyQT::cb_initGL()
{
64
	Utils::GLSLShader::setCurrentOGLVersion(2);
untereiner's avatar
untereiner committed
65

66 67
	// create the render
	m_render = new Algo::Render::GL2::MapRender();
untereiner's avatar
untereiner committed
68

69 70 71
	// create VBO for position
	m_positionVBO = new Utils::VBO();
	m_positionVBO->updateData(position);
untereiner's avatar
untereiner committed
72

73
	m_normalVBO =  new Utils::VBO();
untereiner's avatar
untereiner committed
74

75 76 77
	m_shader = new Utils::ShaderSimpleColor();
	m_shader->setAttributePosition(m_positionVBO);
	m_shader->setColor(Geom::Vec4f(1.,1.,0.,0.));
untereiner's avatar
untereiner committed
78

79 80 81 82 83
	m_lines = new Utils::ShaderVectorPerVertex();
	m_lines->setAttributePosition(m_positionVBO);
	m_lines->setAttributeVector(m_normalVBO);
	m_lines->setScale(0.2f);
	m_lines->setColor(Geom::Vec4f(0.0f, 1.0f, 0.2f, 0.0f));
untereiner's avatar
untereiner committed
84

85
	Algo::Surface::Geometry::computeNormalVertices<PFP>(myMap, position, normal) ;
untereiner's avatar
untereiner committed
86

87 88
	m_render->initPrimitives<PFP>(myMap, Algo::Render::GL2::LINES);
	m_render->initPrimitives<PFP>(myMap, Algo::Render::GL2::POINTS);
untereiner's avatar
untereiner committed
89

90 91
	registerShader(m_shader);
	registerShader(m_lines);
untereiner's avatar
untereiner committed
92

93
	m_normalVBO->updateData(normal);
untereiner's avatar
untereiner committed
94 95 96 97
}

void MyQT::cb_redraw()
{
98 99 100 101 102 103 104 105 106 107 108 109
	m_render->draw(m_shader, Algo::Render::GL2::LINES);
	m_render->draw(m_lines, Algo::Render::GL2::POINTS);
}

void MyQT::squareTiling(int code)
{
	switch(code)
	{
		case 1:
		{
			std::cout << "square grid tiling" << std::endl;
			Algo::Surface::Tilings::Square::Grid<PFP> g(myMap,10,10,true);
untereiner's avatar
untereiner committed
110 111 112 113 114 115
			g.embedIntoGrid(position,1.0,1.0, 0.0);

			Geom::Matrix44f trf;
			trf.identity();
			Geom::translate<float>(0.5,0.5,0.0,trf);
			g.transform(position, trf);
116

untereiner's avatar
untereiner committed
117 118
			g.exportPositions(position, "grid.bs");

119 120 121 122 123 124
			break;
		}
		case 2:
		{
			std::cout << "square grid twisted strip tiling" << std::endl;
			Algo::Surface::Tilings::Square::Grid<PFP> g(myMap,10,10,true);
untereiner's avatar
untereiner committed
125 126 127 128 129 130 131
			g.embedIntoTwistedStrip(position, 0.2, 0.7, 5);

			Geom::Matrix44f trf;
			trf.identity();
			Geom::translate<float>(0.5,0.5,0.0,trf);
			g.transform(position, trf);

132

untereiner's avatar
untereiner committed
133 134
			g.exportPositions(position, "gridtwisted.bs");

135 136 137 138
			break;
		}
		case 3:
		{
untereiner's avatar
untereiner committed
139
			std::cout << "square grid helicoid tiling" << std::endl;
140
			Algo::Surface::Tilings::Square::Grid<PFP> g(myMap,20,20,true);
untereiner's avatar
untereiner committed
141 142 143 144 145 146 147
			g.embedIntoHelicoid(position, 0.2,  0.3, 0.8, 2.0);

//			Geom::Matrix44f trf;
//			trf.identity();
//			Geom::translate<float>(0.5,0.5,0.0,trf);
//			g.transform(position, trf);

148

untereiner's avatar
untereiner committed
149 150
			g.exportPositions(position, "gridhelicoid.bs");

151 152 153 154 155 156
			break;
		}
		case 4:
		{
			std::cout << "square cylinder tiling" << std::endl;
			Algo::Surface::Tilings::Square::Cylinder<PFP> c(myMap,20,20);
untereiner's avatar
untereiner committed
157 158 159 160 161 162 163
			c.embedIntoCylinder(position,0.5,0.7,1.0);

			Geom::Matrix44f trf;
			trf.identity();
			Geom::translate<float>(0.5,0.5,0.0,trf);
			c.transform(position, trf);

164

untereiner's avatar
untereiner committed
165 166
			c.exportPositions(position, "cylinder.bs");

167 168 169 170 171 172 173 174
			break;
		}
		case 5:
		{
			std::cout << "square cylinder tiling" << std::endl;
			Algo::Surface::Tilings::Square::Cylinder<PFP> c(myMap,20,20);
			c.triangleBottom();
			c.triangleTop();
untereiner's avatar
untereiner committed
175
			c.embedIntoCylinder(position,0.5,0.7,1.0);
176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194

			break;
		}
		case 6:
		{
			std::cout << "square cylinder sphere tiling" << std::endl;
			Algo::Surface::Tilings::Square::Cylinder<PFP> c(myMap,20,20);
			c.triangleTop();
			c.triangleBottom();
			c.embedIntoSphere(position,0.5);

			break;
		}
		case 7:
		{
			std::cout << "square cylinder cone tiling" << std::endl;
			Algo::Surface::Tilings::Square::Cylinder<PFP> c(myMap,20,20);
			c.triangleTop();
			c.triangleBottom();
untereiner's avatar
untereiner committed
195
			c.embedIntoCone(position,0.5, 1.0);
196 197 198 199 200

			break;
		}
		case 8:
		{
Lionel Untereiner's avatar
Lionel Untereiner committed
201 202
            std::cout << "square cube tiling" << std::endl;
            Algo::Surface::Tilings::Square::Cube<PFP> c(myMap,1,1,2);
untereiner's avatar
untereiner committed
203 204 205 206 207 208
			c.embedIntoCube(position,1.0, 1.0, 1.0);

			Geom::Matrix44f trf;
			trf.identity();
			Geom::translate<float>(0.5,0.5,0.5,trf);
			c.transform(position, trf);
209

untereiner's avatar
untereiner committed
210 211
			c.exportPositions(position, "cube.bs");

212 213 214 215
			break;
		}
		case 9:
		{
untereiner's avatar
untereiner committed
216
			std::cout << "square tore tiling" << std::endl;
217
			Algo::Surface::Tilings::Square::Tore<PFP> c(myMap,20,10);
untereiner's avatar
untereiner committed
218
			c.embedIntoTore(position, 0.9, 0.5);
219

untereiner's avatar
untereiner committed
220
			c.exportPositions(position, "tore.bs");
221 222 223 224 225 226 227 228
			break;
		}

		default:
		{
			break;
		}
	}
untereiner's avatar
untereiner committed
229 230

	myMap.check();
untereiner's avatar
untereiner committed
231 232 233
}

// mouse picking
234
void MyQT::triangularTiling(int code)
Lionel Untereiner's avatar
Lionel Untereiner committed
235
{
236 237 238 239
	switch(code)
	{
		case 1:
		{
untereiner's avatar
untereiner committed
240
			std::cout << "triangle grid tiling" << std::endl;
241
            Algo::Surface::Tilings::Triangular::Grid<PFP> g(myMap,1,5,true);
untereiner's avatar
untereiner committed
242
			g.embedIntoGrid(position,1,1);
243 244 245 246 247

			break;
		}
		case 2:
		{
untereiner's avatar
untereiner committed
248 249 250
			std::cout << "triangle grid twisted strip tiling" << std::endl;
			Algo::Surface::Tilings::Triangular::Grid<PFP> g(myMap,10,10,true);
			g.embedIntoTwistedStrip(position, 0.3, 0.8, 5);
251 252 253 254 255

			break;
		}
		case 3:
		{
untereiner's avatar
untereiner committed
256 257 258
			std::cout << "triangle grid helocoid tiling" << std::endl;
			Algo::Surface::Tilings::Triangular::Grid<PFP> g(myMap,20,20,true);
			g.embedIntoHelicoid(position, 0.3,  0.8, 5.0, 2.0);
259 260 261 262 263

			break;
		}
		case 4:
		{
untereiner's avatar
untereiner committed
264 265 266
			std::cout << "triangle cylinder tiling" << std::endl;
			Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
			c.embedIntoCylinder(position,0.5,0.7,5.0);
267 268 269 270 271

			break;
		}
		case 5:
		{
untereiner's avatar
untereiner committed
272 273 274 275 276
			std::cout << "triangle cylinder tiling" << std::endl;
			Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
			c.triangleBottom();
			c.triangleTop();
			c.embedIntoCylinder(position,0.5,0.7,5.0);
277 278 279 280 281

			break;
		}
		case 6:
		{
untereiner's avatar
untereiner committed
282 283 284 285 286
			std::cout << "triangle cylinder sphere tiling" << std::endl;
			Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
			c.triangleTop();
			c.triangleBottom();
			c.embedIntoSphere(position,0.5);
287 288 289 290 291

			break;
		}
		case 7:
		{
untereiner's avatar
untereiner committed
292 293 294 295 296
			std::cout << "triangle cylinder cone tiling" << std::endl;
			Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
			c.triangleTop();
			c.triangleBottom();
			c.embedIntoCone(position,0.5, 5.0);
297 298 299 300 301

			break;
		}
		case 8:
		{
Lionel Untereiner's avatar
Lionel Untereiner committed
302 303 304
            std::cout << "triangle cube tiling" << std::endl;
            Algo::Surface::Tilings::Triangular::Cube<PFP> c(myMap,4,4,4);
            c.embedIntoCube(position,5.0,5.0, 5.0);
305 306 307 308 309

			break;
		}
		case 9:
		{
Lionel Untereiner's avatar
Lionel Untereiner committed
310 311 312
            std::cout << "triangle tore tiling" << std::endl;
            Algo::Surface::Tilings::Triangular::Tore<PFP> c(myMap,20,10);
            c.embedIntoTore(position,5.0,2.0);
313 314 315 316 317 318 319 320 321

			break;
		}

		default:
		{
			break;
		}
	}
untereiner's avatar
untereiner committed
322
}
Lionel Untereiner's avatar
Lionel Untereiner committed
323

324 325 326 327 328 329 330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354 355 356 357 358 359 360 361 362 363 364 365 366 367 368 369 370 371 372 373 374 375 376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409 410 411 412 413

void MyQT::hexagonalTiling(int code)
{
    switch(code)
    {
        case 1:
        {
            std::cout << "hexagonal grid tiling" << std::endl;
            Algo::Surface::Tilings::Hexagonal::Grid<PFP> g(myMap,10,6,true);
            g.embedIntoGrid(position,1,1);

            break;
        }
        case 2:
        {
//            std::cout << "hexagonal grid twisted strip tiling" << std::endl;
//            Algo::Surface::Tilings::Triangular::Grid<PFP> g(myMap,10,10,true);
//            g.embedIntoTwistedStrip(position, 0.3, 0.8, 5);

            break;
        }
        case 3:
        {
//            std::cout << "hexagonal grid helocoid tiling" << std::endl;
//            Algo::Surface::Tilings::Triangular::Grid<PFP> g(myMap,20,20,true);
//            g.embedIntoHelicoid(position, 0.3,  0.8, 5.0, 2.0);

            break;
        }
        case 4:
        {
//            std::cout << "hexagonal cylinder tiling" << std::endl;
//            Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
//            c.embedIntoCylinder(position,0.5,0.7,5.0);

            break;
        }
        case 5:
        {
//            std::cout << "hexagonal cylinder tiling" << std::endl;
//            Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
//            c.triangleBottom();
//            c.triangleTop();
//            c.embedIntoCylinder(position,0.5,0.7,5.0);

            break;
        }
        case 6:
        {
//			std::cout << "hexagonal cylinder sphere tiling" << std::endl;
//			Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
//			c.triangleTop();
//			c.triangleBottom();
//			c.embedIntoSphere(position,0.5);

            break;
        }
        case 7:
        {
//			std::cout << "hexagonal cylinder cone tiling" << std::endl;
//			Algo::Surface::Tilings::Triangular::Cylinder<PFP> c(myMap,20,20);
//			c.triangleTop();
//			c.triangleBottom();
//			c.embedIntoCone(position,0.5, 5.0);

            break;
        }
        case 8:
        {
//            std::cout << "hexagonal cube tiling" << std::endl;
//            Algo::Surface::Tilings::Triangular::Cube<PFP> c(myMap,4,4,4);
//            c.embedIntoCube(position,5.0,5.0, 5.0);

            break;
        }
        case 9:
        {
//            std::cout << "hexagonal tore tiling" << std::endl;
//            Algo::Surface::Tilings::Triangular::Tore<PFP> c(myMap,20,10);
//            c.embedIntoTore(position,5.0,2.0);

            break;
        }

        default:
        {
            break;
        }
    }
}
untereiner's avatar
untereiner committed
414 415
int main(int argc, char **argv)
{
416 417 418
	// interface:
	QApplication app(argc, argv);
	MyQT sqt;
untereiner's avatar
untereiner committed
419

420 421
	sqt.statusMsg("Neww to create a sphere or Load for a mesh file");
	CGoGNStream::allToConsole(&sqt);
Lionel Untereiner's avatar
Lionel Untereiner committed
422

423 424
	if (!position.isValid())
		position = myMap.addAttribute<PFP::VEC3, VERTEX>("position");
Lionel Untereiner's avatar
Lionel Untereiner committed
425

426 427
	if (!normal.isValid())
		normal = myMap.addAttribute<PFP::VEC3, VERTEX>("normal");
untereiner's avatar
untereiner committed
428

429 430 431 432 433 434
	if(argc == 3)
	{
		if(argv[1][0] == 'S')
			sqt.squareTiling(atoi(argv[2]));
		else if(argv[1][0] == 'T')
			sqt.triangularTiling(atoi(argv[2]));
435 436 437 438
        else if(argv[1][0] == 'H')
            sqt.hexagonalTiling(atoi(argv[2]));
        else
            return -1;
439
	}
440
	myMap.check();
Lionel Untereiner's avatar
Lionel Untereiner committed
441

442 443 444 445 446 447 448
	//  bounding box
	Geom::BoundingBox<PFP::VEC3> bb = Algo::Geometry::computeBoundingBox<PFP>(myMap, position);
	float lWidthObj = std::max<PFP::REAL>(std::max<PFP::REAL>(bb.size(0), bb.size(1)), bb.size(2));
	Geom::Vec3f lPosObj = (bb.min() +  bb.max()) / PFP::REAL(2);
	CGoGNout << "lPosObj=" << lPosObj << CGoGNendl;
	CGoGNout << "lWidthObj=" << lWidthObj << CGoGNendl;
	sqt.setParamObject(lWidthObj,lPosObj.data());
untereiner's avatar
untereiner committed
449

450 451
	//	myMap.enableQuickTraversal<EDGE>() ;
	//	myMap.enableQuickTraversal<VERTEX>() ;
untereiner's avatar
untereiner committed
452

453
	sqt.show();
untereiner's avatar
untereiner committed
454

455
	return app.exec();
Lionel Untereiner's avatar
Lionel Untereiner committed
456
}