PrimitivesGenerator.cc 33.2 KB
Newer Older
Jan Möbius's avatar
Jan Möbius committed
1
2
3
/*===========================================================================*\
*                                                                            *
*                              OpenFlipper                                   *
4
*      Copyright (C) 2001-2013 by Computer Graphics Group, RWTH Aachen       *
Jan Möbius's avatar
Jan Möbius committed
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
*                           www.openflipper.org                              *
*                                                                            *
*--------------------------------------------------------------------------- *
*  This file is part of OpenFlipper.                                         *
*                                                                            *
*  OpenFlipper 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 3 of            *
*  the License, or (at your option) any later version with the               *
*  following exceptions:                                                     *
*                                                                            *
*  If other files instantiate templates or use macros                        *
*  or inline functions from this file, or you compile this file and          *
*  link it with other files to produce an executable, this file does         *
*  not by itself cause the resulting executable to be covered by the         *
*  GNU Lesser General Public License. This exception does not however        *
*  invalidate any other reasons why the executable file might be             *
*  covered by the GNU Lesser General Public License.                         *
*                                                                            *
*  OpenFlipper 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 LesserGeneral Public           *
*  License along with OpenFlipper. If not,                                   *
*  see <http://www.gnu.org/licenses/>.                                       *
*                                                                            *
\*===========================================================================*/

/*===========================================================================*\
*                                                                            *
*   $Revision$                                                       *
*   $LastChangedBy$                                                *
*   $Date$                     *
*                                                                            *
\*===========================================================================*/

#include "PrimitivesGenerator.hh"

45
46
47
#ifdef ENABLE_OPENVOLUMEMESH_POLYHEDRAL_SUPPORT
#include "TetrahedralCuboidGenerator.hh"
#endif
Jan Möbius's avatar
Jan Möbius committed
48
49

#include <OpenFlipper/BasePlugin/PluginFunctions.hh>
50
#include <ACG/Geometry/Algorithms.hh>
Jan Möbius's avatar
Jan Möbius committed
51

Jan Möbius's avatar
Jan Möbius committed
52
53
PrimitivesGeneratorPlugin::PrimitivesGeneratorPlugin() :
        triMesh_(0),
54
55
        polyMesh_(0),
        slices_(50),
Jan Möbius's avatar
Jan Möbius committed
56
57
        stacks_(50),
        primitivesMenu_(0)
Jan Möbius's avatar
Jan Möbius committed
58
59
60
61
62
63
64
{

}

PrimitivesGeneratorPlugin::~PrimitivesGeneratorPlugin()
{

Jan Möbius's avatar
Jan Möbius committed
65
66
67
68
  if ( OpenFlipper::Options::gui()) {
    delete primitivesMenu_;
  }

Jan Möbius's avatar
Jan Möbius committed
69
70
}

Jan Möbius's avatar
Jan Möbius committed
71
72
void PrimitivesGeneratorPlugin::initializePlugin()
{
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
  emit setSlotDescription("addTetrahedron(Vector, double)",
                          tr("Generates a tetrahedron (ObjectId is returned)"),
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

  emit setSlotDescription("addIcosahedron(Vector, double)",
                          tr("Generates an icosahedron (ObjectId is returned)"),
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

  emit setSlotDescription("addPyramid(Vector, double)",
                          tr("Generates a pyramid (ObjectId is returned)"),
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

  emit setSlotDescription("addOctahedron(Vector, double)",
                          tr("Generates an octahedron (ObjectId is returned)"),
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

  emit setSlotDescription("addDodecahedron(Vector, double)",
                          tr("Generates a dodecahedron (ObjectId is returned)"),
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

  emit setSlotDescription("addSphere(Vector, double)",
                          tr("Generates a sphere (ObjectId is returned)"),
                          QString("Position, Radius").split(","),
                          QString("Center position,Radius").split(","));

  emit setSlotDescription("addTriangulatedCube(Vector, double)",
104
                          tr("Generates a triangular mesh of cube (ObjectId is returned)"),
105
106
107
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

108
109
  emit setSlotDescription("addTriangulatedCylinder(Vector,Vector,double,double)",
                          tr("Generates a triangulated cylinder (ObjectId is returned)")  ,
110
                          QString("Position,Axis,Radius,Height,Top,Bottom").split(","),
111
                          QString("Bottom center vertex position,Center axis,radius,height,add top vertex,add bottom vertex").split(","));
112
113
114
115
116
117
118
119
120
121
122
123
124

#ifdef ENABLE_OPENVOLUMEMESH_POLYHEDRAL_SUPPORT
  emit setSlotDescription("addTetrahedralCube(Vector,double)",
                          tr("Generates a tetrahedral mesh of a cube (ObjectId is returned)"),
                          QString("Position,Length").split(","),
                          QString("Center position,Length of each edge").split(","));

  emit setSlotDescription("addTetrahedralCuboid(Vector,Vector,uint,uint,uint)",
                          tr("Generates a tetrahedral mesh of a cuboid (ObjectId is returned)"),
                          QString("Position,Lengths,Count,Count,Count").split(","),
                          QString("Center position,Length of each side,Number of units in x-axis,Number of units in y-axis,Number of units in z-axis").split(","));
#endif

Jan Möbius's avatar
Jan Möbius committed
125
126
127
128
}

void PrimitivesGeneratorPlugin::pluginsInitialized() {

Jan Möbius's avatar
Jan Möbius committed
129
  if ( OpenFlipper::Options::gui()) {
Jan Möbius's avatar
Jan Möbius committed
130

Jan Möbius's avatar
Jan Möbius committed
131
    emit getMenubarMenu(tr("&Primitives"), primitivesMenu_, true );
Jan Möbius's avatar
Jan Möbius committed
132

Matthias Möller's avatar
Matthias Möller committed
133
    QAction* action;
Jan Möbius's avatar
Jan Möbius committed
134

Jan Möbius's avatar
Jan Möbius committed
135
136
    action = primitivesMenu_->addAction("Cube (Triangle Mesh)"       ,this,SLOT(addTriangulatedCube()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_cube.png"));
137

138
139
    WhatsThisGenerator whatsThisGen("PrimitivesGenerator");
    whatsThisGen.setWhatsThis(action,tr("Create a Cube."),"Cube");
Jan Möbius's avatar
Jan Möbius committed
140

Jan Möbius's avatar
Jan Möbius committed
141
142
    action = primitivesMenu_->addAction("Dodecahedron"               ,this,SLOT(addDodecahedron()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_dodecahedron.png"));
143
    whatsThisGen.setWhatsThis(action,tr("Create a Dodecahedron."), "Dodecahedron");
Jan Möbius's avatar
Jan Möbius committed
144

Jan Möbius's avatar
Jan Möbius committed
145
146
    action = primitivesMenu_->addAction("Icosahedron"                ,this,SLOT(addIcosahedron()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_icosahedron.png"));
147
    whatsThisGen.setWhatsThis(action,tr("Create a Icosahedron.","Icosahedron"));
Jan Möbius's avatar
Jan Möbius committed
148

Jan Möbius's avatar
Jan Möbius committed
149
150
    action = primitivesMenu_->addAction("Octahedron"                  ,this,SLOT(addOctahedron()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_octahedron.png"));
151
    whatsThisGen.setWhatsThis(action,tr("Create an Octahedron."),"Octahedron");
Jan Möbius's avatar
Jan Möbius committed
152

Jan Möbius's avatar
Jan Möbius committed
153
154
    action = primitivesMenu_->addAction("Pyramid"                    ,this,SLOT(addPyramid()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_pyramid.png"));
155
    whatsThisGen.setWhatsThis(action,tr("Create a Pyramid."),"Pyramid");
Jan Möbius's avatar
Jan Möbius committed
156

Jan Möbius's avatar
Jan Möbius committed
157
158
    action = primitivesMenu_->addAction("Cylinder (Triangle Mesh)"       ,this,SLOT(addTriangulatedCylinder()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_cylinder.png"));
Jan Möbius's avatar
Jan Möbius committed
159

Jan Möbius's avatar
Jan Möbius committed
160
161
    action = primitivesMenu_->addAction("Sphere",this,SLOT(addSphere()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_sphere.png"));
162
    whatsThisGen.setWhatsThis(action,tr("Create a Sphere. "),"Sphere");
Jan Möbius's avatar
Jan Möbius committed
163

Jan Möbius's avatar
Jan Möbius committed
164
165
    action = primitivesMenu_->addAction("Tetrahedron",this,SLOT(addTetrahedron()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_tetrahedron.png"));
166
    whatsThisGen.setWhatsThis(action,tr("Create a Tetrahedron."),"Tetrahedron");
Jan Möbius's avatar
Jan Möbius committed
167

168
169
170
171
172
173
174
175
176
177
#ifdef ENABLE_OPENVOLUMEMESH_POLYHEDRAL_SUPPORT
    action = primitivesMenu_->addAction("Cube (Tetrahedral Mesh)"    ,this,SLOT(addTetrahedralCube()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_cube.png"));
    whatsThisGen.setWhatsThis(action,tr("Create a Tetrahedral Cube."), "Cube");

    action = primitivesMenu_->addAction("Cuboid (Tetrahedral Mesh)"  ,this,SLOT(addTetrahedralCuboid()));
    action->setIcon(QIcon(OpenFlipper::Options::iconDirStr()+OpenFlipper::Options::dirSeparator()+"primitive_cube.png"));
    whatsThisGen.setWhatsThis(action,tr("Create a Tetrahedral Cuboid."), "Cuboid");
#endif
  }
Jan Möbius's avatar
Jan Möbius committed
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
}


int PrimitivesGeneratorPlugin::addTriMesh() {
  int objectId = -1;

  emit addEmptyObject( DATA_TRIANGLE_MESH, objectId );

  TriMeshObject* object;
  if ( !PluginFunctions::getObject(objectId,object) ) {
    emit log(LOGERR,"Unable to create new Object");
    return -1;
  }

  return objectId;
}

Jan Möbius's avatar
Jan Möbius committed
195
196
197
198
199
200
201
202
203
204
205
206
207
208
int PrimitivesGeneratorPlugin::addPolyMesh() {
  int objectId = -1;

  emit addEmptyObject( DATA_POLY_MESH, objectId );

  PolyMeshObject* object;
  if ( !PluginFunctions::getObject(objectId,object) ) {
    emit log(LOGERR,"Unable to create new Object");
    return -1;
  }

  return objectId;
}

209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225

#ifdef ENABLE_OPENVOLUMEMESH_POLYHEDRAL_SUPPORT
int PrimitivesGeneratorPlugin::addPolyhedralMesh() {
    int objectId = -1;

    emit addEmptyObject( DATA_POLYHEDRAL_MESH, objectId );

    PolyhedralMeshObject* object;
    if (!PluginFunctions::getObject(objectId, object) ) {
        emit log(LOGERR, "Unable to create new PolyhedralMesh object");
        return -1;
    }

    return objectId;
}
#endif

226
int PrimitivesGeneratorPlugin::addTetrahedron(const Vector& _position, const double _length) {
Jan Möbius's avatar
Jan Möbius committed
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244

  int newObject = addTriMesh();

  TriMeshObject* object;
  if ( !PluginFunctions::getObject(newObject,object) ) {
        emit log(LOGERR,"Unable to create new Object");
        return -1;
  } else {

    object->setName( "Tetrahedron " + QString::number(newObject) );

    triMesh_ =  object->mesh();

    triMesh_->clear();

    // Add 4 vertices
    vhandles_.resize(4);

245
    const double halfSize = 0.5*_length;
Jan Möbius's avatar
Jan Möbius committed
246

247
248
249
250
    vhandles_[0] = triMesh_->add_vertex(TriMesh::Point(-halfSize, -halfSize,  halfSize)+_position);
    vhandles_[1] = triMesh_->add_vertex(TriMesh::Point( halfSize,  halfSize,  halfSize)+_position);
    vhandles_[2] = triMesh_->add_vertex(TriMesh::Point(-halfSize,  halfSize, -halfSize)+_position);
    vhandles_[3] = triMesh_->add_vertex(TriMesh::Point( halfSize, -halfSize, -halfSize)+_position);
Jan Möbius's avatar
Jan Möbius committed
251
252
253
254
255
256
257
258
259
260
261

    // Add 4 faces
    add_face(0,1,2);
    add_face(0,2,3);
    add_face(2,1,3);
    add_face(3,1,0);

    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

262
263
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
264
265
266
267
268
269
    return newObject;
  }

  return -1;
}

270
int PrimitivesGeneratorPlugin::addTriangulatedCube(const Vector& _position,const double _length) {
Jan Möbius's avatar
Jan Möbius committed
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289


  int newObject = addTriMesh();

  TriMeshObject* object;
  if ( !PluginFunctions::getObject(newObject,object) ) {
        emit log(LOGERR,"Unable to create new Object");
        return -1;
  } else {

    object->setName( "Cube " + QString::number(newObject) );

    triMesh_ =  object->mesh();

    triMesh_->clear();

    // Add 8 vertices
    vhandles_.resize(8);

290
    const double halfSize = 0.5*_length;
Jan Möbius's avatar
Jan Möbius committed
291

292
293
294
295
296
297
298
299
    vhandles_[0] = triMesh_->add_vertex(TriMesh::Point( halfSize, -halfSize, halfSize)+_position);
    vhandles_[1] = triMesh_->add_vertex(TriMesh::Point( halfSize,  halfSize, halfSize)+_position);
    vhandles_[2] = triMesh_->add_vertex(TriMesh::Point(-halfSize,  halfSize, halfSize)+_position);
    vhandles_[3] = triMesh_->add_vertex(TriMesh::Point(-halfSize, -halfSize, halfSize)+_position);
    vhandles_[4] = triMesh_->add_vertex(TriMesh::Point( halfSize, -halfSize,-halfSize)+_position);
    vhandles_[5] = triMesh_->add_vertex(TriMesh::Point( halfSize,  halfSize,-halfSize)+_position);
    vhandles_[6] = triMesh_->add_vertex(TriMesh::Point(-halfSize,  halfSize,-halfSize)+_position);
    vhandles_[7] = triMesh_->add_vertex(TriMesh::Point(-halfSize, -halfSize,-halfSize)+_position);
Jan Möbius's avatar
Jan Möbius committed
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320


    // Add 12 faces
    add_face(0,1,2);
    add_face(0,2,3);
    add_face(0,5,1);
    add_face(5,0,4);
    add_face(4,0,7);
    add_face(7,0,3);

    add_face(7,3,6);
    add_face(6,3,2);
    add_face(6,2,5);
    add_face(5,2,1);
    add_face(6,5,4);
    add_face(6,4,7);

    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

321
322
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
323
324
325
326
327
328
    return newObject;
  }

  return -1;
}

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
//========================================================================
// Tetrahedral cube
//========================================================================

#ifdef ENABLE_OPENVOLUMEMESH_POLYHEDRAL_SUPPORT

int PrimitivesGeneratorPlugin::addTetrahedralCube(const Vector& _position, const double _length)
{
    return addTetrahedralCuboid(_position, Vector(_length, _length, _length), 1, 1, 1);
}

#endif

//========================================================================
// Tetrahedral Cuboid
//========================================================================

#ifdef ENABLE_OPENVOLUMEMESH_POLYHEDRAL_SUPPORT

int PrimitivesGeneratorPlugin::addTetrahedralCuboid(const Vector& _position,
        const Vector& _length, const unsigned int  n_x, const unsigned int  n_y, const unsigned int  n_z)
{

    int object_id = addPolyhedralMesh();
    PolyhedralMeshObject* object;

    if (!PluginFunctions::getObject(object_id, object)) {
        return -1;
    }

    object->setName("Cuboid " + QString::number(object_id));

    TetrahedralCuboidGenerator gen;
    gen.generate(*(object->mesh()), _position, _length, n_x, n_y, n_z);

    emit updatedObject(object_id, UPDATE_ALL);

    object->setObjectDrawMode(ACG::SceneGraph::DrawModes::getDrawMode("Cells (flat shaded)"));
    PluginFunctions::viewAll();

    return object_id;
}

#endif

Jan Möbius's avatar
Jan Möbius committed
374
375
376
377
378
//========================================================================
// Cylinder
//========================================================================


379
380
381
382
383
384
ACG::Vec3d PrimitivesGeneratorPlugin::positionOnCylinder(const int _sliceNumber,
                                                         const int _stackNumber,
                                                         const Vector _position,
                                                         const Vector _axis,
                                                         const double _radius,
                                                         const double _height)
Jan Möbius's avatar
Jan Möbius committed
385
386
387
{
  ACG::Vec3d position;

388
389
  const ACG::Vec3d right = (ACG::Geometry::perpendicular(_axis)).normalized();
  const ACG::Vec3d left  = (cross( _axis, right)).normalized();
Jan Möbius's avatar
Jan Möbius committed
390
391
392
393
394
395

  double beta = ((2.0 * M_PI) / double(slices_)) * double(_sliceNumber);

  if ( _sliceNumber == 0 &&  _stackNumber == 0) {
    position[0] = 0.0;
    position[1] = 0.0;
396
    position[2] = _height;
Jan Möbius's avatar
Jan Möbius committed
397
398
399
400
401
  } else  if ( _sliceNumber == slices_ &&  _stackNumber == stacks_ ) {
    position[0] = 0.0;
    position[1] = 0.0;
    position[2] = 0.0;
  } else {
402
403
404
    position[0] = sin(beta) * _radius;
    position[1] = cos(beta) * _radius;
    position[2] = _height * double(stacks_ - _stackNumber -1 ) / double(stacks_-2);
Jan Möbius's avatar
Jan Möbius committed
405
406
  }

407
  position = _position + position[0] * right + position[1] * left + position[2] * _axis ;
408

Jan Möbius's avatar
Jan Möbius committed
409
410
411
  return position;
}

412
int PrimitivesGeneratorPlugin::addTriangulatedCylinder(const Vector& _position,const Vector& _axis,const double _radius,const double _height,const bool _top,const bool   _bottom ) {
Jan Möbius's avatar
Jan Möbius committed
413

Jan Möbius's avatar
Jan Möbius committed
414
  // TODO: Generate texture coordinates for cylinder (Glu compatible)
Jan Möbius's avatar
Jan Möbius committed
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
  int newObject = addTriMesh();

  TriMeshObject* object;
  if (!PluginFunctions::getObject(newObject, object)) {
    emit log(LOGERR, "Unable to create new Object");
    return -1;
  } else {

    object->setName( "Cylinder " + QString::number(newObject) );

    triMesh_ = object->mesh();

    triMesh_->clear();

    //triMesh_->request_vertex_texcoords2D();

    TriMesh::VertexHandle vh;
432
    TriMesh::VertexHandle top = triMesh_->add_vertex(positionOnCylinder(0, 0,_position,_axis,_radius,_height));
Jan Möbius's avatar
Jan Möbius committed
433
434
435
436
    //triMesh_->set_texcoord2D(vh, texCoordOnSphere(0, 0));

    for (int st = 1; st < stacks_; ++st) {
      for (int sl = 0; sl < slices_; ++sl) {
437
        vh = triMesh_->add_vertex(positionOnCylinder(sl, st,_position,_axis,_radius,_height));
Jan Möbius's avatar
Jan Möbius committed
438
439
440
441
        //triMesh_->set_texcoord2D(vh, texCoordOnSphere(sl, st));
      }
    }

442
    TriMesh::VertexHandle bottom = triMesh_->add_vertex(positionOnCylinder(slices_, stacks_,_position,_axis,_radius,_height));
Jan Möbius's avatar
Jan Möbius committed
443
444
445
446
447
    //triMesh_->set_texcoord2D(vh, texCoordOnSphere(slices_, stacks_));

    std::vector<TriMesh::VertexHandle> vhandles;

    // Add top triangle fan ( Vertex index is shifted by one for the first slice )
448
449
    if ( _top ) {
      for (int sl = 1; sl < slices_ + 1; ++sl) {
Jan Möbius's avatar
Jan Möbius committed
450

451
        vhandles.clear();
Jan Möbius's avatar
Jan Möbius committed
452

453
454
455
        vhandles.push_back(triMesh_->vertex_handle(sl));
        vhandles.push_back(triMesh_->vertex_handle(0));
        vhandles.push_back(triMesh_->vertex_handle(1 * 1 + (sl % slices_)));
Jan Möbius's avatar
Jan Möbius committed
456

457
458
459
460
        triMesh_->add_face(vhandles);
      }
    } else {
      triMesh_->delete_vertex(top);
Jan Möbius's avatar
Jan Möbius committed
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
    }

    for (int st = 0; st < stacks_ - 2; ++st) {

      // Move around one slice
      for (int sl = 0; sl < slices_; ++sl) {

        // Offset 1 because of singular vertex
        unsigned int startTop = 1 + slices_ * st;
        unsigned int startBottom = 1 + slices_ * (st + 1);

        vhandles.clear();

        vhandles.push_back(triMesh_->vertex_handle(startTop + sl));
        vhandles.push_back(triMesh_->vertex_handle(startTop + ((sl + 1) % slices_)));
        vhandles.push_back(triMesh_->vertex_handle(startBottom + sl));

        triMesh_->add_face(vhandles);

        vhandles.clear();

        vhandles.push_back(triMesh_->vertex_handle(startBottom + sl));
        vhandles.push_back(triMesh_->vertex_handle(startTop + ((sl + 1) % slices_)));
        vhandles.push_back(triMesh_->vertex_handle(startBottom + ((sl + 1) % slices_)));

        triMesh_->add_face(vhandles);
      }

    }

    const int startTop     = 1 + (stacks_ - 2) * slices_;
    const int bottomVertex = 1 + (stacks_ - 1) * slices_;

    // Add bottom triangle fan
495
496
    if ( _bottom) {
      for (int sl = 0; sl < slices_; ++sl) {
Jan Möbius's avatar
Jan Möbius committed
497

498
        vhandles.clear();
Jan Möbius's avatar
Jan Möbius committed
499

500
501
502
        vhandles.push_back(triMesh_->vertex_handle(bottomVertex));
        vhandles.push_back(triMesh_->vertex_handle(startTop + sl));
        vhandles.push_back(triMesh_->vertex_handle(startTop + ((sl + 1) % slices_)));
Jan Möbius's avatar
Jan Möbius committed
503

504
505
506
507
        triMesh_->add_face(vhandles);
      }
    } else {
      triMesh_->delete_vertex(bottom);
Jan Möbius's avatar
Jan Möbius committed
508
509
    }

510
511
512
    // Cleanup if bottom or top vertex is missing
    triMesh_->garbage_collection();

Jan Möbius's avatar
Jan Möbius committed
513
514
515
516
    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

517
518
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
519
520
521
522
523
524
    return object->id();
  }

}


525
526
527
528
529
//========================================================================
// Sphere
//========================================================================


530
ACG::Vec3d PrimitivesGeneratorPlugin::positionOnSphere(int _sliceNumber, int _stackNumber, double _radius, const Vector& _position)
531
532
533
534
535
536
537
{
  ACG::Vec3d position;

  double alpha = (M_PI / double(stacks_)) * double(_stackNumber);
  double beta = ((2.0 * M_PI) / double(slices_)) * double(_sliceNumber);

  double ringRadius = sin(alpha);
538
539
540
  position[0] = sin(beta) * ringRadius * _radius;
  position[1] = cos(beta) * ringRadius * _radius;
  position[2] = cos(alpha)* _radius;
541

542
  return _position+position;
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
}

//------------------------------------------------------------------------

ACG::Vec2f PrimitivesGeneratorPlugin::texCoordOnSphere(int _sliceNumber, int _stackNumber)
{
  ACG::Vec2f texCoord;

  double alpha = (M_PI / double(stacks_)) * double(_stackNumber);
  texCoord[0] = double(_sliceNumber) / double(slices_);
  texCoord[1] = 0.5 * (cos(alpha) + 1.0);

  return texCoord;
}

Jan Möbius's avatar
Jan Möbius committed
558

559
560
//------------------------------------------------------------------------

561
562
int PrimitivesGeneratorPlugin::addSphere(const Vector& _position, const double _radius)
{
563
564
565
566
567
568
569
570
  int newObject = addTriMesh();

  TriMeshObject* object;
  if (!PluginFunctions::getObject(newObject, object)) {
    emit log(LOGERR, "Unable to create new Object");
    return -1;
  } else {

Jan Möbius's avatar
Jan Möbius committed
571
572
    object->setName( "Sphere " + QString::number(newObject) );

573
    triMesh_ = object->mesh();
Jan Möbius's avatar
Jan Möbius committed
574
575
576

    triMesh_->clear();

577
578
579
580
    triMesh_->request_vertex_texcoords2D();

    TriMesh::VertexHandle vh;

581
    vh = triMesh_->add_vertex(positionOnSphere(0, 0, _radius,_position));
582
583
    triMesh_->set_texcoord2D(vh, texCoordOnSphere(0, 0));

584
    for (int st = 1; st < stacks_; ++st) {
585
      for (int sl = 0; sl < slices_; ++sl) {
586
        vh = triMesh_->add_vertex(positionOnSphere(sl, st, _radius,_position));
587
588
589
590
        triMesh_->set_texcoord2D(vh, texCoordOnSphere(sl, st));
      }
    }

591
    vh = triMesh_->add_vertex(positionOnSphere(slices_, stacks_, _radius,_position));
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
    triMesh_->set_texcoord2D(vh, texCoordOnSphere(slices_, stacks_));

    std::vector<TriMesh::VertexHandle> vhandles;

    // Add top triangle fan ( Vertex index is shifted by one for the first slice )
    for (int sl = 1; sl < slices_ + 1; ++sl) {

      vhandles.clear();

      vhandles.push_back(triMesh_->vertex_handle(sl));
      vhandles.push_back(triMesh_->vertex_handle(0));
      vhandles.push_back(triMesh_->vertex_handle(1 * 1 + (sl % slices_)));

      triMesh_->add_face(vhandles);
    }

608
    for (int st = 0; st < stacks_ - 2; ++st) {
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635

      // Move around one slice
      for (int sl = 0; sl < slices_; ++sl) {

        // Offset 1 because of singular vertex
        unsigned int startTop = 1 + slices_ * st;
        unsigned int startBottom = 1 + slices_ * (st + 1);

        vhandles.clear();

        vhandles.push_back(triMesh_->vertex_handle(startTop + sl));
        vhandles.push_back(triMesh_->vertex_handle(startTop + ((sl + 1) % slices_)));
        vhandles.push_back(triMesh_->vertex_handle(startBottom + sl));

        triMesh_->add_face(vhandles);

        vhandles.clear();

        vhandles.push_back(triMesh_->vertex_handle(startBottom + sl));
        vhandles.push_back(triMesh_->vertex_handle(startTop + ((sl + 1) % slices_)));
        vhandles.push_back(triMesh_->vertex_handle(startBottom + ((sl + 1) % slices_)));

        triMesh_->add_face(vhandles);
      }

    }

636
637
    const int startTop     = 1 + (stacks_ - 2) * slices_;
    const int bottomVertex = 1 + (stacks_ - 1) * slices_;
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654

    // Add bottom triangle fan
    for (int sl = 0; sl < slices_; ++sl) {

      vhandles.clear();

      vhandles.push_back(triMesh_->vertex_handle(bottomVertex));
      vhandles.push_back(triMesh_->vertex_handle(startTop + sl));
      vhandles.push_back(triMesh_->vertex_handle(startTop + ((sl + 1) % slices_)));

      triMesh_->add_face(vhandles);
    }

    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

655
656
    PluginFunctions::viewAll();

657
658
659
660
    return object->id();
  }


Jan Möbius's avatar
Jan Möbius committed
661
662
}

663
664
665
666
667

//========================================================================
// Pyramid
//========================================================================

668
int PrimitivesGeneratorPlugin::addPyramid(const Vector& _position,const double _length) {
Jan Möbius's avatar
Jan Möbius committed
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
  int newObject = addTriMesh();

  TriMeshObject* object;
  if ( !PluginFunctions::getObject(newObject,object) ) {
        emit log(LOGERR,"Unable to create new Object");
        return -1;
  } else {

    object->setName( "Pyramid " + QString::number(newObject) );

    triMesh_ =  object->mesh();

    triMesh_->clear();

    // Add 5 vertices
    vhandles_.resize(5);

686
    const double halfLength   = 0.5*_length;
Jan Möbius's avatar
Jan Möbius committed
687

688
689
690
691
    vhandles_[0] = triMesh_->add_vertex(TriMesh::Point( halfLength, -halfLength, 0.0)+_position);
    vhandles_[1] = triMesh_->add_vertex(TriMesh::Point( halfLength,  halfLength, 0.0)+_position);
    vhandles_[2] = triMesh_->add_vertex(TriMesh::Point(-halfLength,  halfLength, 0.0)+_position);
    vhandles_[3] = triMesh_->add_vertex(TriMesh::Point(-halfLength, -halfLength, 0.0)+_position);
Jan Möbius's avatar
Jan Möbius committed
692

693
    vhandles_[4] = triMesh_->add_vertex(TriMesh::Point(0.0, 0.0, sqrt(2.0)*halfLength));
Jan Möbius's avatar
Jan Möbius committed
694
695
696
697
698
699
700
701
702
703
704
705
706
707

    // Add 6 faces
    add_face(2,1,0);
    add_face(3,2,0);
    add_face(4,0,1);

    add_face(3,0,4);
    add_face(4,2,3);
    add_face(1,2,4);

    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

708
709
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
    return newObject;
  }

  return -1;
}

void PrimitivesGeneratorPlugin::add_face( int _vh1 , int _vh2, int _vh3 ) {
 std::vector<TriMesh::VertexHandle> vhandles;

 vhandles.push_back(vhandles_[_vh1]);
 vhandles.push_back(vhandles_[_vh2]);
 vhandles.push_back(vhandles_[_vh3]);

 triMesh_->add_face(vhandles);

}

Jan Möbius's avatar
Jan Möbius committed
727
728
729
730
731
732
733
734
735
736
737
738
void PrimitivesGeneratorPlugin::add_face( int _vh1 , int _vh2, int _vh3, int _vh4 , int _vh5 ) {
  std::vector<PolyMesh::VertexHandle> vhandles;

  vhandles.push_back(vphandles_[_vh1]);
  vhandles.push_back(vphandles_[_vh2]);
  vhandles.push_back(vphandles_[_vh3]);
  vhandles.push_back(vphandles_[_vh4]);
  vhandles.push_back(vphandles_[_vh5]);

  polyMesh_->add_face(vhandles);
}

739
int PrimitivesGeneratorPlugin::addIcosahedron(const Vector& _position,const double _length) {
Jan Möbius's avatar
Jan Möbius committed
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
  int newObject = addTriMesh();

  TriMeshObject* object;
  if ( !PluginFunctions::getObject(newObject,object) ) {
        emit log(LOGERR,"Unable to create new Object");
        return -1;
  } else {

    object->setName( "Icosahedron " + QString::number(newObject) );

    triMesh_ =  object->mesh();

    triMesh_->clear();

    // Add 12 vertices
    vhandles_.resize(12);

757
    const double phi  = 0.5 * (1.0 + sqrt(5.0));
Jan Möbius's avatar
Jan Möbius committed
758
    //double norm = 1.0 / sqrt(1.0 + phi*phi);
759
760
    const double norm = 1.0;
    const double halfLength = 0.5*_length;
Jan Möbius's avatar
Jan Möbius committed
761

762
763
764
765
    vhandles_[0 ] = triMesh_->add_vertex(norm * TriMesh::Point(  0.0 ,  -halfLength , -phi )+_position);
    vhandles_[1 ] = triMesh_->add_vertex(norm * TriMesh::Point(  0.0 ,   halfLength , -phi )+_position);
    vhandles_[2 ] = triMesh_->add_vertex(norm * TriMesh::Point(  0.0 ,   halfLength ,  phi )+_position);
    vhandles_[3 ] = triMesh_->add_vertex(norm * TriMesh::Point(  0.0 ,  -halfLength ,  phi )+_position);
Jan Möbius's avatar
Jan Möbius committed
766

767
768
769
770
    vhandles_[4 ] = triMesh_->add_vertex(norm * TriMesh::Point( -halfLength , -phi ,  0.0  )+_position);
    vhandles_[5 ] = triMesh_->add_vertex(norm * TriMesh::Point(  halfLength , -phi ,  0.0  )+_position);
    vhandles_[6 ] = triMesh_->add_vertex(norm * TriMesh::Point(  halfLength ,  phi ,  0.0  )+_position);
    vhandles_[7 ] = triMesh_->add_vertex(norm * TriMesh::Point( -halfLength ,  phi ,  0.0  )+_position);
Jan Möbius's avatar
Jan Möbius committed
771

772
773
774
775
    vhandles_[8 ] = triMesh_->add_vertex(norm * TriMesh::Point( -phi ,  0.0 , -halfLength  )+_position);
    vhandles_[9 ] = triMesh_->add_vertex(norm * TriMesh::Point( -phi ,  0.0 ,  halfLength  )+_position);
    vhandles_[10] = triMesh_->add_vertex(norm * TriMesh::Point(  phi ,  0.0 ,  halfLength  )+_position);
    vhandles_[11] = triMesh_->add_vertex(norm * TriMesh::Point(  phi ,  0.0 , -halfLength  )+_position);
Jan Möbius's avatar
Jan Möbius committed
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812


    // Add 20 faces
    add_face(2,6,7);
    add_face(7,6,1);

    add_face(11,0,1);
    add_face(0, 8,1);

    add_face(4,9,8);
    add_face(8,9,7);

    add_face(9,3,2);
    add_face(10,2,3);

    add_face(5,11,10);
    add_face(11,6,10);

    add_face(0,5,4);
    add_face(5,3,4);

    // Upper block
    add_face(6,2,10);
    add_face(6,11,1);
    add_face(1,8,7);
    add_face(9,2,7);

    // Lower block
    add_face(3,5,10);
    add_face(0,11,5);
    add_face(3,9,4);
    add_face(0,4,8);

    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

813
814
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
815
816
817
818
819
820
    return newObject;
  }

  return -1;
}

821
int PrimitivesGeneratorPlugin::addOctahedron(const Vector& _position,const double _length) {
Jan Möbius's avatar
Jan Möbius committed
822
823
824
825
826
827
828
829
830
831
832
833
834
835
  int newObject = addTriMesh();

  TriMeshObject* object;
  if ( !PluginFunctions::getObject(newObject,object) ) {
    emit log(LOGERR,"Unable to create new Object");
    return -1;
  } else {

    object->setName( "Octahedron " + QString::number(newObject) );

    triMesh_ =  object->mesh();

    triMesh_->clear();

Jan Möbius's avatar
Jan Möbius committed
836
    // Add 6 vertices
Jan Möbius's avatar
Jan Möbius committed
837
838
    vhandles_.resize(6);

839
    const double sqrtLength = sqrt(_length);
Jan Möbius's avatar
Jan Möbius committed
840

841
842
843
844
    vhandles_[0 ] = triMesh_->add_vertex(TriMesh::Point(-sqrtLength,   0.0,   0.0)+_position);
    vhandles_[1 ] = triMesh_->add_vertex(TriMesh::Point(  0.0, -sqrtLength,   0.0)+_position);
    vhandles_[2 ] = triMesh_->add_vertex(TriMesh::Point( sqrtLength,   0.0,   0.0)+_position);
    vhandles_[3 ] = triMesh_->add_vertex(TriMesh::Point(  0.0,  sqrtLength,   0.0)+_position);
Jan Möbius's avatar
Jan Möbius committed
845

846
847
    vhandles_[4 ] = triMesh_->add_vertex(TriMesh::Point(  0.0,   0.0,  sqrtLength)+_position);
    vhandles_[5 ] = triMesh_->add_vertex(TriMesh::Point(  0.0,   0.0, -sqrtLength)+_position);
Jan Möbius's avatar
Jan Möbius committed
848
849
850
851
852
853
854
855


    // Add 8 faces
    add_face(0,1,4);
    add_face(1,2,4);

    add_face(2,3,4);
    add_face(0,4,3);
Jan Möbius's avatar
Jan Möbius committed
856

Jan Möbius's avatar
Jan Möbius committed
857
858
859
860
861
862
863
864
865
866
    add_face(5,1,0);
    add_face(5,2,1);

    add_face(5,3,2);
    add_face(5,0,3);

    triMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

867
868
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
869
870
871
872
873
874
    return newObject;
  }

  return -1;
}

875
int PrimitivesGeneratorPlugin::addDodecahedron(const Vector& _position,const double _length) {
Jan Möbius's avatar
Jan Möbius committed
876
877
878
879
880
881
882
883
  int newObject = addPolyMesh();

  PolyMeshObject* object;
  if ( !PluginFunctions::getObject(newObject,object) ) {
    emit log(LOGERR,"Unable to create new Object");
    return -1;
  } else {

Jan Möbius's avatar
Jan Möbius committed
884
    object->setName( "Dodecahedron " + QString::number(newObject) );
Jan Möbius's avatar
Jan Möbius committed
885
886
887
888
889
890
891
892

    polyMesh_ =  object->mesh();

    polyMesh_->clear();

    // Add 20 vertices
    vphandles_.resize(20);

893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
    const double phi = (1.0 + sqrt(5.0)) / 2.0;
    const double halfLength = 0.5*_length;


    vphandles_[0 ] = polyMesh_->add_vertex(TriMesh::Point(  halfLength ,  halfLength , halfLength )+_position);
    vphandles_[1 ] = polyMesh_->add_vertex(TriMesh::Point(  halfLength ,  halfLength ,-halfLength )+_position);
    vphandles_[2 ] = polyMesh_->add_vertex(TriMesh::Point(  halfLength , -halfLength , halfLength )+_position);
    vphandles_[3 ] = polyMesh_->add_vertex(TriMesh::Point(  halfLength , -halfLength ,-halfLength )+_position);
    vphandles_[4 ] = polyMesh_->add_vertex(TriMesh::Point( -halfLength ,  halfLength , halfLength )+_position);
    vphandles_[5 ] = polyMesh_->add_vertex(TriMesh::Point( -halfLength ,  halfLength ,-halfLength )+_position);
    vphandles_[6 ] = polyMesh_->add_vertex(TriMesh::Point( -halfLength , -halfLength , halfLength )+_position);
    vphandles_[7 ] = polyMesh_->add_vertex(TriMesh::Point( -halfLength , -halfLength ,-halfLength )+_position);

    vphandles_[8 ] = polyMesh_->add_vertex(TriMesh::Point( 0.0 ,  halfLength / phi ,  phi )+_position);
    vphandles_[9 ] = polyMesh_->add_vertex(TriMesh::Point( 0.0 ,  halfLength / phi , -phi )+_position);
    vphandles_[10] = polyMesh_->add_vertex(TriMesh::Point( 0.0 , -halfLength / phi ,  phi )+_position);
    vphandles_[11] = polyMesh_->add_vertex(TriMesh::Point( 0.0 , -halfLength / phi , -phi )+_position);

    vphandles_[12] = polyMesh_->add_vertex(TriMesh::Point(  halfLength / phi ,  phi, 0.0)+_position);
    vphandles_[13] = polyMesh_->add_vertex(TriMesh::Point(  halfLength / phi , -phi, 0.0)+_position);
    vphandles_[14] = polyMesh_->add_vertex(TriMesh::Point( -halfLength / phi ,  phi, 0.0)+_position);
    vphandles_[15] = polyMesh_->add_vertex(TriMesh::Point( -halfLength / phi , -phi, 0.0)+_position);

    vphandles_[16] = polyMesh_->add_vertex(TriMesh::Point(  phi , 0.0 , halfLength / phi)+_position);
    vphandles_[17] = polyMesh_->add_vertex(TriMesh::Point(  phi , 0.0 ,-halfLength / phi)+_position);
    vphandles_[18] = polyMesh_->add_vertex(TriMesh::Point( -phi , 0.0 , halfLength / phi)+_position);
    vphandles_[19] = polyMesh_->add_vertex(TriMesh::Point( -phi , 0.0 ,-halfLength / phi)+_position);
Jan Möbius's avatar
Jan Möbius committed
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940

    // Add 12 faces
    add_face(14, 5,19,18, 4);
    add_face( 5, 9,11, 7,19);
    add_face( 6,15,13, 2,10);
    add_face(12, 0,16,17, 1);

    add_face( 0, 8,10, 2,16);
    add_face(16, 2,13, 3,17);
    add_face( 3,13,15, 7,11);
    add_face( 7,15, 6,18,19);

    add_face( 4,18, 6,10, 8);
    add_face( 4, 8, 0,12,14);
    add_face(14,12, 1, 9, 5);
    add_face( 9, 1,17, 3,11);

    polyMesh_->update_normals();

    emit updatedObject(newObject,UPDATE_ALL);

941
942
    PluginFunctions::viewAll();

Jan Möbius's avatar
Jan Möbius committed
943
944
945
946
947
    return newObject;
  }

  return -1;
}
Jan Möbius's avatar
Jan Möbius committed
948
949


Matthias Möller's avatar
Matthias Möller committed
950
951
952
#if QT_VERSION < 0x050000
  Q_EXPORT_PLUGIN2( primitivesgeneratorplugin , PrimitivesGeneratorPlugin );
#endif