]> Creatis software - cpMesh.git/blob - appli/examples/example_SimplexMeshDeleteEdge.cxx
24fc7a6f6aec6842f3168fd0ed174f06d9651606
[cpMesh.git] / appli / examples / example_SimplexMeshDeleteEdge.cxx
1 #include <cmath>
2 #include <cstdlib>
3 #include <iostream>
4 #include <string>
5
6 #include <cpm/DataStructures/QuadEdgeMesh.h>
7 #include <cpm/DataStructures/SimplexMesh.h>
8 #include <cpm/Algorithms/QuadEdge/MeshToDualFilter.h>
9 #include <cpm/IO/MeshReader.h>
10
11 #include <vtkActor.h>
12 #include <vtkCallbackCommand.h>
13 #include <vtkProperty.h>
14 #include <vtkRenderer.h>
15 #include <vtkRenderWindow.h>
16 #include <vtkRenderWindowInteractor.h>
17 #include <vtkSmartPointer.h>
18 #include <cpm/VTK/MeshMapper.h>
19
20 // -------------------------------------------------------------------------
21 const unsigned int Dimension = 3;
22 const unsigned int Order = 2;
23
24 typedef float TScalar;
25 typedef cpm::DataStructures::QuadEdgeMesh< TScalar, Dimension > TMesh;
26 typedef cpm::DataStructures::SimplexMesh< TScalar, Order, Dimension > TSimplex;
27 typedef cpm::VTK::MeshMapper< TMesh > TMeshMapper;
28
29 // -------------------------------------------------------------------------
30 int main( int argc, char* argv[] )
31 {
32   if( argc < 2 )
33   {
34     std::cerr
35       << "Usage: " << argv[ 0 ]
36       << " input_mesh"
37       << std::endl;
38     return( 1 );
39
40   } // fi
41
42   typedef cpm::IO::MeshReader< TMesh > TReader;
43   TReader::Pointer reader = TReader::New( );
44   reader->SetFileName( argv[ 1 ] );
45
46   typedef
47     cpm::Algorithms::QuadEdge::MeshToDualFilter< TMesh, TSimplex >
48     TDualFilter;
49   TDualFilter::Pointer dual_filter = TDualFilter::New( );
50   dual_filter->SetInput( reader->GetOutput( ) );
51   dual_filter->Update( );
52   TSimplex::Pointer simplex = dual_filter->GetOutput( );
53
54   std::cout
55     << " --> Read " << simplex->GetNumberOfPoints( )
56     << " simplex points."
57     << std::endl;
58
59   bool finish = false;
60   while( !finish )
61   {
62     std::cout
63       << " --> Delete an edge between a and b (-1 to finish) : "
64       << std::flush;
65     int a, b;
66     std::cin >> a;
67     if( a == -1 )
68     {
69       finish = true;
70       continue;
71
72     } // fi
73     std::cin >> b;
74     if( simplex->OperatorDeleteEdge( a, b ) )
75       std::cout << "Deleted" << std::endl;
76     else
77       std::cout << "Not present" << std::endl;
78     
79   } // elihw
80
81   // Map mesh
82   vtkSmartPointer< TMeshMapper > mapper =
83     vtkSmartPointer< TMeshMapper >::New( );
84   mapper->SetInputData( reader->GetOutput( ) );
85
86   // Map dual mesh
87   vtkSmartPointer< TMeshMapper > dual_mapper =
88     vtkSmartPointer< TMeshMapper >::New( );
89   dual_mapper->SetInputData( dual_filter->GetOutput( ) );
90
91   // Create actor
92   vtkSmartPointer< vtkActor > actor =
93     vtkSmartPointer< vtkActor >::New( );
94   actor->SetMapper( mapper );
95   actor->GetProperty( )->SetColor( 1, 1, 0 );
96   actor->GetProperty( )->SetOpacity( 0.5 );
97
98   // Create dual actor
99   vtkSmartPointer< vtkActor > dual_actor =
100     vtkSmartPointer< vtkActor >::New( );
101   dual_actor->SetMapper( dual_mapper );
102   dual_actor->GetProperty( )->SetColor( 1, 0, 0 );
103   dual_actor->GetProperty( )->SetOpacity( 1 );
104   dual_actor->GetProperty( )->SetLineWidth( 2 );
105   dual_actor->GetProperty( )->SetRepresentationToWireframe( );
106
107   // Configure visualization objects
108   vtkSmartPointer< vtkRenderer > renderer =
109     vtkSmartPointer< vtkRenderer >::New( );
110   renderer->SetBackground( 0.1, 0.3, 0.5 );
111
112   vtkSmartPointer< vtkRenderWindow > window =
113     vtkSmartPointer< vtkRenderWindow >::New( );
114   window->AddRenderer( renderer );
115   window->SetSize( 800, 800 );
116
117   // Set up the interaction
118   vtkSmartPointer< vtkRenderWindowInteractor > interactor =
119     vtkSmartPointer< vtkRenderWindowInteractor >::New( );
120   window->SetInteractor( interactor );
121
122   // Associate actors
123   //renderer->AddActor( actor );
124   renderer->AddActor( dual_actor );
125
126   // Begin interaction
127   window->Render( );
128   interactor->Start( );
129
130   return( 0 );
131 }
132
133 // eof - $RCSfile$