0024167: Compiler warnings 'unreacheable code' and 'conditional expression is constan...
[occt.git] / samples / qt / VoxelDemo / src / Application.cpp
1 #include "Application.h"
2 #include "ConversionThread.h"
3 #include "Timer.h"
4
5 #include <QPixmap.h>
6 #include <QToolButton.h>
7 #include <QWhatsThis.h>
8 #include <QMenu.h>
9 #include <QMenuBar.h>
10 #include <QStatusBar.h>
11 #include <QApplication.h>
12 #include <QFileDialog.h>
13 #include <QMessageBox.h>
14 #include <QInputDialog.h>
15 #include <QCloseEvent>
16
17 #include <Voxel_BoolDS.hxx>
18 #include <Voxel_ColorDS.hxx>
19 #include <Voxel_FloatDS.hxx>
20 #include <Voxel_OctBoolDS.hxx>
21 #include <Voxel_ROctBoolDS.hxx>
22 #include <Voxel_BooleanOperation.hxx>
23 #include <Voxel_CollisionDetection.hxx>
24 #include <Voxel_FastConverter.hxx>
25 #include <Voxel_Writer.hxx>
26 #include <Voxel_Reader.hxx>
27 #include <VoxelClient_VisDrawer.h>
28
29 #include <BRepTools.hxx>
30 #include <BRepBndLib.hxx>
31 #include <BRep_Builder.hxx>
32 #include <BRepPrimAPI_MakeBox.hxx>
33 #include <BRepPrimAPI_MakeTorus.hxx>
34 #include <BRepPrimAPI_MakeSphere.hxx>
35 #include <BRepPrimAPI_MakeCylinder.hxx>
36 #include <Aspect_ColorScale.hxx>
37
38 #include <Windows.h>
39
40 Application::Application()
41     : QMainWindow( 0 )
42 {
43     // File
44     QMenu * file = menuBar()->addMenu( "&File" );
45
46     QAction* a;
47     // Box
48     a = new QAction("Box", this);
49     connect(a, SIGNAL(triggered()), this, SLOT(box()));
50     file->addAction(a);
51     // Cylinder
52     a = new QAction("Cylinder", this);
53     connect(a, SIGNAL(triggered()), this, SLOT(cylinder()));
54     file->addAction(a);
55     // Torus
56     a = new QAction("Torus", this);
57     connect(a, SIGNAL(triggered()), this, SLOT(torus()));
58     file->addAction(a);
59     // Sphere
60     a = new QAction("Sphere", this);
61     connect(a, SIGNAL(triggered()), this, SLOT(sphere()));
62     file->addAction(a);
63     // Load shape...
64     a = new QAction("Load shape...", this);
65     a->setShortcut(tr("Ctrl+O"));
66     connect(a, SIGNAL(triggered()), this, SLOT(choose()));
67     file->addAction(a);
68
69     file->addSeparator();
70
71     // Open
72     a = new QAction("Open", this);
73     connect(a, SIGNAL(triggered()), this, SLOT(open()));
74     file->addAction(a);
75
76     // Save
77     a = new QAction("Save", this);
78     connect(a, SIGNAL(triggered()), this, SLOT(save()));
79     file->addAction(a);
80
81     file->addSeparator();
82
83     // Quit
84     a = new QAction("&Quit", this);
85     a->setShortcut(tr("Ctrl+Q"));
86     connect(a, SIGNAL(triggered()), qApp, SLOT(closeAllWindows()));
87     file->addAction(a);
88
89     menuBar()->addSeparator();
90
91
92 #ifdef TEST
93     QMenu * test = menuBar()->addMenu( "Test" );
94
95     a = new QAction("Test boolean", this);
96     connect(a, SIGNAL(triggered()), this, SLOT(testBoolDS()));
97     test->addAction(a);
98
99     a = new QAction("Test color", this);
100     connect(a, SIGNAL(triggered()), this, SLOT(testColorDS()));
101     test->addAction(a);
102
103     a = new QAction("Test float", this);
104     connect(a, SIGNAL(triggered()), this, SLOT(testFloatDS()));
105     test->addAction(a);
106
107     a = new QAction("Test boolean / 8", this);
108     connect(a, SIGNAL(triggered()), this, SLOT(testOctBoolDS()));
109     test->addAction(a);
110
111     a = new QAction("Test boolean / 8 / 8..", this);
112     connect(a, SIGNAL(triggered()), this, SLOT(testROctBoolDS()));
113     test->addAction(a);
114
115     test->addSeparator();
116
117     a = new QAction("Test fusion of booleans", this);
118     connect(a, SIGNAL(triggered()), this, SLOT(testFuseBoolDS()));
119     test->addAction(a);
120
121     a = new QAction("Test fusion of colors", this);
122     connect(a, SIGNAL(triggered()), this, SLOT(testFuseColorDS()));
123     test->addAction(a);
124
125     a = new QAction("Test fusion of floating-points", this);
126     connect(a, SIGNAL(triggered()), this, SLOT(testFuseFloatDS()));
127     test->addAction(a);
128
129     a = new QAction("Test cutting of booleans", this);
130     connect(a, SIGNAL(triggered()), this, SLOT(testCutBoolDS()));
131     test->addAction(a);
132
133     a = new QAction("Test cutting of booleans", this);
134     connect(a, SIGNAL(triggered()), this, SLOT(testCutColorDS()));
135     test->addAction(a);
136
137     a = new QAction("Test cutting of floating-points", this);
138     connect(a, SIGNAL(triggered()), this, SLOT(testCutFloatDS()));
139     test->addAction(a);
140
141 #endif // TEST
142
143     QMenu * converter = menuBar()->addMenu( "Converter" );
144
145 #ifdef TEST
146     
147     a = new QAction("Number of splits along X", this);
148     connect(a, SIGNAL(triggered()), this, SLOT(setNbX()));
149     converter->addAction(a);
150
151     a = new QAction("Number of splits along Y", this);
152     connect(a, SIGNAL(triggered()), this, SLOT(setNbY()));
153     converter->addAction(a);
154
155     a = new QAction("Number of splits along Z", this);
156     connect(a, SIGNAL(triggered()), this, SLOT(setNbZ()));
157     converter->addAction(a);
158
159     converter->addSeparator();
160
161     a = new QAction("Side of scanning", this);
162     connect(a, SIGNAL(triggered()), this, SLOT(setScanSide()));
163     converter->addAction(a);
164
165     converter->addSeparator();
166
167     a = new QAction("Volumic value of 1bit voxels", this);
168     connect(a, SIGNAL(triggered()), this, SLOT(setVolumicBoolValue()));
169     converter->addAction(a);
170     
171     a = new QAction("Volumic value of 4bit voxels", this);
172     connect(a, SIGNAL(triggered()), this, SLOT(setVolumicColorValue()));
173     converter->addAction(a);
174
175     converter->addSeparator();
176
177 #endif // TEST
178
179     a = new QAction("Convert to 1bit voxels", this);
180     connect(a, SIGNAL(triggered()), this, SLOT(convert2bool()));
181     converter->addAction(a);
182
183     a = new QAction("Convert to 4bit voxels", this);
184     connect(a, SIGNAL(triggered()), this, SLOT(convert2color()));
185     converter->addAction(a);
186
187     QMenu * vis = menuBar()->addMenu( "&Visualization" );
188
189     a = new QAction("Points", this);
190     connect(a, SIGNAL(triggered()), this, SLOT(displayPoints()));
191     vis->addAction(a);
192
193 #ifdef TEST
194     
195     a = new QAction("Nearest points", this);
196     connect(a, SIGNAL(triggered()), this, SLOT(displayNearestPoints()));
197     vis->addAction(a);
198
199 #endif // TEST
200
201     a = new QAction("Boxes", this);
202     connect(a, SIGNAL(triggered()), this, SLOT(displayBoxes()));
203     vis->addAction(a);
204
205 #ifdef TEST
206
207     a = new QAction("Nearest boxes", this);
208     connect(a, SIGNAL(triggered()), this, SLOT(displayNearestBoxes()));
209     vis->addAction(a);
210
211 #endif // TEST
212     
213     vis->addSeparator();
214
215     a = new QAction("Point size", this);
216     connect(a, SIGNAL(triggered()), this, SLOT(setPointSize()));
217     vis->addAction(a);
218     
219     a = new QAction("Quadrangle size (%)", this);
220     connect(a, SIGNAL(triggered()), this, SLOT(setQuadrangleSize()));
221     vis->addAction(a);
222
223     vis->addSeparator();
224
225     a = new QAction("Color min value", this);
226     connect(a, SIGNAL(triggered()), this, SLOT(setColorMinValue()));
227     vis->addAction(a);
228
229     a = new QAction("Color max value", this);
230     connect(a, SIGNAL(triggered()), this, SLOT(setColorMaxValue()));
231     vis->addAction(a);
232
233 #ifdef TEST
234
235     vis->addSeparator();
236
237     a = new QAction("Use GL lists", this);
238     connect(a, SIGNAL(triggered()), this, SLOT(setUsageOfGLlists()));
239     vis->addAction(a);
240
241 #endif // TEST
242
243     vis->addSeparator();
244
245     a = new QAction("Displayed X min", this);
246     connect(a, SIGNAL(triggered()), this, SLOT(setDisplayedXMin()));
247     vis->addAction(a);
248
249     a = new QAction("Displayed X max", this);
250     connect(a, SIGNAL(triggered()), this, SLOT(setDisplayedXMax()));
251     vis->addAction(a);
252
253     a = new QAction("Displayed Y min", this);
254     connect(a, SIGNAL(triggered()), this, SLOT(setDisplayedYMin()));
255     vis->addAction(a);
256
257     a = new QAction("Displayed Y max", this);
258     connect(a, SIGNAL(triggered()), this, SLOT(setDisplayedYMax()));
259     vis->addAction(a);
260
261     a = new QAction("Displayed Z min", this);
262     connect(a, SIGNAL(triggered()), this, SLOT(setDisplayedZMin()));
263     vis->addAction(a);
264
265     a = new QAction("Displayed Z max", this);
266     connect(a, SIGNAL(triggered()), this, SLOT(setDisplayedZMax()));
267     vis->addAction(a);
268
269
270     QMenu * demo = menuBar()->addMenu( "Demo" );
271
272     a = new QAction("Waves", this);
273     connect(a, SIGNAL(triggered()), this, SLOT(displayWaves()));
274     demo->addAction(a);
275
276     a = new QAction("Cut", this);
277     connect(a, SIGNAL(triggered()), this, SLOT(displayCut()));
278     demo->addAction(a);
279
280     a = new QAction("Collisions", this);
281     connect(a, SIGNAL(triggered()), this, SLOT(displayCollisions()));
282     demo->addAction(a);
283
284
285     QMenu * help = menuBar()->addMenu( "Help" );
286
287     a = new QAction("About", this);
288     a->setShortcut(tr("F1"));
289     connect(a, SIGNAL(triggered()), this, SLOT(about()));
290     help->addAction(a);
291
292
293     myViewer = new Viewer( this );
294     myViewer->setFocus();
295     setCentralWidget( myViewer );
296     statusBar()->showMessage( "Ready", 2000 );
297
298         myNbX = 100;
299         myNbY = 100;
300         myNbZ = 100;
301
302     myScanSide = 7;
303
304     myVolumicBoolValue = false;
305     myVolumicColorValue = 0;
306
307     myQuadrangleSize = 40;
308     
309     myColorMinValue = 1;
310     myColorMaxValue = 15;
311
312         myBoolVoxels = 0;
313         myColorVoxels = 0;
314
315     myDisplayedXMin = -DBL_MAX;
316     myDisplayedXMax =  DBL_MAX;
317     myDisplayedYMin = -DBL_MAX;
318     myDisplayedYMax =  DBL_MAX;
319     myDisplayedZMin = -DBL_MAX;
320     myDisplayedZMax =  DBL_MAX;
321
322         VoxelClient_VisDrawer::Init(myViewer->getGraphicDriver());
323
324     resize( 450, 600 );
325 }
326
327 Application::~Application()
328 {
329         if (myBoolVoxels)
330                 delete myBoolVoxels;
331         if (myColorVoxels)
332                 delete myColorVoxels;
333 }
334
335 void Application::choose()
336 {
337     QString fn = QFileDialog::getOpenFileName( this, QString::null, QString::null, "*.brep");
338     if ( !fn.isEmpty() )
339                 load( fn );
340     else
341                 statusBar()->showMessage( "Loading aborted", 2000 );
342 }
343
344 void Application::load( const QString &fileName )
345 {
346     QFile f( fileName );
347     if ( !f.open( QIODevice::ReadOnly ) )
348                 return;
349
350     // Read shape
351         TopoDS_Shape S;
352         BRep_Builder B;
353         if (!BRepTools::Read(S, (char*) fileName.constData(), B))
354                 statusBar()->showMessage( "Loading failed", 2000 );
355     
356     load(S);
357 }
358
359 void Application::open()
360 {
361     QString fn = QFileDialog::getOpenFileName( this, QString::null, QString::null, "*.vx");
362     if ( fn.isEmpty() || !QFile::exists(fn) )
363     {
364                 statusBar()->showMessage( "Open aborted", 2000 );
365         return;
366     }
367
368     Timer timer;
369     timer.Start();
370
371     // Read the voxels
372     Voxel_Reader reader;
373     if (!reader.Read((char*)fn.constData()))
374     {
375                 statusBar()->showMessage( "Open failed... sorry", 2000 );
376         return;
377     }
378
379     timer.Stop();
380     timer.Print("Open");
381
382     // Release current voxels
383     if (myBoolVoxels)
384     {
385         delete myBoolVoxels;
386         myBoolVoxels = 0;
387     }
388     if (myColorVoxels)
389     {
390         delete myColorVoxels;
391         myColorVoxels = 0;
392     }
393
394     // Take the voxels
395     if (reader.IsBoolVoxels())
396     {
397         myBoolVoxels = (Voxel_BoolDS*) reader.GetBoolVoxels();
398         myViewer->getSelector().SetVoxels(*myBoolVoxels);
399     }
400     else if (reader.IsColorVoxels())
401     {
402         myColorVoxels = (Voxel_ColorDS*) reader.GetColorVoxels();
403         myViewer->getSelector().SetVoxels(*myColorVoxels);
404     }
405     
406     // Display the voxels
407     myViewer->getIC()->EraseAll(false, false);
408     Voxel_DS* ds = myBoolVoxels;
409     if (!ds)
410         ds = myColorVoxels;
411     if (ds)
412     {
413         myDisplayedXMin = ds->GetX() - 10.0 * Precision::Confusion();
414         myDisplayedXMax = ds->GetX() + ds->GetXLen() + 10.0 * Precision::Confusion();
415         myDisplayedYMin = ds->GetY() - 10.0 * Precision::Confusion();
416         myDisplayedYMax = ds->GetY() + ds->GetYLen() + 10.0 * Precision::Confusion();
417         myDisplayedZMin = ds->GetZ() - 10.0 * Precision::Confusion();
418         myDisplayedZMax = ds->GetZ() + ds->GetZLen() + 10.0 * Precision::Confusion();
419     }
420
421     // Init visual data
422     initPrs();
423
424     // Set voxels and display
425     Handle(Poly_Triangulation) empty;
426     myVoxels->SetBoolVoxels(myBoolVoxels);
427         myVoxels->SetColorVoxels(myColorVoxels);
428     myVoxels->SetTriangulation(empty);
429     if (myViewer->getIC()->IsDisplayed(myVoxels))
430         myViewer->getIC()->Redisplay(myVoxels, false);
431     else
432             myViewer->getIC()->Display(myVoxels, false);
433    
434     // Color scale
435     if (myColorVoxels)
436         displayColorScale();
437     else
438         myViewer->getView()->ColorScaleErase();
439
440         myViewer->getView()->FitAll();
441
442     statusBar()->showMessage( "Ready.", 2000 );
443 }
444
445 void Application::save()
446 {
447     QString fn = QFileDialog::getSaveFileName( this, QString::null, QString::null, "*.vx");
448     if ( fn.isEmpty() )
449     {
450                 statusBar()->showMessage( "Storage aborted", 2000 );
451         return;
452     }
453     if (fn.indexOf(".vx", -1, Qt::CaseInsensitive) == -1)
454         fn += ".vx";
455
456     Timer timer;
457     timer.Start();
458
459     // Write the voxels
460     Voxel_Writer writer;
461     writer.SetFormat(Voxel_VFF_BINARY);
462     if (myBoolVoxels)
463         writer.SetVoxels(*myBoolVoxels);
464     else if (myColorVoxels)
465         writer.SetVoxels(*myColorVoxels);
466     else
467     {
468                 statusBar()->showMessage( "Nothing to store", 2000 );
469         return;
470     }
471     if (!writer.Write((char*)fn.constData()))
472     {
473                 statusBar()->showMessage( "Storage failed... sorry", 2000 );
474         return;
475     }
476
477     timer.Stop();
478     timer.Print("Save");
479
480     statusBar()->showMessage( "Saved.", 2000 );
481 }
482
483 void Application::closeEvent( QCloseEvent* ce )
484 {
485     ce->accept();
486 }
487
488 void Application::about()
489 {
490     QMessageBox::about( this, "Voxel demo-application",
491                         "This example demonstrates simple usage of "
492                         "voxel models of Open CASCADE.");
493 }
494
495 void Application::testBoolDS()
496 {
497         Timer timer;
498         int ix, iy, iz;
499         int nbx = 100, nby = 100, nbz = 100;
500
501
502         // 1. BoolDS:
503
504         timer.Start();
505
506         Voxel_BoolDS ds(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
507
508         for (ix = 0; ix < nbx; ix++)
509         {
510                 for (iy = 0; iy < nby; iy++)
511                 {
512                         for (iz = 0; iz < nbz; iz++)
513                         {
514                                 if (ix & 0x01)
515                                         ds.Set(ix, iy, iz, false);
516                                 else
517                                         ds.Set(ix, iy, iz, true);
518                         }
519                 }
520         }
521
522         for (ix = 0; ix < nbx; ix++)
523         {
524                 for (iy = 0; iy < nby; iy++)
525                 {
526                         for (iz = 0; iz < nbz; iz++)
527                         {
528                                 bool value = ds.Get(ix, iy, iz) == Standard_True;
529                                 if (ix & 0x01)
530                                 {
531                                         if (value != false)
532                                                 cout<<"Wrong value!"<<endl;
533                                 }
534                                 else
535                                 {
536                                         if (value != true)
537                                                 cout<<"Wrong value!"<<endl;
538                                 }
539                         }
540                 }
541         }
542
543         timer.Stop();
544         timer.Print("BoolDS");
545 }
546
547 void Application::testColorDS()
548 {
549         Timer timer;
550         int ix, iy, iz;
551         int nbx = 100, nby = 100, nbz = 100;
552
553
554         // 1. ColorDS:
555
556         timer.Start();
557
558         Voxel_ColorDS ds(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
559
560         for (ix = 0; ix < nbx; ix++)
561         {
562                 for (iy = 0; iy < nby; iy++)
563                 {
564                         for (iz = 0; iz < nbz; iz++)
565                         {
566                 if (ix & 0x01)
567                             ds.Set(ix, iy, iz, 8);
568                 else
569                     ds.Set(ix, iy, iz, 7);
570                         }
571                 }
572         }
573
574         for (ix = 0; ix < nbx; ix++)
575         {
576                 for (iy = 0; iy < nby; iy++)
577                 {
578                         for (iz = 0; iz < nbz; iz++)
579                         {
580                                 unsigned char value = ds.Get(ix, iy, iz);
581                 if (ix & 0x01)
582                 {
583                                     if (value != 8)
584                                             cout<<"Wrong value!"<<endl;
585                 }
586                 else
587                 {
588                                     if (value != 7)
589                                             cout<<"Wrong value!"<<endl;
590                 }
591                         }
592                 }
593         }
594
595         timer.Stop();
596         timer.Print("ColorDS");
597 }
598
599 void Application::testFloatDS()
600 {
601         Timer timer;
602         int ix, iy, iz;
603         int nbx = 100, nby = 100, nbz = 100;
604
605
606         // 1. FloatDS:
607
608         timer.Start();
609
610         Voxel_FloatDS ds(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
611
612         for (ix = 0; ix < nbx; ix++)
613         {
614                 for (iy = 0; iy < nby; iy++)
615                 {
616                         for (iz = 0; iz < nbz; iz++)
617                         {
618                 if (ix & 0x01)
619                             ds.Set(ix, iy, iz, 8.8f);
620                 else
621                     ds.Set(ix, iy, iz, 7.7f);
622                         }
623                 }
624         }
625
626         for (ix = 0; ix < nbx; ix++)
627         {
628                 for (iy = 0; iy < nby; iy++)
629                 {
630                         for (iz = 0; iz < nbz; iz++)
631                         {
632                                 float value = ds.Get(ix, iy, iz);
633                 if (ix & 0x01)
634                 {
635                                     if (value != 8.8f)
636                                             cout<<"Wrong value!"<<endl;
637                 }
638                 else
639                 {
640                                     if (value != 7.7f)
641                                             cout<<"Wrong value!"<<endl;
642                 }
643                         }
644                 }
645         }
646
647         timer.Stop();
648         timer.Print("FloatDS");
649 }
650
651 void Application::testOctBoolDS()
652 {
653         Timer timer;
654         int ix, iy, iz;
655         int nbx = 30, nby = 30, nbz = 30;
656
657
658         // 1. OctBoolDS:
659
660         timer.Start();
661
662         Voxel_OctBoolDS ds(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
663
664         for (ix = 0; ix < nbx; ix++)
665         {
666                 for (iy = 0; iy < nby; iy++)
667                 {
668                         for (iz = 0; iz < nbz; iz++)
669                         {
670                                 if (ix & 0x01)
671                 {
672                                         ds.Set(ix, iy, iz, true);
673                 }
674                                 else
675                 {
676                     for (int i = 0; i < 8; i++)
677                     {
678                         if (i & 0x01)
679                                                 ds.Set(ix, iy, iz, i, true);
680                         else
681                             ds.Set(ix, iy, iz, i, false);
682                     }
683                 }
684                         }
685                 }
686         }
687
688         for (ix = 0; ix < nbx; ix++)
689         {
690                 for (iy = 0; iy < nby; iy++)
691                 {
692                         for (iz = 0; iz < nbz; iz++)
693                         {
694                                 if (ix & 0x01)
695                                 {
696                                 bool value = ds.Get(ix, iy, iz) == Standard_True;
697                                         if (value != true)
698                                                 cout<<"Wrong value!"<<endl;
699                                 }
700                                 else
701                                 {
702                     for (int i = 0; i < 8; i++)
703                     {
704                         if (i & 0x01)
705                         {
706                                         bool value = ds.Get(ix, iy, iz, i) == Standard_True;
707                                                 if (value != true)
708                                                         cout<<"Wrong value!"<<endl;
709                         }
710                         else
711                         {
712                                         bool value = ds.Get(ix, iy, iz, i) == Standard_True;
713                                                 if (value != false)
714                                                         cout<<"Wrong value!"<<endl;
715                         }
716                     }
717                                 }
718                         }
719                 }
720         }
721
722         for (ix = 0; ix < nbx; ix++)
723         {
724                 for (iy = 0; iy < nby; iy++)
725                 {
726                         for (iz = 0; iz < nbz; iz++)
727                         {
728                                 if (ix & 0x01)
729                                 {
730                     for (int i = 0; i < 8; i++)
731                     {
732                         ds.Set(ix, iy, iz, i, true);
733                     }
734                                 }
735                 else
736                 {
737                     for (int i = 0; i < 8; i++)
738                     {
739                         ds.Set(ix, iy, iz, i, false);
740                     }
741                 }
742                         }
743                 }
744         }
745
746     ds.OptimizeMemory();
747
748         timer.Stop();
749         timer.Print("OctBoolDS");
750 }
751
752 void Application::testROctBoolDS()
753 {
754         Timer timer;
755         int ix, iy, iz, i, j;
756         int nbx = 30, nby = 30, nbz = 30;
757
758         // 1. ROctBoolDS:
759
760         timer.Start();
761
762         Voxel_ROctBoolDS ds(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
763
764     for (ix = 0; ix < nbx; ix++)
765         {
766                 for (iy = 0; iy < nby; iy++)
767                 {
768                         for (iz = 0; iz < nbz; iz++)
769                         {
770                 ds.Set(ix, iy, iz, true);
771                         }
772                 }
773         }
774
775     for (ix = 0; ix < nbx; ix++)
776         {
777                 for (iy = 0; iy < nby; iy++)
778                 {
779                         for (iz = 0; iz < nbz; iz++)
780                         {
781                 for (i = 0; i < 8; i++)
782                 {
783                     for (j = 0; j < 8; j++)
784                     {
785                         ds.Set(ix, iy, iz, i, j, true);
786                     }
787                 }
788                         }
789                 }
790         }
791
792     ds.OptimizeMemory();
793
794         for (ix = 0; ix < nbx; ix++)
795         {
796                 for (iy = 0; iy < nby; iy++)
797                 {
798                         for (iz = 0; iz < nbz; iz++)
799                         {
800                 if (ds.Deepness(ix, iy, iz) == 0)
801                 {
802                     bool value = ds.Get(ix, iy, iz);
803                     if (value != true)
804                         cout<<"Wrong value..."<<endl;
805                 }
806                 if (ds.Deepness(ix, iy, iz) == 1)
807                 {
808                     for (i = 0; i < 8; i++)
809                     {
810                         bool value = ds.Get(ix, iy, iz, i);
811                         if (value != true)
812                             cout<<"Wrong value..."<<endl;
813                     }
814                 }
815                 if (ds.Deepness(ix, iy, iz) == 2)
816                 {
817                     for (i = 0; i < 8; i++)
818                     {
819                         for (j = 0; j < 8; j++)
820                         {
821                             bool value = ds.Get(ix, iy, iz, i, j);
822                             if (value != true)
823                                 cout<<"Wrong value..."<<endl;
824                         }
825                     }
826                 }
827                         }
828                 }
829         }
830
831         timer.Stop();
832         timer.Print("ROctBoolDS");
833
834
835     // Test converter
836     TopoDS_Shape S = BRepPrimAPI_MakeSphere(100.0);
837
838     timer.Start();
839
840     int progress = 0;
841     Voxel_ROctBoolDS* ds2 = new Voxel_ROctBoolDS;
842     Voxel_FastConverter converter(S, *ds2, 0.1, myNbX, myNbY, myNbZ, 1);
843     converter.Convert(progress);
844     ds2->OptimizeMemory();
845
846         timer.Stop();
847         timer.Print("ROctBoolDS::converter");
848
849
850     // Display 
851     myViewer->getIC()->EraseAll(false, false);
852     initPrs();
853     myVoxels->SetBoolVoxels(0);
854         myVoxels->SetColorVoxels(0);
855     Handle(Poly_Triangulation) empty;
856     myVoxels->SetTriangulation(empty);
857     myVoxels->SetROctBoolVoxels(ds2);
858     myViewer->getIC()->Display(myVoxels, false);
859     myViewer->getView()->ColorScaleErase();
860         myViewer->getView()->FitAll();
861     myViewer->getSelector().SetVoxels(*ds2);
862 }
863
864 void Application::testFuseBoolDS()
865 {
866         Timer timer;
867         int ix, iy, iz;
868         int nbx = 100, nby = 100, nbz = 100;
869
870
871         // 1. Set two BoolDS:
872
873         timer.Start();
874
875         Voxel_BoolDS ds1(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
876         Voxel_BoolDS ds2(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
877
878         for (ix = 0; ix < nbx; ix++)
879         {
880                 for (iy = 0; iy < nby; iy++)
881                 {
882                         for (iz = 0; iz < nbz; iz++)
883                         {
884                                 if (ix & 0x01)
885                                         ds2.Set(ix, iy, iz, false);
886                                 else
887                                         ds2.Set(ix, iy, iz, true);
888                         }
889                 }
890         }
891
892     // 2. Fuse them
893
894     Voxel_BooleanOperation fuser;
895     if (!fuser.Fuse(ds1, ds2))
896         cout<<"The operation failed..."<<endl;
897
898     // 3. Check result
899
900         for (ix = 0; ix < nbx; ix++)
901         {
902                 for (iy = 0; iy < nby; iy++)
903                 {
904                         for (iz = 0; iz < nbz; iz++)
905                         {
906                                 bool value = ds1.Get(ix, iy, iz) == Standard_True;
907                                 if (ix & 0x01)
908                                 {
909                                         if (value != false)
910                                         cout<<"Wrong value!"<<endl;
911                                 }
912                                 else
913                                 {
914                                         if (value != true)
915                                                 cout<<"Wrong value!"<<endl;
916                                 }
917                         }
918                 }
919         }
920
921         timer.Stop();
922         timer.Print("Fusion of BoolDS");
923 }
924
925 void Application::testFuseColorDS()
926 {
927         Timer timer;
928         int ix, iy, iz;
929         int nbx = 100, nby = 100, nbz = 100;
930
931
932         // 1. Set two ColorDS:
933
934         timer.Start();
935
936         Voxel_ColorDS ds1(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
937         Voxel_ColorDS ds2(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
938
939         for (ix = 0; ix < nbx; ix++)
940         {
941                 for (iy = 0; iy < nby; iy++)
942                 {
943                         for (iz = 0; iz < nbz; iz++)
944                         {
945                                 ds1.Set(ix, iy, iz, 11);
946                         }
947                 }
948         }
949         for (ix = 0; ix < nbx; ix++)
950         {
951                 for (iy = 0; iy < nby; iy++)
952                 {
953                         for (iz = 0; iz < nbz; iz++)
954                         {
955                                 if (ix & 0x01)
956                                         ds2.Set(ix, iy, iz, 3);
957                                 else
958                                         ds2.Set(ix, iy, iz, 5);
959                         }
960                 }
961         }
962
963     // 2. Fuse them
964
965     Voxel_BooleanOperation fuser;
966     if (!fuser.Fuse(ds1, ds2))
967         cout<<"The operation failed..."<<endl;
968
969     // 3. Check result
970
971         for (ix = 0; ix < nbx; ix++)
972         {
973                 for (iy = 0; iy < nby; iy++)
974                 {
975                         for (iz = 0; iz < nbz; iz++)
976                         {
977                                 unsigned char value = ds1.Get(ix, iy, iz);
978                                 if (ix & 0x01)
979                                 {
980                                         if (value != 14)
981                                         cout<<"Wrong value!"<<endl;
982                                 }
983                                 else
984                                 {
985                                         if (value != 15)
986                                                 cout<<"Wrong value!"<<endl;
987                                 }
988                         }
989                 }
990         }
991
992         timer.Stop();
993         timer.Print("Fusion of ColorDS");
994 }
995
996 void Application::testFuseFloatDS()
997 {
998         Timer timer;
999         int ix, iy, iz;
1000         int nbx = 100, nby = 100, nbz = 100;
1001
1002
1003         // 1. Set two FloatDS:
1004
1005         timer.Start();
1006
1007         Voxel_FloatDS ds1(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1008         Voxel_FloatDS ds2(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1009
1010         for (ix = 0; ix < nbx; ix++)
1011         {
1012                 for (iy = 0; iy < nby; iy++)
1013                 {
1014                         for (iz = 0; iz < nbz; iz++)
1015                         {
1016                                 ds1.Set(ix, iy, iz, 11.1f);
1017                         }
1018                 }
1019         }
1020         for (ix = 0; ix < nbx; ix++)
1021         {
1022                 for (iy = 0; iy < nby; iy++)
1023                 {
1024                         for (iz = 0; iz < nbz; iz++)
1025                         {
1026                                 if (ix & 0x01)
1027                                         ds2.Set(ix, iy, iz, 3.3f);
1028                                 else
1029                                         ds2.Set(ix, iy, iz, 5.5f);
1030                         }
1031                 }
1032         }
1033
1034     // 2. Fuse them
1035
1036     Voxel_BooleanOperation fuser;
1037     if (!fuser.Fuse(ds1, ds2))
1038         cout<<"The operation failed..."<<endl;
1039
1040     // 3. Check result
1041
1042         for (ix = 0; ix < nbx; ix++)
1043         {
1044                 for (iy = 0; iy < nby; iy++)
1045                 {
1046                         for (iz = 0; iz < nbz; iz++)
1047                         {
1048                                 float value = ds1.Get(ix, iy, iz);
1049                                 if (ix & 0x01)
1050                                 {
1051                     if (fabs(value - 14.4f) > 0.001)
1052                                         cout<<"Wrong value!"<<endl;
1053                                 }
1054                                 else
1055                                 {
1056                     if (fabs(value - 16.6f) > 0.001)
1057                                                 cout<<"Wrong value!"<<endl;
1058                                 }
1059                         }
1060                 }
1061         }
1062
1063         timer.Stop();
1064         timer.Print("Fusion of FloatDS");
1065 }
1066
1067 void Application::testCutBoolDS()
1068 {
1069         Timer timer;
1070         int ix, iy, iz;
1071         int nbx = 100, nby = 100, nbz = 100;
1072
1073
1074         // 1. Set two BoolDS:
1075
1076         timer.Start();
1077
1078         Voxel_BoolDS ds1(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1079         Voxel_BoolDS ds2(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1080
1081         for (ix = 0; ix < nbx; ix++)
1082         {
1083                 for (iy = 0; iy < nby; iy++)
1084                 {
1085                         for (iz = 0; iz < nbz; iz++)
1086                         {
1087                                 ds1.Set(ix, iy, iz, true);
1088                         }
1089                 }
1090         }
1091         for (ix = 0; ix < nbx; ix++)
1092         {
1093                 for (iy = 0; iy < nby; iy++)
1094                 {
1095                         for (iz = 0; iz < nbz; iz++)
1096                         {
1097                                 if (ix & 0x01)
1098                                         ds2.Set(ix, iy, iz, false);
1099                                 else
1100                                         ds2.Set(ix, iy, iz, true);
1101                         }
1102                 }
1103         }
1104
1105     // 2. Cut them
1106
1107     Voxel_BooleanOperation cutter;
1108     if (!cutter.Cut(ds1, ds2))
1109         cout<<"The operation failed..."<<endl;
1110
1111     // 3. Check result
1112
1113         for (ix = 0; ix < nbx; ix++)
1114         {
1115                 for (iy = 0; iy < nby; iy++)
1116                 {
1117                         for (iz = 0; iz < nbz; iz++)
1118                         {
1119                                 bool value = ds1.Get(ix, iy, iz) == Standard_True;
1120                                 if (ix & 0x01)
1121                                 {
1122                                         if (value != true)
1123                                         cout<<"Wrong value!"<<endl;
1124                                 }
1125                                 else
1126                                 {
1127                                         if (value != false)
1128                                                 cout<<"Wrong value!"<<endl;
1129                                 }
1130                         }
1131                 }
1132         }
1133
1134         timer.Stop();
1135         timer.Print("Cut of BoolDS");
1136 }
1137
1138 void Application::testCutColorDS()
1139 {
1140         Timer timer;
1141         int ix, iy, iz;
1142         int nbx = 100, nby = 100, nbz = 100;
1143
1144
1145         // 1. Set two ColorDS:
1146
1147         timer.Start();
1148
1149         Voxel_ColorDS ds1(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1150         Voxel_ColorDS ds2(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1151
1152         for (ix = 0; ix < nbx; ix++)
1153         {
1154                 for (iy = 0; iy < nby; iy++)
1155                 {
1156                         for (iz = 0; iz < nbz; iz++)
1157                         {
1158                                 ds1.Set(ix, iy, iz, 11);
1159                         }
1160                 }
1161         }
1162         for (ix = 0; ix < nbx; ix++)
1163         {
1164                 for (iy = 0; iy < nby; iy++)
1165                 {
1166                         for (iz = 0; iz < nbz; iz++)
1167                         {
1168                                 if (ix & 0x01)
1169                                         ds2.Set(ix, iy, iz, 3);
1170                                 else
1171                                         ds2.Set(ix, iy, iz, 5);
1172                         }
1173                 }
1174         }
1175
1176     // 2. Cut them
1177
1178     Voxel_BooleanOperation cutter;
1179     if (!cutter.Cut(ds1, ds2))
1180         cout<<"The operation failed..."<<endl;
1181
1182     // 3. Check result
1183
1184         for (ix = 0; ix < nbx; ix++)
1185         {
1186                 for (iy = 0; iy < nby; iy++)
1187                 {
1188                         for (iz = 0; iz < nbz; iz++)
1189                         {
1190                                 unsigned char value = ds1.Get(ix, iy, iz);
1191                                 if (ix & 0x01)
1192                                 {
1193                                         if (value != 8)
1194                                         cout<<"Wrong value!"<<endl;
1195                                 }
1196                                 else
1197                                 {
1198                                         if (value != 6)
1199                                                 cout<<"Wrong value!"<<endl;
1200                                 }
1201                         }
1202                 }
1203         }
1204
1205         timer.Stop();
1206         timer.Print("Cut of ColorDS");
1207 }
1208
1209 void Application::testCutFloatDS()
1210 {
1211         Timer timer;
1212         int ix, iy, iz;
1213         int nbx = 100, nby = 100, nbz = 100;
1214
1215
1216         // 1. Set two FloatDS:
1217
1218         timer.Start();
1219
1220         Voxel_FloatDS ds1(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1221         Voxel_FloatDS ds2(0, 0, 0, 1, 1, 1, nbx, nby, nbz);
1222
1223         for (ix = 0; ix < nbx; ix++)
1224         {
1225                 for (iy = 0; iy < nby; iy++)
1226                 {
1227                         for (iz = 0; iz < nbz; iz++)
1228                         {
1229                                 ds1.Set(ix, iy, iz, 11.1f);
1230                         }
1231                 }
1232         }
1233         for (ix = 0; ix < nbx; ix++)
1234         {
1235                 for (iy = 0; iy < nby; iy++)
1236                 {
1237                         for (iz = 0; iz < nbz; iz++)
1238                         {
1239                                 if (ix & 0x01)
1240                                         ds2.Set(ix, iy, iz, 3.3f);
1241                                 else
1242                                         ds2.Set(ix, iy, iz, 5.5f);
1243                         }
1244                 }
1245         }
1246
1247     // 2. Cut them
1248
1249     Voxel_BooleanOperation cutter;
1250     if (!cutter.Cut(ds1, ds2))
1251         cout<<"The operation failed..."<<endl;
1252
1253     // 3. Check result
1254
1255         for (ix = 0; ix < nbx; ix++)
1256         {
1257                 for (iy = 0; iy < nby; iy++)
1258                 {
1259                         for (iz = 0; iz < nbz; iz++)
1260                         {
1261                                 float value = ds1.Get(ix, iy, iz);
1262                                 if (ix & 0x01)
1263                                 {
1264                     if (fabs(value - 7.8f) > 0.001)
1265                                         cout<<"Wrong value!"<<endl;
1266                                 }
1267                                 else
1268                                 {
1269                     if (fabs(value - 5.6f) > 0.001)
1270                                                 cout<<"Wrong value!"<<endl;
1271                                 }
1272                         }
1273                 }
1274         }
1275
1276         timer.Stop();
1277         timer.Print("Cut of FloatDS");
1278 }
1279
1280 void Application::convert2bool()
1281 {
1282     convert(0);
1283 }
1284
1285 void Application::convert2color()
1286 {
1287     convert(1);
1288 }
1289
1290 void Application::convert(const int ivoxel)
1291 {
1292         TopoDS_Shape S;
1293         if (!myShape.IsNull())
1294                 S = myShape->Shape();
1295     if (S.IsNull())
1296         {
1297                 QMessageBox::warning( this, "Voxel demo-application", "No shape for conversion!");
1298                 return;
1299         }
1300
1301     switch (ivoxel)
1302     {
1303         case 0:
1304         {
1305                 if (!myBoolVoxels)
1306                         myBoolVoxels = new Voxel_BoolDS;
1307             if (myColorVoxels)
1308             {
1309                 delete myColorVoxels;
1310                 myColorVoxels = 0;
1311             }
1312             break;
1313         }
1314         case 1:
1315         {
1316                 if (!myColorVoxels)
1317                         myColorVoxels = new Voxel_ColorDS;
1318             if (myBoolVoxels)
1319             {
1320                 delete myBoolVoxels;
1321                 myBoolVoxels = 0;
1322             }
1323             break;
1324         }
1325     }
1326
1327     switch (ivoxel)
1328     {
1329         case 0:
1330         {
1331             Timer timer;
1332             timer.Start();
1333
1334             /*
1335             int progress;
1336             Voxel_Converter converter(S, *myBoolVoxels, myNbX, myNbY, myNbZ);
1337                 if (!converter.Convert(progress, myVolumicBoolValue, myScanSide))
1338                 {
1339                         QMessageBox::warning( this, "Voxel demo-application", "Conversion failed...");
1340                         return;
1341                 }
1342             */
1343
1344             /*
1345             Voxel_Converter converter(S, *myBoolVoxels, myNbX, myNbY, myNbZ, 2);
1346             ConversionThread thread1, thread2;
1347             
1348             thread1.setConverter(&converter);
1349             thread2.setConverter(&converter);
1350
1351             thread1.setVolumicValue(myVolumicBoolValue);
1352             thread2.setVolumicValue(myVolumicBoolValue);
1353
1354             thread1.setScanSide(myScanSide);
1355             thread2.setScanSide(myScanSide);
1356
1357             thread1.setThreadIndex(1);
1358             thread2.setThreadIndex(2);
1359
1360             thread1.start();
1361             thread2.start();
1362
1363             while (thread1.running() || thread2.running())
1364             {
1365                 ::Sleep(100);
1366             }
1367             */
1368
1369             /*
1370             int progress;
1371             Voxel_FastConverter converter(S, *myBoolVoxels, 0.1, myNbX, myNbY, myNbZ, 1);
1372             converter.Convert(progress, 1);
1373             //if (myVolumicBoolValue)
1374             //    converter.FillInVolume(myVolumicBoolValue);
1375             */
1376
1377             Voxel_FastConverter converter(S, *myBoolVoxels, 0.1, myNbX, myNbY, myNbZ, 2);
1378             ConversionThread thread1, thread2;
1379             
1380             thread1.setConverter(&converter);
1381             thread2.setConverter(&converter);
1382
1383             thread1.setThreadIndex(1);
1384             thread2.setThreadIndex(2);
1385
1386             thread1.start();
1387             thread2.start();
1388
1389             while (thread1.isRunning() || thread2.isRunning())
1390             {
1391                 ::Sleep(100);
1392             }
1393
1394             timer.Print("Converter");
1395
1396             myViewer->getSelector().SetVoxels(*myBoolVoxels);
1397             break;
1398         }
1399         case 1:
1400         {
1401
1402             Timer timer;
1403             timer.Start();
1404
1405             /*
1406             int progress;
1407             Voxel_Converter converter(S, *myColorVoxels, myNbX, myNbY, myNbZ);
1408                 if (!converter.Convert(progress, myVolumicColorValue, myScanSide))
1409                 {
1410                         QMessageBox::warning( this, "Voxel demo-application", "Conversion failed...");
1411                         return;
1412                 }
1413             */
1414
1415             /*
1416             Voxel_Converter converter(S, *myColorVoxels, myNbX, myNbY, myNbZ, 2);
1417             ConversionThread thread1, thread2;
1418             
1419             thread1.setConverter(&converter);
1420             thread2.setConverter(&converter);
1421
1422             thread1.setVolumicValue(myVolumicColorValue);
1423             thread2.setVolumicValue(myVolumicColorValue);
1424
1425             thread1.setScanSide(myScanSide);
1426             thread2.setScanSide(myScanSide);
1427
1428             thread1.setThreadIndex(1);
1429             thread2.setThreadIndex(2);
1430
1431             thread1.start();
1432             thread2.start();
1433
1434             while (thread1.running() || thread2.running())
1435             {
1436                 ::Sleep(100);
1437             }
1438             */
1439
1440             /*
1441             int progress;
1442             Voxel_FastConverter converter(S, *myColorVoxels, myNbX, myNbY, myNbZ, 1);
1443             converter.Convert(progress, 1);
1444             if (myVolumicColorValue)
1445                 converter.FillInVolume(myVolumicColorValue);
1446             */
1447
1448             Voxel_FastConverter converter(S, *myColorVoxels, 0.1, myNbX, myNbY, myNbZ, 2);
1449             ConversionThread thread1, thread2;
1450             
1451             thread1.setConverter(&converter);
1452             thread2.setConverter(&converter);
1453
1454             thread1.setThreadIndex(1);
1455             thread2.setThreadIndex(2);
1456
1457             thread1.start();
1458             thread2.start();
1459
1460             while (thread1.isRunning() || thread2.isRunning())
1461             {
1462                 ::Sleep(100);
1463             }
1464
1465             timer.Print("Converter");
1466
1467
1468             // Set color for demonstration
1469             double maxd = 
1470                 fabs(myColorVoxels->GetX()) > fabs(myColorVoxels->GetY()) ? 
1471                 fabs(myColorVoxels->GetX()) : fabs(myColorVoxels->GetY());
1472             maxd = maxd > fabs(myColorVoxels->GetZ()) ? maxd : fabs(myColorVoxels->GetZ());
1473             maxd = maxd > fabs(myColorVoxels->GetX() + myColorVoxels->GetXLen()) ?
1474                 maxd : fabs(myColorVoxels->GetX() + myColorVoxels->GetXLen());
1475             maxd = maxd > fabs(myColorVoxels->GetY() + myColorVoxels->GetYLen()) ?
1476                 maxd : fabs(myColorVoxels->GetY() + myColorVoxels->GetYLen());
1477             maxd = maxd > fabs(myColorVoxels->GetZ() + myColorVoxels->GetZLen()) ?
1478                 maxd : fabs(myColorVoxels->GetZ() + myColorVoxels->GetZLen());
1479             for (int ix = 0; ix < myNbX; ix++)
1480             {
1481                 for (int iy = 0; iy < myNbY; iy++)
1482                 {
1483                     for (int iz = 0; iz < myNbZ; iz++)
1484                     {
1485                         unsigned char value = myColorVoxels->Get(ix, iy, iz);
1486                         if (value)
1487                         {
1488                             double xc, yc, zc, xd, yd, zd;
1489                             myColorVoxels->GetCenter(ix, iy, iz, xc, yc, zc);
1490                             xd = fabs(xc);
1491                             yd = fabs(yc);
1492                             zd = fabs(zc);
1493                             double mind = xd < yd ? xd : yd;
1494                             mind = zd < mind ? zd : mind;
1495                             value = unsigned char(15.0 * (maxd - mind) / maxd);
1496                             if (value <= 0)
1497                                 value = 1;
1498                             myColorVoxels->Set(ix, iy, iz, value);
1499                         }
1500                     }
1501                 }
1502             }
1503
1504             myViewer->getSelector().SetVoxels(*myColorVoxels);
1505             break;
1506         }
1507     }
1508
1509     myViewer->getIC()->EraseAll(false, false);
1510
1511     Voxel_DS* ds = myBoolVoxels;
1512     if (!ds)
1513         ds = myColorVoxels;
1514     if (ds)
1515     {
1516         myDisplayedXMin = ds->GetX() - 10.0 * Precision::Confusion();
1517         myDisplayedXMax = ds->GetX() + ds->GetXLen() + 10.0 * Precision::Confusion();
1518         myDisplayedYMin = ds->GetY() - 10.0 * Precision::Confusion();
1519         myDisplayedYMax = ds->GetY() + ds->GetYLen() + 10.0 * Precision::Confusion();
1520         myDisplayedZMin = ds->GetZ() - 10.0 * Precision::Confusion();
1521         myDisplayedZMax = ds->GetZ() + ds->GetZLen() + 10.0 * Precision::Confusion();
1522     }
1523
1524     // Init visual data
1525     initPrs();
1526
1527     // Set voxels and display
1528     Handle(Poly_Triangulation) empty;
1529     myVoxels->SetBoolVoxels(myBoolVoxels);
1530         myVoxels->SetColorVoxels(myColorVoxels);
1531     myVoxels->SetTriangulation(empty);
1532     if (myViewer->getIC()->IsDisplayed(myVoxels))
1533         myViewer->getIC()->Redisplay(myVoxels, false);
1534     else
1535             myViewer->getIC()->Display(myVoxels, false);
1536     
1537     // Color scale
1538     if (myColorVoxels)
1539         displayColorScale();
1540     else
1541         myViewer->getView()->ColorScaleErase();
1542
1543         myViewer->getView()->FitAll();
1544 }
1545
1546 void Application::setNbX()
1547 {
1548         bool ok;
1549         myNbX = 
1550                 QInputDialog::getInteger(this, "Voxel demo-application", "Number of splits in X-direction:", myNbX,
1551                                                                                                                  1, 100000, 1, &ok);
1552 }
1553
1554 void Application::setNbY()
1555 {
1556         bool ok;
1557         myNbY = 
1558                 QInputDialog::getInteger(this, "Voxel demo-application", "Number of splits in X-direction:", myNbY,
1559                                                                                                                  1, 100000, 1, &ok);
1560 }
1561
1562 void Application::setNbZ()
1563 {
1564         bool ok;
1565         myNbZ = 
1566                 QInputDialog::getInteger(this, "Voxel demo-application", "Number of splits in X-direction:", myNbZ,
1567                                                                                                                  1, 100000, 1, &ok);
1568 }
1569
1570 void Application::setColorMinValue()
1571 {
1572         bool ok;
1573         myColorMinValue = 
1574                 QInputDialog::getInteger(this, "Voxel demo-application", "Minimum value for color [0 .. 15]:", myColorMinValue,
1575                                                                  0, 15, 1, &ok);
1576     if (!myVoxels.IsNull())
1577         myVoxels->SetColorRange(myColorMinValue, myColorMaxValue);
1578 }
1579
1580 void Application::setColorMaxValue()
1581 {
1582         bool ok;
1583         myColorMaxValue = 
1584                 QInputDialog::getInteger(this, "Voxel demo-application", "Maximum value for color [0 .. 15]:", myColorMaxValue,
1585                                                                  0, 15, 1, &ok);
1586     if (!myVoxels.IsNull())
1587         myVoxels->SetColorRange(myColorMinValue, myColorMaxValue);
1588 }
1589
1590 void Application::setUsageOfGLlists()
1591 {
1592     int res = QMessageBox::question( this, "Voxel demo-application", "Press Yes to use GL lists and No not to use them.", QMessageBox::Yes, QMessageBox::No);
1593     if (!myVoxels.IsNull())
1594         myVoxels->SetUsageOfGLlists(res == QMessageBox::Yes);
1595 }
1596
1597 void Application::setDisplayedXMin()
1598 {
1599         myDisplayedXMin = QInputDialog::getDouble(this, "Voxel demo-application", "Minimum X value:", myDisplayedXMin);
1600     if (!myVoxels.IsNull())
1601     {
1602         myVoxels->SetSizeRange(myDisplayedXMin, myDisplayedXMax,
1603                                myDisplayedYMin, myDisplayedYMax,
1604                                myDisplayedZMin, myDisplayedZMax);
1605     }
1606 }
1607
1608 void Application::setDisplayedXMax()
1609 {
1610         myDisplayedXMax = QInputDialog::getDouble(this, "Voxel demo-application", "Maximum X value:", myDisplayedXMax);
1611     if (!myVoxels.IsNull())
1612     {
1613         myVoxels->SetSizeRange(myDisplayedXMin, myDisplayedXMax,
1614                                myDisplayedYMin, myDisplayedYMax,
1615                                myDisplayedZMin, myDisplayedZMax);
1616     }
1617 }
1618
1619 void Application::setDisplayedYMin()
1620 {
1621         myDisplayedYMin = QInputDialog::getDouble(this, "Voxel demo-application", "Minimum Y value:", myDisplayedYMin);
1622     if (!myVoxels.IsNull())
1623     {
1624         myVoxels->SetSizeRange(myDisplayedXMin, myDisplayedXMax,
1625                                myDisplayedYMin, myDisplayedYMax,
1626                                myDisplayedZMin, myDisplayedZMax);
1627     }
1628 }
1629
1630 void Application::setDisplayedYMax()
1631 {
1632         myDisplayedYMax = QInputDialog::getDouble(this, "Voxel demo-application", "Maximum Y value:", myDisplayedYMax);
1633     if (!myVoxels.IsNull())
1634     {
1635         myVoxels->SetSizeRange(myDisplayedXMin, myDisplayedXMax,
1636                                myDisplayedYMin, myDisplayedYMax,
1637                                myDisplayedZMin, myDisplayedZMax);
1638     }
1639 }
1640
1641 void Application::setDisplayedZMin()
1642 {
1643         myDisplayedZMin = QInputDialog::getDouble(this, "Voxel demo-application", "Minimum Z value:", myDisplayedZMin);
1644     if (!myVoxels.IsNull())
1645     {
1646         myVoxels->SetSizeRange(myDisplayedXMin, myDisplayedXMax,
1647                                myDisplayedYMin, myDisplayedYMax,
1648                                myDisplayedZMin, myDisplayedZMax);
1649     }
1650 }
1651
1652 void Application::setDisplayedZMax()
1653 {
1654         myDisplayedZMax = QInputDialog::getDouble(this, "Voxel demo-application", "Maximum Z value:", myDisplayedZMax);
1655     if (!myVoxels.IsNull())
1656     {
1657         myVoxels->SetSizeRange(myDisplayedXMin, myDisplayedXMax,
1658                                myDisplayedYMin, myDisplayedYMax,
1659                                myDisplayedZMin, myDisplayedZMax);
1660     }
1661 }
1662
1663 void Application::setScanSide()
1664 {
1665         myScanSide = 
1666                 QInputDialog::getInteger(this, "Voxel demo-application", "Side of scanning (1: +X side, 2: +Y side, 3: +Z side, 4: +X & +Y sides, .. 7: +X, +Y,& +Z sides):", 
1667                                  myScanSide, 1, 7, 1);
1668 }
1669
1670 void Application::setVolumicBoolValue()
1671 {
1672         myVolumicBoolValue = 
1673                 QInputDialog::getInteger(this, "Voxel demo-application", "Volumic value on voxelization [0 .. 1]:", 
1674                                  myVolumicBoolValue, 0, 1, 1);
1675 }
1676
1677 void Application::setVolumicColorValue()
1678 {
1679         myVolumicColorValue = 
1680                 QInputDialog::getInteger(this, "Voxel demo-application", "Volumic value on voxelization [0 .. 15]:", 
1681                                  myVolumicColorValue, 0, 15, 1);
1682 }
1683
1684 void Application::setQuadrangleSize()
1685 {
1686         myQuadrangleSize = 
1687                 QInputDialog::getInteger(this, "Voxel demo-application", "Size of quadrangles (0% .. 100%):", 
1688                                  myQuadrangleSize, 1, 100, 10);
1689     if (!myVoxels.IsNull())
1690     {
1691         myVoxels->SetQuadrangleSize(myQuadrangleSize);
1692     }
1693 }
1694
1695 void Application::setPointSize()
1696 {
1697         myPointSize = 
1698                 QInputDialog::getInteger(this, "Voxel demo-application", "Size of points (1 .. 10):", 
1699                                  myPointSize, 1, 10, 1);
1700     if (!myVoxels.IsNull())
1701     {
1702         myVoxels->SetPointSize(myPointSize);
1703     }
1704 }
1705
1706 void Application::display(Voxel_VoxelDisplayMode mode)
1707 {
1708     if (myVoxels.IsNull() || !myViewer->getIC()->IsDisplayed(myVoxels))
1709     {
1710                 QMessageBox::warning( this, "Voxel demo-application", "Voxels are not displayed");
1711                 return;
1712         }
1713
1714     myVoxels->SetDisplayMode(mode);
1715
1716     if (myColorVoxels)
1717         displayColorScale();
1718     else
1719         myViewer->getView()->ColorScaleErase();
1720
1721     myViewer->getIC()->Redisplay(myVoxels, true);
1722 }
1723
1724 void Application::displayPoints()
1725 {
1726     display(Voxel_VDM_POINTS);
1727 }
1728
1729 void Application::displayNearestPoints()
1730 {
1731     display(Voxel_VDM_NEARESTPOINTS);
1732 }
1733
1734 void Application::displayBoxes()
1735 {
1736     display(Voxel_VDM_BOXES);
1737 }
1738
1739 void Application::displayNearestBoxes()
1740 {
1741     display(Voxel_VDM_NEARESTBOXES);
1742 }
1743
1744 void Application::displayColorScale()
1745 {
1746         Handle(Aspect_ColorScale) color_scale = myViewer->getView()->ColorScale();
1747         if (!color_scale.IsNull())
1748         {
1749         int nb_colors = 1<<4 /* 4 bits */;
1750         color_scale->SetRange(0, nb_colors - 1);
1751         color_scale->SetNumberOfIntervals(nb_colors);
1752                 color_scale->SetPosition(0.01, 0.5 - 0.01);
1753                 color_scale->SetSize(0.5 - 0.01, 0.5 - 0.01);
1754         }
1755     myViewer->getView()->ColorScaleDisplay();
1756 }
1757
1758 void Application::displayWaves()
1759 {
1760     myViewer->getIC()->EraseAll(false, false);
1761
1762     // Make voxels
1763     if (myBoolVoxels)
1764     {
1765         delete myBoolVoxels;
1766         myBoolVoxels = 0;
1767     }
1768     if (myColorVoxels)
1769         delete myColorVoxels;
1770
1771     int nbx = 500, nby = 50, nbz = 50;
1772     double xlen = 100.0, ylen = 100.0, zlen = 20.0;
1773     double dx = xlen / (double) nbx, dy = ylen / (double) nby, dz = zlen / (double) nbz;
1774     myColorVoxels = new Voxel_ColorDS(0.0, 0.0, 0.0, xlen, ylen, zlen, nbx, nby, nbz);
1775
1776     // Initial state - no colors
1777     int ix, iy, iz;
1778     for (ix = 0; ix < nbx; ix++)
1779     {
1780         for (iy = 0; iy < nby; iy++)
1781         {
1782             for (iz = 0; iz < nbz; iz++)
1783             {
1784                 myColorVoxels->Set(ix, iy, iz, 0);
1785             }
1786         }
1787     }
1788
1789     // Init visual data
1790     initPrs();
1791     myVoxels->SetDisplayMode(Voxel_VDM_POINTS);
1792     myVoxels->SetUsageOfGLlists(false);
1793         myVoxels->SetBoolVoxels(myBoolVoxels);
1794         myVoxels->SetColorVoxels(myColorVoxels);
1795     if (myViewer->getIC()->IsDisplayed(myVoxels))
1796         myViewer->getIC()->Redisplay(myVoxels, false);
1797     else
1798             myViewer->getIC()->Display(myVoxels, false);
1799     myViewer->getView()->FitAll();
1800
1801     // Prepare arrays of values
1802     // X&Z values
1803     int i = 0, di = 5 /* nb waves */;
1804     int* zvalues = new int[nbx];
1805     unsigned char* xvalues = new unsigned char[nbx];
1806     for (ix = 0; ix < nbx; ix++, i += di)
1807     {
1808         if (i > nbx || i < 0)
1809         {
1810             di *= -1;
1811             i += di;
1812         }
1813         double rad = -M_PI / 2.0 + double(i) / (double) nbx * M_PI;
1814         double c = cos(rad);
1815         xvalues[ix] = 15.0 * c;
1816         if (xvalues[ix] == 0)
1817             xvalues[ix] = 1;
1818         zvalues[ix] = (nbz - 2) * c;
1819     }
1820
1821     // Make waves
1822     unsigned char value = 0;
1823     for (i = 0; i <= 100; i++)
1824     {
1825         for (ix = 0; ix < nbx; ix++)
1826         {
1827             int ixi = ix + i;
1828             if (ixi >= nbx)
1829                 ixi -= nbx;
1830             for (iz = 0; iz < nbz; iz++)
1831             {
1832                 value = 0;
1833                 if (iz < zvalues[ixi])
1834                     value = xvalues[ixi];
1835                 for (iy = 0; iy < nby; iy++)
1836                 {
1837                     myColorVoxels->Set(ix, iy, iz, value);
1838                 }
1839             }
1840         }
1841         myViewer->getIC()->Redisplay(myVoxels, true);
1842         qApp->processEvents();
1843     }
1844
1845     delete[] xvalues;
1846     delete[] zvalues;
1847 }
1848
1849 void Application::initPrs()
1850 {
1851         if (myVoxels.IsNull())
1852     {
1853                 myVoxels = new Voxel_Prs;
1854         myVoxels->SetDisplayMode(Voxel_VDM_POINTS);
1855         myVoxels->SetColor(Quantity_NOC_WHITE);
1856         myVoxels->SetPointSize(1.0);
1857         myVoxels->SetSmoothPoints(false);
1858         myVoxels->SetQuadrangleSize(myQuadrangleSize);
1859         myVoxels->SetColorRange(myColorMinValue, myColorMaxValue);
1860         // Colors of ColorDS
1861         int nb_colors = 16 /* 4 bits */;
1862         Handle(Quantity_HArray1OfColor) colors = new Quantity_HArray1OfColor(0, nb_colors - 1);
1863         for (int icolor = 0; icolor < nb_colors; icolor++)
1864         {
1865             Quantity_Color color;
1866             Aspect_ColorScale::FindColor(icolor, 0, nb_colors - 1, nb_colors, color);
1867             colors->SetValue(icolor, color);
1868         }
1869         myVoxels->SetColors(colors);
1870         myViewer->setPrs(myVoxels);
1871     }
1872     else
1873     {
1874         myViewer->getIC()->RecomputePrsOnly(myVoxels, false);
1875     }
1876 }
1877
1878 void Application::box()
1879 {
1880     gp_Ax2 axes(gp_Pnt(0, 0, 0), gp_Dir(0, 0, 1));
1881     TopoDS_Shape S = BRepPrimAPI_MakeBox(axes, 100, 100, 100);
1882     load(S);
1883 }
1884
1885 void Application::cylinder()
1886 {
1887     TopoDS_Shape S = BRepPrimAPI_MakeCylinder(50, 100);
1888     load(S);
1889 }
1890
1891 void Application::torus()
1892 {
1893     TopoDS_Shape S = BRepPrimAPI_MakeTorus(100, 20);
1894     load(S);
1895 }
1896
1897 void Application::sphere()
1898 {
1899     TopoDS_Shape S = BRepPrimAPI_MakeSphere(100);
1900     load(S);
1901 }
1902
1903 void Application::load(const TopoDS_Shape& S)
1904 {
1905         myViewer->getIC()->EraseAll(false, false);
1906
1907     // Delete voxels of previous shape.
1908     if (myBoolVoxels)
1909     {
1910         delete myBoolVoxels;
1911         myBoolVoxels = 0;
1912     }
1913     if (myColorVoxels)
1914     {
1915         delete myColorVoxels;
1916         myColorVoxels = 0;
1917     }
1918     
1919     // Set view size
1920     Bnd_Box box;
1921     double xmin, ymin, zmin, xmax, ymax, zmax, length = 0;
1922     BRepBndLib::Add(S, box);
1923     box.Get(xmin, ymin, zmin, xmax, ymax, zmax);
1924     length = xmax - xmin > ymax - ymin ? xmax - xmin : ymax - ymin;
1925     length = length > zmax - zmin ? length : zmax - zmin;
1926     length *= 2.0;
1927     myViewer->getView()->SetSize(length);
1928     myViewer->getView()->SetZSize(length);
1929
1930     // Display shape
1931         if (myShape.IsNull())
1932         {
1933                 myShape = new AIS_Shape(S);
1934                 myShape->SetDisplayMode(1);
1935         myShape->UnsetSelectionMode();
1936         }
1937         else
1938         {
1939                 myShape->Set(S);
1940         myViewer->getIC()->RecomputePrsOnly(myShape, false);
1941         }
1942     if (myViewer->getIC()->IsDisplayed(myShape))
1943         myViewer->getIC()->Redisplay(myShape, false);
1944     else
1945             myViewer->getIC()->Display(myShape, false);
1946         myViewer->getView()->FitAll();
1947 }
1948
1949 void Application::displayCut()
1950 {
1951     myViewer->getIC()->EraseAll(false, false);
1952
1953     // Make a sphere with a lot of toruses, 
1954     // cut the toruses from the sphere.
1955     TopoDS_Shape sphere = BRepPrimAPI_MakeSphere(100.0);
1956     TopoDS_Shape torus1 = BRepPrimAPI_MakeTorus(gp_Ax2(gp_Pnt( 80,   0, 20), gp::DZ()), 30, 10);
1957     TopoDS_Shape torus2 = BRepPrimAPI_MakeTorus(gp_Ax2(gp_Pnt(  0,  80, 20), gp::DZ()), 30, 10);
1958     TopoDS_Shape torus3 = BRepPrimAPI_MakeTorus(gp_Ax2(gp_Pnt(-80,   0, 20), gp::DZ()), 30, 10);
1959     TopoDS_Shape torus4 = BRepPrimAPI_MakeTorus(gp_Ax2(gp_Pnt(  0, -80, 20), gp::DZ()), 30, 10);
1960
1961     // Compute bounding box of the shapes
1962     Bnd_Box box;
1963     BRepBndLib::Add(sphere, box);
1964     BRepBndLib::Add(torus1, box);
1965     BRepBndLib::Add(torus2, box);
1966     BRepBndLib::Add(torus3, box);
1967     BRepBndLib::Add(torus4, box);
1968
1969     // Nullify voxels
1970     if (myColorVoxels)
1971     {
1972         delete myColorVoxels;
1973         myColorVoxels = 0;
1974     }
1975     if (myBoolVoxels)
1976     {
1977         delete myBoolVoxels;
1978         myBoolVoxels = 0;
1979     }
1980
1981     Timer timer;
1982     timer.Start();
1983
1984     // Create a cube of voxels
1985     int nbx = 100, nby = 100, nbz = 100;
1986     double xmin, ymin, zmin, xmax, ymax, zmax;
1987     box.Get(xmin, ymin, zmin, xmax, ymax, zmax);
1988     myColorVoxels = new Voxel_ColorDS(xmin, ymin, zmin, 
1989                                       xmax - xmin, ymax - ymin, zmax - zmin,
1990                                       nbx, nby, nbz);
1991     Voxel_ColorDS vtorus(xmin, ymin, zmin, 
1992                          xmax - xmin, ymax - ymin, zmax - zmin,
1993                          nbx, nby, nbz);
1994
1995     // Make a cube of voxels for the sphere.
1996     int progress;
1997     Voxel_FastConverter converter(sphere, *myColorVoxels, 0.1, nbx, nby, nbz);
1998     converter.Convert(progress);
1999     converter.FillInVolume(15);
2000    
2001     // Torus 1
2002     Voxel_FastConverter converter1(torus1, vtorus, 0.1, nbx, nby, nbz);
2003     converter1.Convert(progress);
2004     converter1.FillInVolume(3);
2005    
2006     // Torus 2
2007     Voxel_FastConverter converter2(torus2, vtorus, 0.1, nbx, nby, nbz);
2008     converter2.Convert(progress);
2009     converter2.FillInVolume(7);
2010    
2011     // Torus 3
2012     Voxel_FastConverter converter3(torus3, vtorus, 0.1, nbx, nby, nbz);
2013     converter3.Convert(progress);
2014     converter3.FillInVolume(10);
2015    
2016     // Torus 4
2017     Voxel_FastConverter converter4(torus4, vtorus, 0.1, nbx, nby, nbz);
2018     converter4.Convert(progress);
2019     converter4.FillInVolume(12);
2020
2021     // Cut
2022     Voxel_BooleanOperation cutter;
2023     cutter.Cut(*myColorVoxels, vtorus);
2024
2025     // Remove volumic voxels
2026     converter.FillInVolume(0);
2027
2028     timer.Stop();
2029     timer.Print("Cut");
2030
2031     // Display
2032     initPrs();
2033     myVoxels->SetDisplayMode(Voxel_VDM_POINTS);
2034     myVoxels->SetUsageOfGLlists(true);
2035         myVoxels->SetBoolVoxels(myBoolVoxels);
2036         myVoxels->SetColorVoxels(myColorVoxels);
2037     if (myViewer->getIC()->IsDisplayed(myVoxels))
2038         myViewer->getIC()->Redisplay(myVoxels, false);
2039     else
2040             myViewer->getIC()->Display(myVoxels, false);
2041     myViewer->getView()->FitAll();
2042 }
2043
2044 void Application::displayCollisions()
2045 {
2046     myViewer->getIC()->EraseAll(false, false);
2047
2048     // Make a big box with a lot of small spheres inside.
2049     double x = 0.0, y = 0.0, z = 0.0, xlen = 100.0, ylen = 100.0, zlen = 100.0, r = 10.0;
2050     gp_Pnt P1(x, y, z); // center point of moving sphere (S1).
2051     TopoDS_Shape B   = BRepPrimAPI_MakeBox(gp_Pnt(x-r, y-r, z-r), gp_Pnt(xlen+r, ylen+r, zlen+r));
2052     TopoDS_Shape S1  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       y,       z),       gp::DZ()), r / 2.0);
2053     TopoDS_Shape S2  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., y,       z),       gp::DZ()), r);
2054     TopoDS_Shape S3  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    y,       z),       gp::DZ()), r);
2055     TopoDS_Shape S4  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       ylen/2., z),       gp::DZ()), r);
2056     TopoDS_Shape S5  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., ylen/2., z),       gp::DZ()), r);
2057     TopoDS_Shape S6  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    ylen/2., z),       gp::DZ()), r);
2058     TopoDS_Shape S7  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       ylen,    z),       gp::DZ()), r);
2059     TopoDS_Shape S8  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., ylen,    z),       gp::DZ()), r);
2060     TopoDS_Shape S9  = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    ylen,    z),       gp::DZ()), r);
2061     TopoDS_Shape S10 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       y,       zlen/2.), gp::DZ()), r);
2062     TopoDS_Shape S11 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., y,       zlen/2.), gp::DZ()), r);
2063     TopoDS_Shape S12 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    y,       zlen/2.), gp::DZ()), r);
2064     TopoDS_Shape S13 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       ylen/2., zlen/2.), gp::DZ()), r);
2065     TopoDS_Shape S14 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., ylen/2., zlen/2.), gp::DZ()), r);
2066     TopoDS_Shape S15 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    ylen/2., zlen/2.), gp::DZ()), r);
2067     TopoDS_Shape S16 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       ylen,    zlen/2.), gp::DZ()), r);
2068     TopoDS_Shape S17 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., ylen,    zlen/2.), gp::DZ()), r);
2069     TopoDS_Shape S18 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    ylen,    zlen/2.), gp::DZ()), r);
2070     TopoDS_Shape S19 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       y,       zlen),    gp::DZ()), r);
2071     TopoDS_Shape S20 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., y,       zlen),    gp::DZ()), r);
2072     TopoDS_Shape S21 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    y,       zlen),    gp::DZ()), r);
2073     TopoDS_Shape S22 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       ylen/2., zlen),    gp::DZ()), r);
2074     TopoDS_Shape S23 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., ylen/2., zlen),    gp::DZ()), r);
2075     TopoDS_Shape S24 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    ylen/2., zlen),    gp::DZ()), r);
2076     TopoDS_Shape S25 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(x,       ylen,    zlen),    gp::DZ()), r);
2077     TopoDS_Shape S26 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen/2., ylen,    zlen),    gp::DZ()), r);
2078     TopoDS_Shape S27 = BRepPrimAPI_MakeSphere(gp_Ax2(gp_Pnt(xlen,    ylen,    zlen),    gp::DZ()), r);
2079
2080     // Planes of the big box
2081     gp_Ax2 xminusPlane(gp_Pnt(x,    y,    z),    gp::DX());
2082     gp_Ax2 xplusPlane (gp_Pnt(xlen, y,    z),    gp::DX());
2083     gp_Ax2 yminusPlane(gp_Pnt(x,    y,    z),    gp::DY());
2084     gp_Ax2 yplusPlane (gp_Pnt(x,    ylen, z),    gp::DY());
2085     gp_Ax2 zminusPlane(gp_Pnt(x,    y,    z),    gp::DZ());
2086     gp_Ax2 zplusPlane (gp_Pnt(x,    y,    zlen), gp::DZ());
2087
2088     // Nullify voxels
2089     if (myColorVoxels)
2090     {
2091         delete myColorVoxels;
2092         myColorVoxels = 0;
2093     }
2094     if (myBoolVoxels)
2095     {
2096         delete myBoolVoxels;
2097         myBoolVoxels = 0;
2098     }
2099
2100     // Prepare visualization
2101     initPrs();
2102     myVoxels->SetDisplayMode(Voxel_VDM_POINTS);
2103     myVoxels->SetColor(Quantity_NOC_RED);
2104     myVoxels->SetPointSize(4);
2105     myVoxels->SetSmoothPoints(false);
2106     myVoxels->SetUsageOfGLlists(false);
2107         myVoxels->SetColorVoxels(myColorVoxels);
2108
2109     // Display all shapes
2110     double transparency = 0.9;
2111     Handle(AIS_Shape) aisB   = new AIS_Shape(B);
2112     Handle(AIS_Shape) aisS1  = new AIS_Shape(S1);
2113     Handle(AIS_Shape) aisS2  = new AIS_Shape(S2);
2114     Handle(AIS_Shape) aisS3  = new AIS_Shape(S3);
2115     Handle(AIS_Shape) aisS4  = new AIS_Shape(S4);
2116     Handle(AIS_Shape) aisS5  = new AIS_Shape(S5);
2117     Handle(AIS_Shape) aisS6  = new AIS_Shape(S6);
2118     Handle(AIS_Shape) aisS7  = new AIS_Shape(S7);
2119     Handle(AIS_Shape) aisS8  = new AIS_Shape(S8);
2120     Handle(AIS_Shape) aisS9  = new AIS_Shape(S9);
2121     Handle(AIS_Shape) aisS10 = new AIS_Shape(S10);
2122     Handle(AIS_Shape) aisS11 = new AIS_Shape(S11);
2123     Handle(AIS_Shape) aisS12 = new AIS_Shape(S12);
2124     Handle(AIS_Shape) aisS13 = new AIS_Shape(S13);
2125     Handle(AIS_Shape) aisS14 = new AIS_Shape(S14);
2126     Handle(AIS_Shape) aisS15 = new AIS_Shape(S15);
2127     Handle(AIS_Shape) aisS16 = new AIS_Shape(S16);
2128     Handle(AIS_Shape) aisS17 = new AIS_Shape(S17);
2129     Handle(AIS_Shape) aisS18 = new AIS_Shape(S18);
2130     Handle(AIS_Shape) aisS19 = new AIS_Shape(S19);
2131     Handle(AIS_Shape) aisS20 = new AIS_Shape(S20);
2132     Handle(AIS_Shape) aisS21 = new AIS_Shape(S21);
2133     Handle(AIS_Shape) aisS22 = new AIS_Shape(S22);
2134     Handle(AIS_Shape) aisS23 = new AIS_Shape(S23);
2135     Handle(AIS_Shape) aisS24 = new AIS_Shape(S24);
2136     Handle(AIS_Shape) aisS25 = new AIS_Shape(S25);
2137     Handle(AIS_Shape) aisS26 = new AIS_Shape(S26);
2138     Handle(AIS_Shape) aisS27 = new AIS_Shape(S27);
2139     aisS1-> SetDisplayMode(1);
2140     aisS2-> SetDisplayMode(1);
2141     aisS3-> SetDisplayMode(1);
2142     aisS4-> SetDisplayMode(1);
2143     aisS5-> SetDisplayMode(1);
2144     aisS6-> SetDisplayMode(1);
2145     aisS7-> SetDisplayMode(1);
2146     aisS8-> SetDisplayMode(1);
2147     aisS9-> SetDisplayMode(1);
2148     aisS10->SetDisplayMode(1);
2149     aisS11->SetDisplayMode(1);
2150     aisS12->SetDisplayMode(1);
2151     aisS13->SetDisplayMode(1);
2152     aisS14->SetDisplayMode(1);
2153     aisS15->SetDisplayMode(1);
2154     aisS16->SetDisplayMode(1);
2155     aisS17->SetDisplayMode(1);
2156     aisS18->SetDisplayMode(1);
2157     aisS19->SetDisplayMode(1);
2158     aisS20->SetDisplayMode(1);
2159     aisS21->SetDisplayMode(1);
2160     aisS22->SetDisplayMode(1);
2161     aisS23->SetDisplayMode(1);
2162     aisS24->SetDisplayMode(1);
2163     aisS25->SetDisplayMode(1);
2164     aisS26->SetDisplayMode(1);
2165     aisS27->SetDisplayMode(1);
2166     aisB->  UnsetSelectionMode();
2167     aisS1-> UnsetSelectionMode();
2168     aisS2-> UnsetSelectionMode();
2169     aisS3-> UnsetSelectionMode();
2170     aisS4-> UnsetSelectionMode();
2171     aisS5-> UnsetSelectionMode();
2172     aisS6-> UnsetSelectionMode();
2173     aisS7-> UnsetSelectionMode();
2174     aisS8-> UnsetSelectionMode();
2175     aisS9-> UnsetSelectionMode();
2176     aisS10->UnsetSelectionMode();
2177     aisS11->UnsetSelectionMode();
2178     aisS12->UnsetSelectionMode();
2179     aisS13->UnsetSelectionMode();
2180     aisS14->UnsetSelectionMode();
2181     aisS15->UnsetSelectionMode();
2182     aisS16->UnsetSelectionMode();
2183     aisS17->UnsetSelectionMode();
2184     aisS18->UnsetSelectionMode();
2185     aisS19->UnsetSelectionMode();
2186     aisS20->UnsetSelectionMode();
2187     aisS21->UnsetSelectionMode();
2188     aisS22->UnsetSelectionMode();
2189     aisS23->UnsetSelectionMode();
2190     aisS24->UnsetSelectionMode();
2191     aisS25->UnsetSelectionMode();
2192     aisS26->UnsetSelectionMode();
2193     aisS27->UnsetSelectionMode();
2194     aisS1-> SetTransparency(2.0 * transparency / 3.0);
2195     aisS2-> SetTransparency(transparency);
2196     aisS3-> SetTransparency(transparency);
2197     aisS4-> SetTransparency(transparency);
2198     aisS5-> SetTransparency(transparency);
2199     aisS6-> SetTransparency(transparency);
2200     aisS7-> SetTransparency(transparency);
2201     aisS8-> SetTransparency(transparency);
2202     aisS9-> SetTransparency(transparency);
2203     aisS10->SetTransparency(transparency);
2204     aisS11->SetTransparency(transparency);
2205     aisS12->SetTransparency(transparency);
2206     aisS13->SetTransparency(transparency);
2207     aisS14->SetTransparency(transparency);
2208     aisS15->SetTransparency(transparency);
2209     aisS16->SetTransparency(transparency);
2210     aisS17->SetTransparency(transparency);
2211     aisS18->SetTransparency(transparency);
2212     aisS19->SetTransparency(transparency);
2213     aisS20->SetTransparency(transparency);
2214     aisS21->SetTransparency(transparency);
2215     aisS22->SetTransparency(transparency);
2216     aisS23->SetTransparency(transparency);
2217     aisS24->SetTransparency(transparency);
2218     aisS25->SetTransparency(transparency);
2219     aisS26->SetTransparency(transparency);
2220     aisS27->SetTransparency(transparency);
2221     myViewer->getIC()->Display(aisB, false);
2222     myViewer->getIC()->Display(aisS1, false);
2223     myViewer->getIC()->Display(aisS2, false);
2224     myViewer->getIC()->Display(aisS3, false);
2225     myViewer->getIC()->Display(aisS4, false);
2226     myViewer->getIC()->Display(aisS5, false);
2227     myViewer->getIC()->Display(aisS6, false);
2228     myViewer->getIC()->Display(aisS7, false);
2229     myViewer->getIC()->Display(aisS8, false);
2230     myViewer->getIC()->Display(aisS9, false);
2231     myViewer->getIC()->Display(aisS10, false);
2232     myViewer->getIC()->Display(aisS11, false);
2233     myViewer->getIC()->Display(aisS12, false);
2234     myViewer->getIC()->Display(aisS13, false);
2235     myViewer->getIC()->Display(aisS14, false);
2236     myViewer->getIC()->Display(aisS15, false);
2237     myViewer->getIC()->Display(aisS16, false);
2238     myViewer->getIC()->Display(aisS17, false);
2239     myViewer->getIC()->Display(aisS18, false);
2240     myViewer->getIC()->Display(aisS19, false);
2241     myViewer->getIC()->Display(aisS20, false);
2242     myViewer->getIC()->Display(aisS21, false);
2243     myViewer->getIC()->Display(aisS22, false);
2244     myViewer->getIC()->Display(aisS23, false);
2245     myViewer->getIC()->Display(aisS24, false);
2246     myViewer->getIC()->Display(aisS25, false);
2247     myViewer->getIC()->Display(aisS26, false);
2248     myViewer->getIC()->Display(aisS27, false);
2249
2250     // Prepare computer of collisions
2251     double deflection = 0.1;
2252     int nbx = 100, nby = 100, nbz = 100;
2253     
2254     Voxel_CollisionDetection coldet(deflection, nbx, nby, nbz);
2255     coldet.SetUsageOfVolume(false);
2256     coldet.KeepCollisions(false);
2257     coldet.AddShape(S1);
2258     coldet.AddShape(S2);
2259     coldet.AddShape(S3);
2260     coldet.AddShape(S4);
2261     coldet.AddShape(S5);
2262     coldet.AddShape(S6);
2263     coldet.AddShape(S7);
2264     coldet.AddShape(S8);
2265     coldet.AddShape(S9);
2266     coldet.AddShape(S10);
2267     coldet.AddShape(S11);
2268     coldet.AddShape(S12);
2269     coldet.AddShape(S13);
2270     coldet.AddShape(S14);
2271     coldet.AddShape(S15);
2272     coldet.AddShape(S16);
2273     coldet.AddShape(S17);
2274     coldet.AddShape(S18);
2275     coldet.AddShape(S19);
2276     coldet.AddShape(S20);
2277     coldet.AddShape(S21);
2278     coldet.AddShape(S22);
2279     coldet.AddShape(S23);
2280     coldet.AddShape(S24);
2281     coldet.AddShape(S25);
2282     coldet.AddShape(S26);
2283     coldet.AddShape(S27);
2284     //coldet.AddShape(BRepPrimAPI_MakeBox(gp_Pnt(x, y, z), gp_Pnt(xlen, ylen, zlen)));
2285     
2286     Bnd_Box box;
2287     BRepBndLib::Add(B, box);
2288     coldet.SetBoundaryBox(box);
2289
2290     coldet.Voxelize();
2291
2292     // Move one of the spheres inside the box
2293     // and compute collisions
2294     gp_Trsf trsf;
2295     gp_Vec vmove(1, 0.5, 0.25);
2296
2297     int imove = 0, nb_moves = 900;
2298     while (imove < nb_moves)
2299     {
2300         // Move
2301         trsf.SetTranslation(vmove);
2302         TopLoc_Location loc(trsf);
2303         S1.Move(loc);
2304         P1.Translate(vmove);
2305         
2306         // Check whether S1 is inside the big box
2307         // Detect the plane S1 touches to.
2308         if (P1.X() < x)
2309             vmove.Mirror(xminusPlane);
2310         else if (P1.X() > xlen)
2311             vmove.Mirror(xplusPlane);
2312         else if (P1.Y() < y)
2313             vmove.Mirror(yminusPlane);
2314         else if (P1.Y() > ylen)
2315             vmove.Mirror(yplusPlane);
2316         else if (P1.Z() < z)
2317             vmove.Mirror(zminusPlane);
2318         else if (P1.Z() > zlen)
2319             vmove.Mirror(zplusPlane);
2320
2321         // Compute collisions
2322         coldet.ReplaceShape(1, S1);
2323         coldet.Voxelize(1); // only the first sphere (S1)
2324         coldet.Compute();
2325         myBoolVoxels = &((Voxel_BoolDS&) coldet.GetCollisions());
2326
2327         // Redisplay S1
2328         aisS1->Set(S1);
2329         myViewer->getIC()->Redisplay(aisS1, false);
2330
2331         // Display the collisions
2332         myVoxels->SetBoolVoxels(myBoolVoxels);
2333         if (myViewer->getIC()->IsDisplayed(myVoxels))
2334             myViewer->getIC()->Redisplay(myVoxels, true);
2335         else
2336         {
2337                 myViewer->getIC()->Display(myVoxels, false);
2338             myViewer->getView()->FitAll();
2339         }
2340
2341         imove++;
2342         qApp->processEvents();
2343     }
2344
2345     // Copy the result of collision detection
2346     int ix, iy, iz;
2347     myBoolVoxels = new Voxel_BoolDS(coldet.GetCollisions().GetX(), 
2348                                     coldet.GetCollisions().GetY(), 
2349                                     coldet.GetCollisions().GetZ(), 
2350                                     coldet.GetCollisions().GetXLen(), 
2351                                     coldet.GetCollisions().GetYLen(), 
2352                                     coldet.GetCollisions().GetZLen(), 
2353                                     nbx, nby, nbz);
2354     for (ix = 0; ix < nbx; ix++)
2355     {
2356         for (iy = 0; iy < nby; iy++)
2357         {
2358             for (iz = 0; iz < nbz; iz++)
2359             {
2360                 if (coldet.GetCollisions().Get(ix, iy, iz))
2361                     myBoolVoxels->Set(ix, iy, iz, Standard_True);
2362             }
2363         }
2364     }
2365     myVoxels->SetBoolVoxels(myBoolVoxels);
2366 }